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

subtropical

平流层的萝卜 发表于 2013-10-28 20:23
楼主的公式有点问题：算es的时候，exp()括号里的系数用17.67和29.65是不对的，因为仅在求0度以下的水面饱和 ...

因为、楼主计算的是低层850的温度还没达到0度以下

subtropical

徐徐向前 发表于 2014-4-21 13:31
l楼主，我按你的方法只改了文件路径运行出来是这样的
Data Request Error: Invalid grid coordinates
  W ...

我也出现过这种问题，请问解决了吗？？

ymn

好。。。。。

limoxi

谢谢楼主分享！我按照您的思路计算完假相当位温把他们读取到grd文件中，画出来的图却是这样的，不知道为什么……


以下是我的gs，不知道哪里出了问题
'reinit'
'sdfopen  f:\data1\air2014.nc'
'sdfopen  f:\data1\rh2014.nc'
'set lon 60 180'
'set lat 10 70'
'set lev 850'
'set gxout fwrite'
'set fwrite f:/data1/theta.grd'
i=1
while(i<=92)
'set t 'i''
'define t=air.1'
'define rh=rhum.2'
'define prs=lev'
'define es=6.1078*exp(17.27*(t-273.26)/(t-35.86))'
'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))'
'd eqt-273.15
i=i+1
endwhile
'disable fwrite'
'reinit'
;

另外我的ctl文件是这样的
dset f:\data1\theta.grd
undef -9.99E+08
title Daily theta Data during last summer
xdef    144 linear 0 2.5
ydef    29 linear 0 2.5
zdef     1 levels 850
tdef   92 linear 01Jun2014 1dy
vars  1
theta   0  99  hks
endvars

求大神指导一下，我的数据是6月到8月逐日资料，画的是t=1时次的数据
先谢过了！

whmz

收藏了，谢谢！

zhyc

感谢分享，算了下，和ncl的结果很接近

zhyc

发现ncl算的就是假相当位温

superyoung16

感谢楼主！我画出来啦啦啦

dzn

谢谢，不错的帖子。

hugeal10000

原创的真厉害！