! This PyFerret script calculates the values of Table S1 and 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_Pipes_noLand_COU85/tavg.02025.01.01.nc               ! [d=1]	ArtUp_8.5
use /Users/mjuerchott/runs/2021_02_noPipes_noLand_COU85_cont/tavg.02025.01.01.nc        ! [d=2]	REF_8.5



let Change_DIC = O_DIC[d=1] - O_DIC[d=2]
let Change_preDIC = O_preDIC[d=1] - O_preDIC[d=2]
let Change_DICr = O_DICR[d=1] - O_DICR[d=2]



! Budget DIC & DIC>1200m calculations

let budget_DIC = Change_DIC[x=@DIN, y=@DIN, z=@DIN, l=16] * 12 * 10^(-15)
let budget_preDIC = Change_preDIC[x=@DIN, y=@DIN, z=@DIN, l=16] * 12 * 10^(-15)
let budget_remDIC = Change_DICr[x=@DIN, y=@DIN, z=@DIN, l=16] * 12 * 10^(-15)
let budget_CaDIC = budget_DIC - budget_preDIC - budget_remDIC

let budget_DIC_1200 = Change_DIC[x=@DIN, y=@DIN, z=1200:6100@DIN, l=16] * 12 * 10^(-15)
let budget_preDIC_1200 = Change_preDIC[x=@DIN, y=@DIN, z=1200:6100@DIN, l=16] * 12 * 10^(-15)
let budget_remDIC_1200 = Change_DICr[x=@DIN, y=@DIN, z=1200:6100@DIN, l=16] * 12 * 10^(-15)
let budget_CaDIC_1200 = budget_DIC_1200 - budget_preDIC_1200 - budget_remDIC_1200



! Transfer Efficiency calculations

let xEXP_REF_85 = O_DETRREMI[d=2]
let xEXP_AU_85 = O_DETRREMI[d=1]
let xREMI_REF_85 = O_DETRREMI[d=2]
let xREMI_AU_85 = O_DETRREMI[d=1]


! >130m
let EXP_REF_85 = xEXP_REF_85[x=@DIN, y=@DIN, z=130:6100@DIN, l=1:16@AVE] * 6.625 * 12 * 10^(-15) *60*60*24*365
let REMI_REF_85_DOWN = xREMI_REF_85[x=@DIN, y=@DIN, z=1200:6100@DIN, l=1:16@AVE] * 6.625 * 12 * 10^(-15) *60*60*24*365

! >1200m
let EXP_AU_85 = xEXP_AU_85[x=@DIN, y=@DIN, z=130:6100@DIN, l=1:16@AVE] * 6.625 * 12 * 10^(-15) *60*60*24*365
let REMI_AU_85_DOWN = xREMI_AU_85[x=@DIN, y=@DIN, z=1200:6100@DIN, l=1:16@AVE] * 6.625 * 12 * 10^(-15) *60*60*24*365




! Output Numbers

! Table S1

list budget_DIC
list budget_DIC_1200

list budget_preDIC
list budget_preDIC_1200

list budget_remDIC
list budget_remDIC_1200

list budget_CaDIC
list budget_CaDIC_1200


! Table S2

list EXP_REF_85
list REMI_REF_85_DOWN

list EXP_AU_85
list REMI_AU_85_DOWN



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