"reinit"
"sdfopen Ave_Drought_IDX.nc"
IDX.1="PDSI"; IDX.2="scPDSI"; IDX.3="SPI01"; IDX.4="SPI03"; IDX.5="SPI06"; IDX.6="SPI12"
R.1="R1"; R.2="R2"; R.3="R3"; R.4="R4"; R.5="R5"; R.6="R6"; R.7="R7"
L1=1.0; L2=1.5; L3=2.0

NT="1 480"
Nidx=42

x00=1.0; x01=8; dy=1.2

Ndra=6
Nreg=7
id=0

While (id<Ndra)
  id=id+1
  idreg=(id-1)*Nreg+1
  iidx=idreg+Nreg
  "set display color white"
  "c"
  iy=0
  "set string 1 c 7"
  "set strsiz 0.12 0.15"
  "draw string 4.5 10.0 Areal mean of the "IDX.id" index"
  "set strsiz 0.11 0.12"
  "draw string 4.5 9.8 (Top to bottom: R1-R7 sub-regions)"
  "set xlab on"
  ireg=7
  While (iidx > idreg)
    iidx=iidx-1
    iy=iy+1
    if (iy=1);  y1=iy*dy ;  endif
    if (iy>1);  y1=y2;   "set xlab off";  endif
    y2=y1+dy
    "set parea "x00" "x01" "y1 " " y2
    "set grads off"; "set gxout linefill"; "set lfcols 4 2"
    if (id<=2); "set vrange -4 4"; else; "set vrange -3 3"; endif
    "set t "NT; "set y "iidx
    "define a1=0"; "d ave_idx; a1"
     draw_tit(x00+0.3, y2-0.2, R.ireg); ireg=ireg-1
*    if (id<=2); drline (L1*2, L2*2, L3*2); else; drline (L1, L2, L3); endif
  Endwhile
  "printim Pics/"IDX.id".png"
  pull dummy
Endwhile

"close 1"
*******************************************************
function drline (d1, d2, d3)
  "set gxout line"; "set line 1 2 3"; "set cmark 0"
  "define c1="d1; "define c2="d2; "define c3="d3

  "set line 1 2 3"; "set cmark 0"; "d c1"; 
  "set line 1 2 3"; "set cmark 0"; "d c2"; 
  "set line 1 2 3"; "set cmark 0"; "d c3"
  "set line 1 2 3"; "set cmark 0"; "d c1*(-1)"; 
  "set line 1 2 3"; "set cmark 0"; "d c2*(-1)"; 
  "set line 1 2 3"; "set cmark 0"; "d c3*(-1)"
  "set gxout linefill"
return
*******************************************************
function draw_tit(x, y, title)
  "set string 1 c 10"; "set strsiz 0.15 0.15"
  "draw string "x" "y" "title
return
*******************************************************