ERF_InitRayleighDamping.H File Reference

#include <ERF_ParFunctions.H>
#include <AMReX_ParmParse.H>

Include dependency graph for ERF_InitRayleighDamping.H:

Go to the source code of this file.

Functions
void	erf_init_rayleigh (amrex::Vector< amrex::Vector< amrex::Real > > &rayleigh_ptrs, amrex::Geometry const &geom, std::unique_ptr< amrex::MultiFab > &, amrex::Real) override

Function Documentation

erf_init_rayleigh()

void erf_init_rayleigh ( amrex::Vector< amrex::Vector< amrex::Real > > &  rayleigh_ptrs, amrex::Geometry const &  geom, std::unique_ptr< amrex::MultiFab > &  , amrex::Real    )

override

Initialize a Rayleigh damping layer with the same structure as in WRF, based on Durran and Klemp 1983

{
    amrex::ParmParse pp_prob("prob");
    amrex::Real U_0 =   0.0; pp_prob.query("rayleigh_U_0", U_0);
    amrex::Real V_0 =   0.0; pp_prob.query("rayleigh_V_0", V_0);
    amrex::Real W_0 =   0.0; pp_prob.query("rayleigh_W_0", W_0);
    amrex::Real T_0 = 300.0; pp_prob.query("rayleigh_T_0", T_0);
 
    const int khi = geom.Domain().bigEnd()[2];
 
    for (int k = 0; k <= khi; k++)
    {
        rayleigh_ptrs[Rayleigh::ubar][k]     = U_0;
        rayleigh_ptrs[Rayleigh::vbar][k]     = V_0;
        rayleigh_ptrs[Rayleigh::wbar][k]     = W_0;
        rayleigh_ptrs[Rayleigh::thetabar][k] = T_0;
    }
 
    amrex::ParmParse pp("erf");
    std::string prob_name; pp.get("prob_name",prob_name);
    if (prob_name == "Moving Terrain") {
        for (int k = 60; k <= khi; k++)
        {
            rayleigh_ptrs[Rayleigh::ubar][k]     = 2.0;
            rayleigh_ptrs[Rayleigh::vbar][k]     = 1.0;
            rayleigh_ptrs[Rayleigh::wbar][k]     = 0.0;
            rayleigh_ptrs[Rayleigh::thetabar][k] = T_0;
        }
    }
}

Here is the call graph for this function: