nc数据格式文件画出的500hPa高度图不对劲？！

945924058

'reinit'
'sdfopen hgt.1970.nc'
'sdfopen uwnd.1970.nc'
'sdfopen vwnd.1970.nc'
'set vpage 0 11 0 8.5'
'set parea 1 10 1 7.5'
'set mpdset cnworld'
'set lon 80 135'
'set lat 0 50'
*更改时间
'set time 00Z01JAN1970'
*更改高度
'set lev 500'
'set grid off'
'set grads off'
'set xlpos -20'
'set ylpos -20'
*绘制v风场
'set gxout shaded'
'run define_colors.gs'
'set clevs 20 25 30 35 40 45 50'
'set ccols 0 52 54 56 37 35 22 24'
'd vwnd.3'
'cbar 1 0'
*绘制风场
'set z 1'
'set gxout barb'
'set digsize 0.04'
'd uwnd.2;vwnd.3'
*绘制高度场
'set gxout contour'
'set ccolor 1'
'set clab forced'
'set cstyle 1'
'set cthick 6'
'set cint 2'
'set clopts 1 1 0.12'
'd hgt/10'
'run sansiaxis.gs -type b -label on -position o -sinterval 2  -lfont 5 -lsize 0.15 -lthick 0.8 -suffix `3. -interval 15'
'run sansiaxis.gs -type L -label on -position o -sinterval 2  -lfont 5 -lsize 0.15 -lthick 0.8 -suffix `3. -interval 15'
'draw string 5.5 7.8 at 1970/01/01/00(UTC)'
'printim 1970010100.png white'

你好，我想请问我用这个gs画出的500hPa的高度图不对劲，请问是哪里出了问题

占戈示申

set lev * 和set z *，的意义是一样的，都是对高度进行定义。
gs中，最上面有set lev 500，是定义到500hPa层，而下面的set z 1 ，则是把高度定义到第1层了，所以画的等和高度场都是第1层的，应该是10hPa层吧（近地面层）。删除set z 1，应该就好了。

来自苹果APP客户端

river

一楼基本说的是没错的，唯一的问题是对于 set z 1 的解释。这个是指资料的第一层，按照GrADS的资料排列原则是从低到高排列，所以set z 1 指的是这个资料的最下面一层，如果这个资料最底下是1000hPa，则set z 1 等价于 set lev 1000