|
本帖最后由 冬腌菜 于 2014-5-19 14:31 编辑
gfs 0.5 deg +00hrs 是这样处理的
http://nomads.ncep.noaa.gov/pub/ ... 8/gfs.t18z.pgrb2f00
保存为 1118.bin
- for d in 11 12 13 14 15 ; do for h in 00 06 12 18 ; do g2ctl $d$h.bin >$d$h.ctl ; gribmap -i $d$h.ctl ; done ; done
复制代码
然后修改 gfs.ctl 如下:
- dset ^%d2%h2.bin
- options template
- index ^gfs.bin.idx
- undef 9.999E+20
- title 1100.bin
- * produced by g2ctl v0.0.9
- * command line options: 1100.bin
- * griddef=1.1:0:(720 x 361):grid_template=0:winds(N/S): lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48 lat 90.000000 to -90.000000 by 0.500000 lon 0.000000 to 359.500000 by 0.500000 #points=259920:winds(N/S)
- dtype grib2
- ydef 361 linear -90.000000 0.5
- xdef 720 linear 0.000000 0.500000
- tdef 20 linear 00Z11may2014 6hr
- * PROFILE hPa
- zdef 26 levels 100000 97500 95000 92500 90000 85000 80000 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000 15000 10000 7000 5000 3000 2000 1000
- options pascals
- vars 121
- ...
- ENDVARS
复制代码
然后是 gribmap -i gfs.ctl
grads脚本:
- open gfs.ctl
- set t 1 20
- set lev 850
- d ave(aave(UGRDprs, lon=110, lon=120, lat=5, lat=10), t=1, t=20)
复制代码
这么做应该没问题吧?
ps: http://www.sac345.org.cn/upload/userfiles/NHJF.pdf 这个现在还没起用?
更新:JMA第27侯(5.11-15)报告已经完成
10-20N应该是高于235W/m2的,不过CMA似乎不重视这个值
 |
|