[ncl-talk] Help: Errors when trying to specify xyMarkers, Marker Colors, and Marker Labels using gsn_csm_xy
Tyler Herrington
therring at uwaterloo.ca
Thu May 31 12:13:53 MDT 2018
Good afternoon,
I am attempting to create a scatterplot in which each data point has a
specified color and marker style:
; These files are loaded by default in NCL V6.2.0 and newer
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
;*********************
; Create model array
;*********************
mdl = new((/28,1/),"string")
mdl(0,0) = "ACCESS1-0"
mdl(1,0) = "ACCESS1-3"
mdl(2,0) = "bcc-csm1-1"
mdl(3,0) = "bcc-csm1-1-m"
mdl(4,0) = "BNU-ESM"
mdl(5,0) = "CanESM2"
mdl(6,0) = "CCSM4"
mdl(7,0) = "CNRM-CM5"
mdl(8,0) = "CSIRO-Mk3-6-0"
mdl(9,0) = "FGOALS-g2"
mdl(10,0) = "FGOALS-s2"
mdl(11,0) = "GFDL-CM3"
mdl(12,0) = "GFDL-ESM2G"
mdl(13,0) = "GFDL-ESM2M"
mdl(14,0) = "GISS-E2-H"
mdl(15,0) = "GISS-E2-R"
mdl(16,0) = "HadGEM2-ES"
mdl(17,0) = "inmcm4"
mdl(18,0) = "IPSL-CM5A-LR"
mdl(19,0) = "IPSL-CM5A-MR"
mdl(20,0) = "IPSL-CM5B-LR"
mdl(21,0) = "MIROC5"
mdl(22,0) = "MIROC-ESM"
mdl(23,0) = "MPI-ESM-LR"
mdl(24,0) = "MPI-ESM-MR"
mdl(25,0) = "MPI-ESM-P"
mdl(26,0) = "MRI-CGM3"
mdl(27,0) = "NorESM1-M"
;********************************************
; Read in Tundra and Boreal Albedo Feedbacks
;********************************************
TundirMAM = "/praid/users/herringtont/Caldwell_2016/monthly_data/
CLM/by_model/alb_fbk_maps/All.Tundra/fldmean/seasmean/MAM/"
BordirMAM = "/praid/users/herringtont/Caldwell_2016/monthly_data/
CLM/by_model/alb_fbk_maps/NET.Boreal/fldmean/seasmean/MAM/"
TunfilMAM = TundirMAM+"CMIP5_fbks_monres_MAM_AllModels_alb_fbk_maps_
All.Tundra_seasmean.nc
<http://cmip5_fbks_monres_mam_allmodels_alb_fbk_maps_all.tundra_seasmean.nc/>
"
BorfilMAM = BordirMAM+"CMIP5_fbks_monres_MAM_AllModels_alb_fbk_maps_
NET.Boreal_seasmean.nc
<http://cmip5_fbks_monres_mam_allmodels_alb_fbk_maps_net.boreal_seasmean.nc/>
"
TunMAM = addfile(TunfilMAM,"r")
BorMAM = addfile(BorfilMAM,"r")
TunMAM2 = TunMAM->alb_fbk_maps(:,0,0)
BorMAM2 = BorMAM->alb_fbk_maps(:,0,0)
ccr = escorc(BorMAM2,TunMAM2) ; calculate correlation coefficient
ccr2 = decimalPlaces(ccr,3,True)
;******************************************************
; Create MAM Seasonal Mean Albedo Feedback Scatterplot
;******************************************************
wks = gsn_open_wks("png","MAM_Albedo_Feedback_ScatterPlot_corr_labels")
; send graphics to PNG file
res2 = "True" ; plot mods desired
res2 at tmLabelAutoStride = True ; nice tick mark labels
res2 at xyMarkLineMode = "Markers"
res2 at xyMarkers = (/0,0,1,1,2,3,4,5,6,7,7,8,8,8,
9,9,10,11,12,12,12,13,13,14,14,14,15,16/)
res2 at xyMonoMarkerColor = "False" ; Allow different colors for
markers
res2 at xyMarkerColors = (/"dodgerblue2","dodgerblue4",
"firebrick4","tomato3","darkslategrey","forestgreen","
chocolate4","plum3","seagreen","cyan2","cyan4","
mediumorchid2","mediumorchid4","mediumpurple3","orange2","
orange4","antiquewhite4","violetred","salmon","salmon2",
"salmon4","blue","blue4","gold4","goldenrod2","yellow4",
"magenta4","darkolivegreen4"/)
res2 at xyExplicitLabels = (/mdl(0,0),mdl(1,0),mdl(2,0),
mdl(3,0),mdl(4,0),mdl(5,0),mdl(6,0),mdl(7,0),mdl(8,0),
mdl(9,0),mdl(10,0),mdl(11,0),mdl(12,0),mdl(13,0),mdl(14,0),
mdl(15,0),mdl(16,0),mdl(17,0),mdl(18,0),mdl(19,0),mdl(20,0),
mdl(21,0),mdl(22,0),mdl(23,0),mdl(24,0),mdl(25,0),mdl(26,0),mdl(27,0)/)
res2 at xyMarkerSizeF = 0.005 ; Marker size (default
0.01)
res2 at tiYAxisString = "Tundra Albedo Feedback (Wm~S~-2~N~K~S~-1~N~)"
res2 at tiXAxisString = "Boreal Albedo Feedback (Wm~S~-2~N~K~S~-1~N~)"
res2 at tiMainString = "MAM Mean Albedo Feedbacks"
res2 at trYMinF = 0.0
res2 at trYMaxF = 4.0
res2 at trXMinF = -0.75
res2 at trXMaxF = 1.5
res2 at gsnDraw = False ; do not draw
res2 at gsnFrame = False ; do not advance
txres = True
txres at txFontHeightF = 0.015 ; change font height
;---Set some legend resurces
res2 at pmLegendDisplayMode = "Always" ; Turn on the legend
res2 at lgOrientation = "horizontal" ; Default is vertical
res2 at pmLegendWidthF = 0.5 ; Make it wider
res2 at pmLegendOrthogonalPosF = -0.2 ; Move it up slightly
res2 at lgPerimOn = False ; Turn off the perimeter box
plot = gsn_csm_xy(wks,BorMAM2,TunMAM2,res2) ; create plot
draw(plot)
dum = gsn_add_text(wks,plot,"Pearson Correlation: "+ccr2,-0.2,3.8,txres)
frame(wks)
system("mv *.png /home/herringtont/Caldwell_2016/figures/scatterplot")
end
Upon running the script, however, I get the following error:
fatal:Argument type mismatch on argument (3) of (gsn_csm_xy) can not coerce
fatal:["Execute.c":8640]:Execute: Error occurred at or near line 107 in
file Boreal_Tundra_Scatterplot.ncl
Which leads me to suspect that something is wrong with the way that I am
specifying resources for the plot. Any suggestions on what could be causing
this?
Thank You,
Sincerely,
Tyler Herrington, MSc
PhD Student (Climate Science)
Geography and Environmental Management
University of Waterloo
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180531/0e7c3c65/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: CMIP5_fbks_monres_MAM_AllModels_alb_fbk_maps_All.Tundra_seasmean.nc
Type: application/octet-stream
Size: 6296 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180531/0e7c3c65/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: CMIP5_fbks_monres_MAM_AllModels_alb_fbk_maps_NET.Boreal_seasmean.nc
Type: application/octet-stream
Size: 6296 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180531/0e7c3c65/attachment-0001.obj>
More information about the ncl-talk
mailing list