简单介绍如何使用gfs数据跑wrf做预报（仅在下陋识，望高人指点）

andrewsoong

西西 发表于 2015-1-27 21:40
你自己不会修改下么。。。得了。我给你弄个。。。

我的意思是那个网站没有发布0.25的。。。。。。

tornadoarps

全国这么多气象业务单位，跑wrf模式的都用我以前帖子里提到的0.5度分辨率gfs数据来启动的。wrf、arps模式直接支持这种数据。

西西

tornadoarps 发表于 2015-1-30 09:36
全国这么多气象业务单位，跑wrf模式的都用我以前帖子里提到的0.5度分辨率gfs数据来启动的。wrf、arps模式直 ...

额。。。好像不都是0.5吧。0.5的数据太大。有的也用1的吧。

tornadoarps

西西 发表于 2015-1-30 14:48
额。。。好像不都是0.5吧。0.5的数据太大。有的也用1的吧。

0.5度的一个文件才50M，不算大吧？再说0.5度比1度分辨率高一倍，为什么不用呢？

西西

tornadoarps 发表于 2015-1-30 17:53
0.5度的一个文件才50M，不算大吧？再说0.5度比1度分辨率高一倍，为什么不用呢？

额。。。你说的是单个的数据吧。完整的一个大概是4g多。分辨率高，做出来的效果不一定就好的。有些时候也不一定需要那么高的精度。

tornadoarps

西西 发表于 2015-2-1 13:03
额。。。你说的是单个的数据吧。完整的一个大概是4g多。分辨率高，做出来的效果不一定就好的。有些时候也 ...

    高分辨率启动wrf模式效果不一定好？这个说法我还是第一次听说。低分辨率资料不能完整包含中小尺度天气系统，用它来启动wrf，预报效果如何好？科学的东西不是随便说说的。这个网站还是有很多人看的，不确切的话还是不说的好，以免误导初学者。
    0.5度的分辨率不算高，现在国家局下发的ecmwf地面层分辨率达到了0.125度。气压层分辨率0.25度。如果你想转wrf，初始场不用高分辨率资料的话，转wrf根本就没什么意义，因为wrf根本达不到ec的预报效果。我在省级气象台工作，我们单位的wrf每天运行，做96小时预报，但输出产品很少有人看，预报员最信赖的还是ec。wrf最佳的应用领域是初始场加入观测资料，用来做0～24小时的中尺度预报。
    数值预报的预报效果和初始场质量有直接关系，初始场越接近真实大气，预报效果就越好。这就是目前资料同化技术如火如荼的原因。资料同化的目的就是用三维变分、四维变分、集合卡尔曼滤波等技术，将雷达、卫星、探空、风廓线、地面自动站等资料加入初始场，改善初始场质量，使其更加接近真实大气，这样才能提高模式的预报效果。
    0.5度gfs资料分辨率还不够高，我下一步的工作就是用国家局下发的0.25分辨率ec做初始场和边界条件来启动wrf，并实现业务化。

schliezer

西西 发表于 2015-2-1 13:03
额。。。你说的是单个的数据吧。完整的一个大概是4g多。分辨率高，做出来的效果不一定就好的。有些时候也 ...

你好，我也试着用过gfs.t00z.pgrb2.XpXX.fXXX系列数据中的gfs.t00z.pgrb2.1p50.anl，但metgrid总不过，提示：In entry 65 of METGRID.TBL, the wind U-component field must be interpolated to the V staggered grid points。我也不知道这个XpXX（比如2p50,1p00,0p25,0p50）是什么意思。
因此想问下你譬如我要以某天00时开始，用wrf做未来三天预报，究竟应该下GFS预报场网站上哪个系列的哪些数据以及两个namelist（WPS,WRF）里的起止时间又应该如何设定，是否就像ARW里的一样？谢谢你!

西西

schliezer 发表于 2015-2-2 22:54
你好，我也试着用过gfs.t00z.pgrb2.XpXX.fXXX系列数据中的gfs.t00z.pgrb2.1p50.anl，但metgrid总不过，提 ...

2p50 应该是指2.5*2.5°的，1p00应该是1*1的。依次类推吧。因为我原来用gfs的时候，数据名称格式不是这样的。ftp://ftp.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/在这个网址下，你可以下载。根据官网上给出的数据形式，我估计
gfs.t00z.pgrb2.1p00.f00应该是第一个预报数据。3小时为间隔。你要做3天预报，那么至少要72小时。那你可以依次下载到f72甚至更多。具体的还是由你自己多尝试。时间设置和ARW是一样的。你预报哪几天就设置为哪几天。

schliezer

西西 发表于 2015-2-2 23:18
2p50 应该是指2.5*2.5°的，1p00应该是1*1的。依次类推吧。因为我原来用gfs的时候，数据名称格式不是这样 ...

你好，我下载了从f00到f72的数据，geogrid和ungrib都顺利通过，但metgrid总是过不去（我用的是NMM内核），错误显示ERROR: In entry 65 of METGRID.TBL, the wind U-component field must be interpolated to the V staggered grid points.这我在以前用ARW时从没见过的。
METGRID.TBL的第65个词条是：
name=UU
        mandatory=yes    # MUST HAVE THIS FIELD
        interp_option=sixteen_pt+four_pt+average_4pt
        is_u_field=yes
        output_stagger=U
        fill_missing=0.
        fill_lev=200100:const(-1.E30)
但它说的U向风场须内插到V向交叉格点是啥意思呢？

一直胖乎乎

好东西，感谢分享！