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


! Control COU85
        use /Users/mjuerchott/runs/2021_02_noPipes_noLand_COU85_cont/tavg.02025.01.01.nc        ! [d=1]	REF_8.5

! Pipes CUO85
        use /Users/mjuerchott/runs/2021_02_Pipes_noLand_COU85/tavg.02025.01.01.nc               ! [d=2]	ArtUp_8.5



! Temperature

	let Temp_Control = O_TEMP[d=1]
	let Temp_Pipe = O_TEMP[d=2]
	let Temp_Change = Temp_Pipe - Temp_Control


! Plot Routine

	cancel windows/all
	set windows/new
	set text/font='Times New Roman'
	set text/isiz=4

	go portrait1x3.jnl
	define view/xlimits=0,0.48/ylimits=0.2,0.8 lef          ! bot
	define view/xlimits=0.33,0.62/ylimits=0.215,0.785 mid   ! mid
	define view/xlimits=0.47,0.95/ylimits=0.2,0.8 rig       ! top



        set view lef    ! Atlantic N-S
        fill/nolabels/nokey/AXES=1,1,1,0/pal=fDIC_rev/line/ylimits=0:6000:500/levels=(-1.5,3.5,0.5)/xlimits=-60:80:-20/X=30W/Y=80N:60S/l=16 Temp_Change

        set view mid    ! Southern Ocean W-E
        fill/nolabel/nokey/AXES=1,1,0,0/pal=fDIC_rev/line/ylimits=0:6000:500/levels=(-1.5,3.5,0.5)/xlimits=330:510:45/X=31W:149W/Y=60S/l=16 Temp_Change

        set view rig    ! Pacific S-N
        fill/nolabels/key/AXES=1,1,0,0/pal=fDIC_rev/line/ylimits=0:6000:500/levels=(-1.5,3.5,0.5)/xlimits=-60:80:20/X=150W/Y=80N:60S/l=16 Temp_Change


        label/nouser -7.45,2.3,0,90,0.13 "depth [m]"
	label/nouser 5.0,4.8,1,0,0.13 "<#916>T [°C]"
	label/nouser -4.4,-0.8,0,0,0.13 "Atlantic"
	label/nouser -1.1,-0.8,0,0,0.13 "Southern Ocean"
	label/nouser 2.4,-0.8,0,0,0.13 "Pacific"


        ! Dotted Lines
        LABEL -126.8,6000,0,90,0.24 "_"
        LABEL -126.8,5500,0,90,0.24 "_"
        LABEL -126.8,5000,0,90,0.24 "_"
        LABEL -126.8,4500,0,90,0.24 "_"
        LABEL -126.8,4000,0,90,0.24 "_"
        LABEL -126.8,3500,0,90,0.24 "_"
        LABEL -126.8,3000,0,90,0.24 "_"
        LABEL -126.8,2500,0,90,0.24 "_"
        LABEL -126.8,2000,0,90,0.24 "_"
        LABEL -126.8,1500,0,90,0.24 "_"
        LABEL -126.8,1000,0,90,0.24 "_"
        LABEL -126.8,500,0,90,0.24 "_"
        LABEL -126.8,0,0,90,0.24 "_"

        LABEL -62.2,6000,0,90,0.24 "_"
        LABEL -62.2,5500,0,90,0.24 "_"
        LABEL -62.2,5000,0,90,0.24 "_"
        LABEL -62.2,4500,0,90,0.24 "_"
        LABEL -62.2,4000,0,90,0.24 "_"
        LABEL -62.2,3500,0,90,0.24 "_"
        LABEL -62.2,3000,0,90,0.24 "_"
        LABEL -62.2,2500,0,90,0.24 "_"
        LABEL -62.2,2000,0,90,0.24 "_"
        LABEL -62.2,1500,0,90,0.24 "_"
        LABEL -62.2,1000,0,90,0.24 "_"
        LABEL -62.2,500,0,90,0.24 "_"
        LABEL -62.2,0,0,90,0.24 "_"



frame/file = "FIG_S9.png"




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