ERF_InitCustomPertVels_CouettePoiseuille.H File Reference

Go to the source code of this file.

Functions
ParmParse	pp_for_pert_vels ("prob")
pp_for_pert_vels	query ("U_0", U_0)
pp_for_pert_vels	query ("V_0", V_0)
pp_for_pert_vels	query ("W_0", W_0)
pp_for_pert_vels	get ("prob_type", prob_type)
	AMREX_ALWAYS_ASSERT (prob_type==1||prob_type==10||prob_type==11||prob_type==20||prob_type==21)
pp_for_pert_vels	query ("pert_periods_u", pert_periods_u)
pp_for_pert_vels	query ("pert_periods_v", pert_periods_v)
pp_for_pert_vels	query ("pert_delta_u", pert_delta_u)
pp_for_pert_vels	query ("pert_delta_v", pert_delta_v)
	if (prob_type==1)

Variables
const bool	use_z_phys = use_physical_z()
const Real	phys_z_lo = phys_zlo(geomdata, lev)
const Real	phys_z_hi = phys_zhi(geomdata, lev)
Real	U_0 = zero
Real	V_0 = zero
Real	W_0 = zero
int	prob_type
Real	pert_periods_u = amrex::Real(5.0)
Real	pert_periods_v = amrex::Real(5.0)
Real	pert_delta_u = zero
Real	pert_delta_v = zero
auto	problo = geomdata.ProbLo()
auto	probhi = geomdata.ProbHi()
Real	aval = pert_periods_u * two * PI / (probhi[1] - problo[1])
Real	bval = pert_periods_v * two * PI / (probhi[0] - problo[0])

Function Documentation

AMREX_ALWAYS_ASSERT()

AMREX_ALWAYS_ASSERT

(

prob_type

= =1||prob_type==10||prob_type==11||prob_type==20||prob_type==21

)

get()

pp_for_pert_vels get	(	"prob_type"	,
		prob_type
	)

Referenced by ERF::Advance(), TwoStreamRadiation::advance(), ERF::advance_dycore(), ERF::advance_lsm(), ERF::advance_radiation(), SurfaceModel::calculate_weight_average(), SAM::Compute_Coefficients(), SurfaceLayer::compute_fluxes(), ERF::compute_max_buoyancy_gradp_diagnostic(), ERF::compute_max_pressure_gradient_diagnostic(), SurfaceLayer::compute_pblh(), SurfaceLayer::fill_lateral_surface_parameter_ghosts(), ERF::FillPlot2DVars(), ERF::FillPlot3DVars(), eb_aux_::getAreaFrac(), eb_aux_::getFaceCent(), if(), ERF::init_from_input_sounding(), ERF::init_ibseb(), ERF::init_obs_nudging(), ERF::init_stuff(), ERF::InitData_post(), SurfaceModel::initialize_for_level(), ERF::Interp2DArrays(), ERF::Interp2DTimeSlice(), ERF::interp_land_masks_from_coarse(), SurfaceLayer::make_SurfaceLayer_at_level(), ERF::MakeNewLevelFromCoarse(), ERF::MakeNewLevelFromScratch(), ERF::RemakeLevel(), ERF::set_cf_fill_options(), ERF::timeStep(), ERF::turbPert_update_automatic(), ERF::turbPert_update_uniform_wind(), ERF::volWgtSumMF(), and ERF::WriteCheckpointFile().

Here is the caller graph for this function:

if()

if

(

prob_type

= = 1

)

                        {
 
        ParallelFor(xbx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept {
            const auto *const dx       = geomdata.CellSize();
            const Real z = use_z_phys ? erf_prob_heights::z_xface(i,j,k,z_nd) - phys_z_lo : (k + myhalf) * dx[2];
            x_vel_pert(i, j, k) = U_0 * z / (phys_z_hi - phys_z_lo);
        });
 
        ParallelFor(ybx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept {
            const auto *const dx       = geomdata.CellSize();
            const Real z = use_z_phys ? erf_prob_heights::z_yface(i,j,k,z_nd) - phys_z_lo : (k + myhalf) * dx[2];
            y_vel_pert(i, j, k) = V_0 * z / (phys_z_hi - phys_z_lo);
        });
 
        ParallelFor(zbx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept {
            z_vel_pert(i, j, k) = W_0;
        });
 
    // Poiseuille flow
    } else if (prob_type == 10 || prob_type == 11) {

Here is the call graph for this function:

pp_for_pert_vels()

ParmParse pp_for_pert_vels

(

"prob"

)

query() [1/7]

pp_for_pert_vels query	(	"pert_delta_u"	,
		pert_delta_u
	)

query() [2/7]

pp_for_pert_vels query	(	"pert_delta_v"	,
		pert_delta_v
	)

query() [3/7]

pp_for_pert_vels query	(	"pert_periods_u"	,
		pert_periods_u
	)

query() [4/7]

pp_for_pert_vels query	(	"pert_periods_v"	,
		pert_periods_v
	)

query() [5/7]

pp_for_pert_vels query	(	"U_0"	,
		U_0
	)

query() [6/7]

pp_for_pert_vels query	(	"V_0"	,
		V_0
	)

query() [7/7]

pp_for_pert_vels query	(	"W_0"	,
		W_0
	)

Variable Documentation

aval

Real aval = pert_periods_u * two * PI / (probhi[1] - problo[1])

bval

Real bval = pert_periods_v * two * PI / (probhi[0] - problo[0])

pert_delta_u

Real pert_delta_u = zero

pert_delta_v

Real pert_delta_v = zero

pert_periods_u

Real pert_periods_u = amrex::Real(5.0)

pert_periods_v

Real pert_periods_v = amrex::Real(5.0)

phys_z_hi

const Real phys_z_hi = phys_zhi(geomdata, lev)

Referenced by if().

phys_z_lo

const Real phys_z_lo = phys_zlo(geomdata, lev)

Referenced by if().

prob_type

int prob_type

probhi

auto probhi = geomdata.ProbHi()

problo

auto problo = geomdata.ProbLo()

U_0

Real U_0 = zero

Referenced by if().

use_z_phys

const bool use_z_phys = use_physical_z()

Referenced by if().

V_0

Real V_0 = zero

Referenced by if().

W_0

Real W_0 = zero

Referenced by if().