[ncl-talk] Dividing a legend over 2 cells in a matrix plot

Toni Klemm toni-klemm at tamu.edu
Wed Jul 25 17:22:48 MDT 2018


Howdy!

I forgot the script in my earlier email. Apologies!

I wrote an NCL script to plot 16 XY-plots on a 3x6 matrix. Each plot contains 9 lines, so the legend needs to have 9 elements. Because of this large number I wanted to divide the legend up into 2 “columns” (i.e., 2 matrix cells), 5 elements in the first and 4 in the second. I got everything to work, but unfortunately not the legend. The text is too large and the lines are dashed instead of solid. I tried to change the size of the text, but it didn’t do anything. Could someone help me? Could some other commands overrule the settings for the legends?

Thanks a lot for your help! Below is the script, and attached are the netCDF files in a folder and the matrix plot.

Best,
Toni


Toni Klemm, Ph.D.
Postdoctoral Research Associate
Department of Ecosystem Science and Management
Texas A&M University, College Station, TX
Contributor at the Early Career Climate Forum <http://www.eccforum.org/>
www.toni-klemm.de <http://www.toni-klemm.de/> | @toniklemm <http://twitter.com/toniklemm>


Script:


; ---------------------------------------------------

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"

; ---------------------------------------------------

begin

RCP_list = [/"rcp_45","rcp_85"/]
fs_list = [/"fs","nfs"/]
model_list = [/"CCSM4","GFDL-ESM2M","HadGEM2-ES365","IPSL-CM5A-LR","MRI-CGCM3"/]

RCP_title = [/"RCP 4.5","RCP 8.5"/]
fs_title = [/"Fire Suppr. off","Fire Suppr. on"/]
model_title = [/"CCSM4","GFDL-ESM2M","HadGEM2-ES365","IPSL-CM5A-LR","MRI-CGCM3"/]

