E-P通量计算程序

wrf_lanzhou

一个用grads编的e-p通量计算程序，使用ncep资料
'reinit'
'sdfopen uwnd.mon.mean.nc'
'sdfopen vwnd.mon.mean.nc'
'sdfopen air.mon.mean.nc'
i=565
while(i<=575)
'set lon 0 360'
'set lat 0 90'
'set z 1 17'
'set t 'i
'define prs=lev'
'define f=2*7.29e-5*sin(lat*3.1415/180)'
'define a=6.37e6'
'define uu=uwnd.1-ave(uwnd.1,lon=0,lon=360,-b)'
'define vv=vwnd.2-ave(vwnd.2,lon=0,lon=360,-b)'
'define uuvv=uu*vv'
'define uv=ave(uuvv,lon=0,lon=360,-b)'
'define temp=air.3+273.15'
'define theta=temp*pow(1000/prs,0.286)'
'define theaa=ave(theta,lon=0,lon=360,-b)'
'define dtheta=theta-theaa'
'define dthetavv=dtheta*vv'
'define thetav=ave(dthetavv,lon=0,lon=360,-b)'
'define dp=100*(lev(z-1)-lev(z+1))'
"define dtheta1=(air.3(z-1)+273.15)*pow(1000/lev(z-1),0.286)-(air.3(z+1)+273.15)*pow(1000/lev(z+1),0.286)"
'define dthetadp=dtheta1/dp'
'define fx=-1*a*cos(lat*3.1415/180)*uv'
'define fy=f*a*cos(lat*3.1415/180)*thetav/dthetadp'
'set z 2 16'
'define epdiv=hdivg(fx,fy)'
'c'
'set t 'i
'set z 2 16'
'set lat 0 90'
'set lon 60'
'set mpdset cnworld'
'set gxout shaded'
'd uu'
'cbarn'
'draw title 'i-564
'printim 'i-564'_60E.bmp white'  
'c'
'set z 2 16'
'set lat 0 90'
'set lon 60'
'set mpdset cnworld'
'set gxout shaded'
'd uu*uu'
'cbarn'
'draw title 'i-564
'printim 'i-564'_60E.bmp white'  
'c'
'set gxout shaded'
'set cthick 10'
'set z 2 16'
'set lat 0 90'
'set lon 60'
'set mpdset cnworld'
'd vv'
'cbarn'
'draw title 'i-564
'printim 'i-564'_60E.bmp white'  
'c'
'set gxout shaded'
'set z 2 16'
'set lat 0 90'
'set lon 60'
'set mpdset cnworld'
'd (uu*uu+vv*vv)/2'
'cbarn'
'draw title 'i-564
'printim 'i-564'_60E.bmp white'  
'c'
'set gxout shaded'
'set z 2 16'
'set lat 0 90'
'set lon 60'
'set mpdset cnworld'
'd uv'
'cbarn'
'draw title 'i-564
'printim 'i-564'_60E.bmp white'  
'c'
'set gxout contour'
'set ccolor 1'
'set cthick 10'
'set z 2 16'
'set lat 0 90'
'set lon 60'
'set mpdset cnworld'
'd thetav'
'draw title 'i-564
'printim 'i-564'_60E.bmp white'
'c'
'set gxout contour'
'set ccolor 1'
'set cthick 10'
'set z 2 16'
'set lat 0 90'
'set lon 60'
'set mpdset cnworld'
'd uv'
'draw title 'i-564
'printim 'i-564'_60E.bmp white'
'c'
'set gxout contour'
'set ccolor 1'
'set cthick 10'
'set z 2 16'
'set lat 0 90'
'set lon 60'
'set mpdset cnworld'
'd fx'
'draw title 'i-564
'printim 'i-564'_60E.bmp white'
'c'
'set gxout contour'
'set ccolor 1'
'set cthick 10'
'set z 2 16'
'set lat 0 90'
'set lon 60'
'set mpdset cnworld'
'd fy'
'draw title 'i-564
'printim 'i-564'_60E.bmp white'
  'c'
'set zlog on'
'set z 2 16'
'set lat 0 90'
'set lon 60'
'set gxout shaded'
'd ave(uwnd.1,lon=0,lon=360,-b)'
'cbarn'
'set gxout contour'
'set ccolor 1'
'set cthick 10'
'set mpdset cnworld'
'd fx;fy*-100'
'set gxout contour'
'set ccolor 1'
'set cthick 10'
'd epdiv'
'draw title 'i-564
'printim 'i-564'_60E.bmp white'
  i=i+1
endwhile

/xin村儿/

请给出贴图。

liuzf406

你好，正好要用到Ｅ-Ｐflux, 有图吗？谢谢了！

pqman3

算法说明中说：“for top and bottom levels one sided finite difference is used.”

sxq

跟着一起学习学习

sxq

楼主用的是哪个E-P通量的公式呀

xy15036

太好了，谢谢分享

xugk

感谢楼主的分享

★翼★

挺牛啊