;; -*- 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("sfc_toa_fine.nc","r")
Absp= ana->Absp_fine_no_ds1
sfc= ana->sfc_fine_no_ds1
delete(ana)


ana=addfile("Absp3D_fine_no_ds_zero.nc","r")
lev3=ana->lev
Absp3D_fine_no_ds3=ana->Absp3D_fine_no_ds
delete(ana)

Absp3D_fine_no_ds4=dim_avg_n_Wrap(Absp3D_fine_no_ds3,0)
Absp3D_fine_no_ds5=dim_avg_n_Wrap(Absp3D_fine_no_ds4,(/1,2/))*24*3600




Absp3D_summer= Absp3D_fine_no_ds3(5:7,:,:,:)
Absp3D_spring= Absp3D_fine_no_ds3(2:4,:,:,:)
Absp3D_fall= Absp3D_fine_no_ds3(8:10,:,:,:)
Absp3D_winter= new((/3,25,64,128/),typeof(Absp3D_fine_no_ds3))
Absp3D_winter(0:1,:,:,:)= Absp3D_fine_no_ds3(0:1,:,:,:)
Absp3D_winter(2,:,:,:)= Absp3D_fine_no_ds3(11,:,:,:)




Absp3D_summer1= dim_avg_n_Wrap(Absp3D_summer,0)
Absp3D_summer2= dim_avg_n_Wrap(Absp3D_summer1,(/1,2/))*24*3600
Absp3D_spring1= dim_avg_n_Wrap(Absp3D_spring,0)
Absp3D_spring2= dim_avg_n_Wrap(Absp3D_spring1,(/1,2/))*24*3600
Absp3D_fall1= dim_avg_n_Wrap(Absp3D_fall,0)
Absp3D_fall2= dim_avg_n_Wrap(Absp3D_fall1,(/1,2/))*24*3600
Absp3D_winter1= dim_avg_n_Wrap(Absp3D_winter,0)
Absp3D_winter2= dim_avg_n_Wrap(Absp3D_winter1,(/1,2/))*24*3600





wks = gsn_open_wks("eps","atm_sfc_fine")
;gsn_define_colormap(wks,"BlAqGrYeOrRe")
;gsn_define_colormap(wks,"BlWhRe")
 plot = new(3,graphic)   
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
  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@gsnAddCyclic        = False 
  res@gsnDraw             = False           ; don't draw
  res@gsnFrame            = False           ; don't advance frame
  res@cnInfoLabelOn       = False           ; turn off cn info label
  res@cnFillOn            = True            ; turn on color
  res@gsnSpreadColors     = True            ; spread out color table
  res@gsnSpreadColorStart = 2               ; color to startt

  res@lbLabelBarOn        = True          ; turn off individual cb's

;  res@cnLevelSelectionMode =  "AutomaticLevels"   ; AutomaticLevels


  res@cnLevelSelectionMode =  "ManualLevels"   ; AutomaticLevels
  res@cnMinLevelValF       = -40.              ; min level
  res@cnMaxLevelValF       =  40.               ; max level
  res@cnLevelSpacingF      =  4               ; interval

;plot = new(3,graphic)
res@gsnLeftString       = "Solar radiation absorption in the atmosphere"         ; change left string
plot(0) = gsn_csm_contour_map_ce(wks,Absp,res)
res@gsnLeftString       = "Surface net flux"         ; change left string
plot(1) = gsn_csm_contour_map_ce(wks,sfc,res)

res2                   = True                       ; plot mods desired
res2@gsnDraw = False
res2@gsnFrame = False
;res2@tiMainString      = "Vertical profile of atmospheric heating rate"             ; add title
res2@tiYAxisString = "Atmospheric heating rate (K/day)"  ; xaxis string
res2@tiXAxisString = "Pressure (hPa)";Pressure levels (hPa)"

res2@trYMinF     = 100		; set minimum Y-axis value
res2@trYMaxF     =  1000		; set maximum Y-axis value
;res@trXMinF     = -0.5      ; set minimum X-axis value
;res@trXMaxF     = 1.2      ; set maximum X-axis value
 res2@trYReverse = True                     ; reverse Y-axis
res2@xyDashPattern  = 0                         ; Make curves all solid
res2@xyLineColors      = "blue"          ; change line color
res2@xyLineThicknesses = 4.0
 plot(2)  = gsn_csm_xy (wks,Absp3D_summer2,lev3,res) ; create plot

;************************************************
; create panel
;************************************************
  resP                  = True                    ; modify the panel plot

  resP@gsnPanelLabelBar = False                  ; add common colorbar
  resP@gsnPanelFigureStrings = (/"a)","b)","c)"/)      ; add strings to panel
gsn_panel(wks,plot,(/2,2/),resP)