do h = 0, ListCount(RCP_list) - 1 

	do i = 0, ListCount(fs_list) - 1 

		do j = 0, ListCount(model_list) - 1 


  ; ***********************************************
  ; 1 - READ IN MC2 FILE
  ; ***********************************************

	print("")
	print("~/NIFA/MC2/subset_climate_zones/" + RCP_list[h](0) + "/" + fs_list[i](0) + "/" + model_list[j](0))
  print("1 - Read in files")
	
	diri = "~/NIFA/MC2/subset_climate_zones/" + RCP_list[h](0) + "/" + fs_list[i](0) + "/" + model_list[j](0) + "/C_ABOVEGR_GROUPED/"

  zone_01_f = addfile(diri + RCP_list[h](0) + "_" + fs_list[i](0) + "_" + model_list[j](0) + "_C_ABOVEGR_grouped_zone_01.nc","r")
  zone_02_f = addfile(diri + RCP_list[h](0) + "_" + fs_list[i](0) + "_" + model_list[j](0) + "_C_ABOVEGR_grouped_zone_02.nc","r")
  zone_03_f = addfile(diri + RCP_list[h](0) + "_" + fs_list[i](0) + "_" + model_list[j](0) + "_C_ABOVEGR_grouped_zone_03.nc","r")
  zone_04_f = addfile(diri + RCP_list[h](0) + "_" + fs_list[i](0) + "_" + model_list[j](0) + "_C_ABOVEGR_grouped_zone_04.nc","r")
  zone_05_f = addfile(diri + RCP_list[h](0) + "_" + fs_list[i](0) + "_" + model_list[j](0) + "_C_ABOVEGR_grouped_zone_05.nc","r")
  zone_06_f = addfile(diri + RCP_list[h](0) + "_" + fs_list[i](0) + "_" + model_list[j](0) + "_C_ABOVEGR_grouped_zone_06.nc","r")
  zone_07_f = addfile(diri + RCP_list[h](0) + "_" + fs_list[i](0) + "_" + model_list[j](0) + "_C_ABOVEGR_grouped_zone_07.nc","r")
  zone_08_f = addfile(diri + RCP_list[h](0) + "_" + fs_list[i](0) + "_" + model_list[j](0) + "_C_ABOVEGR_grouped_zone_08.nc","r")
  zone_09_f = addfile(diri + RCP_list[h](0) + "_" + fs_list[i](0) + "_" + model_list[j](0) + "_C_ABOVEGR_grouped_zone_09.nc","r")
  zone_10_f = addfile(diri + RCP_list[h](0) + "_" + fs_list[i](0) + "_" + model_list[j](0) + "_C_ABOVEGR_grouped_zone_10.nc","r")
  zone_11_f = addfile(diri + RCP_list[h](0) + "_" + fs_list[i](0) + "_" + model_list[j](0) + "_C_ABOVEGR_grouped_zone_11.nc","r")
  zone_12_f = addfile(diri + RCP_list[h](0) + "_" + fs_list[i](0) + "_" + model_list[j](0) + "_C_ABOVEGR_grouped_zone_12.nc","r")
  zone_13_f = addfile(diri + RCP_list[h](0) + "_" + fs_list[i](0) + "_" + model_list[j](0) + "_C_ABOVEGR_grouped_zone_13.nc","r")
  zone_14_f = addfile(diri + RCP_list[h](0) + "_" + fs_list[i](0) + "_" + model_list[j](0) + "_C_ABOVEGR_grouped_zone_14.nc","r")
  zone_15_f = addfile(diri + RCP_list[h](0) + "_" + fs_list[i](0) + "_" + model_list[j](0) + "_C_ABOVEGR_grouped_zone_15.nc","r")
  zone_16_f = addfile(diri + RCP_list[h](0) + "_" + fs_list[i](0) + "_" + model_list[j](0) + "_C_ABOVEGR_grouped_zone_16.nc","r")


  ; ***********************************************
  ; 2 - CREATE MATRIX PLOT
  ; ***********************************************

  print("2 - Create Plot Matrix")

	plot = new(18,graphic)                ; create 17-plot array
  
  diro = "~/NIFA/MC2/subset_climate_zones/time_series/"
  filo = RCP_list[h](0) + "_" + fs_list[i](0) + "_" + model_list[j](0) + "_C_ABOVEGR_percentage_zone_01-16_group_01-09"

  wks_type = "png"											; plot file type
  wks_type at wkWidth = 3500								; width dimension
  wks_type at wkHeight = 3500 							; height dimension
  wks = gsn_open_wks(wks_type, diro + filo)
  gsn_define_colormap(wks,"wgne15")     ; define colormap


	;*********************************************
	; PLOT MATRIX PARAMETERS
	;********************************************* 

	res               = True
  res at trYMinF       = 0					; min  y axis value
  res at trYMaxF       = 100			 	; max  y axis value
  res at trXMinF       = 2015    	; min  y axis value
  res at trXMaxF       = 2100    	; max  x axis value
  res at tmXBMinorOn   = True 			; turn off minor tickmarks
  res at tmYLMinorOn   = True
  res at tmXTOn        = False 		; turn off top tickmarks
  res at tmYROn        = False 		; turn off right tickmarks
  res at tmXTBorderOn  = False 		; turn off top axis
  res at tmYRBorderOn  = False 		; turn off right axis

  res at tiXAxisString = "Years" ; titles
  res at tiYAxisString = "%"


  res at tmYLLabelFontHeightF       = 0.03 ; set font heights
  res at gsnRightStringFontHeightF  = 0.03
  res at tiXAxisFontHeightF         = 0.03
  res at tiYAxisFontHeightF         = 0.03

  res at vpWidthF       = .8  ; change aspect ratio
  res at vpHeightF      = .4

  res at gsnDraw       = False ; don't draw or advance the frame
  res at gsnFrame      = False 

	res at xyDashPattern  = 0													; Make curves all solid
	res at xyLineColors := (/"blue","red","orange","black","green","pink","brown","purple","grey"/) ; 9 different colors
	res at xyLineThicknesses = (/10,10,10,10,10,10,10,10,10/)		; make second line thicker

