画沿某经度范围的500hpa平均涡度随时间演变的图

cici破晓碎寒

'reinit'
'sdfopen f:\shixi\uvw.nc'
'set t 37 72'
'set lat 15 55'
'set lev 500'
*-----vorticity------
'define vor1=hcurl(u,v)*1e5'
'define vor2=ave(vor1,lon=50,lon=55)'
*--------GRADS设置------------
'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'
'enable print f:\shixi\500vor.gmf'
'set grads off'

'set lon 118'
'set zlog on'
'set ylevs 15 20 25 30 32 35 40 45 50 55'
'set xlabs 06|12|18|00|06|12|18|00|06|12|18|00|06|12|18|00|06|12|18|00|06|12|18|00|06|12|18|00|06|12|18'
'set gxout shaded'
'set cint 1'
'set cmax 0'
'd vor2'
'run cbar.gs'
'set ylpos 0 r'
'set gxout contour'
'set cint 1'
'set cthick 11'
*'set clskip 1'
'd vor2'

'writehz 0.4 4 纬度 1 2 1 1.3 0.6 90 7'

'print'
'disable print'
'printim f:/shixi/500vor.gif white'
;纵坐标显示的一直不对，向请教各位大神，这是为啥？

river

是你自己搞错了。按照你的脚本出来的图纵坐标就是时间，横坐标才是涡度的平均值。而你用set xlabs 强制改变了横坐标。set ylevs 这类命令对时间坐标不起作用。
你把这些设置坐标的东西先都注释掉，按默认出图，先把横纵坐标的真实意义都搞清楚了再说其他的

cici破晓碎寒

river 发表于 2018-4-9 21:28
是你自己搞错了。按照你的脚本出来的图纵坐标就是时间，横坐标才是涡度的平均值。而你用set xlabs 强制改变 ...

谢谢river大神的指导！！！！问题解决了