FV3核心GFS升级后metgrid错误的问题解决

巩崇水 · 发表于 2019-6-17 08:37:22

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由巩崇水于 2019-6-17 08:46 编辑

问题分析
FV3升级后，GFS资料增加了三个层，分别是40,15和0.4hPa。在ungrib之后就会造成层数不匹配，和Metgrid也会冲突
如果直接按if ( level .lt. pmin ) cycle MATCH_LOOP的思路，则会造成real的时候层数小于32而无法正常运行的问题。

问题解决
修改WPS/ungrib/src/rd_grib2.F然后重新编译WPS即可，可以起到略去相关三层GFS读入的作用。
以3.9版本为例，改版本rd_grib2
大约在740行前后搜索如下字段，然后添加红色内容即可
! Level (eg. 10000 mb)
            if(gfld%ipdtmpl(10).eq.100) then
               ! Pressure level (range from 1000mb to 0mb)
               level=gfld%ipdtmpl(12) *
   &                         (10. ** (-1. * gfld%ipdtmpl(11)))
               !if ( level .lt. 5000 ) cycle MATCH_LOOP
               if ( level .eq. 4000 ) cycle MATCH_LOOP
               if ( level .eq. 1500 ) cycle MATCH_LOOP
               if ( level .eq. 40 ) cycle MATCH_LOOP
               !==================================================
            elseif((gfld%ipdtmpl(10).eq.105).or.
   &             (gfld%ipdtmpl(10).eq.118))then

思路原网页，感谢 DCVZ
http://forum.wrfforum.com/viewtopic.php?f=6&t=11314&sid=3280b397e5005fa419704c9476b0d9b2

if(gfld%ipdtmpl(10).eq.100) then
! Pressure level (range from 1000mb to 0mb)
level=gfld%ipdtmpl(12) *
& (10. ** (-1. * gfld%ipdtmpl(11)))
if ( level .lt. pmin ) cycle MATCH_LOOP
elseif((gfld%ipdtmpl(10).eq.105).or.

源网页的做法，是去除掉小于pmin的层数，这样总层数会少于32，而在wrf执行阶段无法运行。主要是num_metgrid_levels默认32与修改后的层数不匹配造成的。
无需修改Vtable之类。

andrewsoong · 发表于 2019-6-17 08:54:46

Another method :update your WPS to v4.0 or later.

巩崇水 · 发表于 2019-6-17 08:39:49

另外，本人还尝试在WPS的namelist里直接修改mod_levs控制输出层数，但是失败了，有没有小伙伴说一下为啥不行？

天下第贰 · 发表于 2019-6-19 18:40:59

请问FV3的数据在哪里下载？https://para.nomads.ncep.noaa.go ... ara/gfs.2019061718/
我在这个网页下载的只有nemsio结尾的数据，这个ungrib不了。。。

巩崇水 · 发表于 2019-6-19 21:44:28

天下第贰发表于 2019-6-19 18:40
请问FV3的数据在哪里下载？https://para.nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/para/gfs.2019061718 ...

ftp://ftpprd.ncep.noaa.gov/pub/data/nccf/com/gfs/prod

天下第贰 · 发表于 2019-6-20 09:43:12

真是太谢谢你了，我昨天找了好久都没找到，请问在哪里可以看到对数据名的描述？有好多文件，不知道都是干什么的，好像FV3 grib2的数据也只有0.25度和0.5度是吗

zhoulm177 · 发表于 2019-6-24 22:12:13

非常感谢楼主，困扰我10来天的问题终于解决了。

Adam-Lee · 发表于 2019-6-27 10:37:23

13号之前的数据都无法下载（自动回复：请不要使用迅雷等下载工具，点我查看下载帮助）了

一颗梅子在南京 · 发表于 2019-7-24 13:24:32

我的WPS是3.5的，根据版主提示修改了rd_grib2.F文件：
! Level (eg. 10000 mb)
            if(gfld%ipdtmpl(10).eq.100) then
               ! Pressure level (range from 1000mb to 0mb)
               level=gfld%ipdtmpl(12) *
   &                         (10. ** (-1. * gfld%ipdtmpl(11)))
if(level.eq.4000)cycle MATCH_LOOP
if(level.eq.1500)cycle MATCH_LOOP
if(level.eq.40)cycle MATCH_LOOP
            elseif(gfld%ipdtmpl(10).eq.105) then
               ! Hybrid level (range from 1 to N)
               level=gfld%ipdtmpl(12)
            elseif(gfld%ipdtmpl(10).eq.104) then
               ! Sigma level (range from 10000 to 0)
               level=gfld%ipdtmpl(12)
            elseif(gfld%ipdtmpl(10).eq.101) then
               ! MSL
               level=201300.
            elseif(gfld%ipdtmpl(10).eq.103) then
         ! Height above ground (m)
   if (gfld%ipdtmpl(12) .eq. 2. .or.
   &             gfld%ipdtmpl(12) .eq. 10. ) then
               level=200100.
   else
               cycle MATCH_LOOP
   endif
然后运行metgrid.exe还是报同样的错误，难道只能把WPS更新到4以上吗？

ZJJJY · 发表于 2019-8-21 17:07:38

一颗梅子在南京发表于 2019-7-24 13:24
我的WPS是3.5的，根据版主提示修改了rd_grib2.F文件：
! Level (eg. 10000 mb)
if(gfld%i ...

请问您这个问题解决了吗？

		自动登录	找回密码
密码			立即注册

FV3核心GFS升级后metgrid错误的问题解决

登录后查看更多精彩内容~

评分