! This PyFerret script calculates the values for Figure 1b and 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 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_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


! DICpre
	let xPREDIC_85 =	 (O_PREDIC[d=1]  - O_PREDIC[d=6])	 *12 *10^(-15)	!unit into Pg C
	let xPREDIC_60 =         (O_PREDIC[d=2]  - O_PREDIC[d=7])        *12 *10^(-15)
	let xPREDIC_45 =         (O_PREDIC[d=3]  - O_PREDIC[d=8])        *12 *10^(-15)
	let xPREDIC_26 =         (O_PREDIC[d=4]  - O_PREDIC[d=9])        *12 *10^(-15)
	let xPREDIC_00 =         (O_PREDIC[d=5]  - O_PREDIC[d=10])       *12 *10^(-15)

! DICrem
	let xDICR_85 =         (O_DICR[d=1]  - O_DICR[d=6])        *12 *10^(-15)
	let xDICR_60 =         (O_DICR[d=2]  - O_DICR[d=7])        *12 *10^(-15)
	let xDICR_45 =         (O_DICR[d=3]  - O_DICR[d=8])        *12 *10^(-15)
	let xDICR_26 =         (O_DICR[d=4]  - O_DICR[d=9])        *12 *10^(-15)
	let xDICR_00 =         (O_DICR[d=5]  - O_DICR[d=10])       *12 *10^(-15)

! DIC
	let xDIC_85 =         (O_DIC[d=1]  - O_DIC[d=6])        *12 *10^(-15)
	let xDIC_60 =         (O_DIC[d=2]  - O_DIC[d=7])        *12 *10^(-15)
	let xDIC_45 =         (O_DIC[d=3]  - O_DIC[d=8])        *12 *10^(-15)
	let xDIC_26 =         (O_DIC[d=4]  - O_DIC[d=9])        *12 *10^(-15)
	let xDIC_00 =         (O_DIC[d=5]  - O_DIC[d=10])       *12 *10^(-15)



! DICpre
	let/title="PREDIC_85"/units="Pg/ocean" PREDIC_85  = xPREDIC_85[I=@DIN,J=@DIN,K=@DIN,L=16]
	let/title="PREDIC_60"/units="Pg/ocean" PREDIC_60  = xPREDIC_60[I=@DIN,J=@DIN,K=@DIN,L=16]
	let/title="PREDIC_45"/units="Pg/ocean" PREDIC_45  = xPREDIC_45[I=@DIN,J=@DIN,K=@DIN,L=16]
	let/title="PREDIC_26"/units="Pg/ocean" PREDIC_26  = xPREDIC_26[I=@DIN,J=@DIN,K=@DIN,L=16]
	let/title="PREDIC_00"/units="Pg/ocean" PREDIC_00  = xPREDIC_00[I=@DIN,J=@DIN,K=@DIN,L=16]

! DICrem
	let/title="DICR_85"/units="Pg/ocean" DICR_85  = xDICR_85[I=@DIN,J=@DIN,K=@DIN,L=16]
	let/title="DICR_60"/units="Pg/ocean" DICR_60  = xDICR_60[I=@DIN,J=@DIN,K=@DIN,L=16]
	let/title="DICR_45"/units="Pg/ocean" DICR_45  = xDICR_45[I=@DIN,J=@DIN,K=@DIN,L=16]
	let/title="DICR_26"/units="Pg/ocean" DICR_26  = xDICR_26[I=@DIN,J=@DIN,K=@DIN,L=16]
	let/title="DICR_00"/units="Pg/ocean" DICR_00  = xDICR_00[I=@DIN,J=@DIN,K=@DIN,L=16]

! DICca
	let/title="CaDIC_85"/units="Pg/ocean" CaDIC_85  = xDIC_85[I=@DIN,J=@DIN,K=@DIN,L=16] - PREDIC_85 - DICR_85
	let/title="CaDIC_60"/units="Pg/ocean" CaDIC_60  = xDIC_60[I=@DIN,J=@DIN,K=@DIN,L=16] - PREDIC_60 - DICR_60
	let/title="CaDIC_45"/units="Pg/ocean" CaDIC_45  = xDIC_45[I=@DIN,J=@DIN,K=@DIN,L=16] - PREDIC_45 - DICR_45
	let/title="CaDIC_26"/units="Pg/ocean" CaDIC_26  = xDIC_26[I=@DIN,J=@DIN,K=@DIN,L=16] - PREDIC_26 - DICR_26
	let/title="CaDIC_00"/units="Pg/ocean" CaDIC_00  = xDIC_00[I=@DIN,J=@DIN,K=@DIN,L=16] - PREDIC_00 - DICR_00


! Calculation 
	list PREDIC_85
	list PREDIC_60
	list PREDIC_45
	list PREDIC_26
	list PREDIC_00

	list DICR_85
	list DICR_60
	list DICR_45
	list DICR_26
	list DICR_00

	list CaDIC_85
	list CaDIC_60
	list CaDIC_45
	list CaDIC_26
	list CaDIC_00


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