整层水汽通量及散度的计算

cafflybird

以下分别计算出各个层次的qu,qv,hd,其中所用的数据是1948年1月到2013年1月的月平均年数据，即共781月的数据
'reinit'
'sdfopen e:\1year\air.mon.mean.nc'
'sdfopen e:\1year\rhum.mon.mean.nc'
'sdfopen e:\1year\uwnd.mon.mean.nc'
'sdfopen e:\1year\vwnd.mon.mean.nc'
'set fwrite e:\1year\1flux.dat'
'set gxout fwrite'
nt=1
'set t 'nt
while(nt<=781)
nz=1
while(nz<=8)
'set z 'nz
'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*q/9.8)'
'd qu'
nz=nz+1
endwhile
nt=nt+1
endwhile
nt=1
'set t 'nt
while(nt<=781)
nz=1
while(nz<=8)
'set z 'nz
'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 qv=(vwnd.4*q/9.8)'
'd qv'
nz=nz+1
endwhile
nt=nt+1
endwhile
nt=1
'set t 'nt
while(nt<=781)
nz=1
while(nz<=8)
'set z 'nz
'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*q/9.8)'
'define qv=(vwnd.4*q/9.8)'
'define hd=hdivg(qu,qv)'
'd hd'
nz=nz+1
endwhile
nt=nt+1
endwhile
'disable fwrite'
'reinit'
;
以下是进行整层计算前的ctl文件
dset e:/1year/1flux.dat
title vertically integrated moisture flux
UNDEF -9.99E33
xdef 144 LINEAR   0 2.5
ydef 73 LINEAR   -90 2.5
zdef 8 LEVELS 1000 925 850 700 600 500 400 300
tdef 781 LINEAR JAN1948 1mo
vars 3
qu 1 99 u_direction
qv 1 99 v_direction
hd 1 99 horizontal divergence
endvars
以下为整层计算
'reinit'
'open e:\1year\1fluxdata.ctl'
'sdfopen e:\1year\pres.sfc.mon.mean.nc'
'set fwrite e:\1year\flux_all.dat'
'set gxout fwrite'
'set z 1'
'define quall=vint(pres.2,qu,300)*9.8'
'define qvall=vint(pres.2,qv,300)*9.8'
'd quall'
'd qvall'
'define hdall=hdivg(quall,qvall)'
'd hdall'
'disable fwrite'
'reinit'
;
以下为整层文件所写的ctl
dset e:/1year/flux_all.dat
title vertically integrated moisture
UNDEF -9.99E33
xdef 144 LINEAR   0 2.5
ydef 73 LINEAR   -90 2.5
zdef 8 LEVELS 1000 925 850 700 600 500 400 300
tdef 781 LINEAR JAN1948 1mo
vars 3
quall 1 99 u_direction
qvall 1 99 v_direction
hdall 1 99 horizontal divergence
endvars
玛法麻烦大家帮我看下上面的程序哪有问题啊？当我最后单独绘制quall,qvall, hdall是发现整个图形一个颜色，期待高手指点啊，在线等待中。。。

cafflybird

顶一下，希望大家帮我看下我的错误啊

cafflybird

自顶一下，希望大家帮我看下我的错误啊

cafflybird

。。。。。。。

taxueyueming

帮顶

cafflybird

为何没人关注我的问题啊，太长了吗？

linjian

老哥，时间应该是最外层循环，你怎么每个量都在循环呀！

yuejun0026

前面提取qu，qv的那些数据你检验过可以出图不没？还有你最后画图的时候的gs文件贴出来一下赛

daodang2008

我也在学习，看不懂你这个

mrguo

帮顶！！！求大神解释