update gwb-grid/main.cc

Author: MFraters

source/gwb-grid/main.cc

@@ -126,7 +126,7 @@ void filter_vtu_mesh(int dim,
                 // The following line checks if we are accessing the velocity data, which is
                 // currently stored as the third property after depth and temperature, and if so,
                 // we need to copy all 3 components.
-                if (d == 2)
+                if (d == 3)
                   {
                     for (unsigned int i=0; i<3; ++i)
                       output_data[d].push_back(input_data[d][src_vid*3+i]);
@@ -659,6 +659,7 @@ int main(int argc, char **argv)
             grid_y.resize(n_p);
 
           grid_depth_wrt_surface.resize(n_p);
+          grid_depth_wrt_reference.resize(n_p);
 
           // compute positions
           size_t counter = 0;
@@ -845,6 +846,7 @@ int main(int argc, char **argv)
           grid_x.resize(n_p);
           grid_z.resize(n_p);
           grid_depth_wrt_surface.resize(n_p);
+          grid_depth_wrt_reference.resize(n_p);
 
           size_t counter = 0;
           for (size_t j = 0; j <= n_cell_z; ++j)