"reinit"
"sdfopen Ave_Drought_IDX_Stn.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"

NT="1 480"
Nidx=42
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"
  "set strsiz 0.12 0.15"
  "draw string 4.5 10.0 Areal mean of the "IDX.id" index (Station data)"
  "set strsiz 0.11 0.12"
  "draw string 4.5 9.8 (Top to bottom: R1-R7 sub-regions)"

  "set strsiz 0.10 0.10"
  While (iidx > idreg)
    "set xlab on"
    "set xaxis 1966 2006 5"
    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 1 8.0 "y1 " " y2

    "set grads off"
    if (iy=1)
    endif
    "set gxout linefill"
    "set lfcols 4 2"
    "set vrange -4 4"
    "set x "iidx
    "d ave_idx; ave_idx-ave_idx"

    "set vpage off"
*say "iidx="iidx
  Endwhile
  "printim Pics/"IDX.id"_Stn.png"
  pull dummy
Endwhile

"close 1"