Commit 645e99d8 authored by David Deepwell's avatar David Deepwell

Fix topography error in wave_reader

The quadrature weights were found to be incorrect for mapped grids
because they were adjusted using the topography which was not
scaled properly. The topography needed to be adjusted to be relative
of the 'minimum z-value'. This is in quation marks because it is
not necessarily the absolute minimum, but rather the reference
This shift also fixed the BPE calculation to have the same reference
as the PE calculation.
parent 4b058551
......@@ -172,7 +172,9 @@ class userControl : public BaseCase {
// Define topography
topo = zg(all,0,0);
// and shift it to match the adjustment made in the PE calculation
// This is also necessary to adjust the voxel sizes properly
topo = zg(all,0,0) - MinZ;
/* Initialize velocities */
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment