帮忙看下这个gs是计算位涡的，错在哪里？

蒹葭

'reinit'
'open f:\2012ncep\720.ctl'
'enable print f:\2012ncep\weiwo.gmf'
'set mpdset cnworld'
'set grads off'
'set ylevs 1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100'
'set lev 995 300'
'set z 1 21'
*'set zlog on'
'set lat 20 60'
'set lon 90 120'
'set lev 1000 100'
'set t 1'
'set clab forced'
'define u=ugrdprs'
'define v=vgrdprs'
'define vor=hcurl(UGRDprs,VGRDprs)*100000'
*'set lat 38'
'define f=2*7.292*sin(lat*3.14159/180.0)'
'define rh=rhprs'
'define T=tmpprs'
'define e1=17.67*(T-273.16)/(T-29.65)'
'define es=6.112*exp(e1)'
'define e0=rh*(es/100.0)'
'define q=(0.622*e0)/(lev-(0.378*e0))'
'define qs=0.62197*es/(lev-0.378*es)'
*'define rh=q/qs*100.'
'define e=lev*q/(0.62197+q)+1e-10'
'define tk=55.0+2840.0/(3.5*log(T)-log(e)-4.805)'
'define pot=(T)*pow((1000/lev),(0.2854*(1.0-0.28*q)))'
'define eqt=pot*exp(((3376./tk)-2.54)*q*(1.0+0.81*q))'
*'d pot'
'set grid off'
'set grads off'
'set cint 2'
'set lon 90 120'
'set lat 30 55'
'set lev 500'
'set t 3'
'draw title eqt'
'd eqt-273.16'
'draw title these'
'print'
'c'
'set lat 38'
'set lon 90 120'
'set lev 900 300'
'define dp=100*(lev(z-1)-lev(z+1))'
'define deqt=eqt(z-1)-eqt(z+1)'
'define du=u(z-1)-u(z+1)'
'define dv=v(z-1)-v(z+1)'
'define dx=2.0*6370949.0*cos(lat*3.14159/180.0)*3.14159/180.0'
'define dy=2.0*6370949.0*3.14159/180.0'
*'define dtx=cdiff(T,x)'
*'define dty=cdiff(T,y)'
'define pv1=-9.8*(vor+f)*deqt/dp'  
'define pv2=9.8*((dv/dp)*(deqt/dx)-(du/dp)*(deqt/dy))'
'define pv=pv1+pv2'
*'set lat 34'
*'set gxout shaded'
*'set cmax 0'
'draw title weiwo'
'd smth9(pv)'
*'set gxout contour'
*'d smth9(pv)'
*'set cint 1'
*'d smth9(pv)'
*'run cbarn.gs'
'print'
'c'
'disable print'
'reinit'

蒹葭

请高手指点一下，不胜感谢

houtian

蒹葭 发表于 2013-1-29 21:05

                               
登录/注册后可看大图

请高手指点一下，不胜感谢

你好， 看到你计算位涡的程序 请问你有没有介绍计算位涡的文档？

1417920146

这句不对哦 'define f=2*7.292*sin(lat*3.14159/180.0)' 少乘1e-5了

jackblak

顶一下！好好学习一下！！

That_gal

我参考一下，谢谢楼主

meizi0874

我参考一下，谢谢楼主

meizi0874

楼主问题解决了吗？要是解决了，放一个完整的程序吧

wuyu

顶一下，学习中

xuehongyang

{:eb303:}{:eb303:}{:eb303:}