From: Ralf Hartmann Date: Thu, 15 Feb 2007 15:27:28 +0000 (+0000) Subject: Generalize GridIn::read_netcdf to multiple bmarker on zero plane. X-Git-Tag: v8.0.0~10528 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=e5fe5f21ce5db362a575df468809330d87fa1b7e;p=dealii.git Generalize GridIn::read_netcdf to multiple bmarker on zero plane. git-svn-id: https://svn.dealii.org/trunk@14475 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/deal.II/source/grid/grid_in.cc b/deal.II/deal.II/source/grid/grid_in.cc index dcfd5d6111..9f5d677396 100644 --- a/deal.II/deal.II/source/grid/grid_in.cc +++ b/deal.II/deal.II/source/grid/grid_in.cc @@ -770,7 +770,7 @@ void GridIn<2>::read_netcdf (const std::string &filename) AssertThrow(false, ExcNeedsNetCDF()); #else const unsigned int dim=2; - const bool output=false; + const bool output=true; Assert (tria != 0, ExcNoTriangulationSelected()); // this function assumes the TAU // grid format. @@ -969,42 +969,59 @@ void GridIn<2>::read_netcdf (const std::string &filename) vertices[i](1)=point_values[y2d][i]; } - // Run over all boundary quads - // until we find a quad in the - // point[coord]=0 plane - unsigned int quad0=0; - for (;quad0 zero_plane_markers; + for (unsigned int quad=0; quad