begin

;************************************************
; variable and file handling
;************************************************
   June = "June@AOD.nc"
     Aero_Ju  = addfile(June,"r")
   aerosol= Aero_Ju->MOD08_M3_6_AOD_550_Dark_Target_Deep_Blue_Combined_Mean_Mean
   anamoly  =  dim_rmvmean_n_Wrap(aerosol,1)
  
  wks = gsn_open_wks("png","Anamoly@June")             ; send graphics to PNG file
  ;gsn_define_colormap(wks, "wgne15")
  
  res            = True                         ; plot mods desired
  res@gsnPolar   = "NH"                         ; specify the hemisphere
  res@mpMinLatF  = 60                           ; minimum lat to plot
  ;res@gsnMaximize          = True             ;  Maximize in frame
  res@gsnDraw              = True            ; Don't draw plots
  ;res@mpFontHeightF   = 0.025
  res@mpFillOn   = False
  res@cnInfoLabelOn        = False        ; Turn off info label
  res@cnLineLabelsOn       = False        ; Turn off line labels
  res@cnLinesOn            = False        ; Turn off contour lines
  
  
  res@cnFillOn          = True                  ; color fill
  ;res@cnLevelSpacingF   = 0.1                   ; interval spacing
  res@cnLinesOn         = False                 ; no contour lines
  res@cnFillPalette     = "ncl_default"
  res@lbOrientation        = "Vertical"   ; vertical label bar
  res@cnLevelSelectionMode = "ExplicitLevels"
  res@lbLabelFontHeightF   = 0.025
  ;res@lbLabelStrings = sprinti("%-5i",res)
  res@cnLevels       = (/-0.2,-0.15,-0.1,0, 0.05,0.1,0.15, 0.2/)
    
  res@lbTitleString  = "AOD"
  plot = gsn_csm_contour_map(wks,anamoly,res)
 
 
   
   
 end