;script to plot the O3 profile of MERRA-2 data at location BBY.
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin
   dir = "/home/wrfuser/jsong/MERRA-2_data/"
   dir = "./"
   fname = "MERRA2_400.inst3_3d_asm_Np.20160601.nc"
   a = addfile(dir + fname,"r")

   type = "x11"
   wks = gsn_open_wks(type, "plt_o3_xy_MERRA2_asm_Jun01_zoom") 

;   o3= a->O3(3,0:22,257,91) ; 18GTM = 11AM PST, every 3-hr data
   o3_p= a->O3(3,4,:,:) ; 18GTM = 11AM PST
   
   o3_p@long_name = "ozone mixing ratio"
   o3_p@units     = "ppbv"

   printVarSummary(o3_p)
   o3_p =(o3_p*1e9)*(28.97/48) 

   printVarSummary(o3_p)
   printMinMax(o3_p,0)

   res = True
   res@gsnMaximize    = True   ; maximize size of plot

   res@cnFillOn       = True
   res@cnLinesOn      = False
   res@cnLineLabelsOn = False
   res@gsnAddCyclic   = True    ; add longitude cyclic point, mostly 
                                ; only useful if plotting global data
   res@mpMinLatF      = 30
   res@mpMaxLatF      = 45
   res@mpMinLonF      = -130
   res@mpMaxLonF      = -110

   res@tiMainString   = "ozone mixing ratio  MERRA2 on June 01, 2016" 
   res@pmTitleZone    = 4   ; move title down
   res@pmTickMarkDisplayMode = "Always"    ; nicer tickmarks

   plot = gsn_csm_contour_map(wks,o3_p,res)
    
end