From: hartmann Date: Thu, 15 Feb 2007 15:27:28 +0000 (+0000) Subject: Generalize GridIn::read_netcdf to multiple bmarker on zero plane. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=d3ad10bf0b1b84d869a0caf8dd6551df486b4297;p=dealii-svn.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