; 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"
;************************************************
begin
;************************************************
; Generate some random data
;************************************************

a = addfile("/home/bliujuss/Plots/g4.timeAvgMap.MOD08_M3_051_Optical_Depth_Land_And_Ocean_Mean_Mean.20150701-20150731.180W_90S_180E_90N.nc","r")
b = addfile("/home/bliujuss/Plots/FIM_Chem_2014112418.nc","r")  

var1 = a->MOD08_M3_051_Optical_Depth_Land_And_Ocean_Mean_Mean
var2 = b->VAR_147_GDS0_EATM
var1_reorder = var1(lat|:,lon|0:358) 
printVarSummary(var1_reorder)
;printVarSummary(var2)

;print(var2(120,199))
;print(var1_reorder(120,199))
 
bias = var1_reorder - var2
wks   = gsn_open_wks ("x11","xy")            ; open workstation
res                   = True                     ; plot mods desired
res@gsnMaximize       = True                     ; maximize plot
res@tiMainString      = "Two curve XY plot"
res@xyLineThickness  = (/1.0,2.0/)
res@xyLineColors   = (/"blue","red"/)      ; add title
res@xyMarkLineMode    = "MarkLines"                ; choose to use markers
res@xyMarkers         =  (/1,2/)                      ; choose type of marker  

print(dimsizes(var1_reorder))
print(dimsizes(var2))
data = new((/2,dimsizes(var1_reorder(120,:))/),float)
;print(data)
;print(dimsizes(data))
data(0,:) = var1_reorder(120,:)
data(1,:) = var2(120,:)
;print(var1_reorder(120,:))
plot  = gsn_csm_y (wks,data,res)                    ; create plot
plot = gsn_csm_y (wks,var2(90,:),res)
end