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

[navidcy]

E.g., here:

https://github.com/simone-silvestri/OceanScalingTests.jl/blob/f034fdb32903c42a09961213308a5b1559b050ed/src/bathymetry_and_forcings.jl#L61-L70

cc @glwagner

[simone-silvestri]

I decided to use a piecewise cubic function for the wind stress (and salinity flux). I guess we could use a parabolic one also

[navidcy]

Is it constant values in between?

[glwagner]

How about

τ₀ * (cos(3π * φ/Δφ) - exp(-φ^2 / 2Δφₑ^2))

?

[navidcy]

I think that should work. Perhaps we should ensure that wind stress goes to zero at the North-South latitudinal boundaries otherwise you get Kelvin waves and circulation at the boundaries.

[glwagner]

Ah nice. Then Δφ is the total latitudinal extent of the domain. Δφₑ is the width of the equatorial bump, typically about 10 degrees latitude. We can also add a parameter to express the relative magnitude of the equatorial bump.