*************************************************
ifile = "PDSI_Drought_Characteristics_Reg.nc"


*figout = "DGeo1_Mid-Future.pdf"
*var = "DGeo1"
*cmip.1 = "MF45"
*cmip.2 = "MF85"

figout = "DGeo1_Far-Future.pdf"
var = "DGeo1"
cmip.1 = "FF45"
cmip.2 = "FF85"


**************************************************
** vrange can be changed at the end of this script
**************************************************


**************************************************
** labels of each region
reg.1 = "R1"
reg.2 = "R2"
reg.3 = "R3"
reg.4 = "R4"
reg.5 = "R5"
reg.6 = "R6"
reg.7 = "R7"
reg.8 = "VN"

**************************************************
** Color and style of marks for each model
cmdl.1 = 2 ; smdl.1 = 2
cmdl.2 = 3 ; smdl.2 = 3
cmdl.3 = 4 ; smdl.3 = 4
cmdl.4 = 5 ; smdl.4 = 5
cmdl.5 = 6 ; smdl.5 = 6
cmdl.6 = 7 ; smdl.6 = 7
cmdl.7 = 8 ; smdl.7 = 8

**************************************************
** outter margin and number of columns and rows
omr = "0.3 0.8 0.5 0 -space 0.05"
nrow = 4
ncol = 4

*************************************************
'reinit'
'sdfopen 'ifile
initialize ()

ir = 1 ; ireg = 1
while (ir <= nrow)
  ic = 1
  while (ic <= ncol)
    im = 1
    'set t 'ireg

    xlb = "off" ; if (ir = nrow) ; xlb = "on" ; endif
    ylb = "off" ; if (ic = 1)    ; ylb = "on" ; endif
    'AIO_plot -nr 'nrow' -nc 'ncol' -ir 'ir' -ic 'ic' -xlb 'xlb' -ylb 'ylb' -omr 'omr
    while (im <= 7)
      'set y 'im
      'set cstyle 0'
      'set ccolor 'cmdl.im
      'set cmark 'smdl.im
      if (ir <= 2)
        'd 'var'_'cmip.1
        if (ic = 1) ; 'nx_title -main 'cmip.1' -pos left -space 0.5 -siz 0.16 -rot 90' ; endif
      else
        'd 'var'_'cmip.2
        if (ic = 1) ; 'nx_title -main 'cmip.2' -pos left -space 0.5 -siz 0.16 -rot 90' ; endif
      endif
      
      'nx_title -main 'reg.ireg' -space -0.18 -siz 0.16'
      im = im + 1
    endwhile

    ic = ic + 1
    if (ireg = 8); ireg = 0 ; endif
    ireg = ireg + 1
  endwhile
  ir = ir + 1
endwhile
'draw string 5.82353 8.2576 'var
'gxprint 'figout


function initialize ()
  'set display color white' ; 'c'
  'set xlint 1' ; 'set ylint 5'
  'set x 0 13'
  'set xlabs |J|F|M|A|M|J|J|A|S|O|N|D|'
  'set grads off'
  'set vrange -18 18'
  'set grid off'
return