'reinit'
'sdfopen uwnd.2020.nc'
'sdfopen vwnd.2020.nc'
'sdfopen uwnd.day.1981-2010.ltm.nc'
'sdfopen vwnd.day.1981-2010.ltm.nc'
*
'set display color white'
'clear'
*
'set lat -50 50'
'set lon -50 75'
*---------------------------
file='dates.txt'
  result=read(file)
  rstr=sublin(result,2)
  ymdb1=subwrd(rstr,1)
  ymde1=subwrd(rstr,2)
  ymdb2=subwrd(rstr,3)
  ymde2=subwrd(rstr,4)
  jdayb=subwrd(rstr,5)
  jdaye=subwrd(rstr,6)
*---------------------------
'set dfile 1'
'set lev 925'
'define u=ave(uwnd,t='jdayb',t='jdaye')'
'set dfile 2'
'set lev 925'
'define v=ave(vwnd,t='jdayb',t='jdaye')'
'set dfile 3'
'set lev 925'
'define uc=ave(uwnd,t='jdayb',t='jdaye')'
'set dfile 4'
'set lev 925'
'define vc=ave(vwnd,t='jdayb',t='jdaye')'
*
'define uan=u-uc'
'define van=v-vc'
*---------------------------
'exec colors'
'set vpage 1 10 0 8.5'
'set mproj scaled'
'set grid off'
'set grads off'
'set map 1 1 6'
'set mpdset hires'
'set gxout shaded'
'set clevs 1 2 4 6 8 10 15 20 25 30'
'set ccols 0 41 43 31 33 35 21 22 23 24 25'
'd mag(u,v)'
'set gxout contour'
'set clevs 1 2 4 6 8 10 15 20 25 30'
'set ccols 1 1 1 1 1 1 1 1 1 1 1'
'd mag(u,v)'
'set gxout vector'
'set arrlab off'
'set ccolor 1'
'd u;v'
'run cbarhov2.gs'
'draw title 925hpa mean wind (m/s)  \ dekad 'ymdb2'-'ymde2''
'printim img/ra1_wind925_dekad_'ymdb2'-'ymde2'.png'
*
*
'set mpdraw off'
'set geotiff img/ra1_wind925_dekad_'ymdb2'-'ymde2''
'set gxout geotiff'
'set mproj latlon'
'd mag(uan,van)'
*
*---------------------------
'clear'
'set mpdraw on'
'exec colors'
'set vpage 1 10 0 8.5'
'set mproj scaled'
'set grid off'
'set grads off'
'set map 1 1 6'
'set mpdset hires'
'set gxout shaded'
'set clevs 1 2 3 4'
'set ccols 0 21 22 23 24'
'd mag(uan,van)'
'set gxout contour'
'set clevs 1 2 3 4'
'set ccols 1 1 1 1 1'
'd mag(uan,van)'
'set gxout vector'
'set arrlab off'
'set ccolor 1'
'd uan;van'
'run cbarhov2.gs'
'draw title 925hpa wind anom (m/s)  \ dekad 'ymdb2'-'ymde2''
'printim img/ra1_wind925_anom_dekad_'ymdb2'-'ymde2'.png'
*
*
'set mpdraw off'
'set geotiff img/ra1_wind925_anom_dekad_'ymdb2'-'ymde2''
'set gxout geotiff'
'set mproj latlon'
'd mag(uan,van)'
*
'quit'