GrADS可以判断图片是否存在吗？

littleboy

GrADS可以判断图片是否存在吗？
GrADS可以对缺失数据的那个时段不画图吗？

缺失一些数据

有数据的

缺失数据的

river

没有数据怎么画图？

littleboy

有数据，就是有时候数据会缺失，想将缺失的数据的不画出图片
ctl文件
*设置数据文件
dset E:/web/gzmb/data/gzcompass/ecmwfthin/ov/20160429000000/ecmwthin_%y4%m2%d2%h2.dat
options template
*定义缺测值或缺值值
undef -999.9
*TITLE string 用字符串string 简略描述数据集内容
title ecmwthindata
*y轴方向上以17.68开始，0.12为间隔的78个对线性映射的网格点
ydef 201 linear 10 0.25
*x轴方向上以106.88开始，0.12为间隔的103个对线性映射的网格点
xdef 361 linear 60.0 0.25
*设置网格值与时间的映射关系，表示共有37 个时次，
*起始时刻2016 年3月8 号12Z 时，增量为1 小时。
tdef 20 linear 00Z29APR2016 3hr
*表示共7 层等压面，其值各为1000，925，850，700,500,200,100等等。
zdef 7 levels  1000 925 850 700 500 200 100
vars 3
hght 7 99 hght
uwnd 7 99 uwnd
vwnd 7 99 vwnd
ENDVARS

ecmwthin.gs文件
'reinit'

'open E:/server/OpenGrADSData/Example/ecwmthin/ecmwthin.ctl'

t=1
while(t<=20)
'set t 't
'set map 1 1 1'
'set font 0'
'set mpdset hires'
'set lev 200'
setcols()
'set xlopts 1 1 0.12'
'set ylopts 1 1 0.12'
'set vpage 0 11 0 8.5'
'set parea 0.5 9.4 0 8.5'
'set clopts 99 1 0.09'
'set clab on'
'set grads off'
'set grid off'
'set ccols 20'
'draw rec 9.8 2.24 10 6.26'
'draw string 9.725 6.4 | V |'
'set clevs 30 40 50 60 70 80 90'
'set ccols 20 21 22 23 24 25 26 27'
'set gxout shaded'
'd mag(uwnd,vwnd)'
'run E:/server/opengrads/Contents/Resources/Scripts/cbar.gs 1 1 '

'set font 0'
'set gxout stream'
'set strmden -5'
'd uwnd;vwnd'
*'set gxout barb'
*'d skip(uwnd*2.5,6);skip(vwnd*2.5,6)'

'set gxout contour'
'set cint 10'
'd hght'
*'set clevs 1244'
*'set cthick 6'
*'d hght'
'set strsiz 0.14'
'set font 5'
'set string 1 tl 3'
*'draw title testsss'
'draw string 0.5 7.8 Wind Speed(m/s) & Wind Stream & Height at 200hPa'
'set string 1 bl 3'
'draw string 0.5 0.5 Ecmwtthin 0.25*0.25'
'draw string 6.8 0.5 '2016030812' +'leadtime(t-1)'h'
*'draw string 8.3 8.0 Height at 200hPa'
'set string 1 tr 2'
'set font 1'
'set strsiz 0.10 0.10'
'draw string 9.5 7.8 'qtime()'(UTC)'
'gxprint E:/server/OpenGrADSData/PicPath/ecmwthin/200pha/div200_2016030812_'leadtime(t-1)'.png x1000 y700 white'
'c'
t=t+1
endwhile

'quit'

function setcols()
'set rgb 20 255 255 255'
'set rgb 21 181 201 255'
'set rgb 22 127 150 255'
'set rgb 23 0 99 255'
'set rgb 24 0 198 51'
'set rgb 25 150 255 0'
'set rgb 26 255 255 0'
'set rgb 27 255 160 0'
'set rgb 28 246 126 49'
'set rgb 29 235 84 41'
'set rgb 30 220 51 40'
'set rgb 31 199 29 37'
'set rgb 99 0 0 0'
'set rgb 100 255 0 0'
'set rbcols 31 30 29 28 27 26 25 24 23 22 21 20'
return

function leadtime(t)
str='00'
if(t<10)
time=substr(str,1,2)''t
else
if(t<100)
time=substr(str,1,1)''t
else
time=t
endif
endif
return time

function qtime()
monthStr='JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC'
monthNum='01 02 03 04 05 06 07 08 09 10 11 12'
'q time'
yy=substr(result,16,4)
MM=substr(result,13,3)
dd=substr(result,11,2)
hh=substr(result,8,2)
loop=1
while(loop<=12)
if(MM=subwrd(monthStr,loop))
MM=subwrd(monthNum,loop)
break
endif
loop=loop+1
endwhile
time=yy''MM''dd''hh
return time