;	res at xyMarkLineMode = "MarkLines"								; Markers *and* lines
;	res at xyMarkers = (/0,2,4,6,7,8,9,14,15/)					; 9 different markers
;	res at xyMarkerColors := (/"blue","red","yellow","black","green","pink","brown","purple","grey"/) ; 9 different colors
;	res at xyMarkerThicknesses = (/10,10,10,10,10,10,10,10,10/)		; make second line thicker


	;*********************************************
	; READ IN PLOT VARIABLES
	;********************************************* 

  print("3 - Read in plot variables")

	year = zone_01_f->year

	units = "%"

	data_01 = new((/9,dimsizes(year)/),float)		; data array for 9 variables and 85 years of record (2015-2099)
	data_02 = new((/9,dimsizes(year)/),float)		; data array for 9 variables and 85 years of record (2015-2099)
	data_03 = new((/9,dimsizes(year)/),float)		; data array for 9 variables and 85 years of record (2015-2099)
	data_04 = new((/9,dimsizes(year)/),float)		; data array for 9 variables and 85 years of record (2015-2099)
	data_05 = new((/9,dimsizes(year)/),float)		; data array for 9 variables and 85 years of record (2015-2099)
	data_06 = new((/9,dimsizes(year)/),float)		; data array for 9 variables and 85 years of record (2015-2099)
	data_07 = new((/9,dimsizes(year)/),float)		; data array for 9 variables and 85 years of record (2015-2099)
	data_08 = new((/9,dimsizes(year)/),float)		; data array for 9 variables and 85 years of record (2015-2099)
	data_09 = new((/9,dimsizes(year)/),float)		; data array for 9 variables and 85 years of record (2015-2099)
	data_10 = new((/9,dimsizes(year)/),float)		; data array for 9 variables and 85 years of record (2015-2099)
	data_11 = new((/9,dimsizes(year)/),float)		; data array for 9 variables and 85 years of record (2015-2099)
	data_12 = new((/9,dimsizes(year)/),float)		; data array for 9 variables and 85 years of record (2015-2099)
	data_13 = new((/9,dimsizes(year)/),float)		; data array for 9 variables and 85 years of record (2015-2099)
	data_14 = new((/9,dimsizes(year)/),float)		; data array for 9 variables and 85 years of record (2015-2099)
	data_15 = new((/9,dimsizes(year)/),float)		; data array for 9 variables and 85 years of record (2015-2099)
	data_16 = new((/9,dimsizes(year)/),float)		; data array for 9 variables and 85 years of record (2015-2099)
	

	print("3a   - Zone 1")

	data_01(0,:) = zone_01_f->C_ABOVEGR_C3G_percentage
	data_01(1,:) = zone_01_f->C_ABOVEGR_C3S_percentage
	data_01(2,:) = zone_01_f->C_ABOVEGR_C4G_percentage
	data_01(3,:) = zone_01_f->C_ABOVEGR_C4S_percentage
	data_01(4,:) = zone_01_f->C_ABOVEGR_OTHER_percentage
	data_01(5,:) = zone_01_f->C_ABOVEGR_SF_percentage
	data_01(6,:) = zone_01_f->C_ABOVEGR_SW_percentage
	data_01(7,:) = zone_01_f->C_ABOVEGR_TF_percentage
	data_01(8,:) = zone_01_f->C_ABOVEGR_TW_percentage
	ALL_data_01	 = zone_01_f->C_ABOVEGR_ALL_annual

	print("3a   - Zone 2")

	data_02(0,:) = zone_02_f->C_ABOVEGR_C3G_percentage
	data_02(1,:) = zone_02_f->C_ABOVEGR_C3S_percentage
	data_02(2,:) = zone_02_f->C_ABOVEGR_C4G_percentage
	data_02(3,:) = zone_02_f->C_ABOVEGR_C4S_percentage
	data_02(4,:) = zone_02_f->C_ABOVEGR_OTHER_percentage
	data_02(5,:) = zone_02_f->C_ABOVEGR_SF_percentage
	data_02(6,:) = zone_02_f->C_ABOVEGR_SW_percentage
	data_02(7,:) = zone_02_f->C_ABOVEGR_TF_percentage
	data_02(8,:) = zone_02_f->C_ABOVEGR_TW_percentage
	ALL_data_02	 = zone_02_f->C_ABOVEGR_ALL_annual

	print("3a   - Zone 3")

	data_03(0,:) = zone_03_f->C_ABOVEGR_C3G_percentage
	data_03(1,:) = zone_03_f->C_ABOVEGR_C3S_percentage
	data_03(2,:) = zone_03_f->C_ABOVEGR_C4G_percentage
	data_03(3,:) = zone_03_f->C_ABOVEGR_C4S_percentage
	data_03(4,:) = zone_03_f->C_ABOVEGR_OTHER_percentage
	data_03(5,:) = zone_03_f->C_ABOVEGR_SF_percentage
	data_03(6,:) = zone_03_f->C_ABOVEGR_SW_percentage
	data_03(7,:) = zone_03_f->C_ABOVEGR_TF_percentage
	data_03(8,:) = zone_03_f->C_ABOVEGR_TW_percentage
	ALL_data_03	 = zone_03_f->C_ABOVEGR_ALL_annual

	print("3a   - Zone 4")

	data_04(0,:) = zone_04_f->C_ABOVEGR_C3G_percentage
	data_04(1,:) = zone_04_f->C_ABOVEGR_C3S_percentage
	data_04(2,:) = zone_04_f->C_ABOVEGR_C4G_percentage
	data_04(3,:) = zone_04_f->C_ABOVEGR_C4S_percentage
	data_04(4,:) = zone_04_f->C_ABOVEGR_OTHER_percentage
	data_04(5,:) = zone_04_f->C_ABOVEGR_SF_percentage
	data_04(6,:) = zone_04_f->C_ABOVEGR_SW_percentage
	data_04(7,:) = zone_04_f->C_ABOVEGR_TF_percentage
	data_04(8,:) = zone_04_f->C_ABOVEGR_TW_percentage
	ALL_data_04	 = zone_04_f->C_ABOVEGR_ALL_annual

	print("3a   - Zone 5")

	data_05(0,:) = zone_05_f->C_ABOVEGR_C3G_percentage
	data_05(1,:) = zone_05_f->C_ABOVEGR_C3S_percentage
	data_05(2,:) = zone_05_f->C_ABOVEGR_C4G_percentage
	data_05(3,:) = zone_05_f->C_ABOVEGR_C4S_percentage
	data_05(4,:) = zone_05_f->C_ABOVEGR_OTHER_percentage
	data_05(5,:) = zone_05_f->C_ABOVEGR_SF_percentage
	data_05(6,:) = zone_05_f->C_ABOVEGR_SW_percentage
	data_05(7,:) = zone_05_f->C_ABOVEGR_TF_percentage
	data_05(8,:) = zone_05_f->C_ABOVEGR_TW_percentage
	ALL_data_05	 = zone_05_f->C_ABOVEGR_ALL_annual

	print("3a   - Zone 6")

	data_06(0,:) = zone_06_f->C_ABOVEGR_C3G_percentage
	data_06(1,:) = zone_06_f->C_ABOVEGR_C3S_percentage
	data_06(2,:) = zone_06_f->C_ABOVEGR_C4G_percentage
	data_06(3,:) = zone_06_f->C_ABOVEGR_C4S_percentage
	data_06(4,:) = zone_06_f->C_ABOVEGR_OTHER_percentage
	data_06(5,:) = zone_06_f->C_ABOVEGR_SF_percentage
	data_06(6,:) = zone_06_f->C_ABOVEGR_SW_percentage
	data_06(7,:) = zone_06_f->C_ABOVEGR_TF_percentage
	data_06(8,:) = zone_06_f->C_ABOVEGR_TW_percentage
	ALL_data_06	 = zone_06_f->C_ABOVEGR_ALL_annual

	print("3a   - Zone 7")

	data_07(0,:) = zone_07_f->C_ABOVEGR_C3G_percentage
	data_07(1,:) = zone_07_f->C_ABOVEGR_C3S_percentage
	data_07(2,:) = zone_07_f->C_ABOVEGR_C4G_percentage
	data_07(3,:) = zone_07_f->C_ABOVEGR_C4S_percentage
	data_07(4,:) = zone_07_f->C_ABOVEGR_OTHER_percentage
	data_07(5,:) = zone_07_f->C_ABOVEGR_SF_percentage
	data_07(6,:) = zone_07_f->C_ABOVEGR_SW_percentage
	data_07(7,:) = zone_07_f->C_ABOVEGR_TF_percentage
	data_07(8,:) = zone_07_f->C_ABOVEGR_TW_percentage
	ALL_data_07	 = zone_07_f->C_ABOVEGR_ALL_annual

	print("3a   - Zone 8")

	data_08(0,:) = zone_08_f->C_ABOVEGR_C3G_percentage
	data_08(1,:) = zone_08_f->C_ABOVEGR_C3S_percentage
	data_08(2,:) = zone_08_f->C_ABOVEGR_C4G_percentage
	data_08(3,:) = zone_08_f->C_ABOVEGR_C4S_percentage
	data_08(4,:) = zone_08_f->C_ABOVEGR_OTHER_percentage
	data_08(5,:) = zone_08_f->C_ABOVEGR_SF_percentage
	data_08(6,:) = zone_08_f->C_ABOVEGR_SW_percentage
	data_08(7,:) = zone_08_f->C_ABOVEGR_TF_percentage
	data_08(8,:) = zone_08_f->C_ABOVEGR_TW_percentage
	ALL_data_08	 = zone_08_f->C_ABOVEGR_ALL_annual

	print("3a   - Zone 9")

	data_09(0,:) = zone_09_f->C_ABOVEGR_C3G_percentage
	data_09(1,:) = zone_09_f->C_ABOVEGR_C3S_percentage
	data_09(2,:) = zone_09_f->C_ABOVEGR_C4G_percentage
	data_09(3,:) = zone_09_f->C_ABOVEGR_C4S_percentage
	data_09(4,:) = zone_09_f->C_ABOVEGR_OTHER_percentage
	data_09(5,:) = zone_09_f->C_ABOVEGR_SF_percentage
	data_09(6,:) = zone_09_f->C_ABOVEGR_SW_percentage
	data_09(7,:) = zone_09_f->C_ABOVEGR_TF_percentage
	data_09(8,:) = zone_09_f->C_ABOVEGR_TW_percentage
	ALL_data_09	 = zone_09_f->C_ABOVEGR_ALL_annual

	print("3a   - Zone 10")

	data_10(0,:) = zone_10_f->C_ABOVEGR_C3G_percentage
	data_10(1,:) = zone_10_f->C_ABOVEGR_C3S_percentage
	data_10(2,:) = zone_10_f->C_ABOVEGR_C4G_percentage
	data_10(3,:) = zone_10_f->C_ABOVEGR_C4S_percentage
	data_10(4,:) = zone_10_f->C_ABOVEGR_OTHER_percentage
	data_10(5,:) = zone_10_f->C_ABOVEGR_SF_percentage
	data_10(6,:) = zone_10_f->C_ABOVEGR_SW_percentage
	data_10(7,:) = zone_10_f->C_ABOVEGR_TF_percentage
	data_10(8,:) = zone_10_f->C_ABOVEGR_TW_percentage
	ALL_data_10	 = zone_10_f->C_ABOVEGR_ALL_annual

	print("3a   - Zone 11")

	data_11(0,:) = zone_11_f->C_ABOVEGR_C3G_percentage
	data_11(1,:) = zone_11_f->C_ABOVEGR_C3S_percentage
	data_11(2,:) = zone_11_f->C_ABOVEGR_C4G_percentage
	data_11(3,:) = zone_11_f->C_ABOVEGR_C4S_percentage
	data_11(4,:) = zone_11_f->C_ABOVEGR_OTHER_percentage
	data_11(5,:) = zone_11_f->C_ABOVEGR_SF_percentage
	data_11(6,:) = zone_11_f->C_ABOVEGR_SW_percentage
	data_11(7,:) = zone_11_f->C_ABOVEGR_TF_percentage
	data_11(8,:) = zone_11_f->C_ABOVEGR_TW_percentage
	ALL_data_11	 = zone_11_f->C_ABOVEGR_ALL_annual

	print("3a   - Zone 12")

	data_12(0,:) = zone_12_f->C_ABOVEGR_C3G_percentage
	data_12(1,:) = zone_12_f->C_ABOVEGR_C3S_percentage
	data_12(2,:) = zone_12_f->C_ABOVEGR_C4G_percentage
	data_12(3,:) = zone_12_f->C_ABOVEGR_C4S_percentage
	data_12(4,:) = zone_12_f->C_ABOVEGR_OTHER_percentage
	data_12(5,:) = zone_12_f->C_ABOVEGR_SF_percentage
	data_12(6,:) = zone_12_f->C_ABOVEGR_SW_percentage
	data_12(7,:) = zone_12_f->C_ABOVEGR_TF_percentage
	data_12(8,:) = zone_12_f->C_ABOVEGR_TW_percentage
	ALL_data_12	 = zone_12_f->C_ABOVEGR_ALL_annual

	print("3a   - Zone 13")

	data_13(0,:) = zone_13_f->C_ABOVEGR_C3G_percentage
	data_13(1,:) = zone_13_f->C_ABOVEGR_C3S_percentage
	data_13(2,:) = zone_13_f->C_ABOVEGR_C4G_percentage
	data_13(3,:) = zone_13_f->C_ABOVEGR_C4S_percentage
	data_13(4,:) = zone_13_f->C_ABOVEGR_OTHER_percentage
	data_13(5,:) = zone_13_f->C_ABOVEGR_SF_percentage
	data_13(6,:) = zone_13_f->C_ABOVEGR_SW_percentage
	data_13(7,:) = zone_13_f->C_ABOVEGR_TF_percentage
	data_13(8,:) = zone_13_f->C_ABOVEGR_TW_percentage
	ALL_data_13	 = zone_13_f->C_ABOVEGR_ALL_annual

	print("3a   - Zone 14")

	data_14(0,:) = zone_14_f->C_ABOVEGR_C3G_percentage
	data_14(1,:) = zone_14_f->C_ABOVEGR_C3S_percentage
	data_14(2,:) = zone_14_f->C_ABOVEGR_C4G_percentage
	data_14(3,:) = zone_14_f->C_ABOVEGR_C4S_percentage
	data_14(4,:) = zone_14_f->C_ABOVEGR_OTHER_percentage
	data_14(5,:) = zone_14_f->C_ABOVEGR_SF_percentage
	data_14(6,:) = zone_14_f->C_ABOVEGR_SW_percentage
	data_14(7,:) = zone_14_f->C_ABOVEGR_TF_percentage
	data_14(8,:) = zone_14_f->C_ABOVEGR_TW_percentage
	ALL_data_14	 = zone_14_f->C_ABOVEGR_ALL_annual

	print("3a   - Zone 15")

	data_15(0,:) = zone_15_f->C_ABOVEGR_C3G_percentage
	data_15(1,:) = zone_15_f->C_ABOVEGR_C3S_percentage
	data_15(2,:) = zone_15_f->C_ABOVEGR_C4G_percentage
	data_15(3,:) = zone_15_f->C_ABOVEGR_C4S_percentage
	data_15(4,:) = zone_15_f->C_ABOVEGR_OTHER_percentage
	data_15(5,:) = zone_15_f->C_ABOVEGR_SF_percentage
	data_15(6,:) = zone_15_f->C_ABOVEGR_SW_percentage
	data_15(7,:) = zone_15_f->C_ABOVEGR_TF_percentage
	data_15(8,:) = zone_15_f->C_ABOVEGR_TW_percentage
	ALL_data_15	 = zone_15_f->C_ABOVEGR_ALL_annual

	print("3a   - Zone 16")

	data_16(0,:) = zone_16_f->C_ABOVEGR_C3G_percentage
	data_16(1,:) = zone_16_f->C_ABOVEGR_C3S_percentage
	data_16(2,:) = zone_16_f->C_ABOVEGR_C4G_percentage
	data_16(3,:) = zone_16_f->C_ABOVEGR_C4S_percentage
	data_16(4,:) = zone_16_f->C_ABOVEGR_OTHER_percentage
	data_16(5,:) = zone_16_f->C_ABOVEGR_SF_percentage
	data_16(6,:) = zone_16_f->C_ABOVEGR_SW_percentage
	data_16(7,:) = zone_16_f->C_ABOVEGR_TF_percentage
	data_16(8,:) = zone_16_f->C_ABOVEGR_TW_percentage
	ALL_data_16	 = zone_16_f->C_ABOVEGR_ALL_annual
	

	;*********************************************************
	; CREATE INDIVIDUAL PLOTS
	;*********************************************************

	print("4 - Create individual plots")

	print("4a   - Zone 1")

  res at gsnRightString = "Zone 1" ; title
  res at gsnLeftString  = "Avg.: "  + decimalPlaces(dim_avg_n(ALL_data_01,0)/1000000,1,True) + " t C m^-2"   ; printed on the top right
	plot(0) = gsn_csm_xy (wks,year,data_01,res) ; create plot
	
	print("4a   - Zone 2")

  res at gsnRightString = "Zone 2" ; title
  res at gsnLeftString  = "Avg.: "  + decimalPlaces(dim_avg_n(ALL_data_02,0)/1000000,1,True) + " t C m^-2"   ; printed on the top right
	plot(1) = gsn_csm_xy (wks,year,data_02,res) ; create plot
	
	print("4a   - Zone 3")

  res at gsnRightString = "Zone 3" ; title
  res at gsnLeftString  = "Avg.: "  + decimalPlaces(dim_avg_n(ALL_data_03,0)/1000000,1,True) + " t C m^-2"   ; printed on the top right
	plot(2) = gsn_csm_xy (wks,year,data_03,res) ; create plot
	
	print("4a   - Zone 4")

  res at gsnRightString = "Zone 4" ; title
  res at gsnLeftString  = "Avg.: "  + decimalPlaces(dim_avg_n(ALL_data_04,0)/1000000,1,True) + " t C m^-2"   ; printed on the top right
	plot(3) = gsn_csm_xy (wks,year,data_04,res) ; create plot
	
	print("4a   - Zone 5")

  res at gsnRightString = "Zone 5" ; title
  res at gsnLeftString  = "Avg.: "  + decimalPlaces(dim_avg_n(ALL_data_05,0)/1000000,1,True) + " t C m^-2"   ; printed on the top right
	plot(4) = gsn_csm_xy (wks,year,data_05,res) ; create plot
	
	print("4a   - Zone 6")

  res at gsnRightString = "Zone 6" ; title
  res at gsnLeftString  = "Avg.: "  + decimalPlaces(dim_avg_n(ALL_data_06,0)/1000000,1,True) + " t C m^-2"   ; printed on the top right
	plot(5) = gsn_csm_xy (wks,year,data_06,res) ; create plot
	
	print("4a   - Zone 7")

  res at gsnRightString = "Zone 7" ; title
  res at gsnLeftString  = "Avg.: "  + decimalPlaces(dim_avg_n(ALL_data_07,0)/1000000,1,True) + " t C m^-2"   ; printed on the top right
	plot(6) = gsn_csm_xy (wks,year,data_07,res) ; create plot
	
	print("4a   - Zone 8")

  res at gsnRightString = "Zone 8" ; title
  res at gsnLeftString  = "Avg.: "  + decimalPlaces(dim_avg_n(ALL_data_08,0)/1000000,1,True) + " t C m^-2"   ; printed on the top right
	plot(7) = gsn_csm_xy (wks,year,data_08,res) ; create plot
	
	print("4a   - Zone 9")

  res at gsnRightString = "Zone 9" ; title
  res at gsnLeftString  = "Avg.: "  + decimalPlaces(dim_avg_n(ALL_data_09,0)/1000000,1,True) + " t C m^-2"   ; printed on the top right
	plot(8) = gsn_csm_xy (wks,year,data_09,res) ; create plot
	
	print("4a   - Zone 10")

  res at gsnRightString = "Zone 10" ; title
  res at gsnLeftString  = "Avg.: "  + decimalPlaces(dim_avg_n(ALL_data_10,0)/1000000,1,True) + " t C m^-2"   ; printed on the top right
	plot(9) = gsn_csm_xy (wks,year,data_10,res) ; create plot
	
	print("4a   - Zone 11")

  res at gsnRightString = "Zone 11" ; title
  res at gsnLeftString  = "Avg.: "  + decimalPlaces(dim_avg_n(ALL_data_11,0)/1000000,1,True) + " t C m^-2"   ; printed on the top right
	plot(10) = gsn_csm_xy (wks,year,data_11,res) ; create plot
	
	print("4a   - Zone 12")

  res at gsnRightString = "Zone 12" ; title
  res at gsnLeftString  = "Avg.: "  + decimalPlaces(dim_avg_n(ALL_data_12,0)/1000000,1,True) + " t C m^-2"   ; printed on the top right
	plot(11) = gsn_csm_xy (wks,year,data_12,res) ; create plot
	
	print("4a   - Zone 13")

  res at gsnRightString = "Zone 13" ; title
  res at gsnLeftString  = "Avg.: "  + decimalPlaces(dim_avg_n(ALL_data_13,0)/1000000,1,True) + " t C m^-2"   ; printed on the top right
	plot(12) = gsn_csm_xy (wks,year,data_13,res) ; create plot
	
	print("4a   - Zone 14")

  res at gsnRightString = "Zone 14" ; title
  res at gsnLeftString  = "Avg.: "  + decimalPlaces(dim_avg_n(ALL_data_14,0)/1000000,1,True) + " t C m^-2"   ; printed on the top right
	plot(13) = gsn_csm_xy (wks,year,data_14,res) ; create plot
	
	print("4a   - Zone 15")

  res at gsnRightString = "Zone 15" ; title
  res at gsnLeftString  = "Avg.: "  + decimalPlaces(dim_avg_n(ALL_data_15,0)/1000000,1,True) + " t C m^-2"   ; printed on the top right
	plot(14) = gsn_csm_xy (wks,year,data_15,res) ; create plot
	
	print("4a   - Zone 16")

  res at gsnRightString = "Zone 16" ; title
  res at gsnLeftString  = "Avg.: "  + decimalPlaces(dim_avg_n(ALL_data_16,0)/1000000,1,True) + " t C m^-2"   ; printed on the top right
	plot(15) = gsn_csm_xy (wks,year,data_16,res) ; create plot
	

