批量读取ctl文件时出错！！！！！

guojiao002

批量读取ctl文件时出错！！！！！
这是ctl. 是在一个ctl基础上修改的。共有8个fnl数据，从fnl0806250000到fnl0806261800。
dset d:/cc/fnl0806%d2%h200
options template
index d:/cc/fnl0806.idx
undef 9.999E+20
title d:/fnl/08/fnl0806010000
*  produced by grib2ctl v0.9.12.5p39c
dtype grib 3
options yrev
ydef 181 linear -90.000000 1
xdef 360 linear 0.000000 1.000000
tdef 8 linear 00Z25jun2008 6hr
*  z has 26 levels, for prs
zdef 26 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
vars 95
no4LFTXsfc  0 132,1,0  ** surface Best (4-layer) lifted index [K]
no5WAVA500mb  0 230,100,500 ** 500 mb 5-wave geopot. height anomaly [gpm]
no5WAVH500mb  0 222,100,500 ** 500 mb 5-wave geopotential height [gpm]
ABSVprs 26 41,100,0 ** (profile) Absolute vorticity
CAPEsfc  0 157,1,0  ** surface Convective Avail. Pot. Energy [J/kg]
CAPE180_0mb  0 157,116,46080 ** 180-0 mb above gnd Convective Avail. Pot. Energy [J/kg]
CINsfc  0 156,1,0  ** surface Convective inhibition [J/kg]
CIN180_0mb  0 156,116,46080 ** 180-0 mb above gnd Convective inhibition [J/kg]
CLWMRprs 21 153,100,0 ** (profile) Cloud water [kg/kg]
CWATclm  0 76,200,0 ** atmos column Cloud water [kg/m^2]
GPAprs 2 27,100,0 ** (profile) Geopotential height anomaly [gpm]
HGTsfc  0 7,1,0  ** surface Geopotential height [gpm]
HGTprs 26 7,100,0 ** (profile) Geopotential height [gpm]
HGTpv2  0 7,117,2000 ** pot vorticity = 2000 units level Geopotential height [gpm]
HGTpvneg2  0 7,117,34768 ** pot vorticity = -2000 units level Geopotential height [gpm]
HGThtfl  0 7,204,0 ** highest trop freezing level Geopotential height [gpm]
HGT0deg  0 7,4,0 ** 0C isotherm level Geopotential height [gpm]
HGTmwl  0 7,6,0 ** max wind level Geopotential height [gpm]
HGTtrp  0 7,7,0 ** tropopause Geopotential height [gpm]
HPBLsfc  0 221,1,0  ** surface Planetary boundary layer height [m]
ICECsfc  0 91,1,0  ** surface Ice concentration (ice=1;no ice=0) [fraction]
LANDsfc  0 81,1,0  ** surface Land cover (land=1;sea=0) [fraction]
LFTXsfc  0 131,1,0  ** surface Surface lifted index [K]
O3MRprs 6 154,100,0 ** (profile) Ozone mixing ratio [kg/kg]
POTsig995   0 13,107,9950 ** sigma=.995  Potential temp. [K]
PRESsfc  0 1,1,0  ** surface Pressure [Pa]
PRESpv2  0 1,117,2000 ** pot vorticity = 2000 units level Pressure [Pa]
PRESpvneg2  0 1,117,34768 ** pot vorticity = -2000 units level Pressure [Pa]
PRESmwl  0 1,6,0 ** max wind level Pressure [Pa]
PREStrp  0 1,7,0 ** tropopause Pressure [Pa]
PRMSLmsl  0 2,102,0  ** unknown level Pressure reduced to MSL [Pa]
PWATclm  0 54,200,0 ** atmos column Precipitable water [kg/m^2]
RHprs 21 52,100,0 ** (profile) Relative humidity [%]
RH2m  0 52,105,2 ** 2 m above ground Relative humidity [%]
RHsig995   0 52,107,9950 ** sigma=.995  Relative humidity [%]
RHsg33_100   0 52,108,8548 ** sigma=0.33-1 layer Relative humidity [%]
RHsg44_72   0 52,108,11336 ** sigma=0.44-0.72 layer Relative humidity [%]
RHsg44_100   0 52,108,11364 ** sigma=0.44-1 layer Relative humidity [%]
RHsg72_94   0 52,108,18526 ** sigma=0.72-0.94 layer Relative humidity [%]
RH30_0mb  0 52,116,7680 ** 30-0 mb above gnd Relative humidity [%]
RHclm  0 52,200,0 ** atmos column Relative humidity [%]
RHhtfl  0 52,204,0 ** highest trop freezing level Relative humidity [%]
RH0deg  0 52,4,0 ** 0C isotherm level Relative humidity [%]
SOILW0_10cm  0 144,112,10 ** 0-10 cm underground Volumetric soil moisture [fraction]
SOILW10_40cm  0 144,112,2600 ** 10-40 cm underground Volumetric soil moisture [fraction]
SOILW40_100cm  0 144,112,10340 ** 40-100 cm underground Volumetric soil moisture [fraction]
SOILW100_200cm  0 144,112,25800 ** 100-200 cm underground Volumetric soil moisture [fraction]
SPFH2m  0 51,105,2 ** 2 m above ground Specific humidity [kg/kg]
SPFH30_0mb  0 51,116,7680 ** 30-0 mb above gnd Specific humidity [kg/kg]
TCDCcvl  0 71,244,0 ** convective cld layer Total cloud cover [%]
TMPsfc  0 11,1,0  ** surface Temp. [K]
TMPprs 26 11,100,0 ** (profile) Temp. [K]
TMP3658m  0 11,103,3658 ** 3658 m above msl Temp. [K]
TMP2743m  0 11,103,2743 ** 2743 m above msl Temp. [K]
TMP1829m  0 11,103,1829 ** 1829 m above msl Temp. [K]
TMP2m  0 11,105,2 ** 2 m above ground Temp. [K]
TMPsig995   0 11,107,9950 ** sigma=.995  Temp. [K]
TMP0_10cm  0 11,112,10 ** 0-10 cm underground Temp. [K]
TMP10_40cm  0 11,112,2600 ** 10-40 cm underground Temp. [K]
TMP40_100cm  0 11,112,10340 ** 40-100 cm underground Temp. [K]
TMP100_200cm  0 11,112,25800 ** 100-200 cm underground Temp. [K]
TMP30_0mb  0 11,116,7680 ** 30-0 mb above gnd Temp. [K]
TMPpv2  0 11,117,2000 ** pot vorticity = 2000 units level Temp. [K]
TMPpvneg2  0 11,117,34768 ** pot vorticity = -2000 units level Temp. [K]
TMPmwl  0 11,6,0 ** max wind level Temp. [K]
TMPtrp  0 11,7,0 ** tropopause Temp. [K]
TOZNEclm  0 10,200,0 ** atmos column Total ozone [Dobson]
UGRDprs 26 33,100,0 ** (profile) u wind [m/s]
UGRD3658m  0 33,103,3658 ** 3658 m above msl u wind [m/s]
UGRD2743m  0 33,103,2743 ** 2743 m above msl u wind [m/s]
UGRD1829m  0 33,103,1829 ** 1829 m above msl u wind [m/s]
UGRD10m  0 33,105,10 ** 10 m above ground u wind [m/s]
UGRDsig995   0 33,107,9950 ** sigma=.995  u wind [m/s]
UGRD30_0mb  0 33,116,7680 ** 30-0 mb above gnd u wind [m/s]
UGRDpv2  0 33,117,2000 ** pot vorticity = 2000 units level u wind [m/s]
UGRDpvneg2  0 33,117,34768 ** pot vorticity = -2000 units level u wind [m/s]
UGRDmwl  0 33,6,0 ** max wind level u wind [m/s]
UGRDtrp  0 33,7,0 ** tropopause u wind [m/s]
VGRDprs 26 34,100,0 ** (profile) v wind [m/s]
VGRD3658m  0 34,103,3658 ** 3658 m above msl v wind [m/s]
VGRD2743m  0 34,103,2743 ** 2743 m above msl v wind [m/s]
VGRD1829m  0 34,103,1829 ** 1829 m above msl v wind [m/s]
VGRD10m  0 34,105,10 ** 10 m above ground v wind [m/s]
VGRDsig995   0 34,107,9950 ** sigma=.995  v wind [m/s]
VGRD30_0mb  0 34,116,7680 ** 30-0 mb above gnd v wind [m/s]
VGRDpv2  0 34,117,2000 ** pot vorticity = 2000 units level v wind [m/s]
VGRDpvneg2  0 34,117,34768 ** pot vorticity = -2000 units level v wind [m/s]
VGRDmwl  0 34,6,0 ** max wind level v wind [m/s]
VGRDtrp  0 34,7,0 ** tropopause v wind [m/s]
VVELprs 21 39,100,0 ** (profile) Pressure vertical velocity [Pa/s]
VVELsig995   0 39,107,9950 ** sigma=.995  Pressure vertical velocity [Pa/s]
VWSHpv2  0 136,117,2000 ** pot vorticity = 2000 units level Vertical speed shear [1/s]
VWSHpvneg2  0 136,117,34768 ** pot vorticity = -2000 units level Vertical speed shear [1/s]
VWSHtrp  0 136,7,0 ** tropopause Vertical speed shear [1/s]
WEASDsfc  0 65,1,0  ** surface Accum. snow [kg/m^2]
ENDVARS

