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

巩崇水

问题分析
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

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

巩崇水

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

天下第贰

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

巩崇水

天下第贰 发表于 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

天下第贰

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

zhoulm177

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

Adam-Lee

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

一颗梅子在南京

我的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-7-24 13:24
我的WPS是3.5的，根据版主提示修改了rd_grib2.F文件：
! Level (eg. 10000 mb)
              if(gfld%i ...

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