假相当位温垂直分布的问题

鄧晓孩

'reinit'
'sdfopen  w:\data\NCEPhighreanalysisdata\air.2013.nc'
'sdfopen  w:\data\NCEPhighreanalysisdata\rhum.2013.nc'
'set lon 95 125'
'set lat 15 40'
'set t 283'
'set lev 1000 100'

'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 grads off'
'set grid off'
'set csmooth on''set map 15 1 9'
'set parea 1 10.2 1.5 8'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set font 0'

'set lon 107.5'
'set lat 20 32.5'
'set zlog on'
'set ylevs 1000 925 850 700 600 500 400 300 250 200 150 100'
'set xlabs 20N|22.5N|25N|27.5N|30N|32.5N'
'set gxout shaded'
'set cint 0.3'
'd eqt'
'cbarn 1 0'
'set gxout contour'
'set cint 0.3'
'set cthick 6'
'd eqt'

'enable print w:\data\outputpicture\jxdwwchuizhi.gmf'
'print'
'disable print'

麻烦各位大大看一下，我想画的是假相当位温的垂直分布，但是出图的时候老是错误，麻烦各位帮忙解决一下。由于第一次发代码，所以排版不是太好，望大家见谅。

river

相对湿度的资料能到100毫巴的高度？

鄧晓孩

river 发表于 2016-5-10 08:23
相对湿度的资料能到100毫巴的高度？

多谢大大的提醒，我忽略了这个问题，只能到300毫巴，已经解决

363358643

这是正常值！

鄧晓孩

363358643 发表于 2017-11-25 22:04
这是正常值！

问题已解决，谢谢

xiadongdong

我也犯了同样的错误，看帖之后改好了，特来发帖感谢！

鄧晓孩

xiadongdong 发表于 2018-5-20 10:20
我也犯了同样的错误，看帖之后改好了，特来发帖感谢！

{:handshake:}

鄧晓孩

xiadongdong 发表于 2018-5-20 10:20
我也犯了同样的错误，看帖之后改好了，特来发帖感谢！

{:handshake:}