请大家看哈这个计算假相当位温的.gs文件正确不

小梅子

'reinit'
'open E:\xiaodat/fnl2011.ctl'
'set lon 105 112'
'set lat 31 40'
'set lev 850'
'set t 19'
'set grid off'


'define prs=lev'
'define es=6.112*exp(17.67*(tmpprs-273.15)/(tmpprs-36.85))'计算饱和水气压
'define rs=0.622*es/(prs-es)'计算饱和混合比
'define q=rhprs*0.62197*es/(prs-0.378*es)'计算比湿
'define theta=tmpprs*pow(1000/(prs-es),(287.05/1004))'计算干空气的位温
'define lv=2.5e6-2323*(tmpprs-273.15)'计算汽化热
'define thetase=theta*exp(rs*lv/1004/tmpprs)'计算假相当位温
'set cint 4'
'set cthick 5'
'set gxout contour'
'd thetase-273.15'
'set mpdset shannx shanxi2_q'
'set map 6 1 6'
'draw map'
'printim E:\xiaodat\2912850.gif white'
;
这些公式是根据大气物理老师课堂上讲的写的，然后我绘制了一张图，可是不知道对不对，也不知道改怎么分析，请大家帮我看哈。在此谢谢啦！

~Catcher~

先注明一下，这也是我从别人的帖子中看到的，对比了好几种方法，发现这种更加严谨：
楼主若使用Tetens公式来计算饱和水汽压更好，下面为gs程序：

if(TMPprs>258.0)
'define a=17.2694'
'define b=35.86'
else
'define a=21.8746'
'define b=7.66'
endif

** Tetens公式***
'define es=6.112*exp(a*(TMPprs-273.15)/(TMPprs-b))'
**************

冬日情愫+

发表以下我个人的见解：1）首先要注意的是单位，比如 rhprs 应该是 % 的吧，比湿的单位我们多用g/g或者 g/kg，如果是前者则要在你定义的基础上除以100，如果是后者则要在你定义的基础上乘以10 ；2）如果初始位置空气已经是饱和湿空气的时候才可以直接用你的定义'define thetase=theta*exp(rs*lv/1004/tmpprs)'否则tmpprs应该用抬升凝结高度出的温度代替。

做个霸气的木头

http://bbs.06climate.com/forum.p ... &extra=page%3D1

huangYstar

完全不懂啊

fen蓝

我运行了一遍，想请教你，垂直剖面图该怎么样画出来？谢谢。

hrsmoto

下载来检验一下看，感谢分享

天边的那边

我发现假相当位温的算法很多

小蝌蚪

楼主你好，最近在研究假相当位温的计算。论坛里大家的算法有些出入，但是都没有人说明每个公式的出处。请问你程序里每个公式具体在大气物理的哪些页呢？非常感谢。
'define es=6.112*exp(17.67*(tmpprs-273.15)/(tmpprs-36.85))'计算饱和水气压
'define rs=0.622*es/(prs-es)'计算饱和混合比
'define q=rhprs*0.62197*es/(prs-0.378*es)'计算比湿
'define theta=tmpprs*pow(1000/(prs-es),(287.05/1004))'计算干空气的位温
'define lv=2.5e6-2323*(tmpprs-273.15)'计算汽化热
'define thetase=theta*exp(rs*lv/1004/tmpprs)'计算假相当位温

徐徐向前

感觉得注意单位

~Catcher~

谢谢楼主分享~