grads里计算位涡的gs，请大家帮忙看看是否正确，同EC自带的位涡不同呢。

会飞的星星星

'open E:\data\ec_data\0.5\high_201802-03.ctl'

path_pic='E:\ensemble\extra_tc\probability\20180228\fig\test'
'enable print 'path_pic'\pv_700_fnl.gmf'

t=98
while(t<=135)

'set gxout shaded'
'set t 't''
'set lat 20 50'
'set lon 110 145'
'set grads off'
'set grid off'

'set z 1 6'
'set lev 1000 200'
'set xlopts 1 6 0.12'
'set ylopts 1 6 0.12'
*------计算位温--------
'define e1=17.67*(T-273.16)/(T-29.65)'
'define es=6.112*exp(e1)'
'define e0=r*(es/100.0)'
'define q=(0.622*e0)/(lev-(0.378*e0))'
'define eqt=(T)*pow((1000/lev),(0.2854*(1.0-0.28*q)))'
*----------pv------------
'define vor=hcurl(u,v)*100000'
'define f=2*7.292*sin(lat*3.14159/180.0)*1e-5'
'set z 4'
'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 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 25 45'
'set lon 115 140'
*'set cint 0.1'
'd pv*10'

tt=(t-1)*6
time_out=time_loop(time_s,tt)

'draw string 4.2 8.0 'time_out' (UTC)'
'cbarn 1 1 10.3 4.1'
'print'

'c'

t=t+1
endwhile
'disable print'

会飞的星星星

大家都没人回复呢，路过的大神可以帮帮忙啦。

和清谈

请问time_out=time_loop(time_s,tt)中，time_loop是什么函数？

会飞的星星星

和清谈 发表于 2019-3-5 14:47
请问time_out=time_loop(time_s,tt)中，time_loop是什么函数？

是计算预报时效对应时间的子程序

mumululu

你计算的结果对么？

mumululu

'set z 1 6'
'set lev 1000 200'
这两个命令不冲突吗？