Remove unnecessary override of init_custom_terrain

Note: The regression tests DynamicRefinement, IsentropicVortex, and ScalarAdvDiff all explicitly initialized z_phys_nd to be k*dz and then called FillBoundary(). However, if use_terrain==true (and init_custom_terrain() was called), then the default behavior should be equivalent: init_zlevels() called without grid stretching, then init_terrain_grid() called to populate z_phys_nd and then FillBoundary()
ewquon · Jan 4, 2024 · 8018f20 · 8018f20
1 parent f9082de
commit 8018f20
Show file tree

Hide file tree

Showing 39 changed files with 0 additions and 628 deletions.
diff --git a/Exec/ABL/prob.H b/Exec/ABL/prob.H
@@ -68,11 +68,6 @@ public:
         amrex::Array4<amrex::Real const> const& mf_v,
         const SolverChoice& sc) override;
 
-    void init_custom_terrain (
-        const amrex::Geometry& geom,
-        amrex::MultiFab& z_phys_nd,
-        const amrex::Real& time) override;
-
 protected:
     std::string name() override { return "ABL"; }
 

diff --git a/Exec/ABL/prob.cpp b/Exec/ABL/prob.cpp
@@ -164,28 +164,3 @@ Problem::init_custom_pert(
     }
   });
 }
-
-void
-Problem::init_custom_terrain(
-    const Geometry& /*geom*/,
-    MultiFab& z_phys_nd,
-    const Real& /*time*/)
-{
-    // Number of ghost cells
-    int ngrow = z_phys_nd.nGrow();
-
-    for ( MFIter mfi(z_phys_nd, TilingIfNotGPU()); mfi.isValid(); ++mfi )
-    {
-        // Grown box with no z range
-        amrex::Box xybx = mfi.growntilebox(ngrow);
-        xybx.setRange(2,0);
-
-        Array4<Real> z_arr = z_phys_nd.array(mfi);
-
-        ParallelFor(xybx, [=] AMREX_GPU_DEVICE (int i, int j, int) {
-
-            // Flat terrain with z = 0 at k = 0
-            z_arr(i,j,0) = 0.;
-        });
-    }
-}
diff --git a/Exec/ABL_input_sounding/prob.H b/Exec/ABL_input_sounding/prob.H
@@ -67,11 +67,6 @@ public:
         amrex::Array4<amrex::Real const> const& mf_v,
         const SolverChoice& sc) override;
 
-    void init_custom_terrain (
-        const amrex::Geometry& geom,
-        amrex::MultiFab& z_phys_nd,
-        const amrex::Real& time) override;
-
 protected:
     std::string name() override { return "ABL with input sounding"; }
 

diff --git a/Exec/ABL_input_sounding/prob.cpp b/Exec/ABL_input_sounding/prob.cpp
@@ -157,28 +157,3 @@ Problem::init_custom_pert(
     }
   });
 }
-
-void
-Problem::init_custom_terrain(
-    const Geometry& /*geom*/,
-    MultiFab& z_phys_nd,
-    const Real& /*time*/)
-{
-    // Number of ghost cells
-    int ngrow = z_phys_nd.nGrow();
-
-    for ( MFIter mfi(z_phys_nd, TilingIfNotGPU()); mfi.isValid(); ++mfi )
-    {
-        // Grown box with no z range
-        amrex::Box xybx = mfi.growntilebox(ngrow);
-        xybx.setRange(2,0);
-
-        Array4<Real> z_arr = z_phys_nd.array(mfi);
-
-        ParallelFor(xybx, [=] AMREX_GPU_DEVICE (int i, int j, int) {
-
-            // Flat terrain with z = 0 at k = 0
-            z_arr(i,j,0) = 0.;
-        });
-    }
-}
diff --git a/Exec/DevTests/MiguelDev/prob.H b/Exec/DevTests/MiguelDev/prob.H
@@ -57,11 +57,6 @@ public:
         amrex::Array4<amrex::Real const> const& mf_v,
         const SolverChoice& sc) override;
 
