求救 关于fnl资料批处理的

wuwei2163

在论坛里学习了关于关于fnl资料批处理的方法，对fnl资料批处理：先生成初始时刻的的ctl文件，然后修改ctl文件后经过gribmap处理，但是出来的每一个时刻的图都是一样的，而且都是第一个时刻的图，什么原因啊？
昨天时间有点仓促，没有贴上，ctl修改后，也用gribmap处理了；部分ctl如下：
dset c:/fnl_201303%d2_%h2_00
options template
index c:/fnl_20130308_12_00.idx
undef 9.999E+20
title c:/fnl_20130308_12_00
* produced by g2ctl v0.0.8.2
* command line options: c:/fnl_20130308_12_00
dtype grib2
ydef 181 linear -90.000000 1
xdef 360 linear 0.000000 1.000000
tdef 8 linear 12Z08mar2013 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 138
gs如下：
'open c:/fnl/fnl_201303%d2_%h2_00'
'set mpdset cnworld'
'set lat 30 40'
'set lon 110 120'
'set lev 850'
'define t850=TMPprs-273'
'define td850=(TMPprs-273)-((14.55+0.114*(TMPprs-273))*(1-0.01*RHprs) + pow((2.5+0.007*(TMPprs-273))*(1-0.01*RHprs),3) + (15.9+0.37*(TMPprs-273))*pow((1-0.01*RHprs),14))'
'set lev 700'
'define t700=TMPprs-273'
'define td700=(TMPprs-273)-((14.55+0.114*(TMPprs-273))*(1-0.01*RHprs) + pow((2.5+0.007*(TMPprs-273))*(1-0.01*RHprs),3) + (15.9+0.37*(TMPprs-273))*pow((1-0.01*RHprs),14))'
'set lev 500'
'define t500=TMPprs-273'
'set grid off'
'set grads off'
'set gxout contour'
'set lat 30 40'
'set lon 110 120'
'set t 1'
'd t850-t500+td850-(t700-td700)'
'enable print c:/55.gmf '
'print'
'disable print'
;把gs文件中的时间改之后，出的图都是一样的

liutaoz21

ctl和gs都贴出来吧，不然空口无凭。我做过几次fnl批处理，就是用的兰溪整合的方法，都没问题的。

尽头的尽头

你不贴CTL和GS，我只能说原因就是你做错了

空谷幽竹

错误原因有多种，需要仔细分析啊

wuwei2163

liutaoz21 发表于 2013-6-22 03:54

                               
登录/注册后可看大图

ctl和gs都贴出来吧，不然空口无凭。我做过几次fnl批处理，就是用的兰溪整合的方法，都没问题的。

昨天时间有点仓促，没有贴上，ctl修改后，也用gribmap处理了；部分ctl如下：
dset c:/fnl_201303%d2_%h2_00
options template
index c:/fnl_20130308_12_00.idx
undef 9.999E+20
title c:/fnl_20130308_12_00
* produced by g2ctl v0.0.8.2
* command line options: c:/fnl_20130308_12_00
dtype grib2
ydef 181 linear -90.000000 1
xdef 360 linear 0.000000 1.000000
tdef 8 linear 12Z08mar2013 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 138
gs如下：
'open c:/fnl/fnl_201303%d2_%h2_00'
'set mpdset cnworld'
'set lat 30 40'
'set lon 110 120'
'set lev 850'
'define t850=TMPprs-273'
'define td850=(TMPprs-273)-((14.55+0.114*(TMPprs-273))*(1-0.01*RHprs) + pow((2.5+0.007*(TMPprs-273))*(1-0.01*RHprs),3) + (15.9+0.37*(TMPprs-273))*pow((1-0.01*RHprs),14))'
'set lev 700'
'define t700=TMPprs-273'
'define td700=(TMPprs-273)-((14.55+0.114*(TMPprs-273))*(1-0.01*RHprs) + pow((2.5+0.007*(TMPprs-273))*(1-0.01*RHprs),3) + (15.9+0.37*(TMPprs-273))*pow((1-0.01*RHprs),14))'
'set lev 500'
'define t500=TMPprs-273'
'set grid off'
'set grads off'
'set gxout contour'
'set lat 30 40'
'set lon 110 120'
'set t 1'
'd t850-t500+td850-(t700-td700)'
'enable print c:/55.gmf '
'print'
'disable print'
;把gs文件中的时间改之后，出的图都是一样的

wuwei2163

尽头的尽头 发表于 2013-6-22 08:13

                               
登录/注册后可看大图

你不贴CTL和GS，我只能说原因就是你做错了

昨天时间有点仓促，没有贴上，ctl修改后，也用gribmap处理了；部分ctl如下：
dset c:/fnl_201303%d2_%h2_00
options template
index c:/fnl_20130308_12_00.idx
undef 9.999E+20
title c:/fnl_20130308_12_00
* produced by g2ctl v0.0.8.2
* command line options: c:/fnl_20130308_12_00
dtype grib2
ydef 181 linear -90.000000 1
xdef 360 linear 0.000000 1.000000
tdef 8 linear 12Z08mar2013 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 138
gs如下：
'open c:/fnl/fnl_201303%d2_%h2_00'
'set mpdset cnworld'
'set lat 30 40'
'set lon 110 120'
'set lev 850'
'define t850=TMPprs-273'
'define td850=(TMPprs-273)-((14.55+0.114*(TMPprs-273))*(1-0.01*RHprs) + pow((2.5+0.007*(TMPprs-273))*(1-0.01*RHprs),3) + (15.9+0.37*(TMPprs-273))*pow((1-0.01*RHprs),14))'
'set lev 700'
'define t700=TMPprs-273'
'define td700=(TMPprs-273)-((14.55+0.114*(TMPprs-273))*(1-0.01*RHprs) + pow((2.5+0.007*(TMPprs-273))*(1-0.01*RHprs),3) + (15.9+0.37*(TMPprs-273))*pow((1-0.01*RHprs),14))'
'set lev 500'
'define t500=TMPprs-273'
'set grid off'
'set grads off'
'set gxout contour'
'set lat 30 40'
'set lon 110 120'
'set t 1'
'd t850-t500+td850-(t700-td700)'
'enable print c:/55.gmf '
'print'
'disable print'
;把gs文件中的时间改之后，出的图都是一样的

尽头的尽头

wuwei2163 发表于 2013-6-22 09:59

                               
登录/注册后可看大图

昨天时间有点仓促，没有贴上，ctl修改后，也用gribmap处理了；部分ctl如下：
dset c:/fnl_201303%d2_%h2 ...

你的高度场是个什么情况

liutaoz21

wuwei2163 发表于 2013-6-22 09:59

                               
登录/注册后可看大图

昨天时间有点仓促，没有贴上，ctl修改后，也用gribmap处理了；部分ctl如下：
dset c:/fnl_201303%d2_%h2 ...

要不要试试把set t 1放到set level之前？

blnuist

gs文件里open的不应该是ctl文件吗？

林宝宝

楼上正解。。。。。。。