Prescribed surface heat and evaporative fluxes

[milankl]

Introduces PrescribedSurfaceHeatFlux and PrescribedSurfaceEvaporation to read those in from the diagnostic variabels for coupling.

[milankl]

Starting with 0 humidity and no surface heat/humidity flux but then suddenly imposing 5e-5 kg/s/m^2 on the northern hemisphere seems to work in order to prescribe fluxes

out.mp4

done via

using SpeedyWeather
spectral_grid = SpectralGrid()

surface_heat_flux = PrescribedSurfaceHeatFlux()
surface_evaporation = PrescribedSurfaceEvaporation()
model = PrimitiveWetModel(spectral_grid; surface_heat_flux, surface_evaporation)
simulation = initialize!(model);

# no fluxes no humidity
set!(simulation, humid=0)
run!(simulation, period=Day(10), output=true)

# constant northern hemisphere fluxes
set!(simulation.diagnostic_variables.physics.evaporative_flux, (λ, φ) -> φ > 0 ? 5e-5 : 0, model.geometry)
run!(simulation, period=Day(10), output=true)

[milankl]

and lowermost layer temperature after switching on a +100W/m^2 flux on the northern but -100W/m^2 flux on the southern hemisphere

out.mp4

via

set!(simulation.diagnostic_variables.physics.sensible_heat_flux, (λ, φ) -> sign(φ)*100, model.geometry)
run!(simulation, period=Day(10), output=true)