* This script is used to plot spatial and timeseries of the seasonal distribution
"reinit"
Var.1=ZPalmer_1
Var.2=ZPalmer_2
Var.3=ZPalmer_3
Var.4=VplmBL_1
Var.5=VplmBL_2
Var.6=VplmBL_3

T.1="Spacial - Temporal Distribution of ZPalmer (1981-1990)"
T.2="Spacial - Temporal Distribution of ZPalmer (1991-2000)"
T.3="Spacial - Temporal Distribution of ZPalmer (2001-2011)"
T.4="Spacial - Temporal Distribution of Vplm (1981-1990)"
T.5="Spacial - Temporal Distribution of Vplm (1991-2000)"
T.6="Spacial - Temporal Distribution of Vplm (2001-2011)"
*_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"
_Xaxis=" 1981 1983 1985 1987 1989 1991 1993 1995 1997 1999 2001 2003 2005 2007 2009 2011"
**********************************************
  nvar=6
  iv=0
While (iv<nvar)
  iv=iv+1
  "open "Var.iv"_BL.ctl"
  "set display color white"
  "c"
  "set vpage 0.3 11 0.2 8.4"
   Set_Plt(iv)
*   Set_Color(iv)
   Set_My_Color() 
  "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"
    
    "set xlint 10"
*    "set ylint 5"
*    "set xlabs "_Xaxis
*    "set ylopts 1 4 0.08"
*    "set ylabs "_Yaxis.t
return
********************************************
function Set_My_Color()
lev1="  -4 -3 -2 -1 -0.5"
lev2=" 0.5 1 2 3 4"
levels = lev1%lev2
"run My_Color.gs -gxout grfill -kind red->white->blue -levs" levels
"set grads off"
return
********************************************