;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.20160530.nc","r")  

;   a = addfile("/home/wrfuser/jsong/MERRA-2_data/MERRA2_400.inst6_3d_ana_Np.20160622.nc","r")  
;   a = addfile("/home/wrfuser/jsong/MERRA-2_data/MERRA2_400.inst3_3d_asm_Np.20160622.nc","r")  

   type = "ps"
   wks = gsn_open_wks(type, "plt_o3_xy_MERRA2_asm_May30") 

;   o3= a->O3(3,0:22,257,91) ; 18GTM = 11AM PST, every 3-hr data
   o3= a->O3(3,4,:,:) ; 18GTM = 11AM PST
printVarSummary(o3)
;   o3_p =(o3*1e9)/2.14  ; ozone rho=2.14kg/m3
   o3_p =(o3*1e9)*(28.97/48) 

   res = True
   res@tiMainString = "ozone mixing ratio horizontal distribution MERRA2 on May 30, 2016" 
;   res@trXMinF   =0
;   res@trXMaxF   =100

   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.
;plot = gsn_contour_map(wks,o3_p(240:260, 80:100),res)     ; Draw contours over a map.
    
end