Why do we need `cubic_profile` in the forcing functions?

Add Julia to path export PATH=/TO/JULIA:${PATH}
EXPERIMENT:
- "Quiescent" no flow, no bathymetry and zero flux boundary conditions
- "DoubleDrake" setup similar to https://doi.org/10.1175/2009JCLI3197.1, default
- "RealisticOcean" a realistic ocean setup with initial conditions and boundary forcing from ECCO2 Version 4 climatological data. It assumes data is available to load in the data folder
PROFILE: 0 (run simulation), 1 (profile 100 time steps marked with NVTX, default)
NZ: vertical levels (default 120)
WITHFLUXES: 0, 1 (only relevant for RealisticOcean, default 1)
RESTART: restart file iterations (only relevant for RealisticOcean, default "")

Satori

EXPERIMENT="DoubleDrake" RESOLUTION=6 sbatch -N1 satori_job.sh

Julia 1.9.0-rc2
Instantiate Julia environment JULIA_LOAD_PATH=:$(pwd)/perlmutter julia --project=. -e "import Pkg; Pkg.instantiate()"
Warmup GPUCompiler cache JULIA_LOAD_PATH=:$(pwd)/perlmutter julia --project=. experiments/run.jl

vorticity (left) and free surface (right) for the "DoubleDrake" experiment at resolution = 3 ran on 2 mpi processors (9 years integration)

double_drake_LQ.mp4

	@inline function wind_stress_coefficients(south_north_limit)
	below_45_coeffs = cubic_profile(-south_north_limit, -45.0, 0.0, 0.2, 0.0, 0.0) ./ 1000
	below_15_coeffs = cubic_profile(-45.0, -15.0, 0.2, -0.1, 0.0, 0.0) ./ 1000
	below_00_coeffs = cubic_profile(-15.0, 0.0, -0.1, -0.02, 0.0, 0.0) ./ 1000
	above_00_coeffs = cubic_profile(0.0, 15.0, -0.02, -0.1, 0.0, 0.0) ./ 1000
	above_15_coeffs = cubic_profile(15.0, 45.0, -0.1, 0.1, 0.0, 0.0) ./ 1000
	above_45_coeffs = cubic_profile(45.0, south_north_limit, 0.1, 0.0, 0.0, 0.0) ./ 1000

	return (below_45_coeffs, below_15_coeffs, below_00_coeffs, above_00_coeffs, above_15_coeffs, above_45_coeffs)
	end