! This PyFerret script plots FIG_S2 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 (long time series)
	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 (long time series)


! A_CO2 CTRL Overview

	let A_CO2_00 = A_CO2[d=1]
	let A_CO2_26 = A_CO2[d=2]
	let A_CO2_45 = A_CO2[d=3]
	let A_CO2_60 = A_CO2[d=4]
	let A_CO2_85 = A_CO2[d=5]

! A_SAT CTRL Overview
	let A_SAT_00 = A_SAT[d=1]
	let A_SAT_26 = A_SAT[d=2]
	let A_SAT_45 = A_SAT[d=3]
	let A_SAT_60 = A_SAT[d=4]
	let A_SAT_85 = A_SAT[d=5]



! 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=250:1000:50/HLIMITS=1750:2101:50/color=(100,68,0)/thick=2 A_CO2_00
	plot/nolabel/overlay/HLIMITS=1750:2101:50/color=green/thick=2 A_CO2_26
	plot/nolabel/overlay/HLIMITS=1750:2101:50/color=red/thick=2 A_CO2_45
	plot/nolabel/overlay/HLIMITS=1750:2101:50/color=blue/thick=2 A_CO2_60
	plot/nolabel/overlay/HLIMITS=1750:2101:50/color=black/thick=2 A_CO2_85
	label/nouser -0.8,-0.8,0,0,0.2 "(a)"
	label/nouser -1.3,1.7,0,90,0.13 "CO<sub>2</sub><sup>atm</sup> [ppm]"
	label/nouser 2.5,-0.8,0,0,0.13 "[years]"


	set view ll4
	plot/nolabel/VLIMIT=12.5:17.5:0.5/HLIMITS=1750:2101:50/color=(100,68,0)/thick=2 A_SAT_00
	plot/nolabel/overlay/HLIMITS=1750:2101:50/color=green/thick=2 A_SAT_26
	plot/nolabel/overlay/HLIMITS=1750:2101:50/color=red/thick=2 A_SAT_45
	plot/nolabel/overlay/HLIMITS=1750:2101:50/color=blue/thick=2 A_SAT_60
	plot/nolabel/overlay/HLIMITS=1750:2101:50/color=black/thick=2 A_SAT_85
	label/nouser -0.8,-0.8,0,0,0.2 "(b)"
	label/nouser -1.3,1.7,0,90,0.13 "SAT [°C]"
	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_S2.png"



! Cancel everything

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

! Restore current state

        set data/restore
        set grid/restore