grads画图求助

木雨朵

'reinit'
'c'
*----------------input--------------------------
'sdfopen e:\new\p.nc'
tt=708
'set t 'tt''
while(tt<=1286)
'define pre1=3*ave(precip,t='tt',t='tt+2')'
tt=tt+12
endwhile
'define pre=ave(pre1,t=1,t=41)'
'set mproj latlon'
'set font 1'
'set lon 60 140'
'set lat 2 60'
'set xlint 10'
'set ylint 10'
'set grads off'
'set ccolor rainbow'
'set clab on'
'set cthick 2'
'set cstyle 1'
'set gxout shaded'
'set clevs  10 25 50 100 150 200 250 300 400 800 1000 2000'
'set mpdset hires cnworld'
'd pre'
'e:\new\station.gs'
'cbarn.gs'

*------画风场------------------------------
'open e:\new\ugrd.ctl'
'define uwnd=ave(uwnd.2,t=1,t=41)'
'open e:\new\vgrd.ctl'
'define vwnd=ave(vwnd.3,t=1,t=41)'
'set lon 60 140'
'set lat 2 60'
'set mpdset hires cnworld'
'set xlint 10'
'set ylint 10'
'set grads off'
'd uwnd;vwnd'

*----------------------------tu-----------------------------------------
'printim e:\new\winter.gif white x1000 y800'
'reinit'
上边是我的gs文件，为什么降水和风场不叠加呢？

木雨朵

lqouc 发表于 2015-7-28 14:10
网格不一样吧，应该贴上来ctl，可以用close关闭已经打开的文件。

dset e:\new\vmean.grd
title 1950-2011 may-sept mean v wind from the NCEP Reanalysis
undef -9.99e+08
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1 linear 0 1
tdef 49 linear 01may1960 12mo
vars 1
vwnd 0 t,y,x Monthly mean v wind
endvars
这是我做的ctl文件

木雨朵

E:\new\6-8.gif

木雨朵

但是我如果单独运行降水量，就可以作图，或者单独运行风场，也可以

木雨朵

file:///C:\Users\zhangxw\AppData\Roaming\Tencent\Users\466202071\QQ\WinTemp\RichOle\90I(ZUNW(LX{VYRWE54$8$D.png

lqouc

网格不一样吧，应该贴上来ctl，可以用close关闭已经打开的文件。

木雨朵

我用close把第一个文件给关掉以后就可以，谢谢啊！

中国电信6

谢谢楼主真的很棒