;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
   a = addfile("/home/wrfuser/jsong/MERRA-2_data/MERRA2_400.inst3_3d_asm_Np.20160601.nc","r")  
;   a = addfile("/home/wrfuser/jsong/MERRA-2_data/MERRA2_400.inst3_3d_asm_Np.20160602.nc","r")  

   type = "ps"
   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)

   res = True
   res@tiMainString = "ozone mixing ratio  MERRA2 on June 01, 2016" 
;   res@trXMinF   =0
;   res@trXMaxF   =100
   res@trGridType = "TriangularMesh"

   res@cnFillOn = True

   res@mpMinLatF               = 30
   res@mpMaxLatF               = 45
   res@mpMinLonF               = -130
   res@mpMaxLonF               = -110

;   res@gsnPaperOrientation = "landscape"
;   res@vpWidthF    = 0.8        ; Make plots wider, default: 0.6
;   res@vpHeightF   = 0.4

 plot = gsn_contour_map(wks,o3_p,res)     ; Draw contours over a map.
    
end