画NCEP2.5*2.5假相当位温

stone · 发表于 2013-11-2 19:56:44

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

看了很多论坛的帖子和大气物理书，附上假相当位温gs'reinit'
'sdfopen g:\air.2012.nc'
'sdfopen g:\rhum.2012.nc'
'set display color white'
'c'
'set xlopts 1 8 0.25'
'set ylopts 1 8 0.25'
'set clopts 1 8 0.25'
'set lon 105 130'
'set lat 25 50'
'set xlevs 110 120 130'
'set ylevs 30 40 50'
'set t 257'
'set lev 850'
'define t=air.1'
'define rh=rhum.2'
'define prs=lev'
'define es=(6.112*exp(17.67*(t-273.15)/(t-29.65)))'
'define q=rh*(0.62197*es/(prs-es))/100.'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(t)-log(e)-4.805)'
'define theta=t*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'
'set ccolor 1'
'set cthick 16'
'set grid off'
'set grads off'
'd eqt'
'enable print g:\tse_0508.gmf'
'print'
'disable print'

'reinit'
'sdfopen g:\air.2012.nc'
'sdfopen g:\rhum.2012.nc'
'set display color white'
'c'
'set xlopts 1 8 0.25'
'set ylopts 1 8 0.25'
'set clopts 1 8 0.25'
'set lon 105 130'
'set lat 25 50'
'set xlevs 110 120 130'
'set ylevs 30 40 50'
'set t 257'
'set lev 850'
'define t=air.1'
'define rh=rhum.2'
'define prs=lev'
'define es=(6.112*exp(17.67*(t-273.15)/(t-29.65)))'
'define q=rh*(0.62197*es/(prs-es))/100.'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(t)-log(e)-4.805)'
'define theta=t*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'
'set ccolor 1'
'set cthick 16'
'set grid off'
'set grads off'
'd eqt'
'enable print g:\tse_0508.gmf'
'print'
'disable print'
http://bbs.06climate.com/forum.php?mod=attachment&aid=MjEwOTl8NzNhODQwYWJiNTlkM2U1ZDNlMGU0ZmRjYjE3YWU0YTR8MTczMjczNjk2Ng%3D%3D&request=yes&_f=.gs

kongfeng0824 · 发表于 2013-11-2 20:29:01

谢谢分享

meizi0874 · 发表于 2013-11-3 16:41:37

好东西，但是我却无法下载（自动回复：请不要使用迅雷等下载工具，点我查看下载帮助）

polointer · 发表于 2013-11-11 22:57:35

感谢分享刚好用上

杨柳松 · 发表于 2013-12-4 18:10:30

楼主，谢谢你的gs,能不能麻烦你把计算公式的每一小步说明一下？？？（万分感谢）

liusx2007 · 发表于 2013-12-5 19:14:21

楼主，请把计算公式的出处详细说明一下。拜托了！！！
'define tlcl=55.0+2840.0/(3.5*log(t)-log(e)-4.805)'
'define theta=t*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'

shayingying · 发表于 2013-12-6 16:06:35

同问最后这四个公式的出处啊
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(t)-log(e)-4.805)'
'define theta=t*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'

风卷云 · 发表于 2014-5-23 17:28:12

最好是有图，这样才能更好地说明问题

白菜 · 发表于 2014-5-23 21:14:30

{:5_235:}

want_data · 发表于 2014-10-24 17:14:52

后面的公式从哪里来的

		自动登录	找回密码
密码			立即注册

[分享资料] 画NCEP2.5*2.5假相当位温

登录后查看更多精彩内容~