ERF_SEBTurbulentFlux.cpp File Reference

#include <ERF_SEBTurbulentFlux.H>
#include <ERF_Plotfile2DFill.H>
#include <ERF_Constants.H>
#include <AMReX_Gpu.H>

Include dependency graph for ERF_SEBTurbulentFlux.cpp:

Functions
void	fill_seb_field_from_land_surface (MultiFab &seb_field, const MultiFab *field, Real fallback, Real scale, const MultiFab *add_field)
SEBFluxOrigin	fill_seb_turbulent_flux (MultiFab &seb_flux, SEBTurbulentFlux kind, const MultiFab *land_surface_field, const MultiFab *surface_layer_flux, int surface_k, SEBTurbulentFluxSource source, Real fallback)

Function Documentation

fill_seb_field_from_land_surface()

void fill_seb_field_from_land_surface	(	MultiFab &	seb_field,
		const MultiFab *	field,
		Real	fallback,
		Real	scale,
		const MultiFab *	add_field
	)

{
    seb_field.setVal(fallback);
    if (field != nullptr) {
        overwrite_valid_land_surface_cells(seb_field, *field, scale, add_field);
    }
}

fill_seb_turbulent_flux()

SEBFluxOrigin fill_seb_turbulent_flux	(	MultiFab &	seb_flux,
		SEBTurbulentFlux	kind,
		const MultiFab *	land_surface_field,
		const MultiFab *	surface_layer_flux,
		int	surface_k,
		SEBTurbulentFluxSource	source,
		Real	fallback
	)

{
    // The fallback fills every cell, halo included; the surface layer's flux then
    // overwrites the valid cells when it is the source, and the land model's valid
    // values overwrite those in turn.
    seb_flux.setVal(fallback);
    const SEBFluxOrigin below = select_seb_flux_origin(false, source,
                                                       surface_layer_flux != nullptr);
    if (below == SEBFluxOrigin::SurfaceLayer) {
        // The 2D-output conversions, so the balance's H and LE are exactly the
        // sensible_heat_flux and latent_heat_flux a plotfile reports.
        if (kind == SEBTurbulentFlux::Sensible) {
            plotfile2d::fill_sensible_heat_flux_from_klevel_or_missing(
                seb_flux, 0, surface_layer_flux, surface_k, fallback);
        } else {
            plotfile2d::fill_latent_heat_flux_from_klevel_or_missing(
                seb_flux, 0, surface_layer_flux, surface_k, fallback);
        }
    }
    if (land_surface_field != nullptr) {
        overwrite_valid_land_surface_cells(seb_flux, *land_surface_field, Real(1.0), nullptr);
    }
    return select_seb_flux_origin(land_surface_field != nullptr, source,
                                  surface_layer_flux != nullptr);
}

Referenced by TwoStreamRadiation::advance().

Here is the call graph for this function:

Here is the caller graph for this function: