'reinit'
'sdfopen rhum.2020.nc'
'sdfopen rhum.day.1981-2010.ltm.nc'
'set display color white'
'clear'
*
'set lat -50 50'
'set lon -50 75'
* *'set lev 850'
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 850'
'define rh=ave(rhum,t='jdayb',t='jdaye')'
'set dfile 2'
'set lev 850'
'define rhc=ave(rhum,t='jdayb',t='jdaye')'
*
'define rhan=rh-rhc'
*--------------------------------------------
'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 10 20 30 40 50 60 70 80 90'
'set ccols 0 32 34 21 22 23 24 25 26 27'
'd rh'
'set gxout contour'
'set clevs 10 20 30 40 50 60 70 80 90'
'set ccols 1 1 1 1 1 1 1 1 1 1 1'
'd rh'
'run cbarhov2.gs'
'draw title 850hpa Mean Relative Humidity (in%) \ dekad 'ymdb2'-'ymde2''
'printim img/ra1_rh850_dekad_'ymdb2'-'ymde2'.png'
*
'set mpdraw off'
'set geotiff img/ra1_rh850_dekad_'ymdb2'-'ymde2''
'set gxout geotiff'
'set mproj latlon'
'd rh'
'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 -40 -30 -20 -10 -5 5 10 20 30 40'
'set ccols 36 35 34 33 32 0 21 22 23 24 25 26'
'd rhan'
'set gxout contour'
'set clevs -40 -30 -20 -10 -5 5 10 20 30 40'
'set ccols 1 1 1 1 1 1 1 1 1 1 1'
'd rhan'
'run cbarhov2.gs'
'draw title 850hpa Relative Humidity Anom (in%) \ dekad 'ymdb2'-'ymde2''
'printim img/ra1_rh850_anom_dekad_'ymdb2'-'ymde2'.png'
*
'set mpdraw off'
'set geotiff img/ra1_rh850_anom_dekad_'ymdb2'-'ymde2''
'set gxout geotiff'
'set mproj latlon'
'd rhan'
*
'quit'