;***********Load the ncl internal libraries*********;
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"

begin

;***get data***
regcm_ctr_slp= addfile("/media/Data/REGCM_SST/atm_files/2001/control/daymean/CTR-ATM-2001-mslp-7-daymean.nc","r")
regcm_sens_slp_plus= addfile("/media/Data/REGCM_SST/atm_files/2001/sens+2/daymean/CTR-ATM-2001-mslp-7-daymean.nc","r")
regcm_ctr_u = addfile("/home/denizbozkurt/Desktop/DATA/REGCM_SST/2001/daymean/CTR-SRF-2001-uas-7-daymean.nc","r")
regcm_ctr_v = addfile("/home/denizbozkurt/Desktop/DATA/REGCM_SST/2001/daymean/CTR-SRF-2001-vas-7-daymean.nc","r")
regcm_sens_u_plus = addfile("/home/denizbozkurt/Desktop/DATA/REGCM_SST/sens+2/2001/daymean/CTR-SRF-2001-uas-7-daymean.nc","r")
regcm_sens_v_plus = addfile("/home/denizbozkurt/Desktop/DATA/REGCM_SST/sens+2/2001/daymean/CTR-SRF-2001-vas-7-daymean.nc","r")
regcm_domain = addfile("/home/denizbozkurt/Desktop/DATA/REGCM_SST/CTR_DOMAIN000.nc","r")

;***select domain***
latsouth = -58
latnorth = -4
lonwest = -122
loneast = -48

;***select variables***
selday=18
regcm_ctr_slp_day = regcm_ctr_slp->mslp(time|selday,iy|:,jx|:)
regcm_sens_slp_plus_day = regcm_sens_slp_plus->mslp(time|selday,iy|:,jx|:)
regcm_sens_slp_plus_diff = regcm_sens_slp_plus_day - regcm_ctr_slp_day
regcm_ctr_u_day = regcm_ctr_u->uas(time|selday,iy|:,jx|:)
regcm_ctr_v_day = regcm_ctr_v->vas(time|selday,iy|:,jx|:)
regcm_sens_u_plus_day = regcm_sens_u_plus->uas(time|selday,iy|:,jx|:)
regcm_sens_v_plus_day = regcm_sens_v_plus->vas(time|selday,iy|:,jx|:)
regcm_sens_u_plus_diff = regcm_sens_u_plus_day - regcm_ctr_u_day
regcm_sens_v_plus_diff = regcm_sens_v_plus_day - regcm_ctr_v_day

;***lat&lon***
lat2d=regcm_domain->xlat(iy|:,jx|:)
lon2d=regcm_domain->xlon(iy|:,jx|:)
nlat = dimsizes(lat2d(:,0))
nlon = dimsizes(lon2d(0,:))


wks=gsn_open_wks("ps","synoptic_slp_wind_diff_NCL-5.2.1")
gsn_define_colormap (wks,"BlueDarkRed18")
res                        = True               
res@cnFillOn               = True               
res@cnLinesOn             = False
res@cnLineLabelsOn         = False              
res@cnInfoLabelOn          = False              
res@lbLabelBarOn         = False
res@mpFillOn = False
res@gsnRightString = ""
res@cnLevelSelectionMode = "ExplicitLevels"
res@cnLevels = (/-5,-4.5,-4,-3.5,-3,-2.5,-2,-1.5,-1,-0.5,0,0.5,1,1.5,2,2.5,3,3.5,4,4.5,5/)
res@gsnSpreadColors = True
res@tmXBLabelFontHeightF  = 0.014
res@tmYLLabelFontHeightF  = 0.014
;res@gsnCenterStringFontHeightF = 0.024
res@gsnLeftString = ""
res@tmXTOn = False
res@tmYROn = False
res@tiMainString           = ""
res@tiMainFontHeightF      = 0.016           
res@gsnDraw              = False             
res@gsnFrame             = False
res@gsnAddCyclic           = False           
res@mpPerimOn = True
res@mpDataBaseVersion = "MediumRes"
res@pmTickMarkDisplayMode  = "Always"        


res@vcRefMagnitudeF         = 3.0            
res@vcRefLengthF            = 0.03           
res@vcGlyphStyle            = "CurlyVector"
res@vcMinDistanceF          = 0.025
res@vcLineArrowColor        = "Black"
res@vcLineArrowThicknessF    = 1.0
res@vcRefAnnoOrthogonalPosF  = -1.0
res@vcVectorDrawOrder        = "Postdraw"
res@gsnScalarContour = True

res@mpMinLatF = latsouth
res@mpMaxLatF = latnorth
res@mpMinLonF = lonwest
res@mpMaxLonF = loneast

res@tfDoNDCOverlay = True

res@tfDoNDCOverlay = True
res@mpLimitMode       = "Corners"          
res@mpLeftCornerLatF  = lat2d(0,0)
res@mpLeftCornerLonF  = lon2d(0,0)
res@mpRightCornerLatF = lat2d(nlon-1,nlat-1)
res@mpRightCornerLonF = lon2d(nlon-1,nlat-1)
prj = regcm_domain@projection
clon = regcm_domain@longitude_of_projection_origin
clat = regcm_domain@latitude_of_projection_origin
res@mpProjection        = "Mercator"
res@mpCenterLonF = clon
res@mpCenterLatF = clat

plot=gsn_csm_vector_scalar_map(wks,regcm_sens_u_plus_diff,regcm_sens_v_plus_diff,regcm_sens_slp_plus_diff,res)


;box
ypts = (/ -30,-30,-38,-38,-30/)
xpts = (/ -74,-70,-70,-74,-74/)
resl = True
resl@gsLineColor = "black"
resl@gsLineThicknessF = 2.0
resl@gsLineLabelString = ""


dum = new(4,graphic)
do i = 0 , 3
dum(i)=gsn_add_polyline(wks,plot,xpts(i:i+1),ypts(i:i+1),resl)
end do


resP = True
resP@gsnPanelLabelBar = True
resP@lbOrientation = "Vertical"
resP@lbLabelAutoStride = True
resP@pmLabelBarWidthF = 0.065
resP@pmLabelBarHeightF = 0.7
resP@lbTitleOn        =  True                
resP@lbTitleString    = "hPa"                
resP@lbTitlePosition  = "Bottom"             
resP@lbTitleFontHeightF= .011              
resP@lbTitleDirection = "Across"             
resP@lbLabelStride        = 2
resP@pmLabelBarParallelPosF = -0.02
resP@lbLabelFontHeightF  = 0.012
resP@pmLabelBarOrthogonalPosF = 0.02
resP@txString   = ""
gsn_panel(wks,plot,(/1,1/),resP)

end