'reinit'
'sdfopen hgt.2020.nc'
'sdfopen hgt.day.1981-2010.ltm.nc'
'set display color white'
'clear'
*
'set lat -50 50'
'set lon -50 75'
* 20200912 20200921 12Sep2020 21Sep2020 256 265
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 500'
'define h=ave(hgt,t='jdayb',t='jdaye')'
'set dfile 2'
'set lev 500'
'define hc=ave(hgt,t='jdayb',t='jdaye')'
*
'define han=h-hc'
*---------------------------
'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 5200 5300 5400 5500 5600 5700 5800'
'set ccols 0 34 32 21 22 23 24 25'
'd h'
'set gxout contour'
'set clevs 5200 5300 5400 5500 5600 5700 5800'
'set ccols 1 1 1 1 1 1 1 1'
'd h'
'run cbarhov2.gs'
'draw title 500hpa Mean Geop.Hgt (dyn mtrs) \ dekad 'ymdb2'-'ymde2''
'printim img/ra1_hgt500_dekad_'ymdb2'-'ymde2'.png'
*
'set mpdraw off'
'set geotiff img/ra1_hgt500_dekad_'ymdb2'-'ymde2''
'set gxout geotiff'
'set mproj latlon'
'd h'
*
*---------------------------
'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 -100 -80 -60 -40 -20 20 40 60 80 100'
'set ccols 36 35 34 33 32 0 21 22 23 24 25 26'
'd han'
'set gxout contour'
'set clevs -100 -80 -60 -40 -20 20 40 60 80 100'
'set ccols 1 1 1 1 1 1 1 1 1 1 1'
'd han'
'run cbarhov2.gs'
'draw title 500hpa GeopHgt Anom \ dekad 'ymdb2'-'ymde2''
'printim img/ra1_hgt500_anom_dekad_'ymdb2'-'ymde2'.png'
*
'set mpdraw off'
'set geotiff img/ra1_hgt500_anom_dekad_'ymdb2'-'ymde2''
'set gxout geotiff'
'set mproj latlon'
'd han'
*
'quit'