多年平均夏季整层水汽通量的计算

胡菡

这是自己根据手册和论坛里编写的脚本写的，但是无法运行。请教高手该如何修改！！！
资料采用NECP再分析月分析数据做出1986-2003年中国夏季（6-8月）平均整层水汽通量的图。
是时间编写出了问题吗？

'reinit'
'set grid off'
'set grads off'
'set poli on'
'set parea 1 10 1 8'
'sdfopen D:\GrADS\shum.mon.mean.nc'  
'sdfopen D:\GrADS\uwnd.mon.mean.nc'
'sdfopen D:\GrADS\vwnd.mon.mean.nc'
'sdfopen D:\GrADS\pres.mon.mean.nc'
'set lat 4 53'
'set lon 73 135'
'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6'
'set mpdset cnworld hires'
'set map 1 1 1'
'draw map'
'set ylopts 1 4 0.12'
'set xlopts 1 4 0.12'
'set cthick 4'
'set z 1 8'
'define qu=uwnd*shum/9.8'
'define qv=vwnd*shum/9.8'
'set z 1'
'set t 456 660'
'define uqs=vint(pres(z=1),qu,300)/1000'
'define vqs=vint(pres(z=1),qv,300)/1000'
'define sq=mag(uqs,vqs)'
'set t 6 8'
'define sq1=ave(sq,t=456,t=660,12)'
'modify sq1 seasonal'
'define sq2=sum(sq1,t=462,t=464)'
'd sq2'
'run cbarn'
'draw title 1986-2003 year summer average vapor flux china'
'printim D:\GrADS.gmf x1000 y800 white'

胡菡

C:\Users\Huhan\Desktop

ric_bob

打酱油滴飘过。。

Crystal_qiong

打开多个文件，要在变量后面标注.1,.2.3的呀，像你的，shum.1,uwnd.2,vwnd.3,pres.4

龙行天下

请问你的整层水汽通量算出来了吗？

baggio

楼主算出来了么？有的话可以分享下么？

玄枫

楼主做出来了没？分享下吧！另外，红色的部分是什么意思，哪位大神能解释下吗？非常感谢

D.B.L

与@Crystal_qiong说的一样，要加打开顺序

LNYZ

想问楼主解决问题没

tcj-mlh

打酱油滴飘过。。