Forcing

apply_global_salt_correction!(FS, ocean_mask, area)

Adjust the virtual salinity flux field so its globally integrated value is zero. This preserves total ocean salinity (Appendix B4 of Willeit et al. 2022).

Arguments

• FS : virtual salinity flux (nlon, nlat), psu·m·s⁻¹ — modified in place
• ocean_mask : Bool (nlon, nlat)
• area : cell area (nlon, nlat), m²

heat_flux_bc(QH, H_top, rho0, cp)

Convert surface heat flux QH (W m⁻²) to a temperature tendency (K s⁻¹) for the top layer. Positive QH = into ocean.

Returns a (nlon, nlat) array of ∂T/∂t forcing due to surface heating.

virtual_salinity_flux!(FS, FW, S_surface, H_top, ocean_mask)

Compute the virtual salinity flux FS (psu m s⁻¹) at the ocean surface.

Arguments

• FS : output virtual salinity flux (nlon, nlat), psu·m·s⁻¹
• FW : net freshwater flux = E − P − R (nlon, nlat), m s⁻¹ (positive = into ocean)
• S_surface : surface salinity (nlon, nlat), psu
• H_top : top layer thickness (m), scalar
• ocean_mask : Bool (nlon, nlat)

geothermal_tendency(Q_geo, H_bottom, ocean_mask; rho0=1025.0, cp=3994.0)

Convert geothermal heat flux (W m⁻²) to a temperature tendency (K s⁻¹) for the bottom ocean layer.

Arguments

• Q_geo : geothermal heat flux (nlon, nlat), W m⁻²
• H_bottom : bottom layer thickness (m)
• ocean_mask: Bool (nlon, nlat)
• rho0 : reference density, kg m⁻³
• cp : specific heat, J kg⁻¹ K⁻¹

uniform_geothermal(nlon, nlat)

Return a uniform geothermal flux field using the default value QGEODEFAULT.