-    void init_custom_terrain (
-        const amrex::Geometry& geom,
-        amrex::MultiFab& z_phys_nd,
-        const amrex::Real& time) override;
-
 protected:
     std::string name() override { return "ABL test"; }
 

diff --git a/Exec/DevTests/MiguelDev/prob.cpp b/Exec/DevTests/MiguelDev/prob.cpp
@@ -140,31 +140,3 @@ Problem::init_custom_pert(
   });
 
 }
-
-void
-Problem::init_custom_terrain(
-    const Geometry& /*geom*/,
-    MultiFab& z_phys_nd,
-    const Real& /*time*/)
-{
-    // Number of ghost cells
-    int ngrow = z_phys_nd.nGrow();
-
-    // Bottom of domain
-    int k0 = 0;
-
-    for ( MFIter mfi(z_phys_nd, TilingIfNotGPU()); mfi.isValid(); ++mfi )
-    {
-        // Grown box with no z range
-        Box xybx = mfi.growntilebox(ngrow);
-        xybx.setRange(2,0);
-
-        Array4<Real> const& z_arr = z_phys_nd.array(mfi);
-
-        ParallelFor(xybx, [=] AMREX_GPU_DEVICE (int i, int j, int) {
-
-            // Flat terrain with z = 0 at k = 0
-            z_arr(i,j,k0) = 0.0;
-        });
-    }
-}
diff --git a/Exec/DevTests/MultiBlock/prob.H b/Exec/DevTests/MultiBlock/prob.H
@@ -46,11 +46,6 @@ public:
         amrex::Array4<amrex::Real const> const& mf_v,
         const SolverChoice& sc) override;
 
-    void init_custom_terrain (
-        const amrex::Geometry& geom,
-        amrex::MultiFab& z_phys_nd,
-        const amrex::Real& time) override;
-
 protected:
     std::string name() override { return "MultiBlock"; }
 

diff --git a/Exec/DevTests/MultiBlock/prob.cpp b/Exec/DevTests/MultiBlock/prob.cpp
@@ -111,31 +111,3 @@ Problem::init_custom_pert(
 
   amrex::Gpu::streamSynchronize();
 }
-
-void
-Problem::init_custom_terrain(
-    const Geometry& /*geom*/,
-    MultiFab& z_phys_nd,
-    const Real& /*time*/)
-{
-    // Number of ghost cells
-    int ngrow = z_phys_nd.nGrow();
-
-    // Bottom of domain
-    int k0 = 0;
-
-    for ( MFIter mfi(z_phys_nd, TilingIfNotGPU()); mfi.isValid(); ++mfi )
-    {
-        // Grown box with no z range
-        amrex::Box xybx = mfi.growntilebox(ngrow);
-        xybx.setRange(2,0);
-
-        Array4<Real> const& z_arr = z_phys_nd.array(mfi);
-
-        ParallelFor(xybx, [=] AMREX_GPU_DEVICE (int i, int j, int) {
-
-            // Flat terrain with z = 0 at k = 0
-            z_arr(i,j,k0) = 0.0;
-        });
-    }
-}
diff --git a/Exec/LLJ/prob.H b/Exec/LLJ/prob.H
@@ -39,11 +39,6 @@ public:
         amrex::Array4<amrex::Real const> const& mf_v,
         const SolverChoice& sc) override;
 
-    void init_custom_terrain (
-        const amrex::Geometry& geom,
-        amrex::MultiFab& z_phys_nd,
-        const amrex::Real& time) override;
-
 protected:
     std::string name() override { return "Low-Level Jet"; }
 

