赤道环流剖面图时出错

海盗船长

在画冬半年的赤道环流时感觉有问题，烦劳大家看一下。谢谢了！
'reinit'
'sdfopen d:\omega.mon.mean.nc'
'sdfopen d:\uwnd.mon.mean.nc'
'enable print  d:\equator.gmf white'

'set dfile 1'
'set lon 20 280'
'set lat 0'
'set lev 1000 100'
a1=11
while(a1<=801)
'define om=ave(omega,t='a1',t='a1'+2)'
a1=a1+12
endwhile

'set dfile 2'
'set lon 20 280'
'set lat 0'
'set lev 1000 100'
a2=11
while(a2<=801)
'define um=ave(uwnd,t='a2',t='a2'+2)'
a2=a2+12
endwhile

'set gxout stream'
'set strmden 6'
'd um;(-om*100)'
'print'
'disable print'
;

river

出错了错误提示呢

对雪吟歌

请附上错误提示……

学而不厌

我认为冬季是指12,1 2月，你改一下试试。a1=a2=12

海盗船长

学而不厌 发表于 2014-10-14 20:57
我认为冬季是指12,1 2月，你改一下试试。a1=a2=12

修改过，但图差别不太大的。

海盗船长

对雪吟歌 发表于 2014-10-14 20:54
请附上错误提示……

抱歉，我没传图，现在你看看

海盗船长

river 发表于 2014-10-14 20:13
出错了错误提示呢

抱歉，我没传图，现在你看看

学而不厌

'reinit'
'sdfopen d:\omega.mon.mean.nc'
'sdfopen d:\uwnd.mon.mean.nc'
'enable print  d:\equator.gmf white'

'set dfile 1'
'set lon 20 280'
'set lat 0'
'set lev 1000 100'
a1=12
n=0
'define q=0'
while(a1<=801)
'define om=ave(omega,t='a1',t='a1'+2)'
a1=a1+12
n=n+1
''q'='q'+om'
endwhile


'set dfile 2'
'set lon 20 280'
'set lat 0'
'set lev 1000 100'
a2=12
m=0
'define p=0'
while(a2<=801)
'define um=ave(uwnd,t='a2',t='a2'+2)'
a2=a2+12
m=m+1
''p'='p'+um'
endwhile


'set gxout stream'
'set strmden 6'
'd 'p'/'m';(-'q'/'n'*100)'
'print'
'disable print'
;

求教一下如何传图片上去，我貌似传不上

学而不厌

                               
登录/注册后可看大图

学而不厌

                               
登录/注册后可看大图