- 积分
- 4987
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-9-9
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 鄧晓孩 于 2016-5-9 23:34 编辑
'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'
麻烦各位大大看一下,我想画的是假相当位温的垂直分布,但是出图的时候老是错误,麻烦各位帮忙解决一下。由于第一次发代码,所以排版不是太好,望大家见谅。
|
-
错误显示
|