假相当位温的剖面图

panghuaji

'reinit'
'open f:\data\grib.ctl'
'set grads off'
'set grid off'
'set font 1'
'set mpdset hires'
'set xlint 1'
'set ylint 100'
'set xlopts 1 4 0.18'
'set ylopts 1 4 0.18'
'set clopts 1 -1 0.15'
'set cint 0.2'

'set lon 60 180'
'set lat 20 70'
'set lev 1000 100'

***气温
'define tmp=TMPprs'
***相对湿度
'define rh=RHprs'
***比湿
'define q=SPFHprs'
'define prs=lev'

'define e1=17.67*(tmp-273.16)/(tmp-29.65)'
'define es=6.112*exp(e1)'
'define qs=0.62197*es/(prs-0.378*es)'
*'define rh=q/qs*100.'
'define e=prs*q/(0.62197+q)+1e-10'
'define tk=55.0+2840.0/(3.5*log(tmp)-log(e)-4.805)'
'define pot=tmp*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt=pot*exp(((3376./tk)-2.54)*q*(1.0+0.81*q))'

'set lat 50'
'set clab on'
'set zlog on'
'set clskip 2'
'set cint 0.2'
'set ccolor 1'
'set cthick 6'
'set lon 145 151'
'set lev 975 300'

*'d pot'
'd eqt'

mofangbao

楼主多少写点文字介绍下哈  比如用的数据源等等

river

说实话，坛子里已经有好多这样的例子了当然多多益善啊，只是楼主也没说用的那些资料什么的

panghuaji

用的是NCEP的再分析资料，格点为0.5*0.5，6小时

nuistyy

多谢LZ分享，试试能不能用

huangYstar

怎么用nc的文件  做假相当位温啊

quentin011099

panghuaji 发表于 2013-3-18 20:44

                               
登录/注册后可看大图

用的是NCEP的再分析资料，格点为0.5*0.5，6小时

请问楼主这个0.5*0.5的资料在哪里下载啊？我只找得到1*1的。还有，楼主我作图的时候一直说 'define es=6.112*exp(e1)'出错了。我实在找不出原因啊。请指点啊谢谢，希望能知道您是怎么读数据什么的。

jingworld

ncep 1*1资料里是没有spfhprs的。。。比湿都是自己算的 楼主你的0.5度资料里有吗

刘振

quentin011099 发表于 2013-5-14 16:12

                               
登录/注册后可看大图

请问楼主这个0.5*0.5的资料在哪里下载啊？我只找得到1*1的。还有，楼主我作图的时候一直说 'define es=6. ...

你可能是数值超出范围了，太大了在终端输入的话。

心语