我想画一张湿位涡时间高度图 可惜 求指正

heyufei

[img]file:///C:\Users\Administrator\AppData\Roaming\Tencent\Users\863705527\QQ\WinTemp\RichOle\G`7YI7(G[A{$)912ZU$9A[F.jpg[/img]'reinit'
'open E:\aa\fnll\fnl_201307.ctl'
'enable print E:\aa\fnll\0706850.jpg'


'set z 1 20'
'set zlog on'
'set lon 80 130'

'set lat 10 60'
'define tc=tmpprs-273.15'
'define rh=rhprs'
'define prs=lev'  
'define es=(6.112*exp((17.67*tc)/(tc+243.5)))'
'define qs=0.62197*es/(prs-0.378*es)'
'define q=rh*qs/100'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(tc+273.16)-log(e)-4.805)'
'define theta=(tc+273.16)*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'

'define u=UGRDprs'
'define v=VGRDprs'
'define vo=hcurl(u,v)'
'define f=2*7.292*sin(lat*3.14159/180.0)*0.00001'
'define g=9.8'
'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=u(z-1)-u(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(eqt,x)'
'define dty=cdiff(eqt,y)'
'define pv1=-g*(vo+f)*deqt/dp'  
'define pv2=g*((dv/dp)*(dtx/dx)-(du/dp)*(dty/dy))'
'define pv=pv1+pv2'
'set xyrev on'
'set parea 2 7 3 7'
'set lon 114'
'set lat 30'
'set t 1 10'
'set grid off'
'set grads off'
'set zlog off'
'set gxout contour'
'set csmooth on'
'set cthick 7'
'set cint 0.3'
'd pv*10e5'
*'d pv2*10e6'
*'d pv*10e5'
'printim E:\aa\fnll\0706850.jpg white'

'reinit'
到底是哪里不对啊

heyufei

到底哪里有问题 应该是横纵标是时间 纵坐标是高度的

萍水啸天

好像不对呵

lqouc

把set t 1 10 放在计算前面。
grads默认的时间是t=1，所以define的pv只是t=1。多了自然没有。

heyufei

lqouc 发表于 2014-6-9 22:07
把set t 1 10 放在计算前面。
grads默认的时间是t=1，所以define的pv只是t=1。多了自然没有。

这样改了之后 变成了无限循环的状态

xry_2005

计算时t是不是应该是循环？然后再做时间剖面？

一水天下

四楼说的对，无限循环是什么意思？能上个错误截图吗

雨竹晓风

没有时间序列

huanggg

学习下，搞不懂

That_gal

比较想知道湿位涡和位涡的区别