diff --git a/Exec/LLJ/prob.cpp b/Exec/LLJ/prob.cpp
@@ -42,28 +42,3 @@ Problem::init_custom_pert(
 {
   amrex::Print() << "Dummy function..Needed for linking" << std::endl;
 }
-
-void
-Problem::init_custom_terrain(
-    const Geometry& /*geom*/,
-    MultiFab& z_phys_nd,
-    const Real& /*time*/)
-{
-    // Number of ghost cells
-    int ngrow = z_phys_nd.nGrow();
-
-    for ( MFIter mfi(z_phys_nd, TilingIfNotGPU()); mfi.isValid(); ++mfi )
-    {
-        // Grown box with no z range
-        amrex::Box xybx = mfi.growntilebox(ngrow);
-        xybx.setRange(2,0);
-
-        Array4<Real> z_arr = z_phys_nd.array(mfi);
-
-        ParallelFor(xybx, [=] AMREX_GPU_DEVICE (int i, int j, int) {
-
-            // Flat terrain with z = 0 at k = 0
-            z_arr(i,j,0) = 0.;
-        });
-    }
-}
diff --git a/Exec/Radiation/prob.H b/Exec/Radiation/prob.H
@@ -46,11 +46,6 @@ public:
         amrex::Array4<amrex::Real const> const& mf_v,
         const SolverChoice& sc) override;
 
-    void init_custom_terrain (
-        const amrex::Geometry& geom,
-        amrex::MultiFab& z_phys_nd,
-        const amrex::Real& time) override;
-
     void erf_init_rayleigh (
         amrex::Vector<amrex::Real>& tau,
         amrex::Vector<amrex::Real>& ubar,

diff --git a/Exec/Radiation/prob.cpp b/Exec/Radiation/prob.cpp
@@ -212,34 +212,6 @@ Problem::init_custom_pert(
   amrex::Gpu::streamSynchronize();
 }
 
-void
-Problem::init_custom_terrain(
-    const Geometry& /*geom*/,
-    MultiFab& z_phys_nd,
-    const Real& /*time*/)
-{
-    // Number of ghost cells
-    int ngrow = z_phys_nd.nGrow();
-
-    // Bottom of domain
-    int k0 = 0;
-
-    for ( MFIter mfi(z_phys_nd, TilingIfNotGPU()); mfi.isValid(); ++mfi )
-    {
-        // Grown box with no z range
-        amrex::Box xybx = mfi.growntilebox(ngrow);
-        xybx.setRange(2,0);
-
-        Array4<Real> const& z_arr = z_phys_nd.array(mfi);
-
-        ParallelFor(xybx, [=] AMREX_GPU_DEVICE (int i, int j, int) {
-
-            // Flat terrain with z = 0 at k = 0
-            z_arr(i,j,k0) = 0.0;
-        });
-    }
-}
-
 void
 Problem::erf_init_rayleigh(
     amrex::Vector<amrex::Real>& tau,

diff --git a/Exec/RegTests/Bubble/prob.H b/Exec/RegTests/Bubble/prob.H
@@ -62,11 +62,6 @@ public:
         amrex::Array4<amrex::Real const> const& mf_v,
         const SolverChoice& sc) override;
 
-    void init_custom_terrain (
-        const amrex::Geometry& geom,
-        amrex::MultiFab& z_phys_nd,
-        const amrex::Real& time) override;
-
 #include "Prob/init_rayleigh_damping.H"
 
 protected:

diff --git a/Exec/RegTests/Bubble/prob.cpp b/Exec/RegTests/Bubble/prob.cpp
@@ -209,31 +209,3 @@ Problem::init_custom_pert(
 
     amrex::Gpu::streamSynchronize();
 }
-
-void
-Problem::init_custom_terrain(
-    const Geometry& /*geom*/,
-    MultiFab& z_phys_nd,
-    const Real& /*time*/)
-{
-    // Number of ghost cells
-    int ngrow = z_phys_nd.nGrow();
-
-    // Bottom of domain
-    int k0 = 0;
-
-    for ( MFIter mfi(z_phys_nd, TilingIfNotGPU()); mfi.isValid(); ++mfi )
-    {
-        // Grown box with no z range
-        amrex::Box xybx = mfi.growntilebox(ngrow);
-        xybx.setRange(2,0);
-
-        Array4<Real> const& z_arr = z_phys_nd.array(mfi);
-
-        ParallelFor(xybx, [=] AMREX_GPU_DEVICE (int i, int j, int) {
-
-            // Flat terrain with z = 0 at k = 0
-            z_arr(i,j,k0) = 0.0;
-        });
-    }
-}
diff --git a/Exec/RegTests/CouetteFlow/prob.H b/Exec/RegTests/CouetteFlow/prob.H
@@ -41,11 +41,6 @@ public:
         amrex::Array4<amrex::Real const> const& mf_v,
         const SolverChoice& sc) override;
 
