【求助】FY2E PRE的处理

凌寒

ctl：
dset REPLACEDSET
title REPLACETITLE
options yrev
undef -999
xdef 1001 linear 55 0.1
ydef 501 linear 0 0.1
zdef 1 linear 300 1
tdef 1 linear 23z26jun2013 1hr
vars 1
pre 0 -1,40,1 PRE(MM/HOUR)
endvars
说明：之所以写REPLACEDSET，是因为要实时处理，处理时，REPLACEDSET会被替换成其他的内容

gs：
'reinit'
'open /home/cfg/FY2E_PRE.ctl'
'set gxout shaded'
'set grads off'
'set mpdset cnworld'
'set map 1'
'run mycolor.gs'
'set clevs 0 5 10 20 50 100 200 300'
'set ccols 0 21 22 23 24 25 26 27 28'
'd pre'
'cbarn 0.9 0'
'q file'
t=sublin(result,3)
FulFile=subwrd(t,2)
length=strlen(FulFile)
PostFileName=substr(FulFile,35,length)
length2=strlen(PostFileName)
PreFileName=substr(PostFileName,1,length2-5)
'q ctlinfo'
TitleString=sublin(result,2)
Title=subwrd(TitleString,2)
'draw title 'Title''
'printim /home/pic/'PreFileName'.png white'
quit

绘图：

把这个图放大看会发现格点值不连续，但是查数据的话没有看到有不连续的地方。
即在一片区域中会有取值为0的格点
求解？

njzqxt

用什么数据啊？

凌寒

njzqxt 发表于 2013-7-22 18:25

                               
登录/注册后可看大图

用什么数据啊？

附件的文件

MeteoInfo

参加此帖：MeteoInfo读取AWX数据的改进
http://bbs.06climate.com/forum.p ... 900&fromuid=106