不同时段夏季水汽通量差值计算

cafflybird

数据及说明;
所用的数据下载地址
http://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.derived.pressure.html
air.mon.mean.nc
rhum.mon.mean.nc
uwnd.mon.mean.nc
vwnd.mon.mean.nc
我所求取的是1960-1980；1981-2011年夏季（6,7,8月）水汽通量的差值图（后段时期减去前短时期）
刚开始nt=145(代表1960年1月的时次)
nt=768(代表2011年12月的时次)
求取水汽通量的公式利用以下：
es=6.112*exp(17.269*(air-273.16)/(air-29.65))  （求饱和水汽压）
qs=qs=rhum.2*es/100'                           (求比湿)
v=mag(uwnd.3, vwnd.4)'                      （风场）
uvq=v*qs/9.8                              （uvq为最后求取的水汽通量）
对1960-2011年每一个月的水汽通量求完后，我最后计算了各时段夏季的平均值
前半期:   1960 -1980
    a1   t=150, t=390  (分别是1960年6月，1980年6月的时次)
   a2   t=151,t= 391  (同上7月)
   a3   t=152, t=392  (同上8月)
后半期： 1981-2011
    a4  t6=402,t=762(分别是1981年6月，2011年6月的时次)
    a5  t7=403,t=63 (同上7月)
   a6  t8=404,t=764 (同上7月)
所用到的程序：
'reinit'
'sdfopen e:\data0220\air.mon.mean.nc'
'sdfopen e:\data0220\rhum.mon.mean.nc'
'sdfopen e:\data0220\uwnd.mon.mean.nc'
'sdfopen e:\data0220\vwnd.mon.mean.nc'
nt=145
while(nt<=768)
'set t 'nt''
'set lev 500'
'define es=6.112*exp(17.269*(air-273.16)/(air-29.65))'
'define qs=rhum.2*es/100'
'define v=mag(uwnd.3,vwnd.4)'
'define uvq=v*qs/9.8'
'set grads off'
'set grid off'
'set mpdset cnworld'
'set lon 40 170'
'set lat 10 70'
'set gxout shaded'
'd uvq'
'run cbarn'
'print'
'c'
nt=nt+1
endwhile
'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'
cbarn
'd sua'
'printim e:\data0220\sua.bmp white'
'print'
'disable print'

残阳飞血

谢谢楼主分享

小船

楼主这个程序用过么，我用的时时候出来的结果是一片红色，分别求a1+a2+a3的值做的图和a4+a5+a6的值做的图是一样的，这是怎么回事

cafflybird

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

                               
登录/注册后可看大图

楼主这个程序用过么，我用的时时候出来的结果是一片红色，分别求a1+a2+a3的值做的图和a4+a5+a6的值做的图是 ...

这个程序有问题啊。。。

龙行天下

非常感谢楼主，介绍得非常详细，适合初学者

龙行天下

你的东西对我太有用了，最近正想做这个，菜鸟一只，以后还望楼主多多指教。给你加点分，支持一下

龙行天下

才发现，你这个程序有问题，请问楼主修改过了吗

shanshanhai

谢谢楼主分享

baggio

{:5_275:}{:5_275:}{:5_275:}{:5_275:}