查看: 3370|回复: 1

求大神帮忙！！关于使用FNL批处理，GS文件不出图的问题

发表于 2014-4-29 21:21:51 | 显示全部楼层 |阅读模式

GrADS
系统平台:	win7
问题截图:
问题概况:	数据资料为FNL资料，16个时次，26个层次。 CTL文件如下： dset D:/FNL/fnl_201309%d2_%h2_00_c options template index D:/FNL/fnl_201309_c.idx undef 9.999E+20 title D:/FNL/fnl_201309_c * 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 16 linear 00Z21sep2013 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 20 no5WAVA500mb 0 230,100,500 500 mb 5-wave geopot. height anomaly [gpm]500百帕等压面位势高度距平 no5WAVH500mb 0 222,100,500 500 mb 5-wave geopotential height [gpm]500百帕位势高度 ABSVprs 26 41,100,0 (profile) Absolute vorticity [/s]绝对涡度 APCPsfc 0 61,1,0 surface Total precipitation [kg/m^2]地面总降水量 GPAprs 2 27,100,0 (profile) Geopotential height anomaly [gpm]位势高度距平 HGTprs 26 7,100,0 (profile) Geopotential height [gpm]位势高度 LHTFLsfc 0 121,1,0 surface Latent heat flux [W/m^2]地表潜热通量 POTsig995 0 13,107,9950 sigma=.995 Potential temp. [K]位温 PRATEsfc 0 59,1,0 surface Precipitation rate [kg/m^2/s]地表降水率 PRESsfc 0 1,1,0 surface Pressure [Pa]地表气压 RHprs 21 52,100,0 (profile) Relative humidity [%]相对湿度 RH2m 0 52,105,2 2 m above ground Relative humidity [%]地面2米以上相对湿度 SHTFLsfc 0 122,1,0 surface Sensible heat flux [W/m^2]地表显热通量 SPFH2m 0 51,105,2 2 m above ground Specific humidity [kg/kg]地面2米以上比湿 TMPsfc 0 11,1,0 surface Temp. [K]地表温度 TMPprs 26 11,100,0 (profile) Temp. [K]温度 UGRDprs 26 33,100,0 (profile) u wind [m/s]U风速 VGRDprs 26 34,100,0 (profile) v wind [m/s]V风速 VVELprs 21 39,100,0 (profile) Pressure vertical velocity [Pa/s]垂直速度 VWSHtrp 0 136,7,0 tropopause Vertical speed shear [1/s]垂直风切变 ENDVARS GS文件如下： 'reinit' 'open D:/FNL/fnl_2013.ctl' tt=1 while(tt<=16) 'set t'tt zz=1 while(zz<=26) 'set z'zz 'set lat 0 60' 'set lon 60 150' 'enable print D:/FNL/picture/HGT/hgt.gmf' 'set grads off' 'set grid off' 'set cterp on' 'set csmooth on' 'set mpdraw on' 'set poli on' 'set gxout contour' 'd HGTprs' 'q dims' it=sublin(result,5) it=subwrd(it,6) iz=sublin(result,4) iz=subwrd(iz,6) 'draw title 'it' HGTprs for 'iz'hPa' 'print' 'c' zz=zz+1 endwhile tt=tt+1 endwhile 'disable print' ;
我看过提问的智慧:	看过
自己思考时长(天):	1