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/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "/WRF_myfunc.ncl"
load "/Shapefile/shapefile_utils.ncl"
begin


a=addfile("../../BEM_April","r")
b=addfile("../../../April_2016_without_SST/NEW_BEM_SST301_April_2016","r")


ta_mean_BEM=a->TA_MEAN
ta=ta_mean_BEM(24:695,0,:,:)
tb_mean_BEM=b->TA_MEAN
tb=tb_mean_BEM(24:695,0,:,:)
LH= ta - tb

time = a->Times
times = time(24:695,:)
ntimes = dimsizes(times)         ; number of times in the file
ntimes_plot = ntimes(0)


  lh3d = new((/28,129,210/), float)
  lh2d = new((/129,210/), float)
    c = 0
    do i=time_serie, ntimes_plot-2,24
       lh2d(:,:) = LH(i,:,:)
       lh3d(c,:,:) = lh2d(:,:)
       c = c + 1

      end do
    lh_avg = dim_avg_n(lh3d,0)


 type="png"
 wks = gsn_open_wks(type,"testplot")

 res                          = True
 res@gsnMaximize              = True
 res@gsnDraw                  = False
 res@gsnFrame                 = False
 res@cnFillOn                 = True
 res@cnLinesOn                = False
 res@tmYROn                   = False
 res@lbOrientation            = "Vertical"
 res@lbLabelFontHeightF       = 0.010
 res@lbLabelFontThicknessF    =1.5
 ;res@lbBoxMajorExtentF        = 0.75
 res@lbBoxEndCapStyle         = "TriangleBothEnds"
 cmap_data=read_colormap_file("wgne15")
 res@cnFillPalette            = cmap_data
 contour_lh                   = gsn_csm_contour(wks,lh_avg,res)
 draw(contour_lh)
 frame(wks)

end