Extend ClimaLand turbulent_fluxes! methods for CoupledAtmosphere

[juliasloan25]

When coupling the full land model, we'll need to compute turbulent fluxes for each component of the land model. The calculations needed differ slightly between components, so when computing fluxes the coupler will call ClimaLand functions that dispatch on the component model type. We currently have these functions defined for PrescribedAtmosphere, but need to extend them for CoupledAtmosphere

Since these functions require ClimaCoupler get_field functions to retrieve the atmos inputs, I think they need to be written in ClimaCoupler rather than ClimaLand. They can still call the ClimaLand (soil/canopy_)turbulent_fluxes_at_a_point functions

To do

• turbulent_fluxes! default (used for snow)
  • get ts_in, u, h from atmos model instead of from driver
  • returns LHF, SHF, E~ (vapor flux), r_ae
• turbulent_fluxes! EnergyHydrology
  • get ts_in, u, h from atmos model instead of from driver
  • returns LHF, SHF, E~_liq (liquid vapor flux), E~_liq (ice vapor flux), r_ae, dlhfdT, dshfdT
• turbulent_fluxes! CanopyModel
  • get ts_in, u, h from atmos model instead of from driver
  • returns LHF, SHF, E~, r_ae, dlhfdqc, dshfdTc

Questions

• Soil and Canopy turbulent_fluxes! functions use atmosphere gustiness (here and here). What is this, and can we get it from atmos?
• The return values of these functions don't match what the coupler/atmos expect. How can we map from LHF, SHF, E~, r_ae to turbulent_energy_flux, turbulent_moisture_flux?
  • For the bucket, we just provide SHF and E~