! This PyFerret script plots FIG_S4 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_cont/tavg.02025.01.01.nc 	![d=1]	REF_0.0

	let biological = O_DICR / O_PO4R
	let solubility = O_preDIC / O_prePO4


! 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
	shade/nolabel/levels=(104.5,107.5,1)/z=1000/l=16 biological
	label/nouser -0.8,-0.8,0,0,0.2 "(a)"
	label/nouser 5.5,3.8,1,0,0.13 "[C<sup>rem</sup> / PO<sub>4</sub><sup>rem</sup>]"

	set view ll4
	shade/nolabel/levels=(1000,5000,500)/z=1000/l=16 solubility
	label/nouser -0.8,-0.8,0,0,0.2 "(b)"
	label/nouser 5.5,3.8,1,0,0.13 "[C<sup>pre</sup> / PO<sub>4</sub><sup>pre</sup>]"



frame/file = "FIG_S4.png"



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