! This PyFerret script plots FIG_S7 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 -nojnl 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

	use /Users/mjuerchott/runs/2021_02_noPipes_noLand_CTRL/tsi.01766.01.01.nc		! [d=1]	REF_0.0
	use /Users/mjuerchott/runs/2021_11_noPipes_noLand_COU26_cont/tsi.02021.01.01.nc		! [d=2]	REF_2.6
	use /Users/mjuerchott/runs/2021_11_noPipes_noLand_COU45_cont/tsi.02021.01.01.nc		! [d=3]	REF_4.5
	use /Users/mjuerchott/runs/2021_11_noPipes_noLand_COU60_cont/tsi.02021.01.01.nc		! [d=4]	REF_6.0
	use /Users/mjuerchott/runs/2021_02_noPipes_noLand_COU85/tsi.01766.01.01.nc		! [d=5]	REF_8.5

	use /Users/mjuerchott/runs/2021_02_Pipes_noLand_CTRL/tsi.02021.01.01.nc			! [d=6]	ArtUp_0.0
	use /Users/mjuerchott/runs/2021_11_Pipes_noLand_COU26/tsi.02021.01.01.nc		! [d=7]	ArtUp_2.6
	use /Users/mjuerchott/runs/2021_11_Pipes_noLand_COU45/tsi.02021.01.01.nc		! [d=8]	ArtUp_4.5
	use /Users/mjuerchott/runs/2021_11_Pipes_noLand_COU60/tsi.02021.01.01.nc		! [d=9]	ArtUp_6.0
	use /Users/mjuerchott/runs/2021_02_Pipes_noLand_COU85/tsi.02021.01.01.nc		! [d=10]ArtUp_8.5


! F_DIC REF Overview

	let F_DIC_00 = F_DIC[d=1] *60*60*24*365
	let F_DIC_26 = F_DIC[d=2] *60*60*24*365
	let F_DIC_45 = F_DIC[d=3] *60*60*24*365
	let F_DIC_60 = F_DIC[d=4] *60*60*24*365
	let F_DIC_85 = F_DIC[d=5] *60*60*24*365


! A_SAT ArtUp Overview

	let F_DIC_AU_00 = F_DIC[d=6] *60*60*24*365
	let F_DIC_AU_26 = F_DIC[d=7] *60*60*24*365
	let F_DIC_AU_45 = F_DIC[d=8] *60*60*24*365
	let F_DIC_AU_60 = F_DIC[d=9] *60*60*24*365
	let F_DIC_AU_85 = F_DIC[d=10] *60*60*24*365

	let F_DIC_change_00 = F_DIC_AU_00 - F_DIC_00
	let F_DIC_change_26 = F_DIC_AU_26 - F_DIC_26
	let F_DIC_change_45 = F_DIC_AU_45 - F_DIC_45
	let F_DIC_change_60 = F_DIC_AU_60 - F_DIC_60
	let F_DIC_change_85 = F_DIC_AU_85 - F_DIC_85



! 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/VLIMIT=-0.2:2.6:0.2/HLIMITS=2020:2101:10/color=(100,68,0)/thick=2 F_DIC_00
	plot/nolabel/overlay/HLIMITS=2020:2101:10/color=green/thick=2 F_DIC_26
	plot/nolabel/overlay/HLIMITS=2020:2101:10/color=red/thick=2 F_DIC_45
	plot/nolabel/overlay/HLIMITS=2020:2101:10/color=blue/thick=2 F_DIC_60
	plot/nolabel/overlay/HLIMITS=2020:2101:10/color=black/thick=2 F_DIC_85
	label/nouser -0.8,-0.8,0,0,0.2 "(a)"
	label/nouser -1.3,1.7,0,90,0.13 "fDIC [mol/m²/year]"
	label/nouser 2.5,-0.8,0,0,0.13 "[years]"


	set view ll4
	plot/nolabel/VLIMIT=-0.2:0.6:0.1/HLIMITS=2020:2101:10/color=(100,68,0)/thick=2 F_DIC_change_00
	plot/nolabel/overlay/HLIMITS=2020:2101:10/color=green/thick=2 F_DIC_change_26
	plot/nolabel/overlay/HLIMITS=2020:2101:10/color=red/thick=2 F_DIC_change_45
	plot/nolabel/overlay/HLIMITS=2020:2101:10/color=blue/thick=2 F_DIC_change_60
	plot/nolabel/overlay/HLIMITS=2020:2101:10/color=black/thick=2 F_DIC_change_85
	label/nouser -0.8,-0.8,0,0,0.2 "(b)"
	label/nouser -1.3,1.7,0,90,0.13 "<#916>fDIC [mol/m²/year]"
	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,8.5,-1,0,0.11 - RCP 8.5

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

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

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

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

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



frame/file = "FIG_S7.png"



! Cancel everything

        cancel data/all
        cancel variable/all
        cancel region/all

! Restore current state

        set data/restore
        set grid/restore