求取500hpa 1960-1980夏季水汽通量的变化

cafflybird

以下是我的程序，我所用的资料都是NCEP/NCAR数据，程序后面a1,a2,a3分别为1960-1980 6月，7月，8月的水汽通量平均值，a4,a5,a6分别为1981-2011 6月，7月，8月的水汽通量平均值，程序运行时，没有错误，但是没有图出现？期待高手解答啊。。。
'reinit'
'sdfopen e:\data0308\air.mon.mean.nc'
'sdfopen e:\data0308\rhum.mon.mean.nc'
'sdfopen e:\data0308\uwnd.mon.mean.nc'
'sdfopen e:\data0308\vwnd.mon.mean.nc'
'set t 145 768'
'set lev 500'
'define p=lev'
'define es=6.112*exp(17.67*air/(air+273.15-29.65))'
'define qs=0.622*es/(p-0.378*es)'
'define q=rhum.2*qs/100'
'define qu=(uwnd.3*10*q/9.8)'
'define qv=(vwnd.4*10*q/9.8)'
'define uvq=1000*mag(qu,qv)'
'set grads off'
'set grid off'
'set mpdset cnworld'
'set lon 40 170'
'set lat 10 70
'set gxout shaded'
'd uvq'
cbarn
'define a1=ave(uvq,t=150,t=390,12)'
'define a2=ave(uvq,t=151,t=391,12)'
'define a3=ave(uvq,t=152,t=392,12)'
'define a4=ave(uvq,t=402,t=762,12)'
'define a5=ave(uvq,t=403,t=763,12)'
'define a6=ave(uvq,t=404,t=764,12)'
'define sua=(a4+a5+a6)/3-(a1+a2+a3)/3'
'set lon 40 170'
'set lat 10 70'
'set gxout shaded'
'd a1'
cbarn
'printim e:\data0308\sua.bmp white'

jacci_zou

pull dummy 貌似是少了这个吧，最后一行，如果你的GS是正确的话

wet510

试试在'd a1'前加上'set t 1'将时间维固定后，仅有空间维，看看能否出图

cafflybird

jacci_zou 发表于 2013-3-9 08:42

                               
登录/注册后可看大图

pull dummy 貌似是少了这个吧，最后一行，如果你的GS是正确的话

我试了一下，还是不行

cafflybird

wet510 发表于 2013-3-9 10:06

                               
登录/注册后可看大图

试试在'd a1'前加上'set t 1'将时间维固定后，仅有空间维，看看能否出图

还是没图啊

wet510

cafflybird 发表于 2013-3-9 15:06

                               
登录/注册后可看大图

还是没图啊

你指图形显示窗口没图？还是输出路径处没图？，如果是后者，建议把输出文件的后缀名改成.png或.gif。不知道printim支不支持.bmp格式的输出。。。。。

小船

我用的程序和你的一样，我的出来了图不过是一片红色，同样纠结中

cafflybird

小船 发表于 2013-6-19 14:33

                               
登录/注册后可看大图

我用的程序和你的一样，我的出来了图不过是一片红色，同样纠结中

我也不知道这样错误在哪里，但是当我把uvq单独提取出来存储成*.dat格式的数据，在以此数据为基础进行后半部分计算时是正确的。