求假相当位温500-850hpa的计算公式和gs

蒹葭

求那位大神做过假相当位温500-850hpa的计算公式和gs，请发一下gs，谢谢，本人正在写论文可加作者

蒹葭

QQ459191259
电话18091239383

Masterpiece

核心代码这样，设置等等的自己写吧...


'set lev 850'
****计算假相当位温
'define prs=lev'
'set ccolor 1'
'define  l1=Tprs'
if(l1>273.15);
'define a1=0.622*6.11*EXP(17.26*(Tprs-273.16)/(Tprs-35.86))'
'define b1=prs-0.278*EXP(17.26*(Tprs-273.16)/(Tprs-35.86))'
'define  qs1=a1/b1'
endif
if(l1<=273.15);
'define a1=0.622*6.11*EXP(21.87*(Tprs-273.16)/(Tprs-7.66))'
'define b1=lev-0.278*EXP(21.87*(Tprs-273.16)/(Tprs-7.66))'
'define  qs1=a1/b1'
endif
'define  q1=qs1*Rprs'
'define e1=lev*q1/100.0/(0.62197+q1/100.0)'
'define tk1=55.0+2840.0/(3.5*log(Tprs)-log(e1)-4.805)'
'define pot1=Tprs*pow((1000/lev),(0.2854*(1.0-0.28*q1/100.0)))'
'define ept850=pot1*exp(((3376./tk1)-2.54)*q1/100.0*(1.0+0.81*q1/100.0))'


'set lev 500'
****计算假相当位温
'define prs=lev'
'set ccolor 1'
'define  l1=Tprs'
if(l1>273.15);
'define a1=0.622*6.11*EXP(17.26*(Tprs-273.16)/(Tprs-35.86))'
'define b1=prs-0.278*EXP(17.26*(Tprs-273.16)/(Tprs-35.86))'
'define  qs1=a1/b1'
endif
if(l1<=273.15);
'define a1=0.622*6.11*EXP(21.87*(Tprs-273.16)/(Tprs-7.66))'
'define b1=lev-0.278*EXP(21.87*(Tprs-273.16)/(Tprs-7.66))'
'define  qs1=a1/b1'
endif
'define  q1=qs1*Rprs'
'define e1=lev*q1/100.0/(0.62197+q1/100.0)'
'define tk1=55.0+2840.0/(3.5*log(Tprs)-log(e1)-4.805)'
'define pot1=Tprs*pow((1000/lev),(0.2854*(1.0-0.28*q1/100.0)))'
'define ept500=pot1*exp(((3376./tk1)-2.54)*q1/100.0*(1.0+0.81*q1/100.0))'

'define dt=ept500-ept850'

蒹葭

好象不行，能分别算出850和500的假相当位温，但是好象不计算差值，画不出来图？

蒹葭

Masterpiece 发表于 2017-10-29 15:53
核心代码这样，设置等等的自己写吧...

请教一下，算出来的ept850和ept500只能直接画图，没有存储吗,为什么再次读取就没有了？

Masterpiece

蒹葭 发表于 2017-10-29 18:43
请教一下，算出来的ept850和ept500只能直接画图，没有存储吗,为什么再次读取就没有了？

这个怎么存的看Grads手册吧

Masterpiece

蒹葭 发表于 2017-10-29 18:39
好象不行，能分别算出850和500的假相当位温，但是好象不计算差值，画不出来图？

计算了呀
'define dt=ept500-ept850'
dt就是你要的东西

蒹葭

做好了，谢谢