求助 如何用FNL资料绘制湿位涡的垂直剖面图

jodia

gs如下，u（z+1）与u（z-1）的值相同，没法控制表示不同层次的变量值，求帮助，万分感谢！
'reinit'
'open d:/data/0506/fnl_20050609_00_00.ctl'
'enable print d:/data/vel_pv_1.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 15 55'
'set lon 105 135'
*'set lev 1000 100'
'set t 1'
'set clab forced'
'define u=ugrdprs'
'define v=vgrdprs'
'define t0=TMPprs'
'define vor=hcurl(UGRDprs,VGRDprs)*100000'
'define f=2*7.292*sin(lat*3.14159/180.0)'
'define g=9.8'

*'define prs=1.333*760*pow((1-hgtprs/44300),5.256)'
'define prs=lev'
'define es=(6.112*exp((17.67*(t0-273.16))/(t0-29.65)))'
'define qs=0.62197*es/(prs-0.378*es)'
'define q=qs*RHprs/100'
'define e=hgtprs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(t0)-log(e)-4.805)'
'define theta=t0*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 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(t0,x)'
'define dty=cdiff(t0,y)'

'define pv1=-g*(vor+f)*deqt/dp'  
'define pv2=g*((dv/dp)*(dtx/dx)-(du/dp)*(dty/dy))'
'define pv=pv1+pv2'
'set lat 34'
'set gxout shaded'
'set cmax 0'
'd smth9(pv)'
'set gxout contour'
'd smth9(pv)'
'set cint 1'
'd smth9(pv)'
'run cbarn.gs'
'print'
'disable print'
;

mofangbao

最近怎么都是这种提问方式，什么信息都没有，就贴个一段程序，自己也不先找找哪里错了，是不是没有错误提示，不知道哪里错了就一行一行调试

jodia

啊？我第一行不是说了吗？画剖面图的时候u(z-1)和u(z+1)的值相同，'du=u(z-1)-u(z+1)'一直是等于0 ，没有办法表示不同高度的值，只能算某一固定高度的湿位涡

mofangbao

分别画一下各层次，看看数值是不是相同的

传说中的谁

不知道什么错误，但是有一点，用这种差分方式，肯定不能从第一层算起，也不能算到最后一层，就是说，最多的层数也只能是'set z 2 20'

jodia

传说中的谁 发表于 2012-4-21 16:51

                               
登录/注册后可看大图

不知道什么错误，但是有一点，用这种差分方式，肯定不能从第一层算起，也不能算到最后一层，就是说，最多的 ...

是的，我后来也发现了，改成了'set z 2 20'，不过还是做不出来

jodia

mofangbao 发表于 2012-4-21 16:01

                               
登录/注册后可看大图

分别画一下各层次，看看数值是不是相同的

嗯，我试了，固定了层次，单独绘制某一高度场可以画出图来，每层高度的图都不同，但是层次取一个范围之后，就没法计算差分，剖面图没法画，不知道什么原因

littlelock

我也遇到过这样的问题，后来我先把每一层的假相当位温算出来，再计算差分就可以了。
只是计算出来的mpv比一般文献中的大，有达6PVU的，不知道楼主有没有这样的问题？

acmilan

楼主，这个问题解决了没啊？

Free__Bird

调好了的gs文件传上来哟！呵呵