- 积分
- 947
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-4-12
- 最后登录
- 1970-1-1
|
发表于 2019-4-4 22:10:29
|
显示全部楼层
请问楼主,为什么我先画的阴影图后画等值线图,等值线还是没有显示呢
'reinit'
'sdfopen D:\20081228\NC\air.2008.nc'
'sdfopen D:\20081228\NC\rhum.2008.nc'
'set grads off'
'set grid off'
'set t 1447'
'set xlint 4'
'set ylint 100'
'set ylopts 1 6 0.15'
'set xlopts 1 6 0.15'
'set lat 20 40'
'set lon -88'
'set lev 1000 300'
'define t=air.1'
'define r=rhum.2'
'define prs=lev'
'define es=(6.112*exp(17.67*(t-273.15)/(t-29.65)))'
'define q=r*(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 2'
'set rgb 21 230 53 50'
'set rgb 22 233 55 58'
'set rgb 23 237 71 78'
'set rgb 24 239 95 102'
'set rgb 25 246 115 120'
'set rgb 26 247 180 183'
'set rgb 27 248 207 208'
'set rgb 28 252 233 233'
'set rgb 29 244 155 160'
'set rgb 30 244 247 249'
'set rgb 31 163 40 46'
'set rgb 32 197 29 24'
'set clevs 318 320 322 324 326 328 330 332 334 336 338'
*'set ccols 21 22 23 24 25 29 26 27 28'
'set ccols 30 28 27 26 29 25 24 23 22 21 32 31'
'set gxout shaded'
*'set cthick 6'
'd eqt'
'c'
*'run cbar_interp.gs'
'cbarn 1 0 6 0.3'
'set gxout contour'
'set cint 2'
*'set csmooth on'
*'set cterp on'
*'set cint 2'
'set cthick 6'
'set cstyle 1'
'set clopts 1 3 0.09'
'd eqt' |
|