'open data/soilw.14mon.mean.ctl'
'sdfopen data/soilw.clim.ctl'
*
'set display color white'
'clear'
*
'set lat -50 50'
'set lon -50 75'
* *'set lev 500'
file='monyear'
  result=read(file)
  rstr=sublin(result,2)
  tm=subwrd(rstr,1)
  mon=subwrd(rstr,2)
  year=subwrd(rstr,3)
*---------------------------
'set dfile 1'
'define sw=soilw(time='mon''year')'
'set dfile 2'
'define swc=soilw(t='tm')'
'define swan=sw-swc'
*--------------------------------------------
'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 200 300 400 500 600 700 800'
'set ccols 0 33 34 36 38 43 45 46 48'
'd sw'
'run cbarhov2.gs'
'draw title Mean Soil Moisture (mm)   'mon''year''
'printim cpc_soilmoisture_'mon''year'.png'
*
'set mpdraw off'
'set geotiff cpc_soilw_'mon''year''
'set gxout geotiff'
'set mproj latlon'
'd sw'
*
*---------------------------
'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 -500 -400 -300 -200 -100 -50 50 100 200 300 400 500'
'set ccols 78 77 76 75 74 72 0 32 33 35 36 37 39'
'd swan'
'run cbarhov2.gs'
'draw title Soil Moisture Anom  (mm)  'mon''year''
'printim cpc_soilmoisture_anom_'mon''year'.png'
*
*
'set mpdraw off'
'set geotiff cpc_soilw_anom_'mon''year''
'set gxout geotiff'
'set mproj latlon'
'd swan'
*
'quit'