* This script is used to plot spatial and timeseries of the seasonal distribution
"reinit"
Var.1=mop
Var.2=mow
Var.3=moz
Var.4=msp
Var.5=msw
Var.6=msz
Var.7=w1p
Var.8=w1w
Var.9=w1z
Var.10=w2p
Var.11=w2w
Var.12=w2z


T.1="Spacial - Temporal Distribution of ZPalmer "
T.2="Spacial - Temporal Distribution of Vplm "
T.3="Spacial - Temporal Distribution of SPI_01 "
T.4="Spacial - Temporal Distribution of SPI_03 "
T.5="Spacial - Temporal Distribution of SPI_06 "
T.6="Spacial - Temporal Distribution of SPI_12 "
_Xaxis=" 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012"
_Xaxis=" 1981 1983 1985 1987 1989 1991 1993 1995 1997 1999 2001 2003 2005 2007 2009 2011"
**********************************************
nvar=12
nvar2=nvar/2
nvar4=nvar/4
iv=0
While (iv<nvar)
  iv=iv+1
  if (iv = 1)
    "open "All_IDX_M.ctl
  else 
    if (iv = nvar2+1)
       "close 1" 
       "open "All_IDX_W1.ctl
    else
       if (iv = nvar2+nvar4+1)
          "close 1"
          "open "All_IDX_W2.ctl
       endif
    endif
  endif
    "set display color white"
    "c"
     "set vpage 0.3 11 0.2 8.4"
     Set_Plt(iv)
*     if (iv<=2)
         Set_Cr_PDSI()
*     else
*         Set_Cr_SPI()
*     endif
     "d "Var.iv
     "cbarn"
     "set vpage off"
     "set font 4"
     "set strsiz 0.2 0.2"
*     "draw string 5.5 8.3 "T.iv "" Period.iperi
*    "draw string 5.5 7.8 (Numbers in Y-axis: Stations in S->N direction)"
     "set font 0"
     "set strsiz 0.12 0.12"
*    "draw ylab No. of Station from S to N"
     "printim " Var.iv ".gif" 
say "Var=" iv " " Var.iv 
pull dum
*  "close 1"
Endwhile
"close 1"
*******************************************
function Set_Plt(t)
    "set grads off"
    "set font 4"
    "set mproj off"
    "set gxout shaded"
*    "set gxout shaded"
    
    "set xlint 2"
*    "set ylint 5"
    "set xlabs "_Xaxis
*    "set ylopts 1 4 0.08"
*    "set ylabs "_Yaxis.t
return
********************************************
function Set_Cr_PDSI()
lev1="  -4 -3 -2 -1 "
lev2="  1 2 3 4"
levels = lev1%lev2
"run My_Color.gs -gxout grfill -kind red->white->blue -levs" levels
"set grads off"
return
********************************************
function Set_Cr_SPI()
*lev1="  -2 -1.5 -1 "
*lev2=" 1 1.5 2 "
*levels = lev1%lev2
levels =" -2. -1.6 -1.3 -0.8 -0.5 0.5 0.8 1.3 1.6 2."
"run My_Color.gs -gxout grfill -kind red->white->blue -levs" levels
"set grads off"
return
********************************************