;***********************************************
; CREATING A LEGEND
;***********************************************

	print("5 - Create legend")
	
;	var_long_name = (/"C3 grassland","C3 shrubland","C4 grassland","C4 shrubland","Other","Subtropical forest","Subtropical woodland","Temperate forest","Temperate woodland"/)

	var_long_name_01 = (/"C3 grassland","C3 shrubland","C4 grassland","C4 shrubland","Other"/)
	var_long_name_02 = (/"Subtropical forest","Subtropical woodland","Temperate forest","Temperate woodland"/)

; line colors: "blue","red","orange","black","green","pink","brown","purple","grey"

	; Legend elements 1-5 	
  lgres                    = True
  lgres at lgAutoManage			 = False
  lgres at lgLineColors     	 = (/"blue","red","orange","black","green"/)
  lgres at lgLineThicknessF   = res at xyLineThicknesses  ; legend line thickness = contour line thickness
  lgres at lgItemType         = "Lines"        				; show lines only (default)
  lgres at lgItemCount				 = 5											; Number of items
  lgres at lgLabelFontHeightF = .04            				; set the legend label font thickness
  lgres at vpWidthF           = 1		           				; width of legend (NDC)
  lgres at vpHeightF          = 1    	        				; height of legend (NDC)
  lgres at lgPerimColor       = "white"	       				; draw the box perimeter in orange
  lgres at lgPerimThicknessF  = 1.0            				; thicken the box perimeter
  plot(16) = gsn_create_legend(wks,5,var_long_name_01,lgres)         ; create legend
  
  delete(lgres)
  
	; Legend elements 6-9 
  lgres                    = True
  lgres at lgAutoManage			 = False
  lgres at lgLineColors     	 = (/"pink","brown","purple","grey"/)
  lgres at lgLineThicknessF   = res at xyLineThicknesses  ; legend line thickness = 1/2 contour line thickness
  lgres at lgItemType         = "Lines"        				; show lines only (default)
  lgres at lgItemCount				 = 4											; Number of items
  lgres at lgLabelFontHeightF = .04            				; set the legend label font thickness
  lgres at vpWidthF           = 1		           				; width of legend (NDC)
  lgres at vpHeightF          = 1    	        				; height of legend (NDC)
  lgres at lgPerimColor       = "white"	       				; draw the box perimeter in orange
  lgres at lgPerimThicknessF  = 1.0            				; thicken the box perimeter
  plot(17) = gsn_create_legend(wks,4,var_long_name_02,lgres)         ; create legend

