脚本

妮妮

我用的是160站的资料，想要画降水分布图。之前的站点资料转格点资料，以及两个CTL编辑都是正确的。但是最终画图时总是提示如图所示的错位，gs文件检查了很多遍还是不知道错在哪，请大家帮帮忙。现附上站点的ctl文件                                                                       

dset c:\rain_data\data2\grid.grddset c:\rain_data\data2\rain6.grd
dtype station
stnmap c:\rain_data\data2\rain6.map
undef -999.0
title rain
tdef 63 linear Jan1951 1yr
vars 1
p 0 99 rainfall data
ENDVARS
站点ctl文件
dset c:\rain_data\data2\grid.grd
undef -999.0
title Sample GRID Data
xdef 71 linear 70 1
ydef 41 linear 15 1
zdef 1 linear 1000 1
tdef 63 linear Jan1951 1yr
vars 1
g 0 99 grid data
endvars


以及gs文件
'open c:\rain_data\data2\grid.ctl'
'open c:\rain_data\data2\rain.ctl'
'enable print c:\rain_data\data2\rain6.gmf'
'set lon 73 135.5'
'set lat 17 55'
'set t 2'
'set mpdset cn cnriver'
'define a=oacres(g,p.2,10,7,4,2,1)'
'define a1=maskout(a,g-o.5)'
'define aa=smth9(a1)'
'set xlopts 1 10 0.18'
'set ylopts 1 10 0.18'
'set gxout shaded'
'c:\rain_data\data2\rgbset.gs'
'set clevs 0 10 50 100 200 300 400 500 600'
'set ccolor rainbow'
'd aa'
'cbarn.gs'
'set cthick 8'
'set clopts 1 6 0.1'
'set gxout contour'
'set cint 50'
'd aa'
'print'
'disable print'
;

四叶草

敢问楼主'define a1=maskout(a,g-o.5)'这一句中的o.5？还是0.5？

妮妮

谢谢啦，检查出问题了，就是那个0.5的问题，把0打成o了