; scatter_2.ncl
;
; Concepts illustrated:
;   - Drawing a scatter plot
;   - Changing the markers in an XY plot
;   - Changing the marker color in an XY plot
;   - Changing the marker size in an XY plot
;   - Rotating markers in an XY plot
;   - Creating your own markers for an XY plot
;   - Generating dummy data using "random_normal"
;   - Clipping markers outside the viewport
;
; ***********************************************
;load "/home/bliujuss/ncl/lib/ncarg/nclscripts/csm/gsn_code.ncl"
;load "/home/bliujuss/ncl/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
;load "/home/bliujuss/ncl/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;load "/home/bliujuss/ncl/lib/ncarg/nclscripts/csm/shea_util.ncl"
;load "/home/bliujuss/ncl/lib/ncarg/nclscripts/csm/contributed.ncl" 
;************************************************
begin

diri ="";
fils = systemfunc("ls "+diri+"*.nc");

f = addfiles(fils,"r")
ListSetType (f, "join") 

var1 = f[:]->MOD08_M3_6_Aerosol_Optical_Depth_Land_Ocean_Mean_Mean
printVarSummary(var1)
var1AVG = dim_avg_n_Wrap(var1,0)
printVarSummary(var1AVG)

wks = gsn_open_wks("png","Aerosol_OD_avg")

sres = True
sres@tiMainString = "Verification Data (Deep Blue)"
sres@cnFillOn = True
sres@gsnSpreadColors = True
sres@gsnAddCyclic = False
sres@gsnLeftString = ""
sres@gsnRightString = ""
sres@gsnCenterString = ""
sres@tiXAxisString = "Aerosol Optical Depth 550nm"
sres@cnLineLabelsOn = False
sres@cnLineLabelBackgroundColor = "white"

plot = gsn_csm_contour_map_ce(wks,var1AVG,sres)

end