;; -*- mode: emacs-lisp; coding: emacs-mule; -*-;  -*- mode: ncl;-*-
;***************************************************************************
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"



ana=addfile("aerosol_RF.nc","r")

Absp_fine_no_ds=ana->Absp_fine_no_ds
Absp3D_fine_no_ds=ana->Absp3D_fine_no_ds
sfc_fine_no_ds=ana->sfc_fine_no_ds

lat=ana->lat
lon=ana->lon
lev=ana->lev

delete(ana)


Absp_fine_no_ds1=dim_avg_n_Wrap(Absp_fine_no_ds,0)
Absp3D_fine_no_ds1=dim_avg_n_Wrap(Absp3D_fine_no_ds,0)
sfc_fine_no_ds1=dim_avg_n_Wrap(sfc_fine_no_ds,0)



wks = gsn_open_wks("eps","atmospheric_heating_rate_1")

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
  res                     = True               ; plot mods desired
  res@cnFillOn            = True              ; turn on color   
  res@gsnSpreadColors     = True               ; use full color map
  res@cnLinesOn           = False              ; no contour lines
  res@cnLineLabelsOn      = False              ; no line labels

  res@mpMaxLatF =  87.8638                         ; specify the plot domain
  res@mpMinLatF = -87.8638                         ;                         
  res@mpMinLonF = 0                          ;
  res@mpMaxLonF = 357.1875     
  res@mpCenterLonF= 0;178.594 
  res@gsnAddCyclic        = False 

  res@cnLevelSelectionMode =  "AutomaticLevels"   ; AutomaticLevels


  ;res@cnLevelSelectionMode =  "ManualLevels"   ; AutomaticLevels
 ; res@cnMinLevelValF       = 0.              ; min level
 ; res@cnMaxLevelValF       =  2.              ; max level
 ; res@cnLevelSpacingF      =  0.1               ; interval

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
res@tiMainString        = " Atmospheric heating rate "

plot= gsn_csm_contour_map(wks,Absp_fine_no_ds1,res)
delete(wks)