; 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 ="/home/bliujuss/Plots/work/";
fils = systemfunc("ls "+diri+"*.nc");

f = addfiles(fils,"r")
print (f)
var1 = addfiles_GetVar(f,fils,"MOD08_M3_6_Aerosol_Optical_Depth_Land_Ocean_Mean_Mean")
printVarSummary(var1)
var1AVG = dim_avg_Wrap(var1)
printVarSummary(var1AVG)

a = addfile("/home/bliujuss/Plots/work/g4.timeAvgMap.MOD08_M3_6_Aerosol_Optical_Depth_Land_Ocean_Mean_Mean.20080701-20080731.180W_90S_180E_90N.nc","r")
b = addfile("/home/bliujuss/Plots/work/g4.timeAvgMap.MOD08_M3_6_Aerosol_Optical_Depth_Land_Ocean_Mean_Mean.20090701-20090731.180W_90S_180E_90N.nc","r")

var3 = a->MOD08_M3_6_Aerosol_Optical_Depth_Land_Ocean_Mean_Mean(:,:)
printVarSummary(var3)

var4 = b->MOD08_M3_6_Aerosol_Optical_Depth_Land_Ocean_Mean_Mean(:,:)
printVarSummary(var4)

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