画逐层的距平值的垂直剖面

爱喝柠檬汁

自己弄了很久，就是画不出文献中的效果，下面是文献中的做法和图

爱喝柠檬汁

然后是我自己的脚本，希望路过的大神帮我看看，感激不尽！
'reinit'
'sdfopen e:\123.nc'
'set gxout fwrite'
'set fwrite e:\123.dat'
'set lon 123'
'set lat 35 60'
'set t 109'


i=1
while(i<=37)
'set lev 'i

'define avez=ave(z,lat=35,lat=60)'
'define avet=ave(t,lat=35,lat=60)'

'd z/9.8-avez/9.8'
'd t-avet'

i=i+1
endwhile

'disable fwrite'
'set gxout contour'
；
这个gs是将距平值转存到一个.dat文件中
然后是对应的ctl
dset e:\123.dat
undef -9.99E+33
title juping
xdef   1 linear  123  0.5
ydef  51 linear   35  0.5
zdef  37 levels  1000 975 950 925 900 875 850 825 800 775 750 700 650 600 550 500 450 400 350 300 250 225 200 175 150 125 100 70 50 30 20 10 7 5 3 2 1
tdef   1 linear  08z28aug2017 6hr
vars  2
zz     37  99 ** <z>     hpa
tt     37  99 ** <t>     k
endvars
最后是画图的gs
*****过冷涡中心温度、位势高度过123E距平
'reinit'
'open E:/08/tzjuping/123.ctl'
'set grads off'
'set mpdset hires'
'set map 15 1 1'

*'set t 1'
*'set x 1'
'set lat 35 60'
'set lev 1000 100'
'set zlog on'
*'set cmax 0'
*'set cint 10'
'set gxout shaded'
'd zz'
'run cbarn.gs'

'set gxout contour'
'set ccolor 15'
'set cstyle 1'
'set cthick 2'
'set cint 20'
'd tt'

'printim E:/tz123456.png white'
;
最后画出来的图是这样的，我觉得是错的，可是改来改去就是不对啊，

爱喝柠檬汁

然后是我自己的脚本，希望路过的大神帮我看看，感激不尽！
'reinit'
'sdfopen e:\123.nc'
'set gxout fwrite'
'set fwrite e:\123.dat'
'set lon 123'
'set lat 35 60'
'set t 109'


i=1
while(i<=37)
'set lev 'i

'define avez=ave(z,lat=35,lat=60)'
'define avet=ave(t,lat=35,lat=60)'

'd z/9.8-avez/9.8'
'd t-avet'

i=i+1
endwhile

'disable fwrite'
'set gxout contour'
；
这个gs是将距平值转存到一个.dat文件中
然后是对应的ctl
dset e:\123.dat
undef -9.99E+33
title juping
xdef   1 linear  123  0.5
ydef  51 linear   35  0.5
zdef  37 levels  1000 975 950 925 900 875 850 825 800 775 750 700 650 600 550 500 450 400 350 300 250 225 200 175 150 125 100 70 50 30 20 10 7 5 3 2 1
tdef   1 linear  08z28aug2017 6hr
vars  2
z      37  99 ** <z>     hpa
t      37  99 ** <t>     k
endvars
最后是画图的gs
'reinit'
'open E:/123.ctl'
'set grads off'
'set mpdset hires'
'set map 15 1 1'

'set t 1'
'set x 1'
'set lat 35 60'
'set lev 1000 100'
'set zlog on'
'set cmax 0'
'set gxout shaded'
'd z'
'run cbarn.gs'

'set gxout contour'
'set ccolor 15'
'set cstyle 1'
'set cthick 2'
'set cint 20'
'd t'

'printim E:/tz123456.png white'
;
最后画出来的图是这样的，我觉得是错的，可是改来改去就是不对啊，

爱喝柠檬汁

自己想了很久，百思不得其解啊！求各位帮我看看啊！

爱喝柠檬汁

下面另一种方法
'reinit'
'sdfopen e:\123.nc'
'set gxout fwrite'
'set fwrite e:\123.dat'
'set lon 123'
'set lat 45 60'
'set t 109'

i=1000
while(i>=100)
'set lev 'i

'define avez=ave(z,lat=45,lat=60)'
'define avet=ave(t,lat=45,lat=60)'

'd z/9.8-avez/9.8'
'd t-avet'

i=i-100
endwhile

'disable fwrite'
'set gxout contour'
；
这个gs是将距平值输出到一个.dat文件里下面是对应的ctl

dset e:\123.dat
undef -9.99E+33
title juping
xdef   1 linear  123  0.5
ydef  31 linear   45  0.5
zdef  10 levels  1000 900 800 700 600 500 400 300 200 100
tdef   1 linear  08z28aug2017 6hr
vars  2
z     10  99 ** <z>     hpa
t     10  99 ** <t>     k
endvars

下面是画图的gs
'reinit'
'open E:/123.ctl'
'set grads off'
'set mpdset hires'
'set map 15 1 1'

'set t 1'
'set x 1'
'set lat 45 60'
'set lev 1000 100'
'set zlog on'
'set cmax 0'
'set cint 10'
'set gxout shaded'
'd z'
'run cbarn.gs'

'set gxout contour'
'set ccolor 15'
'set cstyle 1'
'set cthick 2'
'set cint 20'
'd t'

'printim E:/11.png white'
;

爱喝柠檬汁

我看文献里的图只有一个核心，而我画的图有好几层的核心，我觉得是不对的，但就是画不出来文献中的效果

爱喝柠檬汁

有没有人理理我呀！画了半个多月了，

爱喝柠檬汁

谁理理我呀，求助！

爱喝柠檬汁

帮帮我呀！我实在是画不出来了

river

'reinit'
'sdfopen e:\123.nc'
'set t 109'
'set z 1 37'

'define avez=ave(z,lat=35,lat=60)'
'define avet=ave(t,lat=35,lat=60)'

'set lev 1000 100'
'set lon 123'
'set lat 35 60'
'set zlog on'
'set gxout shade2''set cmax 0'
'd z/9.8-avez/9.8'
'set gxout contour'
'd t-avet'

'printim ``````` '
；
试试看吧，我没太看懂文献上说的计算方法。