;*********************************************************
; Calculate area average of eke and eke for 30 years
;*********************************************************

  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"

  begin

  in                   = addfile("/data2/eke/control/ctrl.cam2.h0.1997-11-01-00000.nc","r")
  temp                 = in->T(0:29,1,:,:)
  temp_FillValue       = -999
  dimz                 = dimsizes(temp)
  nlat                 = 128
  mlon                 = 256
  nens                 = 2
  printVarSummary(temp)
  ta = temp(0,:,:)

;***********************************************************
; Read eke from  ascii files and store seasonal mean into 
; an array
;***********************************************************
  eke_ltw                 = new((/nens,nlat,mlon/),"float")
  eke_ltw@_FillValue      = -999
  eke_ltw!0               = "year"
  eke_ltw!1               = "lat"
  eke_ltw!2               = "lon" 

  t=0 
  do i= 1,2 
     eke_ltw(t,:,:)       = asciiread("/data2/6-hourly/ensemble_run/tracks/eke_test"+i,(/nlat,mlon/),"float")
     t=t+1
  end do

  n=0                             
  if (any(ismissing(eke_ltw)))then
  print("Your data is missing")
  n=n+1
  print(n)
  end if

  printVarSummary(eke_ltw)
  printMinMax(eke_ltw,True)
  
;***********************************************
; Calculating ensemble mean of u
;***********************************************

  eke_ens_av = dim_avg_n_Wrap(eke_ltw, 0)
  copy_VarCoords(ta,eke_ens_av)
  printVarSummary(eke_ens_av)
  printMinMax(eke_ens_av, True)

;*************************************************
;Plot
;*************************************************
  res                              = True    
  res@gsnMaximize                  = True
  res@gsnSpreadColors              = True     

  res@cnFillOn                     = True               
  res@cnLinesOn                    = False  
  res@cnLevelSelectionMode         = "ManualLevels"           
  res@cnMinLevelValF               = 0
  res@cnMaxLevelValF               = 1000
  res@cnLevelSpacingF              = 10
   
  res@lbLabelStride                = 10   
  res@lbOrientation                = "vertical"       ;vertical label bar

  res@tmLabelAutoStride            = True
  res@tmXTOn                       = False
  res@tmYROn                       = False
  
  res@mpLimitMode                  = "LatLon" 
  res@mpMinLatF                    = 10.
  res@mpMaxLatF                    = 80.
  res@mpMinLonF                    = 190.
  res@mpMaxLonF                    = 310.
  res@mpGeophysicalLineThicknessF  = 1.5

;*****************************************************
;               Plot
;*****************************************************

  wks = gsn_open_wks("pdf","eke_mean")        
  gsn_define_colormap(wks,"BlAqGrYeOrReVi200")  

  plot = gsn_csm_contour_map_ce(wks,eke_ens_av(:,:),res)

  end