已经生成了idx文件（总文件），用grads打开上述的ctl后，
set t 6   
set lev 300
d vvelprs
出现：cannot contour grid ----all undefined values
数据没有错误，用单个的ctl 都可打开画图。
求救啊！！！！

guojiao002

@mofangbao
求救啊！

mofangbao

我没有做过批量描述fnl资料，手机上的，没法去查手册，只能猜一下，比如你用了两个options，在写idx文件的时候需要批量描述吗，也就是一个idx文件够用吗？你在试试

天目神眉

呃，没这样批量处理过，纯帮顶

topmad

感觉好复杂的说，帮顶一下吧，等大牛来解惑

@天目神眉 帮顶也有卡啊，木有天理啊

天目神眉

topmad 发表于 2011-10-31 13:25

                               
登录/注册后可看大图

感觉好复杂的说，帮顶一下吧，等大牛来解惑

@天目神眉 帮顶也有卡啊，木有天理啊

人品好，再来一张就去买彩票

善人/jw

没处理过fnl
@传说中的谁来看看吧

guojiao002

mofangbao 发表于 2011-10-31 10:41

                               
登录/注册后可看大图

我没有做过批量描述fnl资料，手机上的，没法去查手册，只能猜一下，比如你用了两个options，在写idx文件的时 ...

谢谢啊。我再试试。

guojiao002

感谢各位管理员、版主捧场。
问题已解决。

mofangbao

guojiao002 发表于 2011-10-31 19:03

                               
登录/注册后可看大图

感谢各位管理员、版主捧场。
问题已解决。

解决的问题请附上方法吧，要不然不觉得大家白讨论了么？谢谢了