ERF_EBStruct.H

Go to the documentation of this file.

#ifndef ERF_EB_STRUCT_H_
#define ERF_EB_STRUCT_H_
 
#include <string>
#include <iostream>
 
#include <AMReX_ParmParse.H>
#include <AMReX_Print.H>
#include <AMReX_Gpu.H>
 
enum struct EBBoundaryType {
    SlipWall, NoSlipWall, SurfaceLayer
};
 
enum EBTauType {
    tau_eb13, tau_eb23
};
 
enum EBGridType {
    xface, yface, zface
};
 
/**
 * Container holding EB-related choices
 */
 
struct EBChoice {
  public:
    void init_params(std::string pp_prefix)
    {
        amrex::ParmParse pp(pp_prefix);
 
        static std::string eb_boundary_type_string = "NoSlipWall";
        pp.query("eb_boundary_type",eb_boundary_type_string);
 
        if (!eb_boundary_type_string.compare("SlipWall")) {
            eb_boundary_type = EBBoundaryType::SlipWall;
        } else if (!eb_boundary_type_string.compare("NoSlipWall")) {
            eb_boundary_type = EBBoundaryType::NoSlipWall;
        } else if (!eb_boundary_type_string.compare("SurfaceLayer")) {
            eb_boundary_type = EBBoundaryType::SurfaceLayer;
        } else {
            amrex::Error("Don't know this eb_boundary_type");
        }
    }
 
    void display()
    {
        amrex::Print() << "EB choices: " << std::endl;
        if (eb_boundary_type == EBBoundaryType::SlipWall) {
            amrex::Print() << "    eb_boundary_type        : SlipWall" << std::endl;
        } else if (eb_boundary_type == EBBoundaryType::NoSlipWall) {
            amrex::Print() << "    eb_boundary_type        : NoSlipWall" << std::endl;
        } else if (eb_boundary_type == EBBoundaryType::SurfaceLayer) {
            amrex::Print() << "    eb_boundary_type        : SurfaceLayer" << std::endl;
        }
    }
 
    EBBoundaryType eb_boundary_type = EBBoundaryType::SlipWall;
};
#endif