-    void init_custom_terrain (
-        const amrex::Geometry& geom,
-        amrex::MultiFab& z_phys_nd,
-        const amrex::Real& time) override;
-
 protected:
     std::string name() override { return "Couette Flow"; }
 

diff --git a/Exec/RegTests/CouetteFlow/prob.cpp b/Exec/RegTests/CouetteFlow/prob.cpp
@@ -75,28 +75,3 @@ Problem::init_custom_pert(
     z_vel(i, j, k) = parms.w_0;
   });
 }
-
-void
-Problem::init_custom_terrain(
-    const Geometry& /*geom*/,
-    MultiFab& z_phys_nd,
-    const Real& /*time*/)
-{
-    // Number of ghost cells
-    int ngrow = z_phys_nd.nGrow();
-
-    for ( MFIter mfi(z_phys_nd, TilingIfNotGPU()); mfi.isValid(); ++mfi )
-    {
-        // Grown box with no z range
-        amrex::Box xybx = mfi.growntilebox(ngrow);
-        xybx.setRange(2,0);
-
-        Array4<Real> z_arr = z_phys_nd.array(mfi);
-
-        ParallelFor(xybx, [=] AMREX_GPU_DEVICE (int i, int j, int) {
-
-            // Flat terrain with z = 0 at k = 0
-            z_arr(i,j,0) = 0.;
-        });
-    }
-}
diff --git a/Exec/RegTests/DensityCurrent/prob.H b/Exec/RegTests/DensityCurrent/prob.H
@@ -44,11 +44,6 @@ public:
         amrex::Array4<amrex::Real const> const& mf_v,
         const SolverChoice& sc) override;
 
-    void init_custom_terrain (
-        const amrex::Geometry& geom,
-        amrex::MultiFab& z_phys_nd,
-        const amrex::Real& time) override;
-
 protected:
     std::string name() override { return "Density Current"; }
 

diff --git a/Exec/RegTests/DensityCurrent/prob.cpp b/Exec/RegTests/DensityCurrent/prob.cpp
@@ -148,32 +148,3 @@ Problem::init_custom_pert(
 
   amrex::Gpu::streamSynchronize();
 }
-
-void
-Problem::init_custom_terrain(
-    const Geometry& /*geom*/,
-    MultiFab& z_phys_nd,
-    const Real& /*time*/)
-{
-    // Number of ghost cells
-    int ngrow = z_phys_nd.nGrow();
-
-    // Bottom of domain
-    int k0 = 0;
-
-    for ( MFIter mfi(z_phys_nd, TilingIfNotGPU()); mfi.isValid(); ++mfi )
-    {
-        // Grown box with no z range
-        amrex::Box xybx = mfi.growntilebox(ngrow);
-        xybx.setRange(2,0);
-
-        Array4<Real> const& z_arr = z_phys_nd.array(mfi);
-
-        ParallelFor(xybx, [=] AMREX_GPU_DEVICE (int i, int j, int) {
-
-            // Flat terrain with z = 0 at k = 0
-            z_arr(i,j,k0) = 0.0;
-        });
-    }
-}
-
diff --git a/Exec/RegTests/DynamicRefinement/prob.H b/Exec/RegTests/DynamicRefinement/prob.H
@@ -51,11 +51,6 @@ public:
         amrex::Array4<amrex::Real const> const& mf_v,
         const SolverChoice& sc) override;
 
-    void init_custom_terrain (
-        const amrex::Geometry& geom,
-        amrex::MultiFab& z_phys_nd,
-        const amrex::Real& time) override;
-
 protected:
     std::string name() override { return "Dynamic Refinement"; }