Update grid_stretching_ratio check for realistic input

ewquon · Nov 13, 2023 · 4cae5f5 · 4cae5f5
1 parent e79c8e1
commit 4cae5f5
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 4 deletions.
diff --git a/Docs/sphinx_doc/Inputs.rst b/Docs/sphinx_doc/Inputs.rst
@@ -183,7 +183,7 @@ List of Parameters
 | Parameter                     | Definition      | Acceptable      | Default     |
 |                               |                 | Values          |             |
 +===============================+=================+=================+=============+
-| **erf.grid_stretching_ratio** | scaling factor  | Real > 0        | 0 (no grid  |
+| **erf.grid_stretching_ratio** | scaling factor  | Real > 1        | 0 (no grid  |
 |                               | applied to      |                 | stretching) |
 |                               | delta z at each |                 |             |
 |                               | level           |                 |             |

diff --git a/Source/DataStructs/DataStruct.H b/Source/DataStructs/DataStruct.H
@@ -46,9 +46,11 @@ struct SolverChoice {
         pp.query("use_terrain", use_terrain);
 
         pp.query("grid_stretching_ratio", grid_stretching_ratio);
-        AMREX_ASSERT_WITH_MESSAGE((grid_stretching_ratio >= 0.),
-                                  "The grid stretching ratio must be greater than 0");
-        if (grid_stretching_ratio > 0) {
+        if (grid_stretching_ratio != 0) {
+            AMREX_ASSERT_WITH_MESSAGE((grid_stretching_ratio >= 1.),
+                                      "The grid stretching ratio must be greater than 1");
+        }
+        if (grid_stretching_ratio >= 1) {
             if (!use_terrain) {
                 amrex::Print() << "Turning terrain on to enable grid stretching" << std::endl;
                 use_terrain = true;