Merge pull request #672 from MFraters/add_gwb-annulus_test_coverage

Add tester for gwb-grid annulus.
GeodynamicWorldBuilder · Feb 23, 2024 · 7b227fc · 7b227fc
2 parents 704e06d + 31bbc56
commit 7b227fc
Show file tree

Hide file tree

Showing 4 changed files with 95 additions and 0 deletions.
diff --git a/source/gwb-grid/main.cc b/source/gwb-grid/main.cc
@@ -853,6 +853,7 @@ int main(int argc, char **argv)
                   grid_x[counter] = std::cos(theta) * (inner_radius + zi);
                   grid_z[counter] = std::sin(theta) * (inner_radius + zi);
                   grid_depth[counter] = outer_radius - std::sqrt(grid_x[counter] * grid_x[counter] + grid_z[counter] * grid_z [counter]);
+                  grid_depth[counter] = (std::fabs(grid_depth[counter]) < 1e-8 ? 0 : grid_depth[counter]);
                   counter++;
                 }
             }

diff --git a/tests/gwb-grid/3d_annulus_plume.grid b/tests/gwb-grid/3d_annulus_plume.grid
@@ -0,0 +1,18 @@
+# output variables
+grid_type = annulus
+dim = 2
+compositions = 3
+vtu_output_format = ASCII
+
+# domain of the grid
+x_max = 25
+x_min = -25
+y_max = 25
+y_min = -25
+z_min = 4371000
+z_max = 6371000
+
+# grid properties
+n_cell_x = 8
+n_cell_y = 8
+n_cell_z = 8