批量画图加地形怎么画

小甜甜

我想在每张图上都画高于1500的地形，但是总画不出，是因为有t循环，但是高度数据只有一个时次的缘故吗？
下面是我的GS文件，和GRADS提示



'reinit'
'sdfopen j:\qx\data\pres.slp\slp.1963.nc'
'sdfopen j:\qx\data\u.sig995\uwnd.sig995.1963.nc'
'sdfopen j:\qx\data\v.sig995\vwnd.sig995.1963.nc'
'sdfopen j:\qx\data\hgt.sfc\hgt.sfc.nc'
'enable print j:\qx\cold-high\1963.12.19.20-22.20.gmf'

i=1411
while(i<=1423)
*'set parea 0.5 10.8 2.0 7.0'

'set t 'i
'set lat 5 85'
'set lon 10 180'
'set grid on'
'set grads off'
'set xlpos -20'
'set ylpos -20'
*'set xlopts 1 5 0.13'
*'set ylopts 1 5 0.13'
*'set font 0'
'set cterp on'
*'set xlint 10'
*'set ylint 10'

'set gxout contour'
'set cthick 8'
'set ccolor 1'
'set cint 2.5'
*'set cint 4'
'set clskip 2'
'set clopts 1 5 0.1'
'd slp.1/100'

'set ccolor 4'
*'set clskip 4'


'set gxout vector'
'set t 'i
'set grads off'
'set cterp on'
'set cthick 6'
'set ccolor 4'
'set arrscl 0.5 20'
'd uwnd.2;vwnd.3'

'set t 1'
'set xlint 10'
'set ylint 10'
'set grads off'
'set gxout shaded'
'set cmin 1500'
'd hgt.4'
'cbar 1.0 0'

a=math_int(i/4)
b=math_mod(i,4)
*it=a+1-180

*第1时次（4月1日）是第361时次

if b=1
itt=08;ittt=00
endif
if b=2
itt=14;ittt=06
endif
if b=3
itt=20;ittt=12
endif
if b=0
itt=02;ittt=18
endif
it=a+1
'draw title  'i' 'it' 'itt' ('ittt')'

*'q w2xy 100 26'
'q w2xy 80 25'
x1=subwrd(result,3)
y1=subwrd(result,6)
*'q w2xy 108 26'
'q w2xy 140 25'
x2=subwrd(result,3)
y2=subwrd(result,6)
*'q w2xy 100 33'
'q w2xy 80 60'
x3=subwrd(result,3)
y3=subwrd(result,6)
*'q w2xy 108 33'
'q w2xy 140 60'
x4=subwrd(result,3)
y4=subwrd(result,6)
*'draw poly x1 y1 x2 y2 x3 y3 x4 y4'还是不对
'draw rec 'x1' 'y1' 'x4' 'y4''
'run axis.gs -type b  -label on -asuffix off -position o -interval 10 -sinterval 5 '
'run axis.gs -type T  -label on -asuffix off -position o -interval 10 -sinterval 5 '
'run axis.gs -type l  -label on -asuffix off -position o -interval 10 -sinterval 5 '
'run axis.gs -type R  -label on -position o -interval 10 -sinterval 5  '
'print'
'c'
i=i+1
endwhile
'disable print'
'close 4'
'close 3'
'close 2'
'close 1'

传说中的谁

大概是的，画地形的时候，在地形变量后面加上(t=1)

梦之蓝

挺好用的啊

Clover

'sdfopen j:\qx\data\hgt.sfc\hgt.sfc.nc'
这里hgt.sfc.nc在哪里下？