假相当位温Θse计算--GS

zsmnantong · 发表于 2011-8-24 17:05:25

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

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

x

本帖最后由 zsmnantong 于 2011-8-26 12:01 编辑

NCEP2.5*2.5格点资料http://www.esrl.noaa.gov/psd/data/gridded/reanalysis/计算假相当位温Θse--脚本
'reinit'
'sdfopen f:\ncep\air.2011.nc'
'sdfopen f:\ncep\rhum.2011.nc'
'set lon 80 130'
'set lat 15 45'
'set lev 850'
'set mpdset cnworld'
'set grid on'
i=765
while(i<=779)
'set t 'i
'set ylopts 1 4 0.15'
'set clopts 1 4 0.09'
'set gxout contour'
'set xlint 2'
*'set arrowhead 0.05'
*'set arrscl 0.3'
'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 cint 4'
'set gxout contour'
'set cthick 5'
'd eqt-273.15'
'printim f:\picture\theta07-'i'.png white'
'c'
i=i+1
endwhile

rhythm9806 · 发表于 2011-8-25 11:09:50

有点错误的
'define tlcl=55.0+2840.0/(3.5*log(t+273.16)-log(e)-4.805)'
这里的log(t+273.16)-不用加273.16了

zsmnantong · 发表于 2011-8-25 11:06:16

自己顶一下，太靠后不好找

善人/jw · 发表于 2011-8-24 17:12:02

做个沙发，顶~~

mofangbao · 发表于 2011-8-24 17:13:52

板凳也要顶起的原创东东难得啊

Mouse · 发表于 2011-8-24 17:26:30

原创贴子，必须狠狠加分

尽头的尽头 · 发表于 2011-8-24 17:32:33

感谢感谢

zsmnantong · 发表于 2011-8-24 17:55:13

各位版主多发发原创啊

artmunich · 发表于 2011-8-24 19:01:57

先留名，感谢一下~
待我回去自己检验一下，可能还得麻烦lz呢

黄河沙 · 发表于 2011-8-24 19:04:25

收藏了，谢谢！

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

[分享资料] 假相当位温Θse计算--GS

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

评分