'sdfopen e:/data/grads/2011/uwnd.2011.nc'
'sdfopen e:/data/grads/2011/vwnd.2011.nc'
'set grads off'
'set lev 500'
'set lon 60 180'
'set lat 0 90'
'set gxout shaded'
'define var=mag(uwnd.1,vwnd.2)'
'set cmin 30'
'd var'
'set gxout vector'
'd skip(uwnd.1,2,2);skip(vwnd.2,2,2)'

'q dims'
i=sublin(result,2)
lon_s=subwrd(i,11)
lon_e=subwrd(i,13)
j=sublin(result,3)
lat_s=subwrd(j,11)
lat_e=subwrd(j,13)
say lon_s' 'lon_e' 'lat_s' 'lat_e;

***当然可以自己设定画轴线的区间****

lon1=lon_s 
while(lon1<=lon_e)
if(lon1=lon_s)
   'set x 'lon1''
   'd maxloc(var,y='lat_s',y='lat_e')'
   a=sublin(result,2)
   lat1=subwrd(a,4)
   'q gr2xy 'lon1' 'lat1''
   x0=subwrd(result,3)
   y0=subwrd(result,6)

else
   'set x 'lon1''
   'd maxloc(var,y='lat_s',y='lat_e')'
   a=sublin(result,2)
   lat1=subwrd(a,4)
   'q gr2xy 'lon1' 'lat1''
   x1=subwrd(result,3)
   y1=subwrd(result,6)
   'set line 3 1 7'
   'draw line 'x0' 'y0' 'x1' 'y1''
   x0=x1
   y0=y1

endif

lon1=lon1+1
endwhile

'printim e:/data/test/ceshi/exam1/windmax.png white'
;