集合预报数据如何写ctl文件

fzmmarco

我现在手头有一个集合预报得到的数据，假如集合成员有14个，想用grads绘图，不知道怎么在ctl文件中体现集合成员数？
不知道版上各位有没有遇到写集合预报数据ctl文件的情况？还望不吝赐教

兰溪之水

chenshengjie 发表于 2016-3-3 21:40
兰溪版主你好，你提供的http://www.iges.org/grads/gadoc/ensembles.html这个网址打开总是域名解析，我最 ...

grads的网址换了，http://cola.gmu.edu/grads/gadoc/ensembles.html

兰溪之水

http://www.iges.org/grads/gadoc/ensembles.html

fzmmarco

兰溪之水 发表于 2015-3-21 13:00
http://www.iges.org/grads/gadoc/ensembles.html

我仔细看完了，收货很大，非常感谢！！

寒枫沁雪

好好
好好

chenshengjie

兰溪之水 发表于 2015-3-21 13:00
http://www.iges.org/grads/gadoc/ensembles.html

兰溪版主你好，你提供的http://www.iges.org/grads/gadoc/ensembles.html这个网址打开总是域名解析，我最近着手做集合预报，有不少要学习的。目前下载好2011-2015年夏季3个月的EC50个成员的数据，希望用grads处理GRIB，但是G2CTL总是出错，正在论坛里搜帖子学习。请赐教，谢谢。

baggio

{:lxm_2:}{:lxm_2:}{:lxm_2:}{:lxm_2:}{:lxm_2:}{:lxm_2:}

chenshengjie

兰溪之水 发表于 2016-3-3 23:34
grads的网址换了，http://cola.gmu.edu/grads/gadoc/ensembles.html

你好，我在http://apps.ecmwf.int/datasets/data/tigge/levtype=sfc/type=cf/下载了2011年6月08时起报，时效为12小时的降水敏感性预报，应该是4维数据吧，X,Y,z,T，也用g2ctl生成了ctl,dset h:/jihedata1/c6712.grib
index h:/jihedata1/c6712.grib.idx
undef 9.999E+20
title h:/jihedata1/c6712.grib
*  produced by g2ctl v0.0.3d
* griddef=1:0:(321 x 141):grid_template=0: lat-lon grid:(321 x 141) units 1e-006 input WE:NS output WE:SN res 48 lat 70.000000 to 0.000000 by 0.500000 lon 0.000000 to 160.000000 by 0.500000 #points=45261

dtype grib2
ydef 141 linear 0.000000 0.5
xdef 321 linear 0.000000 0.500000
tdef 61 linear 12Z01jun2011 1dy
zdef 1 linear 1 1
vars 1
var0152sfc 0,1   0,1,52,1 ** surface acc desc [unit]
ENDVARS

但不懂为何画图，出现cannot grid-all undefined values 出不来图！！！！相同的补充用在forcast上就可以出图。还没有加入多个成员。求助

chenshengjie

兰溪之水 发表于 2016-3-3 23:34
grads的网址换了，http://cola.gmu.edu/grads/gadoc/ensembles.html

相同的步骤处理type为forecast的（http://apps.ecmwf.int/datasets/data/tigge/levtype=sfc/type=fc/）数据就可以正常画图。是g2ctl版本旧？还是怎么回事？下载的数据中还没有加入多个成员。同样的Perturbed forecast数据，只下载了一个成员，照理也是4维数据，同样出现cannot contour的错误。求帮助！谢谢