T639预报数据，如何用grads批处理

兰北

    如题，想用T639的预报数据，数据格式为“年月日+时+3位数的预报时效”，如2010050100003.grb1，003表示3h的预报，然后是006,009,012等。我用grid2ctl生成ctl文件后，想进行批处理，但不知描述文件tdef该怎么写了。
    如果写成 tdef 24 linear 00Z01may2010 3hr，就自动寻找2010050103了，提示找不到文件。
    求大神指导！感谢！

兰北

如果写成分，也不太合适，因为预报时效可以达到240小时，附一张文件名的图

edwardli

兰北 发表于 2013-7-28 22:18

                               
登录/注册后可看大图

如果写成分，也不太合适，因为预报时效可以达到240小时，附一张文件名的图

你这个是T639的高分辨率产品。
先用一个，使用grib2ctl，产生一个时效文件的ctl文件，好好琢磨；
然后，copy /b *.grb1 2010050100.grib1；
接着，用你自己的办法，生成这个时间起报的ctl文件，也就是将第一步的获得的ctl改写一下，比如tdef，其他的基本不变；
继续，用gribmap -verf 2010050100.ctl 生成grib1文件的idx文件
最后，你就能用了，open 2020050100.ctl。。。。。。。。
。gs？手动控制？
。up to you
。

mofangbao

tdef跟你说的什么GrADS找什么是个什么意细思，tdef只是表示时间的信息，跟你文件名没有关系的，是跟你文件里面的时间有关

兰北

mofangbao 发表于 2013-7-29 08:20

                               
登录/注册后可看大图

tdef跟你说的什么GrADS找什么是个什么意细思，tdef只是表示时间的信息，跟你文件名没有关系的，是跟你文件里 ...

谢谢！使用grib2ctl，产生一个预报时效003的文件的ctl文件中，tdef是 1 linear 00Z01may2010 1mo，预报时效006的ctl中，tdef描述也是1 linear  00Z01may2010 1mo
但我还是不太明白，就是，我想进行批处理，如果写tdef 1，就没有时间的循环啊，默认的是gmf.639.2010050100000.grb1

该怎么改？
ctl文件如下：
dset F:/201005-09/T639GSI2GRIB_HR_201005/T639GSI2GRIB_HR_2010050100/gmf.639.2010050100%f3.grb1
options template
index F:/201005-09/T639GSI2GRIB_HR_201005/T639GSI2GRIB_HR_2010050100/gmf.639.2010050100.grb1.idx
undef 9.999E+20
title gmf.639.2010050100003.grb1
*  produced by grib2ctl v0.9.12.5p39c
dtype grib 255
options yrev
ydef 181 linear -89.999000 1
xdef 360 linear 0.000000 1.000000

**问题在这里，就是时间的描述方法？预报资料间隔是3h，003,006,009,012等，tdef该怎么写
***tdef 24 linear 00Z01may2010 1mo
tdef 1 linear 00Z01may2010 1mo

*  z has 32 levels, for prs
zdef 32 levels
1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10 7 5 4 3 2 1
vars 77
ACPCPsfc  0 63,1,0  ** surface Convective precipitation [kg/m^2]
APCPsfc  0 61,1,0  ** surface Total precipitation [kg/m^2]
EVPsfc  0 57,1,0  ** surface Evaporation [kg/m^2]
GFLUXsfc  0 155,1,0  ** surface Ground heat flux [W/m^2]
GRADsfc  0 117,1,0  ** surface Global radiation [W/m^2]
HCDCclb  0 75,2,0 ** cloud base High level cloud cover [%]
HGTclm33025  0 7,10,33025 **level=33025 Geopotential height [gpm]
……
endvars

兰北

edwardli 发表于 2013-7-29 07:36

                               
登录/注册后可看大图

你这个是T639的高分辨率产品。
先用一个，使用grib2ctl，产生一个时效文件的ctl文件，好好琢磨；
然后， ...

嗯，是的，我之前做过ncep资料中6h间隔的再分析资料的循环。现在就是针对预报资料，发现不会描述tdef了，ctl文件在上面，请大力帮助指导下！谢谢！

wuwei2163

兰北 发表于 2013-7-29 14:35

                               
登录/注册后可看大图

嗯，是的，我之前做过ncep资料中6h间隔的再分析资料的循环。现在就是针对预报资料，发现不会描述tdef了， ...

楼主 我想问你做fnl在分析资料 进行批处理的时候 出现过这样的情况吗  如果那样可以的话 你可以把639的文件名改成在分析资料的名字格式啊

兰北

wuwei2163 发表于 2013-7-29 23:20

                               
登录/注册后可看大图

楼主 我想问你做fnl在分析资料 进行批处理的时候 出现过这样的情况吗  如果那样可以的话 你可以把639的文 ...

这倒可以试一下啊，谢谢

wuwei2163

兰北 发表于 2013-7-30 10:43

                               
登录/注册后可看大图

这倒可以试一下啊，谢谢

不客气 我最近也在做这方面的研究  分别作fnl和t639资料批处理 进展不太顺利 问一下 你的fnl批处理能够出来吗？

兰北

wuwei2163 发表于 2013-7-30 11:10

                               
登录/注册后可看大图

不客气 我最近也在做这方面的研究  分别作fnl和t639资料批处理 进展不太顺利 问一下 你的fnl批处理能够出 ...

fnl没问题的