* This script is used to plot spatial and timeseries of the seasonal distribution
"reinit"
Var.1=PDSI_S_Freq
Var.2=PDSI_T_Freq
Var.3=SPI_S_Freq_01
Var.4=SPI_T_Freq_01
Var.5=SPI_S_Freq_03
Var.6=SPI_T_Freq_03
Var.7=SPI_S_Freq_06
Var.8=SPI_T_Freq_06
Var.9=SPI_S_Freq_12
Var.10=SPI_T_Freq_12

T.1="Spacial Distribution of Drought Freq (PDSI) "
T.2="Temporal Distribution of Drought Freq (PDSI) "
T.3="Spacial Distribution of Drought Freq (SPI-01) "
T.4="Temporal Distribution of Drought Freq (SPI-01)"
T.5="Spacial Distribution of Drought Freq (SPI-03) "
T.6="Temporal Distribution of Drought Freq (SPI-03)"
T.7="Spacial Distribution of Drought Freq (SPI-06) "
T.8="Temporal Distribution of Drought Freq (SPI-06)"
T.9="Spacial Distribution of Drought Freq (SPI-12) "
T.10="Temporal Distribution of Drought Freq (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"
_Yaxis1="Ext_Dro|Sev_Dro|Mod_Dro|Mil_Dro|In_DSpe|Nea_Nor|In_WSpe|Sli_Wet|Mod_Wet|Ver_Wet|Ext_Wet"
_Yaxis2="Exc_Dry|Ext_Dry|Sev_Dry|Mod_Dry|Abn_Dry|Nea_Nor|Abn_Moi|Mod_Moi|Ver_Moi|Ext_Wet|Exc_Moi"
**********************************************
nvar=10
iv=0
While (iv<nvar)
  iv=iv+1
  "open "Var.iv"_Freq.ctl"
  "set display color white"
  "c"
  "set vpage 0.8 11 0.2 8.4"
  Set_Plt(iv)
  Set_Cr_Freq()
  "d "Var.iv
  "cbarn"
  "set vpage off"
  "set font 4"
  "set strsiz 0.2 0.2"
  "draw string 5.5 8.3 "T.iv 
* "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
*******************************************
function Set_Plt(t)
    "set grads off"
    "set font 4"
    "set mproj off"
    "set gxout grfill"
*    "set gxout shaded"
    if (t<=2)
       "set ylabs "_Yaxis1
    else
       "set ylabs "_Yaxis2
    endif
    if (t=2|t=4|t=6|t=8|t=10)
       "set xlabs "_Xaxis
    else 
       "set xlint 10"
    endif
*    "set ylint 5"
*    "set xlabs "_Xaxis
    "set ylopts 1 3 0.11"
*    "set ylabs "_Yaxis.t
return
********************************************
function Set_Cr_Freq()
lev1=" 0 0.05 0.1 0.15 0.2"
lev2=" 0.3 0.4 0.5 0.6 0.7 "
levels = lev1%lev2
"run My_Color.gs -gxout grfill -kind white->blue -levs" levels
"set grads off"
return
********************************************