! This PyFerret script plots FIG_S5 of the Paper "Artificial Upwelling - A New Narrative" 
! and was written by M. Juerchott (mjuerchott@geomar.de) on the 12th September 2022.
! In order to run this script without error you have to change the path to the model output files to your personal folder structure.
 
 
! ***  Only run it in -nodisplay mode ***
 
 
! Save current state
 
	define REGION/DEFAULT "save"
	set data/save
	set grid/save
 
! Clean everything
 
	cancel data/all
	cancel variable/all
	cancel region/all
 
 
! Script function
 
 
! Pipe Process Elimination

	use /Users/mjuerchott/runs/2021_02_Pipes_noLand_COU85/tavg.02025.01.01.nc		![d=1]	ArtUp_8.5
	use /Users/mjuerchott/runs/2021_11_Pipes_noLand_COU60/tavg.02025.01.01.nc		![d=2]	ArtUp_6.0
	use /Users/mjuerchott/runs/2021_11_Pipes_noLand_COU45/tavg.02025.01.01.nc		![d=3]	ArtUp_4.5
	use /Users/mjuerchott/runs/2021_11_Pipes_noLand_COU26/tavg.02025.01.01.nc		![d=4]	ArtUp_2.6
	use /Users/mjuerchott/runs/2021_02_Pipes_noLand_CTRL/tavg.02025.01.01.nc		![d=5]	ArtUp_0.0

	use /Users/mjuerchott/runs/2021_02_noPipes_noLand_COU85_cont/tavg.02025.01.01.nc	![d=6]	REF_8.5
	use /Users/mjuerchott/runs/2021_11_noPipes_noLand_COU60_cont/tavg.02025.01.01.nc	![d=7]	REF_6.0
	use /Users/mjuerchott/runs/2021_11_noPipes_noLand_COU45_cont/tavg.02025.01.01.nc	![d=8]	REF_4.5
	use /Users/mjuerchott/runs/2021_11_noPipes_noLand_COU26_cont/tavg.02025.01.01.nc	![d=9]	REF_2.6
	use /Users/mjuerchott/runs/2021_02_noPipes_noLand_CTRL_cont/tavg.02025.01.01.nc		![d=10]	REF_0.0


! Global alkalinity / depth profile

	let ALK_REF_85 = O_ALK[d=6] *1000  ! unit from mol/m3 to mmol/m3
	let ALK_REF_60 = O_ALK[d=7] *1000
	let ALK_REF_45 = O_ALK[d=8] *1000
	let ALK_REF_26 = O_ALK[d=9] *1000
	let ALK_REF_00 = O_ALK[d=10] *1000

	let ALK_REF_ave_depth_85 = ALK_REF_85[i=@AVE,j=@AVE,l=16]
	let ALK_REF_ave_depth_60 = ALK_REF_60[i=@AVE,j=@AVE,l=16]
	let ALK_REF_ave_depth_45 = ALK_REF_45[i=@AVE,j=@AVE,l=16]
	let ALK_REF_ave_depth_26 = ALK_REF_26[i=@AVE,j=@AVE,l=16]
	let ALK_REF_ave_depth_00 = ALK_REF_00[i=@AVE,j=@AVE,l=16]


! Alkalinity change at surface layer via AU

	let xpreALK_85 =	(O_ALK[d=1]  - O_ALK[d=6]) *1000  ! unit from mol/m3 to mmol/m3
	let xpreALK_60 =	(O_ALK[d=2]  - O_ALK[d=7]) *1000  
	let xpreALK_45 =	(O_ALK[d=3]  - O_ALK[d=8]) *1000  
	let xpreALK_26 =	(O_ALK[d=4]  - O_ALK[d=9]) *1000 
	let xpreALK_00 =	(O_ALK[d=5]  - O_ALK[d=10]) *1000

	let/title="RPC 8.5"/units="TA mmol m-3" preALK_85  = xpreALK_85[I=@AVE,J=@AVE,K=1@AVE]
	let/title="RCP 6.0"/units="TA mmol m-3" preALK_60  = xpreALK_60[I=@AVE,J=@AVE,K=1@AVE]
	let/title="RPC 4.5"/units="TA mmol m-3" preALK_45  = xpreALK_45[I=@AVE,J=@AVE,K=1@AVE]
	let/title="RPC 2.6"/units="TA mmol m-3" preALK_26  = xpreALK_26[I=@AVE,J=@AVE,K=1@AVE]
	let/title="RPC 0.0"/units="TA mmol m-3" preALK_00  = xpreALK_00[I=@AVE,J=@AVE,K=1@AVE]



! plot routine

	cancel windows/all
	set windows/new
	set text/font='Times New Roman'
	set text/isiz=4

	go landscape2x2.jnl
	define view/xlimits=0.2,0.7/ylimits=0.05,0.55 ll4
	define view/xlimits=0.2,0.7/ylimits=0.5,1 ul4


	set view ul4
	plot/nolabel/color=black/thick=2 ALK_REF_ave_depth_85
	plot/nolabel/overlay/color=blue/thick=2 ALK_REF_ave_depth_60
	plot/nolabel/overlay/color=red/thick=2 ALK_REF_ave_depth_45
	plot/nolabel/overlay/color=green/thick=2 ALK_REF_ave_depth_26
	plot/nolabel/overlay/color=(100,68,0)/thick=2 ALK_REF_ave_depth_00
	label/nouser -0.8,-0.8,0,0,0.2 "(a)"
	label/nouser -1.3,1.7,0,90,0.13 "depth [m]"
	label/nouser 2.5,-0.8,0,0,0.13 "TA [mmol/m³]"


	set view ll4
	plot/nolabel/hlimit=2020:2101:10/vlimit=-3:5:0.5/color=black/thick=2 preALK_85
	plot/nolabel/overlay/color=blue/thick=2 preALK_60
	plot/nolabel/overlay/color=red/thick=2 preALK_45
	plot/nolabel/overlay/color=green/thick=2 preALK_26
	plot/nolabel/overlay/color=(100,68,0)/thick=2 preALK_00
	label/nouser -0.8,-0.8,0,0,0.2 "(b)"
	label/nouser -1.3,1.7,0,90,0.13 "<#916>TA [mmol/m³]"
	label/nouser 2.5,-0.8,0,0,0.13 "[years]"


! plot legend

	set text/font='Times New Roman'/color=black
	label/nouser 0.2,7.5,-1,0,0.11 - RCP 8.5

	set text/font='Times New Roman'/color=blue
	label/nouser 0.2,7.2,-1,0,0.11 - RCP 6.0

	set text/font='Times New Roman'/color=red
	label/nouser 0.2,6.9,-1,0,0.11 - RCP 4.5

	set text/font='Times New Roman'/color=green
	label/nouser 0.2,6.6,-1,0,0.11 - RCP 2.6

	set text/font='Times New Roman'/color=(100,68,0)
	label/nouser 0.2,6.3,-1,0,0.11 - no emission

	set text/font='Times New Roman'/color=black



frame/file = "FIG_S5.png"


! Cancel everything
 
	cancel data/all
	cancel variable/all
	cancel region/all
 
! Restore current state
 
	set data/restore
	set grid/restore