/;
	amres = True
  amres at amParallelPosF   = 0		  								; move legend to the right
  amres at amOrthogonalPosF = 0         			  	  ; move the legend down
  annoid1 = gsn_add_annotation(plot(16),lbid,amres)   	; add legend to plot
;/

;***********************************************
; panel plots together
;***********************************************

	print("6 - Create plot matrix")

;	print((/decimalPlaces(min(data),0,True)/))
;	print((/decimalPlaces(max(data),0,True)/))

  pres                 = True
  pres at gsnMaximize     = False
  pres at gsnPanelRowSpec = True
  pres at gsnPanelCenter  = False
  pres at gsnFrame        = False
;  pres at gsnPanelXWhiteSpacePercent	= 10	; horizontal white space between columns in percent
;  pres at gsnPanelYWhiteSpacePercent	= 5		; vertical white space between rows in percent
  pres at gsnPanelMainString = RCP_title[h](0) + " | " + fs_title[i](0) + " | " + model_title[j](0) + " | Above ground C per func. group"
  pres at tiMainFontHeightF  = 0.04
  pres at tiMainFont         = 21
;  pres at gsnPanelXF  = (/-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,0.6/)

  gsn_panel(wks,plot,(/3,3,3,3,3,3/),pres)
  
  frame(wks)
  
  print("")

	print("**** DONE ****")
	
	delete([/diri,plot,diro,filo,wks,res,data_01,data_02,data_03,data_04,data_05,data_06,data_07,data_08,data_09,data_10,data_11,data_12,data_13,data_14,data_15,data_16,zone_01_f,zone_02_f,zone_03_f,zone_04_f,zone_05_f,zone_06_f,zone_07_f,zone_08_f,zone_09_f,zone_10_f,zone_11_f,zone_12_f,zone_13_f,zone_14_f,zone_15_f,zone_16_f,pres,lgres/])

	sleep(10)
	
		end do	; end model list loop

	end do	; end fs list loop

end do 	; end RCP list loop

print("**** ALL DONE ****")

end




-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180725/d37f302f/attachment-0003.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: C_ABOVEGR_GROUPED.zip
Type: application/zip
Size: 107697 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180725/d37f302f/attachment-0001.zip>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180725/d37f302f/attachment-0004.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: rcp_45_fs_CCSM4_C_ABOVEGR_percentage_zone_01-16_group_01-09.png
Type: image/png
Size: 1275823 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180725/d37f302f/attachment-0001.png>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180725/d37f302f/attachment-0005.html>


More information about the ncl-talk mailing list