;*********************************
; eosdis_4.ncl
;*********************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;*********************************
begin 
;*********************************
; read in HDF file
;*********************************
  f = addfile("eos.hdf","r")
;*************************************************************
; read in first channel in EV_250_Aggr1km_RefSB (0.65 um)
;*************************************************************
  ch1_int  = f->EV_250_Aggr1km_RefSB(0,2::5,2::5)
; orig is integer, make float
  ch1      = ch1_int*ch1_int@reflectance_scales(0) 

  lat2d    = f->Latitude
  lon2d    = f->Longitude

  ch1@lon2d  = lon2d
  ch1@lat2d  = lat2d
;*************************************************************
; assign named dimensions
;*************************************************************
  ch1@long_name = "VISIBLE REFLECTANCE (0.6 um)"
  ch1@units     = "unitless" 
;*************************************************************
; create plot
;*************************************************************
  wks = gsn_open_wks("ps", "eosdis") 
  gsn_define_colormap(wks,"gsdtol")         ; choose colormap
  
  res                  = True
  res@cnFillOn         = True               ; color Fill 
  res@cnFillMode       = "RasterFill"       ; Raster Mode
  res@cnLinesOn        = False              ; Turn off contour lines
  
  res@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels
  res@cnMinLevelValF   =  0.0               ; set min contour level
  res@cnMaxLevelValF   =  1.0               ; set max contour level
  res@cnLevelSpacingF  =  0.04              ; set contour spacing
  
  res@gsnSpreadColors     = True            ; use full colormap
  res@gsnSpreadColorStart = 11              ; start at color 11
  res@gsnSpreadColorEnd   = 32              ; end at color 32
  res@lbLabelStride       = 3               ; every third label bar label
  
  
  res@mpMinLatF           = min(lat2d)-5      ; choose map limits
  res@mpMaxLatF           = max(lat2d)+5
  res@mpMinLonF           = min(lon2d)-5
  res@mpMaxLonF           = max(lon2d)+5
  
  res@mpCenterLonF        = 180             ; change map center
  res@mpDataBaseVersion   = "Ncarg4_1"      ; higher res coastline
  res@gsnAddCyclic        = False           ; needed b/c data not global
  
  res@tiMainString        = "eos-dis data on a ce projection"
  
  plot = gsn_csm_contour_map_ce(wks,ch1, res) ; create plot
  
end