;************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"  
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"  
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
;************************************************

  a = addfile("/home/swati/work/analysis/2_cities/Patna/14aug2014/wrf_out_data/3rd_domain/patna_data.nc","r")

  wks = gsn_open_wks("png","partial_time_avg_plot")       ; Create a plot workstation

 t   = a->T2(0:7,:,:)
 ter = dim_avg_n( ter1, 0 ) 
 ter = ter - 273 
printVarSummary(ter)
;exit
 
  res = True                                 ; Use basic options for this field
  res@cnFillOn = True                        ; Create a color fill plot
 ; res@ContourParameters = (/ 1., 30., 2. /) ; Set the levels

  contour = wrf_contour(a,wks,ter,res)

  pltres = True                              ; Set plot options
  mpres = True                               ; Set map options
  plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)  ; Plot the data over a map background