Fix project_boundary values for small values of function.

git-svn-id: https://svn.dealii.org/trunk@32024 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/doc/news/changes.h b/deal.II/doc/news/changes.h
index afa759f9093cf46d78d5b596e017a649152e47be..63dd4ad503e3441df5649d9b407d38209ed363af 100644 (file)
--- a/deal.II/doc/news/changes.h
+++ b/deal.II/doc/news/changes.h
@@ -59,6 +59,12 @@ inconvenience this causes.
 <h3>Specific improvements</h3>
 
 <ol>
+  <li> Fixed: VectorTools::project_boundary_values could not deal with
+  function values close to (but not exactly equal to) zero. This is now fixed.
+  <br>
+  (Martin Kronbichler, 2013/12/16)
+  </li>
+
   <li> New: It is now possible to select between different smoothers and coarse
   solvers in the Trilinos AMG preconditioners by a string to the smoother's name.
   <br>

diff --git a/deal.II/include/deal.II/numerics/vector_tools.templates.h b/deal.II/include/deal.II/numerics/vector_tools.templates.h
index b4cdccc0153ca696977a62e3e52c49c93657222e..757c2a35fe1eabd90bcba54736478d836ae13b95 100644 (file)
--- a/deal.II/include/deal.II/numerics/vector_tools.templates.h
+++ b/deal.II/include/deal.II/numerics/vector_tools.templates.h
@@ -2111,23 +2111,27 @@ namespace VectorTools
 
       Vector<double> boundary_projection (rhs.size());
 
-      // Allow for a maximum of 5*n
-      // steps to reduce the residual by
-      // 10^-12. n steps may not be
-      // sufficient, since roundoff
-      // errors may accumulate for badly
-      // conditioned matrices
-      ReductionControl        control(5*rhs.size(), 0., 1.e-12, false, false);
-      GrowingVectorMemory<> memory;
-      SolverCG<>              cg(control,memory);
-
-      PreconditionSSOR<> prec;
-      prec.initialize(mass_matrix, 1.2);
-      filtered_precondition.initialize(prec, true);
-      // solve
-      cg.solve (filtered_mass_matrix, boundary_projection, rhs, filtered_precondition);
-      filtered_precondition.apply_constraints(boundary_projection, true);
-      filtered_precondition.clear();
+      // cannot reduce residual in a useful way if we are close to the square
+      // root of the minimal double value
+      if (rhs.norm_sqr() < 1e28 * std::numeric_limits<double>::min())
+        boundary_projection = 0;
+      else
+        {
+          // Allow for a maximum of 5*n steps to reduce the residual by 10^-12. n
+          // steps may not be sufficient, since roundoff errors may accumulate for
+          // badly conditioned matrices
+          ReductionControl        control(5*rhs.size(), 0., 1.e-12, false, false);
+          GrowingVectorMemory<> memory;
+          SolverCG<>              cg(control,memory);
+
+          PreconditionSSOR<> prec;
+          prec.initialize(mass_matrix, 1.2);
+          filtered_precondition.initialize(prec, true);
+          // solve
+          cg.solve (filtered_mass_matrix, boundary_projection, rhs, filtered_precondition);
+          filtered_precondition.apply_constraints(boundary_projection, true);
+          filtered_precondition.clear();
+        }
       // fill in boundary values
       for (unsigned int i=0; i<dof_to_boundary_mapping.size(); ++i)
         if (dof_to_boundary_mapping[i] != DoFHandler<dim,spacedim>::invalid_dof_index

diff --git a/tests/deal.II/assemble_matrix_parallel_01.cc b/tests/deal.II/assemble_matrix_parallel_01.cc
index 6d842edf1bd971c79857662257ad898916895f41..69f8fefa6dc9324413cbc171d97c3b513379b4d2 100644 (file)
--- a/tests/deal.II/assemble_matrix_parallel_01.cc
+++ b/tests/deal.II/assemble_matrix_parallel_01.cc
@@ -319,6 +319,9 @@ LaplaceProblem<dim>::local_assemble (const typename hp::DoFHandler<dim>::active_
 
   data.local_dof_indices.resize (dofs_per_cell);
   cell->get_dof_indices (data.local_dof_indices);
+  constraints.distribute_local_to_global(data.local_matrix, data.local_rhs,
+                                         data.local_dof_indices,
+                                         test_matrix, test_rhs);
 }
 
 
@@ -327,9 +330,6 @@ template <int dim>
 void
 LaplaceProblem<dim>::copy_local_to_global (const Assembly::Copy::Data &data)
 {
-  constraints.distribute_local_to_global(data.local_matrix, data.local_rhs,
-                                         data.local_dof_indices,
-                                         test_matrix, test_rhs);
 }
 
 

diff --git a/tests/deal.II/project_boundary_01.cc b/tests/deal.II/project_boundary_01.cc
new file mode 100644 (file)
index 0000000..9ffc067
--- /dev/null
+++ b/tests/deal.II/project_boundary_01.cc
@@ -0,0 +1,194 @@
+// ---------------------------------------------------------------------
+// $Id$
+//
+// Copyright (C) 2000 - 2013 by the deal.II authors
+//
+// This file is part of the deal.II library.
+//
+// The deal.II library is free software; you can use it, redistribute
+// it, and/or modify it under the terms of the GNU Lesser General
+// Public License as published by the Free Software Foundation; either
+// version 2.1 of the License, or (at your option) any later version.
+// The full text of the license can be found in the file LICENSE at
+// the top level of the deal.II distribution.
+//
+// ---------------------------------------------------------------------
+
+
+
+/* Purpose: check projection of boundary values for really small functions,
+   otherwise similar to boundaries.cc. */
+
+
+
+#include "../tests.h"
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/function_lib.h>
+#include <deal.II/base/quadrature_lib.h>
+#include <deal.II/lac/vector.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/tria_iterator.h>
+#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/dofs/dof_accessor.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/dofs/dof_tools.h>
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/fe_system.h>
+#include <deal.II/fe/mapping_q.h>
+#include <deal.II/numerics/vector_tools.h>
+
+#include <fstream>
+
+
+template<int dim>
+class MySquareFunction : public Function<dim>
+{
+public:
+  MySquareFunction (const unsigned int n_components) :
+    Function<dim>(n_components),
+    scaling(1)
+  {}
+
+  void set_scaling (const double scaling)
+  {
+    this->scaling = scaling;
+  }
+
+  virtual double value (const Point<dim>   &p,
+                        const unsigned int  component) const
+  {
+    return 100*(component+1)*p.square()*std::sin(p.square()) * scaling;
+  }
+
+  virtual void   vector_value (const Point<dim>   &p,
+                               Vector<double>     &values) const
+  {
+    for (unsigned int d=0; d<this->n_components; ++d) values(d) = value(p,d);
+  }
+
+private:
+  double scaling;
+};
+
+
+template <int dim>
+const Quadrature<dim-1> &
+boundary_q (const DoFHandler<dim> &)
+{
+  static const QGauss<dim-1> q(4);
+  return q;
+}
+
+
+const Quadrature<0> &
+boundary_q (const DoFHandler<1> &)
+{
+  static const Quadrature<0> *q = 0;
+  return *q;
+}
+
+
+void write_map (const std::map<types::global_dof_index,double> &bv)
+{
+  for (std::map<types::global_dof_index,double>::const_iterator
+       i=bv.begin(); i!=bv.end(); ++i)
+    // also output log of value to also display small numbers
+    deallog << i->first << ' ' << i->second << " "
+            << (std::abs(i->second)>0 ? std::log(std::abs(i->second)) : -10000 )
+            << std::endl;
+}
+
+
+
+
+template <int dim>
+void
+check ()
+{
+  Triangulation<dim> tr;
+  if (dim==2)
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    }
+  else
+    GridGenerator::hyper_cube(tr, -1./std::sqrt(static_cast<double>(dim)),1./std::sqrt(static_cast<double>(dim)));
+  static const HyperBallBoundary<dim> boundary;
+  if (dim != 1)
+    {
+      tr.set_boundary (0, boundary);
+    }
+  tr.refine_global (1);
+
+  MappingQ<dim> mapping(3);
+
+
+  // list of finite elements for which we want check, and associated list of
+  // boundary value functions
+  std::vector<const FiniteElement<dim>*> fe_list;
+  std::vector<MySquareFunction<dim>*> function_list;
+  
+  // FE1: a system of a quadratic and a linear element
+  fe_list.push_back (new FESystem<dim> (FE_Q<dim>(2), 1, FE_Q<dim>(1), 1));
+  function_list.push_back (new MySquareFunction<dim>(2));
+
+  // FE2: a linear element, to make things simple
+  fe_list.push_back (new FE_Q<dim> (1));
+  function_list.push_back (new MySquareFunction<dim>(1));
+
+  // test four different cases for the parameter: 1, 1e-20, 1e-170, 1e-800 (= 0)
+  double factors [] = {1., 1e-40, 1e-170, 1e-800};
+  for (unsigned int it=0; it<4; ++it)
+    {
+      deallog.push(Utilities::int_to_string(it,1));
+
+      // check all of them
+      for (unsigned int i=0; i<fe_list.size(); ++i)
+        {
+          function_list[i]->set_scaling(factors[it]);
+          const FiniteElement<dim> &fe = *fe_list[i];
+
+          DoFHandler<dim> dof(tr);
+          dof.distribute_dofs(fe);
+
+          typename FunctionMap<dim>::type function_map;
+          function_map[0] = function_list[i];
+
+          // interpolate boundary values
+          deallog << "Interpolated boundary values" << std::endl;
+          std::map<types::global_dof_index,double> interpolated_bv;
+          VectorTools::interpolate_boundary_values (mapping, dof, function_map,
+                                                    interpolated_bv, std::vector<bool>());
+          write_map (interpolated_bv);
+
+          deallog << "Projected boundary values" << std::endl;
+          std::map<types::global_dof_index,double> projected_bv;
+          VectorTools::project_boundary_values (mapping, dof, function_map,
+                                                boundary_q(dof), projected_bv);
+          write_map (projected_bv);
+        };
+      deallog.pop();
+    }
+
+  // delete objects now no more needed
+  for (unsigned int i=0; i<fe_list.size(); ++i)
+    {
+      delete fe_list[i];
+      delete function_list[i];
+    };
+}
+
+
+int main ()
+{
+  initlog();
+  deallog << std::setprecision (3);
+  deallog << std::fixed;
+
+  deallog.push ("2d");
+  check<2> ();
+  deallog.pop ();
+  deallog.push ("3d");
+  check<3> ();
+  deallog.pop ();
+}

diff --git a/tests/deal.II/project_boundary_01.output b/tests/deal.II/project_boundary_01.output
new file mode 100644 (file)
index 0000000..5b37f60
--- /dev/null
+++ b/tests/deal.II/project_boundary_01.output
@@ -0,0 +1,1489 @@
+
+DEAL:2d:0::Interpolated boundary values
+DEAL:2d:0::0 84.147 4.433
+DEAL:2d:0::1 168.294 5.126
+DEAL:2d:0::2 84.147 4.433
+DEAL:2d:0::3 168.294 5.126
+DEAL:2d:0::10 84.093 4.432
+DEAL:2d:0::13 84.147 4.433
+DEAL:2d:0::14 168.294 5.126
+DEAL:2d:0::18 84.093 4.432
+DEAL:2d:0::34 84.147 4.433
+DEAL:2d:0::35 168.294 5.126
+DEAL:2d:0::38 84.093 4.432
+DEAL:2d:0::47 84.147 4.433
+DEAL:2d:0::48 168.294 5.126
+DEAL:2d:0::51 84.093 4.432
+DEAL:2d:0::80 84.147 4.433
+DEAL:2d:0::81 168.294 5.126
+DEAL:2d:0::85 84.093 4.432
+DEAL:2d:0::88 84.147 4.433
+DEAL:2d:0::89 168.294 5.126
+DEAL:2d:0::93 84.093 4.432
+DEAL:2d:0::100 84.147 4.433
+DEAL:2d:0::101 168.294 5.126
+DEAL:2d:0::104 84.093 4.432
+DEAL:2d:0::110 84.093 4.432
+DEAL:2d:0::Projected boundary values
+DEAL:2d:0::0 84.207 4.433
+DEAL:2d:0::1 168.293 5.126
+DEAL:2d:0::2 84.207 4.433
+DEAL:2d:0::3 168.293 5.126
+DEAL:2d:0::10 84.116 4.432
+DEAL:2d:0::13 84.207 4.433
+DEAL:2d:0::14 168.293 5.126
+DEAL:2d:0::18 84.116 4.432
+DEAL:2d:0::34 84.207 4.433
+DEAL:2d:0::35 168.293 5.126
+DEAL:2d:0::38 84.116 4.432
+DEAL:2d:0::47 84.207 4.433
+DEAL:2d:0::48 168.293 5.126
+DEAL:2d:0::51 84.116 4.432
+DEAL:2d:0::80 84.207 4.433
+DEAL:2d:0::81 168.293 5.126
+DEAL:2d:0::85 84.116 4.432
+DEAL:2d:0::88 84.207 4.433
+DEAL:2d:0::89 168.293 5.126
+DEAL:2d:0::93 84.116 4.432
+DEAL:2d:0::100 84.207 4.433
+DEAL:2d:0::101 168.293 5.126
+DEAL:2d:0::104 84.116 4.432
+DEAL:2d:0::110 84.116 4.432
+DEAL:2d:0::Interpolated boundary values
+DEAL:2d:0::0 84.147 4.433
+DEAL:2d:0::1 84.147 4.433
+DEAL:2d:0::4 84.147 4.433
+DEAL:2d:0::9 84.147 4.433
+DEAL:2d:0::12 84.147 4.433
+DEAL:2d:0::19 84.147 4.433
+DEAL:2d:0::21 84.147 4.433
+DEAL:2d:0::23 84.147 4.433
+DEAL:2d:0::Projected boundary values
+DEAL:2d:0::0 84.147 4.433
+DEAL:2d:0::1 84.147 4.433
+DEAL:2d:0::4 84.147 4.433
+DEAL:2d:0::9 84.147 4.433
+DEAL:2d:0::12 84.147 4.433
+DEAL:2d:0::19 84.147 4.433
+DEAL:2d:0::21 84.147 4.433
+DEAL:2d:0::23 84.147 4.433
+DEAL:2d:1::Interpolated boundary values
+DEAL:2d:1::0 0.000 -87.671
+DEAL:2d:1::1 0.000 -86.978
+DEAL:2d:1::2 0.000 -87.671
+DEAL:2d:1::3 0.000 -86.978
+DEAL:2d:1::10 0.000 -87.671
+DEAL:2d:1::13 0.000 -87.671
+DEAL:2d:1::14 0.000 -86.978
+DEAL:2d:1::18 0.000 -87.671
+DEAL:2d:1::34 0.000 -87.671
+DEAL:2d:1::35 0.000 -86.978
+DEAL:2d:1::38 0.000 -87.671
+DEAL:2d:1::47 0.000 -87.671
+DEAL:2d:1::48 0.000 -86.978
+DEAL:2d:1::51 0.000 -87.671
+DEAL:2d:1::80 0.000 -87.671
+DEAL:2d:1::81 0.000 -86.978
+DEAL:2d:1::85 0.000 -87.671
+DEAL:2d:1::88 0.000 -87.671
+DEAL:2d:1::89 0.000 -86.978
+DEAL:2d:1::93 0.000 -87.671
+DEAL:2d:1::100 0.000 -87.671
+DEAL:2d:1::101 0.000 -86.978
+DEAL:2d:1::104 0.000 -87.671
+DEAL:2d:1::110 0.000 -87.671
+DEAL:2d:1::Projected boundary values
+DEAL:2d:1::0 0.000 -87.670
+DEAL:2d:1::1 0.000 -86.978
+DEAL:2d:1::2 0.000 -87.670
+DEAL:2d:1::3 0.000 -86.978
+DEAL:2d:1::10 0.000 -87.671
+DEAL:2d:1::13 0.000 -87.670
+DEAL:2d:1::14 0.000 -86.978
+DEAL:2d:1::18 0.000 -87.671
+DEAL:2d:1::34 0.000 -87.670
+DEAL:2d:1::35 0.000 -86.978
+DEAL:2d:1::38 0.000 -87.671
+DEAL:2d:1::47 0.000 -87.670
+DEAL:2d:1::48 0.000 -86.978
+DEAL:2d:1::51 0.000 -87.671
+DEAL:2d:1::80 0.000 -87.670
+DEAL:2d:1::81 0.000 -86.978
+DEAL:2d:1::85 0.000 -87.671
+DEAL:2d:1::88 0.000 -87.670
+DEAL:2d:1::89 0.000 -86.978
+DEAL:2d:1::93 0.000 -87.671
+DEAL:2d:1::100 0.000 -87.670
+DEAL:2d:1::101 0.000 -86.978
+DEAL:2d:1::104 0.000 -87.671
+DEAL:2d:1::110 0.000 -87.671
+DEAL:2d:1::Interpolated boundary values
+DEAL:2d:1::0 0.000 -87.671
+DEAL:2d:1::1 0.000 -87.671
+DEAL:2d:1::4 0.000 -87.671
+DEAL:2d:1::9 0.000 -87.671
+DEAL:2d:1::12 0.000 -87.671
+DEAL:2d:1::19 0.000 -87.671
+DEAL:2d:1::21 0.000 -87.671
+DEAL:2d:1::23 0.000 -87.671
+DEAL:2d:1::Projected boundary values
+DEAL:2d:1::0 0.000 -87.671
+DEAL:2d:1::1 0.000 -87.671
+DEAL:2d:1::4 0.000 -87.671
+DEAL:2d:1::9 0.000 -87.671
+DEAL:2d:1::12 0.000 -87.671
+DEAL:2d:1::19 0.000 -87.671
+DEAL:2d:1::21 0.000 -87.671
+DEAL:2d:1::23 0.000 -87.671
+DEAL:2d:2::Interpolated boundary values
+DEAL:2d:2::0 0.000 -387.007
+DEAL:2d:2::1 0.000 -386.314
+DEAL:2d:2::2 0.000 -387.007
+DEAL:2d:2::3 0.000 -386.314
+DEAL:2d:2::10 0.000 -387.008
+DEAL:2d:2::13 0.000 -387.007
+DEAL:2d:2::14 0.000 -386.314
+DEAL:2d:2::18 0.000 -387.008
+DEAL:2d:2::34 0.000 -387.007
+DEAL:2d:2::35 0.000 -386.314
+DEAL:2d:2::38 0.000 -387.008
+DEAL:2d:2::47 0.000 -387.007
+DEAL:2d:2::48 0.000 -386.314
+DEAL:2d:2::51 0.000 -387.008
+DEAL:2d:2::80 0.000 -387.007
+DEAL:2d:2::81 0.000 -386.314
+DEAL:2d:2::85 0.000 -387.008
+DEAL:2d:2::88 0.000 -387.007
+DEAL:2d:2::89 0.000 -386.314
+DEAL:2d:2::93 0.000 -387.008
+DEAL:2d:2::100 0.000 -387.007
+DEAL:2d:2::101 0.000 -386.314
+DEAL:2d:2::104 0.000 -387.008
+DEAL:2d:2::110 0.000 -387.008
+DEAL:2d:2::Projected boundary values
+DEAL:2d:2::0 0.000 -10000.000
+DEAL:2d:2::1 0.000 -10000.000
+DEAL:2d:2::2 0.000 -10000.000
+DEAL:2d:2::3 0.000 -10000.000
+DEAL:2d:2::10 0.000 -10000.000
+DEAL:2d:2::13 0.000 -10000.000
+DEAL:2d:2::14 0.000 -10000.000
+DEAL:2d:2::18 0.000 -10000.000
+DEAL:2d:2::34 0.000 -10000.000
+DEAL:2d:2::35 0.000 -10000.000
+DEAL:2d:2::38 0.000 -10000.000
+DEAL:2d:2::47 0.000 -10000.000
+DEAL:2d:2::48 0.000 -10000.000
+DEAL:2d:2::51 0.000 -10000.000
+DEAL:2d:2::80 0.000 -10000.000
+DEAL:2d:2::81 0.000 -10000.000
+DEAL:2d:2::85 0.000 -10000.000
+DEAL:2d:2::88 0.000 -10000.000
+DEAL:2d:2::89 0.000 -10000.000
+DEAL:2d:2::93 0.000 -10000.000
+DEAL:2d:2::100 0.000 -10000.000
+DEAL:2d:2::101 0.000 -10000.000
+DEAL:2d:2::104 0.000 -10000.000
+DEAL:2d:2::110 0.000 -10000.000
+DEAL:2d:2::Interpolated boundary values
+DEAL:2d:2::0 0.000 -387.007
+DEAL:2d:2::1 0.000 -387.007
+DEAL:2d:2::4 0.000 -387.007
+DEAL:2d:2::9 0.000 -387.007
+DEAL:2d:2::12 0.000 -387.007
+DEAL:2d:2::19 0.000 -387.007
+DEAL:2d:2::21 0.000 -387.007
+DEAL:2d:2::23 0.000 -387.007
+DEAL:2d:2::Projected boundary values
+DEAL:2d:2::0 0.000 -10000.000
+DEAL:2d:2::1 0.000 -10000.000
+DEAL:2d:2::4 0.000 -10000.000
+DEAL:2d:2::9 0.000 -10000.000
+DEAL:2d:2::12 0.000 -10000.000
+DEAL:2d:2::19 0.000 -10000.000
+DEAL:2d:2::21 0.000 -10000.000
+DEAL:2d:2::23 0.000 -10000.000
+DEAL:2d:3::Interpolated boundary values
+DEAL:2d:3::0 0.000 -10000.000
+DEAL:2d:3::1 0.000 -10000.000
+DEAL:2d:3::2 0.000 -10000.000
+DEAL:2d:3::3 0.000 -10000.000
+DEAL:2d:3::10 0.000 -10000.000
+DEAL:2d:3::13 0.000 -10000.000
+DEAL:2d:3::14 0.000 -10000.000
+DEAL:2d:3::18 0.000 -10000.000
+DEAL:2d:3::34 0.000 -10000.000
+DEAL:2d:3::35 0.000 -10000.000
+DEAL:2d:3::38 0.000 -10000.000
+DEAL:2d:3::47 0.000 -10000.000
+DEAL:2d:3::48 0.000 -10000.000
+DEAL:2d:3::51 0.000 -10000.000
+DEAL:2d:3::80 0.000 -10000.000
+DEAL:2d:3::81 0.000 -10000.000
+DEAL:2d:3::85 0.000 -10000.000
+DEAL:2d:3::88 0.000 -10000.000
+DEAL:2d:3::89 0.000 -10000.000
+DEAL:2d:3::93 0.000 -10000.000
+DEAL:2d:3::100 0.000 -10000.000
+DEAL:2d:3::101 0.000 -10000.000
+DEAL:2d:3::104 0.000 -10000.000
+DEAL:2d:3::110 0.000 -10000.000
+DEAL:2d:3::Projected boundary values
+DEAL:2d:3::0 0.000 -10000.000
+DEAL:2d:3::1 0.000 -10000.000
+DEAL:2d:3::2 0.000 -10000.000
+DEAL:2d:3::3 0.000 -10000.000
+DEAL:2d:3::10 0.000 -10000.000
+DEAL:2d:3::13 0.000 -10000.000
+DEAL:2d:3::14 0.000 -10000.000
+DEAL:2d:3::18 0.000 -10000.000
+DEAL:2d:3::34 0.000 -10000.000
+DEAL:2d:3::35 0.000 -10000.000
+DEAL:2d:3::38 0.000 -10000.000
+DEAL:2d:3::47 0.000 -10000.000
+DEAL:2d:3::48 0.000 -10000.000
+DEAL:2d:3::51 0.000 -10000.000
+DEAL:2d:3::80 0.000 -10000.000
+DEAL:2d:3::81 0.000 -10000.000
+DEAL:2d:3::85 0.000 -10000.000
+DEAL:2d:3::88 0.000 -10000.000
+DEAL:2d:3::89 0.000 -10000.000
+DEAL:2d:3::93 0.000 -10000.000
+DEAL:2d:3::100 0.000 -10000.000
+DEAL:2d:3::101 0.000 -10000.000
+DEAL:2d:3::104 0.000 -10000.000
+DEAL:2d:3::110 0.000 -10000.000
+DEAL:2d:3::Interpolated boundary values
+DEAL:2d:3::0 0.000 -10000.000
+DEAL:2d:3::1 0.000 -10000.000
+DEAL:2d:3::4 0.000 -10000.000
+DEAL:2d:3::9 0.000 -10000.000
+DEAL:2d:3::12 0.000 -10000.000
+DEAL:2d:3::19 0.000 -10000.000
+DEAL:2d:3::21 0.000 -10000.000
+DEAL:2d:3::23 0.000 -10000.000
+DEAL:2d:3::Projected boundary values
+DEAL:2d:3::0 0.000 -10000.000
+DEAL:2d:3::1 0.000 -10000.000
+DEAL:2d:3::4 0.000 -10000.000
+DEAL:2d:3::9 0.000 -10000.000
+DEAL:2d:3::12 0.000 -10000.000
+DEAL:2d:3::19 0.000 -10000.000
+DEAL:2d:3::21 0.000 -10000.000
+DEAL:2d:3::23 0.000 -10000.000
+DEAL:3d:0::Interpolated boundary values
+DEAL:3d:0::0 84.147 4.433
+DEAL:3d:0::1 168.294 5.126
+DEAL:3d:0::2 84.147 4.433
+DEAL:3d:0::3 168.294 5.126
+DEAL:3d:0::4 84.147 4.433
+DEAL:3d:0::5 168.294 5.126
+DEAL:3d:0::6 84.147 4.433
+DEAL:3d:0::7 168.294 5.126
+DEAL:3d:0::8 84.147 4.433
+DEAL:3d:0::9 168.294 5.126
+DEAL:3d:0::10 84.147 4.433
+DEAL:3d:0::11 168.294 5.126
+DEAL:3d:0::12 84.147 4.433
+DEAL:3d:0::13 168.294 5.126
+DEAL:3d:0::16 84.127 4.432
+DEAL:3d:0::17 84.093 4.432
+DEAL:3d:0::18 84.127 4.432
+DEAL:3d:0::19 84.093 4.432
+DEAL:3d:0::20 84.093 4.432
+DEAL:3d:0::22 84.093 4.432
+DEAL:3d:0::24 84.127 4.432
+DEAL:3d:0::25 84.093 4.432
+DEAL:3d:0::26 84.093 4.432
+DEAL:3d:0::28 84.027 4.431
+DEAL:3d:0::30 84.027 4.431
+DEAL:3d:0::32 84.027 4.431
+DEAL:3d:0::35 84.147 4.433
+DEAL:3d:0::36 168.294 5.126
+DEAL:3d:0::37 84.147 4.433
+DEAL:3d:0::38 168.294 5.126
+DEAL:3d:0::39 84.147 4.433
+DEAL:3d:0::40 168.294 5.126
+DEAL:3d:0::41 84.147 4.433
+DEAL:3d:0::42 168.294 5.126
+DEAL:3d:0::43 84.127 4.432
+DEAL:3d:0::44 84.127 4.432
+DEAL:3d:0::45 84.093 4.432
+DEAL:3d:0::46 84.093 4.432
+DEAL:3d:0::47 84.093 4.432
+DEAL:3d:0::49 84.127 4.432
+DEAL:3d:0::50 84.093 4.432
+DEAL:3d:0::51 84.027 4.431
+DEAL:3d:0::52 84.027 4.431
+DEAL:3d:0::54 84.027 4.431
+DEAL:3d:0::57 84.147 4.433
+DEAL:3d:0::58 168.294 5.126
+DEAL:3d:0::59 84.147 4.433
+DEAL:3d:0::60 168.294 5.126
+DEAL:3d:0::61 84.147 4.433
+DEAL:3d:0::62 168.294 5.126
+DEAL:3d:0::63 84.147 4.433
+DEAL:3d:0::64 168.294 5.126
+DEAL:3d:0::65 84.127 4.432
+DEAL:3d:0::66 84.093 4.432
+DEAL:3d:0::67 84.127 4.432
+DEAL:3d:0::68 84.093 4.432
+DEAL:3d:0::70 84.093 4.432
+DEAL:3d:0::71 84.127 4.432
+DEAL:3d:0::72 84.093 4.432
+DEAL:3d:0::73 84.027 4.431
+DEAL:3d:0::75 84.027 4.431
+DEAL:3d:0::76 84.027 4.431
+DEAL:3d:0::79 84.147 4.433
+DEAL:3d:0::80 168.294 5.126
+DEAL:3d:0::81 84.147 4.433
+DEAL:3d:0::82 168.294 5.126
+DEAL:3d:0::83 84.127 4.432
+DEAL:3d:0::84 84.127 4.432
+DEAL:3d:0::85 84.093 4.432
+DEAL:3d:0::86 84.093 4.432
+DEAL:3d:0::87 84.127 4.432
+DEAL:3d:0::88 84.027 4.431
+DEAL:3d:0::89 84.027 4.431
+DEAL:3d:0::90 84.027 4.431
+DEAL:3d:0::93 84.147 4.433
+DEAL:3d:0::94 168.294 5.126
+DEAL:3d:0::95 84.147 4.433
+DEAL:3d:0::96 168.294 5.126
+DEAL:3d:0::97 84.147 4.433
+DEAL:3d:0::98 168.294 5.126
+DEAL:3d:0::99 84.147 4.433
+DEAL:3d:0::100 168.294 5.126
+DEAL:3d:0::101 84.127 4.432
+DEAL:3d:0::102 84.093 4.432
+DEAL:3d:0::103 84.127 4.432
+DEAL:3d:0::104 84.093 4.432
+DEAL:3d:0::105 84.127 4.432
+DEAL:3d:0::106 84.093 4.432
+DEAL:3d:0::107 84.093 4.432
+DEAL:3d:0::109 84.027 4.431
+DEAL:3d:0::111 84.027 4.431
+DEAL:3d:0::113 84.027 4.431
+DEAL:3d:0::115 84.147 4.433
+DEAL:3d:0::116 168.294 5.126
+DEAL:3d:0::117 84.147 4.433
+DEAL:3d:0::118 168.294 5.126
+DEAL:3d:0::119 84.127 4.432
+DEAL:3d:0::120 84.127 4.432
+DEAL:3d:0::121 84.093 4.432
+DEAL:3d:0::122 84.127 4.432
+DEAL:3d:0::123 84.093 4.432
+DEAL:3d:0::124 84.027 4.431
+DEAL:3d:0::125 84.027 4.431
+DEAL:3d:0::127 84.027 4.431
+DEAL:3d:0::129 84.147 4.433
+DEAL:3d:0::130 168.294 5.126
+DEAL:3d:0::131 84.147 4.433
+DEAL:3d:0::132 168.294 5.126
+DEAL:3d:0::133 84.127 4.432
+DEAL:3d:0::134 84.093 4.432
+DEAL:3d:0::135 84.127 4.432
+DEAL:3d:0::136 84.127 4.432
+DEAL:3d:0::137 84.093 4.432
+DEAL:3d:0::138 84.027 4.431
+DEAL:3d:0::140 84.027 4.431
+DEAL:3d:0::141 84.027 4.431
+DEAL:3d:0::143 84.147 4.433
+DEAL:3d:0::144 168.294 5.126
+DEAL:3d:0::145 84.127 4.432
+DEAL:3d:0::146 84.127 4.432
+DEAL:3d:0::147 84.127 4.432
+DEAL:3d:0::148 84.027 4.431
+DEAL:3d:0::149 84.027 4.431
+DEAL:3d:0::150 84.027 4.431
+DEAL:3d:0::Projected boundary values
+DEAL:3d:0::0 84.189 4.433
+DEAL:3d:0::1 168.282 5.126
+DEAL:3d:0::2 84.230 4.434
+DEAL:3d:0::3 168.288 5.126
+DEAL:3d:0::4 84.230 4.434
+DEAL:3d:0::5 168.288 5.126
+DEAL:3d:0::6 84.270 4.434
+DEAL:3d:0::7 168.301 5.126
+DEAL:3d:0::8 84.230 4.434
+DEAL:3d:0::9 168.288 5.126
+DEAL:3d:0::10 84.270 4.434
+DEAL:3d:0::11 168.301 5.126
+DEAL:3d:0::12 84.270 4.434
+DEAL:3d:0::13 168.301 5.126
+DEAL:3d:0::16 84.204 4.433
+DEAL:3d:0::17 84.169 4.433
+DEAL:3d:0::18 84.204 4.433
+DEAL:3d:0::19 84.192 4.433
+DEAL:3d:0::20 84.192 4.433
+DEAL:3d:0::22 84.169 4.433
+DEAL:3d:0::24 84.204 4.433
+DEAL:3d:0::25 84.192 4.433
+DEAL:3d:0::26 84.169 4.433
+DEAL:3d:0::28 84.080 4.432
+DEAL:3d:0::30 84.080 4.432
+DEAL:3d:0::32 84.080 4.432
+DEAL:3d:0::35 84.189 4.433
+DEAL:3d:0::36 168.282 5.126
+DEAL:3d:0::37 84.230 4.434
+DEAL:3d:0::38 168.288 5.126
+DEAL:3d:0::39 84.230 4.434
+DEAL:3d:0::40 168.288 5.126
+DEAL:3d:0::41 84.270 4.434
+DEAL:3d:0::42 168.301 5.126
+DEAL:3d:0::43 84.204 4.433
+DEAL:3d:0::44 84.204 4.433
+DEAL:3d:0::45 84.192 4.433
+DEAL:3d:0::46 84.192 4.433
+DEAL:3d:0::47 84.169 4.433
+DEAL:3d:0::49 84.204 4.433
+DEAL:3d:0::50 84.169 4.433
+DEAL:3d:0::51 84.080 4.432
+DEAL:3d:0::52 84.080 4.432
+DEAL:3d:0::54 84.080 4.432
+DEAL:3d:0::57 84.189 4.433
+DEAL:3d:0::58 168.282 5.126
+DEAL:3d:0::59 84.230 4.434
+DEAL:3d:0::60 168.288 5.126
+DEAL:3d:0::61 84.230 4.434
+DEAL:3d:0::62 168.288 5.126
+DEAL:3d:0::63 84.270 4.434
+DEAL:3d:0::64 168.301 5.126
+DEAL:3d:0::65 84.204 4.433
+DEAL:3d:0::66 84.169 4.433
+DEAL:3d:0::67 84.204 4.433
+DEAL:3d:0::68 84.192 4.433
+DEAL:3d:0::70 84.169 4.433
+DEAL:3d:0::71 84.204 4.433
+DEAL:3d:0::72 84.192 4.433
+DEAL:3d:0::73 84.080 4.432
+DEAL:3d:0::75 84.080 4.432
+DEAL:3d:0::76 84.080 4.432
+DEAL:3d:0::79 84.189 4.433
+DEAL:3d:0::80 168.282 5.126
+DEAL:3d:0::81 84.230 4.434
+DEAL:3d:0::82 168.288 5.126
+DEAL:3d:0::83 84.204 4.433
+DEAL:3d:0::84 84.204 4.433
+DEAL:3d:0::85 84.192 4.433
+DEAL:3d:0::86 84.169 4.433
+DEAL:3d:0::87 84.204 4.433
+DEAL:3d:0::88 84.080 4.432
+DEAL:3d:0::89 84.080 4.432
+DEAL:3d:0::90 84.080 4.432
+DEAL:3d:0::93 84.189 4.433
+DEAL:3d:0::94 168.282 5.126
+DEAL:3d:0::95 84.230 4.434
+DEAL:3d:0::96 168.288 5.126
+DEAL:3d:0::97 84.230 4.434
+DEAL:3d:0::98 168.288 5.126
+DEAL:3d:0::99 84.270 4.434
+DEAL:3d:0::100 168.301 5.126
+DEAL:3d:0::101 84.204 4.433
+DEAL:3d:0::102 84.169 4.433
+DEAL:3d:0::103 84.204 4.433
+DEAL:3d:0::104 84.192 4.433
+DEAL:3d:0::105 84.204 4.433
+DEAL:3d:0::106 84.192 4.433
+DEAL:3d:0::107 84.169 4.433
+DEAL:3d:0::109 84.080 4.432
+DEAL:3d:0::111 84.080 4.432
+DEAL:3d:0::113 84.080 4.432
+DEAL:3d:0::115 84.189 4.433
+DEAL:3d:0::116 168.282 5.126
+DEAL:3d:0::117 84.230 4.434
+DEAL:3d:0::118 168.288 5.126
+DEAL:3d:0::119 84.204 4.433
+DEAL:3d:0::120 84.204 4.433
+DEAL:3d:0::121 84.192 4.433
+DEAL:3d:0::122 84.204 4.433
+DEAL:3d:0::123 84.169 4.433
+DEAL:3d:0::124 84.080 4.432
+DEAL:3d:0::125 84.080 4.432
+DEAL:3d:0::127 84.080 4.432
+DEAL:3d:0::129 84.189 4.433
+DEAL:3d:0::130 168.282 5.126
+DEAL:3d:0::131 84.230 4.434
+DEAL:3d:0::132 168.288 5.126
+DEAL:3d:0::133 84.204 4.433
+DEAL:3d:0::134 84.169 4.433
+DEAL:3d:0::135 84.204 4.433
+DEAL:3d:0::136 84.204 4.433
+DEAL:3d:0::137 84.192 4.433
+DEAL:3d:0::138 84.080 4.432
+DEAL:3d:0::140 84.080 4.432
+DEAL:3d:0::141 84.080 4.432
+DEAL:3d:0::143 84.189 4.433
+DEAL:3d:0::144 168.282 5.126
+DEAL:3d:0::145 84.204 4.433
+DEAL:3d:0::146 84.204 4.433
+DEAL:3d:0::147 84.204 4.433
+DEAL:3d:0::148 84.080 4.432
+DEAL:3d:0::149 84.080 4.432
+DEAL:3d:0::150 84.080 4.432
+DEAL:3d:0::Interpolated boundary values
+DEAL:3d:0::0 84.147 4.433
+DEAL:3d:0::1 84.147 4.433
+DEAL:3d:0::2 84.147 4.433
+DEAL:3d:0::3 84.147 4.433
+DEAL:3d:0::4 84.147 4.433
+DEAL:3d:0::5 84.147 4.433
+DEAL:3d:0::6 84.147 4.433
+DEAL:3d:0::8 84.147 4.433
+DEAL:3d:0::9 84.147 4.433
+DEAL:3d:0::10 84.147 4.433
+DEAL:3d:0::11 84.147 4.433
+DEAL:3d:0::12 84.147 4.433
+DEAL:3d:0::13 84.147 4.433
+DEAL:3d:0::14 84.147 4.433
+DEAL:3d:0::15 84.147 4.433
+DEAL:3d:0::16 84.147 4.433
+DEAL:3d:0::17 84.147 4.433
+DEAL:3d:0::18 84.147 4.433
+DEAL:3d:0::19 84.147 4.433
+DEAL:3d:0::20 84.147 4.433
+DEAL:3d:0::21 84.147 4.433
+DEAL:3d:0::22 84.147 4.433
+DEAL:3d:0::23 84.147 4.433
+DEAL:3d:0::24 84.147 4.433
+DEAL:3d:0::25 84.147 4.433
+DEAL:3d:0::26 84.147 4.433
+DEAL:3d:0::Projected boundary values
+DEAL:3d:0::0 84.141 4.432
+DEAL:3d:0::1 84.144 4.433
+DEAL:3d:0::2 84.144 4.433
+DEAL:3d:0::3 84.151 4.433
+DEAL:3d:0::4 84.144 4.433
+DEAL:3d:0::5 84.151 4.433
+DEAL:3d:0::6 84.151 4.433
+DEAL:3d:0::8 84.141 4.432
+DEAL:3d:0::9 84.144 4.433
+DEAL:3d:0::10 84.144 4.433
+DEAL:3d:0::11 84.151 4.433
+DEAL:3d:0::12 84.141 4.432
+DEAL:3d:0::13 84.144 4.433
+DEAL:3d:0::14 84.144 4.433
+DEAL:3d:0::15 84.151 4.433
+DEAL:3d:0::16 84.141 4.432
+DEAL:3d:0::17 84.144 4.433
+DEAL:3d:0::18 84.141 4.432
+DEAL:3d:0::19 84.144 4.433
+DEAL:3d:0::20 84.144 4.433
+DEAL:3d:0::21 84.151 4.433
+DEAL:3d:0::22 84.141 4.432
+DEAL:3d:0::23 84.144 4.433
+DEAL:3d:0::24 84.141 4.432
+DEAL:3d:0::25 84.144 4.433
+DEAL:3d:0::26 84.141 4.432
+DEAL:3d:1::Interpolated boundary values
+DEAL:3d:1::0 0.000 -87.671
+DEAL:3d:1::1 0.000 -86.978
+DEAL:3d:1::2 0.000 -87.671
+DEAL:3d:1::3 0.000 -86.978
+DEAL:3d:1::4 0.000 -87.671
+DEAL:3d:1::5 0.000 -86.978
+DEAL:3d:1::6 0.000 -87.671
+DEAL:3d:1::7 0.000 -86.978
+DEAL:3d:1::8 0.000 -87.671
+DEAL:3d:1::9 0.000 -86.978
+DEAL:3d:1::10 0.000 -87.671
+DEAL:3d:1::11 0.000 -86.978
+DEAL:3d:1::12 0.000 -87.671
+DEAL:3d:1::13 0.000 -86.978
+DEAL:3d:1::16 0.000 -87.671
+DEAL:3d:1::17 0.000 -87.671
+DEAL:3d:1::18 0.000 -87.671
+DEAL:3d:1::19 0.000 -87.671
+DEAL:3d:1::20 0.000 -87.671
+DEAL:3d:1::22 0.000 -87.671
+DEAL:3d:1::24 0.000 -87.671
+DEAL:3d:1::25 0.000 -87.671
+DEAL:3d:1::26 0.000 -87.671
+DEAL:3d:1::28 0.000 -87.672
+DEAL:3d:1::30 0.000 -87.672
+DEAL:3d:1::32 0.000 -87.672
+DEAL:3d:1::35 0.000 -87.671
+DEAL:3d:1::36 0.000 -86.978
+DEAL:3d:1::37 0.000 -87.671
+DEAL:3d:1::38 0.000 -86.978
+DEAL:3d:1::39 0.000 -87.671
+DEAL:3d:1::40 0.000 -86.978
+DEAL:3d:1::41 0.000 -87.671
+DEAL:3d:1::42 0.000 -86.978
+DEAL:3d:1::43 0.000 -87.671
+DEAL:3d:1::44 0.000 -87.671
+DEAL:3d:1::45 0.000 -87.671
+DEAL:3d:1::46 0.000 -87.671
+DEAL:3d:1::47 0.000 -87.671
+DEAL:3d:1::49 0.000 -87.671
+DEAL:3d:1::50 0.000 -87.671
+DEAL:3d:1::51 0.000 -87.672
+DEAL:3d:1::52 0.000 -87.672
+DEAL:3d:1::54 0.000 -87.672
+DEAL:3d:1::57 0.000 -87.671
+DEAL:3d:1::58 0.000 -86.978
+DEAL:3d:1::59 0.000 -87.671
+DEAL:3d:1::60 0.000 -86.978
+DEAL:3d:1::61 0.000 -87.671
+DEAL:3d:1::62 0.000 -86.978
+DEAL:3d:1::63 0.000 -87.671
+DEAL:3d:1::64 0.000 -86.978
+DEAL:3d:1::65 0.000 -87.671
+DEAL:3d:1::66 0.000 -87.671
+DEAL:3d:1::67 0.000 -87.671
+DEAL:3d:1::68 0.000 -87.671
+DEAL:3d:1::70 0.000 -87.671
+DEAL:3d:1::71 0.000 -87.671
+DEAL:3d:1::72 0.000 -87.671
+DEAL:3d:1::73 0.000 -87.672
+DEAL:3d:1::75 0.000 -87.672
+DEAL:3d:1::76 0.000 -87.672
+DEAL:3d:1::79 0.000 -87.671
+DEAL:3d:1::80 0.000 -86.978
+DEAL:3d:1::81 0.000 -87.671
+DEAL:3d:1::82 0.000 -86.978
+DEAL:3d:1::83 0.000 -87.671
+DEAL:3d:1::84 0.000 -87.671
+DEAL:3d:1::85 0.000 -87.671
+DEAL:3d:1::86 0.000 -87.671
+DEAL:3d:1::87 0.000 -87.671
+DEAL:3d:1::88 0.000 -87.672
+DEAL:3d:1::89 0.000 -87.672
+DEAL:3d:1::90 0.000 -87.672
+DEAL:3d:1::93 0.000 -87.671
+DEAL:3d:1::94 0.000 -86.978
+DEAL:3d:1::95 0.000 -87.671
+DEAL:3d:1::96 0.000 -86.978
+DEAL:3d:1::97 0.000 -87.671
+DEAL:3d:1::98 0.000 -86.978
+DEAL:3d:1::99 0.000 -87.671
+DEAL:3d:1::100 0.000 -86.978
+DEAL:3d:1::101 0.000 -87.671
+DEAL:3d:1::102 0.000 -87.671
+DEAL:3d:1::103 0.000 -87.671
+DEAL:3d:1::104 0.000 -87.671
+DEAL:3d:1::105 0.000 -87.671
+DEAL:3d:1::106 0.000 -87.671
+DEAL:3d:1::107 0.000 -87.671
+DEAL:3d:1::109 0.000 -87.672
+DEAL:3d:1::111 0.000 -87.672
+DEAL:3d:1::113 0.000 -87.672
+DEAL:3d:1::115 0.000 -87.671
+DEAL:3d:1::116 0.000 -86.978
+DEAL:3d:1::117 0.000 -87.671
+DEAL:3d:1::118 0.000 -86.978
+DEAL:3d:1::119 0.000 -87.671
+DEAL:3d:1::120 0.000 -87.671
+DEAL:3d:1::121 0.000 -87.671
+DEAL:3d:1::122 0.000 -87.671
+DEAL:3d:1::123 0.000 -87.671
+DEAL:3d:1::124 0.000 -87.672
+DEAL:3d:1::125 0.000 -87.672
+DEAL:3d:1::127 0.000 -87.672
+DEAL:3d:1::129 0.000 -87.671
+DEAL:3d:1::130 0.000 -86.978
+DEAL:3d:1::131 0.000 -87.671
+DEAL:3d:1::132 0.000 -86.978
+DEAL:3d:1::133 0.000 -87.671
+DEAL:3d:1::134 0.000 -87.671
+DEAL:3d:1::135 0.000 -87.671
+DEAL:3d:1::136 0.000 -87.671
+DEAL:3d:1::137 0.000 -87.671
+DEAL:3d:1::138 0.000 -87.672
+DEAL:3d:1::140 0.000 -87.672
+DEAL:3d:1::141 0.000 -87.672
+DEAL:3d:1::143 0.000 -87.671
+DEAL:3d:1::144 0.000 -86.978
+DEAL:3d:1::145 0.000 -87.671
+DEAL:3d:1::146 0.000 -87.671
+DEAL:3d:1::147 0.000 -87.671
+DEAL:3d:1::148 0.000 -87.672
+DEAL:3d:1::149 0.000 -87.672
+DEAL:3d:1::150 0.000 -87.672
+DEAL:3d:1::Projected boundary values
+DEAL:3d:1::0 0.000 -87.670
+DEAL:3d:1::1 0.000 -86.978
+DEAL:3d:1::2 0.000 -87.670
+DEAL:3d:1::3 0.000 -86.978
+DEAL:3d:1::4 0.000 -87.670
+DEAL:3d:1::5 0.000 -86.978
+DEAL:3d:1::6 0.000 -87.669
+DEAL:3d:1::7 0.000 -86.978
+DEAL:3d:1::8 0.000 -87.670
+DEAL:3d:1::9 0.000 -86.978
+DEAL:3d:1::10 0.000 -87.669
+DEAL:3d:1::11 0.000 -86.978
+DEAL:3d:1::12 0.000 -87.669
+DEAL:3d:1::13 0.000 -86.978
+DEAL:3d:1::16 0.000 -87.670
+DEAL:3d:1::17 0.000 -87.671
+DEAL:3d:1::18 0.000 -87.670
+DEAL:3d:1::19 0.000 -87.670
+DEAL:3d:1::20 0.000 -87.670
+DEAL:3d:1::22 0.000 -87.671
+DEAL:3d:1::24 0.000 -87.670
+DEAL:3d:1::25 0.000 -87.670
+DEAL:3d:1::26 0.000 -87.671
+DEAL:3d:1::28 0.000 -87.672
+DEAL:3d:1::30 0.000 -87.672
+DEAL:3d:1::32 0.000 -87.672
+DEAL:3d:1::35 0.000 -87.670
+DEAL:3d:1::36 0.000 -86.978
+DEAL:3d:1::37 0.000 -87.670
+DEAL:3d:1::38 0.000 -86.978
+DEAL:3d:1::39 0.000 -87.670
+DEAL:3d:1::40 0.000 -86.978
+DEAL:3d:1::41 0.000 -87.669
+DEAL:3d:1::42 0.000 -86.978
+DEAL:3d:1::43 0.000 -87.670
+DEAL:3d:1::44 0.000 -87.670
+DEAL:3d:1::45 0.000 -87.670
+DEAL:3d:1::46 0.000 -87.670
+DEAL:3d:1::47 0.000 -87.671
+DEAL:3d:1::49 0.000 -87.670
+DEAL:3d:1::50 0.000 -87.671
+DEAL:3d:1::51 0.000 -87.672
+DEAL:3d:1::52 0.000 -87.672
+DEAL:3d:1::54 0.000 -87.672
+DEAL:3d:1::57 0.000 -87.670
+DEAL:3d:1::58 0.000 -86.978
+DEAL:3d:1::59 0.000 -87.670
+DEAL:3d:1::60 0.000 -86.978
+DEAL:3d:1::61 0.000 -87.670
+DEAL:3d:1::62 0.000 -86.978
+DEAL:3d:1::63 0.000 -87.669
+DEAL:3d:1::64 0.000 -86.978
+DEAL:3d:1::65 0.000 -87.670
+DEAL:3d:1::66 0.000 -87.671
+DEAL:3d:1::67 0.000 -87.670
+DEAL:3d:1::68 0.000 -87.670
+DEAL:3d:1::70 0.000 -87.671
+DEAL:3d:1::71 0.000 -87.670
+DEAL:3d:1::72 0.000 -87.670
+DEAL:3d:1::73 0.000 -87.672
+DEAL:3d:1::75 0.000 -87.672
+DEAL:3d:1::76 0.000 -87.672
+DEAL:3d:1::79 0.000 -87.670
+DEAL:3d:1::80 0.000 -86.978
+DEAL:3d:1::81 0.000 -87.670
+DEAL:3d:1::82 0.000 -86.978
+DEAL:3d:1::83 0.000 -87.670
+DEAL:3d:1::84 0.000 -87.670
+DEAL:3d:1::85 0.000 -87.670
+DEAL:3d:1::86 0.000 -87.671
+DEAL:3d:1::87 0.000 -87.670
+DEAL:3d:1::88 0.000 -87.672
+DEAL:3d:1::89 0.000 -87.672
+DEAL:3d:1::90 0.000 -87.672
+DEAL:3d:1::93 0.000 -87.670
+DEAL:3d:1::94 0.000 -86.978
+DEAL:3d:1::95 0.000 -87.670
+DEAL:3d:1::96 0.000 -86.978
+DEAL:3d:1::97 0.000 -87.670
+DEAL:3d:1::98 0.000 -86.978
+DEAL:3d:1::99 0.000 -87.669
+DEAL:3d:1::100 0.000 -86.978
+DEAL:3d:1::101 0.000 -87.670
+DEAL:3d:1::102 0.000 -87.671
+DEAL:3d:1::103 0.000 -87.670
+DEAL:3d:1::104 0.000 -87.670
+DEAL:3d:1::105 0.000 -87.670
+DEAL:3d:1::106 0.000 -87.670
+DEAL:3d:1::107 0.000 -87.671
+DEAL:3d:1::109 0.000 -87.672
+DEAL:3d:1::111 0.000 -87.672
+DEAL:3d:1::113 0.000 -87.672
+DEAL:3d:1::115 0.000 -87.670
+DEAL:3d:1::116 0.000 -86.978
+DEAL:3d:1::117 0.000 -87.670
+DEAL:3d:1::118 0.000 -86.978
+DEAL:3d:1::119 0.000 -87.670
+DEAL:3d:1::120 0.000 -87.670
+DEAL:3d:1::121 0.000 -87.670
+DEAL:3d:1::122 0.000 -87.670
+DEAL:3d:1::123 0.000 -87.671
+DEAL:3d:1::124 0.000 -87.672
+DEAL:3d:1::125 0.000 -87.672
+DEAL:3d:1::127 0.000 -87.672
+DEAL:3d:1::129 0.000 -87.670
+DEAL:3d:1::130 0.000 -86.978
+DEAL:3d:1::131 0.000 -87.670
+DEAL:3d:1::132 0.000 -86.978
+DEAL:3d:1::133 0.000 -87.670
+DEAL:3d:1::134 0.000 -87.671
+DEAL:3d:1::135 0.000 -87.670
+DEAL:3d:1::136 0.000 -87.670
+DEAL:3d:1::137 0.000 -87.670
+DEAL:3d:1::138 0.000 -87.672
+DEAL:3d:1::140 0.000 -87.672
+DEAL:3d:1::141 0.000 -87.672
+DEAL:3d:1::143 0.000 -87.670
+DEAL:3d:1::144 0.000 -86.978
+DEAL:3d:1::145 0.000 -87.670
+DEAL:3d:1::146 0.000 -87.670
+DEAL:3d:1::147 0.000 -87.670
+DEAL:3d:1::148 0.000 -87.672
+DEAL:3d:1::149 0.000 -87.672
+DEAL:3d:1::150 0.000 -87.672
+DEAL:3d:1::Interpolated boundary values
+DEAL:3d:1::0 0.000 -87.671
+DEAL:3d:1::1 0.000 -87.671
+DEAL:3d:1::2 0.000 -87.671
+DEAL:3d:1::3 0.000 -87.671
+DEAL:3d:1::4 0.000 -87.671
+DEAL:3d:1::5 0.000 -87.671
+DEAL:3d:1::6 0.000 -87.671
+DEAL:3d:1::8 0.000 -87.671
+DEAL:3d:1::9 0.000 -87.671
+DEAL:3d:1::10 0.000 -87.671
+DEAL:3d:1::11 0.000 -87.671
+DEAL:3d:1::12 0.000 -87.671
+DEAL:3d:1::13 0.000 -87.671
+DEAL:3d:1::14 0.000 -87.671
+DEAL:3d:1::15 0.000 -87.671
+DEAL:3d:1::16 0.000 -87.671
+DEAL:3d:1::17 0.000 -87.671
+DEAL:3d:1::18 0.000 -87.671
+DEAL:3d:1::19 0.000 -87.671
+DEAL:3d:1::20 0.000 -87.671
+DEAL:3d:1::21 0.000 -87.671
+DEAL:3d:1::22 0.000 -87.671
+DEAL:3d:1::23 0.000 -87.671
+DEAL:3d:1::24 0.000 -87.671
+DEAL:3d:1::25 0.000 -87.671
+DEAL:3d:1::26 0.000 -87.671
+DEAL:3d:1::Projected boundary values
+DEAL:3d:1::0 0.000 -87.671
+DEAL:3d:1::1 0.000 -87.671
+DEAL:3d:1::2 0.000 -87.671
+DEAL:3d:1::3 0.000 -87.671
+DEAL:3d:1::4 0.000 -87.671
+DEAL:3d:1::5 0.000 -87.671
+DEAL:3d:1::6 0.000 -87.671
+DEAL:3d:1::8 0.000 -87.671
+DEAL:3d:1::9 0.000 -87.671
+DEAL:3d:1::10 0.000 -87.671
+DEAL:3d:1::11 0.000 -87.671
+DEAL:3d:1::12 0.000 -87.671
+DEAL:3d:1::13 0.000 -87.671
+DEAL:3d:1::14 0.000 -87.671
+DEAL:3d:1::15 0.000 -87.671
+DEAL:3d:1::16 0.000 -87.671
+DEAL:3d:1::17 0.000 -87.671
+DEAL:3d:1::18 0.000 -87.671
+DEAL:3d:1::19 0.000 -87.671
+DEAL:3d:1::20 0.000 -87.671
+DEAL:3d:1::21 0.000 -87.671
+DEAL:3d:1::22 0.000 -87.671
+DEAL:3d:1::23 0.000 -87.671
+DEAL:3d:1::24 0.000 -87.671
+DEAL:3d:1::25 0.000 -87.671
+DEAL:3d:1::26 0.000 -87.671
+DEAL:3d:2::Interpolated boundary values
+DEAL:3d:2::0 0.000 -387.007
+DEAL:3d:2::1 0.000 -386.314
+DEAL:3d:2::2 0.000 -387.007
+DEAL:3d:2::3 0.000 -386.314
+DEAL:3d:2::4 0.000 -387.007
+DEAL:3d:2::5 0.000 -386.314
+DEAL:3d:2::6 0.000 -387.007
+DEAL:3d:2::7 0.000 -386.314
+DEAL:3d:2::8 0.000 -387.007
+DEAL:3d:2::9 0.000 -386.314
+DEAL:3d:2::10 0.000 -387.007
+DEAL:3d:2::11 0.000 -386.314
+DEAL:3d:2::12 0.000 -387.007
+DEAL:3d:2::13 0.000 -386.314
+DEAL:3d:2::16 0.000 -387.007
+DEAL:3d:2::17 0.000 -387.008
+DEAL:3d:2::18 0.000 -387.007
+DEAL:3d:2::19 0.000 -387.008
+DEAL:3d:2::20 0.000 -387.008
+DEAL:3d:2::22 0.000 -387.008
+DEAL:3d:2::24 0.000 -387.007
+DEAL:3d:2::25 0.000 -387.008
+DEAL:3d:2::26 0.000 -387.008
+DEAL:3d:2::28 0.000 -387.008
+DEAL:3d:2::30 0.000 -387.008
+DEAL:3d:2::32 0.000 -387.008
+DEAL:3d:2::35 0.000 -387.007
+DEAL:3d:2::36 0.000 -386.314
+DEAL:3d:2::37 0.000 -387.007
+DEAL:3d:2::38 0.000 -386.314
+DEAL:3d:2::39 0.000 -387.007
+DEAL:3d:2::40 0.000 -386.314
+DEAL:3d:2::41 0.000 -387.007
+DEAL:3d:2::42 0.000 -386.314
+DEAL:3d:2::43 0.000 -387.007
+DEAL:3d:2::44 0.000 -387.007
+DEAL:3d:2::45 0.000 -387.008
+DEAL:3d:2::46 0.000 -387.008
+DEAL:3d:2::47 0.000 -387.008
+DEAL:3d:2::49 0.000 -387.007
+DEAL:3d:2::50 0.000 -387.008
+DEAL:3d:2::51 0.000 -387.008
+DEAL:3d:2::52 0.000 -387.008
+DEAL:3d:2::54 0.000 -387.008
+DEAL:3d:2::57 0.000 -387.007
+DEAL:3d:2::58 0.000 -386.314
+DEAL:3d:2::59 0.000 -387.007
+DEAL:3d:2::60 0.000 -386.314
+DEAL:3d:2::61 0.000 -387.007
+DEAL:3d:2::62 0.000 -386.314
+DEAL:3d:2::63 0.000 -387.007
+DEAL:3d:2::64 0.000 -386.314
+DEAL:3d:2::65 0.000 -387.007
+DEAL:3d:2::66 0.000 -387.008
+DEAL:3d:2::67 0.000 -387.007
+DEAL:3d:2::68 0.000 -387.008
+DEAL:3d:2::70 0.000 -387.008
+DEAL:3d:2::71 0.000 -387.007
+DEAL:3d:2::72 0.000 -387.008
+DEAL:3d:2::73 0.000 -387.008
+DEAL:3d:2::75 0.000 -387.008
+DEAL:3d:2::76 0.000 -387.008
+DEAL:3d:2::79 0.000 -387.007
+DEAL:3d:2::80 0.000 -386.314
+DEAL:3d:2::81 0.000 -387.007
+DEAL:3d:2::82 0.000 -386.314
+DEAL:3d:2::83 0.000 -387.007
+DEAL:3d:2::84 0.000 -387.007
+DEAL:3d:2::85 0.000 -387.008
+DEAL:3d:2::86 0.000 -387.008
+DEAL:3d:2::87 0.000 -387.007
+DEAL:3d:2::88 0.000 -387.008
+DEAL:3d:2::89 0.000 -387.008
+DEAL:3d:2::90 0.000 -387.008
+DEAL:3d:2::93 0.000 -387.007
+DEAL:3d:2::94 0.000 -386.314
+DEAL:3d:2::95 0.000 -387.007
+DEAL:3d:2::96 0.000 -386.314
+DEAL:3d:2::97 0.000 -387.007
+DEAL:3d:2::98 0.000 -386.314
+DEAL:3d:2::99 0.000 -387.007
+DEAL:3d:2::100 0.000 -386.314
+DEAL:3d:2::101 0.000 -387.007
+DEAL:3d:2::102 0.000 -387.008
+DEAL:3d:2::103 0.000 -387.007
+DEAL:3d:2::104 0.000 -387.008
+DEAL:3d:2::105 0.000 -387.007
+DEAL:3d:2::106 0.000 -387.008
+DEAL:3d:2::107 0.000 -387.008
+DEAL:3d:2::109 0.000 -387.008
+DEAL:3d:2::111 0.000 -387.008
+DEAL:3d:2::113 0.000 -387.008
+DEAL:3d:2::115 0.000 -387.007
+DEAL:3d:2::116 0.000 -386.314
+DEAL:3d:2::117 0.000 -387.007
+DEAL:3d:2::118 0.000 -386.314
+DEAL:3d:2::119 0.000 -387.007
+DEAL:3d:2::120 0.000 -387.007
+DEAL:3d:2::121 0.000 -387.008
+DEAL:3d:2::122 0.000 -387.007
+DEAL:3d:2::123 0.000 -387.008
+DEAL:3d:2::124 0.000 -387.008
+DEAL:3d:2::125 0.000 -387.008
+DEAL:3d:2::127 0.000 -387.008
+DEAL:3d:2::129 0.000 -387.007
+DEAL:3d:2::130 0.000 -386.314
+DEAL:3d:2::131 0.000 -387.007
+DEAL:3d:2::132 0.000 -386.314
+DEAL:3d:2::133 0.000 -387.007
+DEAL:3d:2::134 0.000 -387.008
+DEAL:3d:2::135 0.000 -387.007
+DEAL:3d:2::136 0.000 -387.007
+DEAL:3d:2::137 0.000 -387.008
+DEAL:3d:2::138 0.000 -387.008
+DEAL:3d:2::140 0.000 -387.008
+DEAL:3d:2::141 0.000 -387.008
+DEAL:3d:2::143 0.000 -387.007
+DEAL:3d:2::144 0.000 -386.314
+DEAL:3d:2::145 0.000 -387.007
+DEAL:3d:2::146 0.000 -387.007
+DEAL:3d:2::147 0.000 -387.007
+DEAL:3d:2::148 0.000 -387.008
+DEAL:3d:2::149 0.000 -387.008
+DEAL:3d:2::150 0.000 -387.008
+DEAL:3d:2::Projected boundary values
+DEAL:3d:2::0 0.000 -10000.000
+DEAL:3d:2::1 0.000 -10000.000
+DEAL:3d:2::2 0.000 -10000.000
+DEAL:3d:2::3 0.000 -10000.000
+DEAL:3d:2::4 0.000 -10000.000
+DEAL:3d:2::5 0.000 -10000.000
+DEAL:3d:2::6 0.000 -10000.000
+DEAL:3d:2::7 0.000 -10000.000
+DEAL:3d:2::8 0.000 -10000.000
+DEAL:3d:2::9 0.000 -10000.000
+DEAL:3d:2::10 0.000 -10000.000
+DEAL:3d:2::11 0.000 -10000.000
+DEAL:3d:2::12 0.000 -10000.000
+DEAL:3d:2::13 0.000 -10000.000
+DEAL:3d:2::16 0.000 -10000.000
+DEAL:3d:2::17 0.000 -10000.000
+DEAL:3d:2::18 0.000 -10000.000
+DEAL:3d:2::19 0.000 -10000.000
+DEAL:3d:2::20 0.000 -10000.000
+DEAL:3d:2::22 0.000 -10000.000
+DEAL:3d:2::24 0.000 -10000.000
+DEAL:3d:2::25 0.000 -10000.000
+DEAL:3d:2::26 0.000 -10000.000
+DEAL:3d:2::28 0.000 -10000.000
+DEAL:3d:2::30 0.000 -10000.000
+DEAL:3d:2::32 0.000 -10000.000
+DEAL:3d:2::35 0.000 -10000.000
+DEAL:3d:2::36 0.000 -10000.000
+DEAL:3d:2::37 0.000 -10000.000
+DEAL:3d:2::38 0.000 -10000.000
+DEAL:3d:2::39 0.000 -10000.000
+DEAL:3d:2::40 0.000 -10000.000
+DEAL:3d:2::41 0.000 -10000.000
+DEAL:3d:2::42 0.000 -10000.000
+DEAL:3d:2::43 0.000 -10000.000
+DEAL:3d:2::44 0.000 -10000.000
+DEAL:3d:2::45 0.000 -10000.000
+DEAL:3d:2::46 0.000 -10000.000
+DEAL:3d:2::47 0.000 -10000.000
+DEAL:3d:2::49 0.000 -10000.000
+DEAL:3d:2::50 0.000 -10000.000
+DEAL:3d:2::51 0.000 -10000.000
+DEAL:3d:2::52 0.000 -10000.000
+DEAL:3d:2::54 0.000 -10000.000
+DEAL:3d:2::57 0.000 -10000.000
+DEAL:3d:2::58 0.000 -10000.000
+DEAL:3d:2::59 0.000 -10000.000
+DEAL:3d:2::60 0.000 -10000.000
+DEAL:3d:2::61 0.000 -10000.000
+DEAL:3d:2::62 0.000 -10000.000
+DEAL:3d:2::63 0.000 -10000.000
+DEAL:3d:2::64 0.000 -10000.000
+DEAL:3d:2::65 0.000 -10000.000
+DEAL:3d:2::66 0.000 -10000.000
+DEAL:3d:2::67 0.000 -10000.000
+DEAL:3d:2::68 0.000 -10000.000
+DEAL:3d:2::70 0.000 -10000.000
+DEAL:3d:2::71 0.000 -10000.000
+DEAL:3d:2::72 0.000 -10000.000
+DEAL:3d:2::73 0.000 -10000.000
+DEAL:3d:2::75 0.000 -10000.000
+DEAL:3d:2::76 0.000 -10000.000
+DEAL:3d:2::79 0.000 -10000.000
+DEAL:3d:2::80 0.000 -10000.000
+DEAL:3d:2::81 0.000 -10000.000
+DEAL:3d:2::82 0.000 -10000.000
+DEAL:3d:2::83 0.000 -10000.000
+DEAL:3d:2::84 0.000 -10000.000
+DEAL:3d:2::85 0.000 -10000.000
+DEAL:3d:2::86 0.000 -10000.000
+DEAL:3d:2::87 0.000 -10000.000
+DEAL:3d:2::88 0.000 -10000.000
+DEAL:3d:2::89 0.000 -10000.000
+DEAL:3d:2::90 0.000 -10000.000
+DEAL:3d:2::93 0.000 -10000.000
+DEAL:3d:2::94 0.000 -10000.000
+DEAL:3d:2::95 0.000 -10000.000
+DEAL:3d:2::96 0.000 -10000.000
+DEAL:3d:2::97 0.000 -10000.000
+DEAL:3d:2::98 0.000 -10000.000
+DEAL:3d:2::99 0.000 -10000.000
+DEAL:3d:2::100 0.000 -10000.000
+DEAL:3d:2::101 0.000 -10000.000
+DEAL:3d:2::102 0.000 -10000.000
+DEAL:3d:2::103 0.000 -10000.000
+DEAL:3d:2::104 0.000 -10000.000
+DEAL:3d:2::105 0.000 -10000.000
+DEAL:3d:2::106 0.000 -10000.000
+DEAL:3d:2::107 0.000 -10000.000
+DEAL:3d:2::109 0.000 -10000.000
+DEAL:3d:2::111 0.000 -10000.000
+DEAL:3d:2::113 0.000 -10000.000
+DEAL:3d:2::115 0.000 -10000.000
+DEAL:3d:2::116 0.000 -10000.000
+DEAL:3d:2::117 0.000 -10000.000
+DEAL:3d:2::118 0.000 -10000.000
+DEAL:3d:2::119 0.000 -10000.000
+DEAL:3d:2::120 0.000 -10000.000
+DEAL:3d:2::121 0.000 -10000.000
+DEAL:3d:2::122 0.000 -10000.000
+DEAL:3d:2::123 0.000 -10000.000
+DEAL:3d:2::124 0.000 -10000.000
+DEAL:3d:2::125 0.000 -10000.000
+DEAL:3d:2::127 0.000 -10000.000
+DEAL:3d:2::129 0.000 -10000.000
+DEAL:3d:2::130 0.000 -10000.000
+DEAL:3d:2::131 0.000 -10000.000
+DEAL:3d:2::132 0.000 -10000.000
+DEAL:3d:2::133 0.000 -10000.000
+DEAL:3d:2::134 0.000 -10000.000
+DEAL:3d:2::135 0.000 -10000.000
+DEAL:3d:2::136 0.000 -10000.000
+DEAL:3d:2::137 0.000 -10000.000
+DEAL:3d:2::138 0.000 -10000.000
+DEAL:3d:2::140 0.000 -10000.000
+DEAL:3d:2::141 0.000 -10000.000
+DEAL:3d:2::143 0.000 -10000.000
+DEAL:3d:2::144 0.000 -10000.000
+DEAL:3d:2::145 0.000 -10000.000
+DEAL:3d:2::146 0.000 -10000.000
+DEAL:3d:2::147 0.000 -10000.000
+DEAL:3d:2::148 0.000 -10000.000
+DEAL:3d:2::149 0.000 -10000.000
+DEAL:3d:2::150 0.000 -10000.000
+DEAL:3d:2::Interpolated boundary values
+DEAL:3d:2::0 0.000 -387.007
+DEAL:3d:2::1 0.000 -387.007
+DEAL:3d:2::2 0.000 -387.007
+DEAL:3d:2::3 0.000 -387.007
+DEAL:3d:2::4 0.000 -387.007
+DEAL:3d:2::5 0.000 -387.007
+DEAL:3d:2::6 0.000 -387.007
+DEAL:3d:2::8 0.000 -387.007
+DEAL:3d:2::9 0.000 -387.007
+DEAL:3d:2::10 0.000 -387.007
+DEAL:3d:2::11 0.000 -387.007
+DEAL:3d:2::12 0.000 -387.007
+DEAL:3d:2::13 0.000 -387.007
+DEAL:3d:2::14 0.000 -387.007
+DEAL:3d:2::15 0.000 -387.007
+DEAL:3d:2::16 0.000 -387.007
+DEAL:3d:2::17 0.000 -387.007
+DEAL:3d:2::18 0.000 -387.007
+DEAL:3d:2::19 0.000 -387.007
+DEAL:3d:2::20 0.000 -387.007
+DEAL:3d:2::21 0.000 -387.007
+DEAL:3d:2::22 0.000 -387.007
+DEAL:3d:2::23 0.000 -387.007
+DEAL:3d:2::24 0.000 -387.007
+DEAL:3d:2::25 0.000 -387.007
+DEAL:3d:2::26 0.000 -387.007
+DEAL:3d:2::Projected boundary values
+DEAL:3d:2::0 0.000 -10000.000
+DEAL:3d:2::1 0.000 -10000.000
+DEAL:3d:2::2 0.000 -10000.000
+DEAL:3d:2::3 0.000 -10000.000
+DEAL:3d:2::4 0.000 -10000.000
+DEAL:3d:2::5 0.000 -10000.000
+DEAL:3d:2::6 0.000 -10000.000
+DEAL:3d:2::8 0.000 -10000.000
+DEAL:3d:2::9 0.000 -10000.000
+DEAL:3d:2::10 0.000 -10000.000
+DEAL:3d:2::11 0.000 -10000.000
+DEAL:3d:2::12 0.000 -10000.000
+DEAL:3d:2::13 0.000 -10000.000
+DEAL:3d:2::14 0.000 -10000.000
+DEAL:3d:2::15 0.000 -10000.000
+DEAL:3d:2::16 0.000 -10000.000
+DEAL:3d:2::17 0.000 -10000.000
+DEAL:3d:2::18 0.000 -10000.000
+DEAL:3d:2::19 0.000 -10000.000
+DEAL:3d:2::20 0.000 -10000.000
+DEAL:3d:2::21 0.000 -10000.000
+DEAL:3d:2::22 0.000 -10000.000
+DEAL:3d:2::23 0.000 -10000.000
+DEAL:3d:2::24 0.000 -10000.000
+DEAL:3d:2::25 0.000 -10000.000
+DEAL:3d:2::26 0.000 -10000.000
+DEAL:3d:3::Interpolated boundary values
+DEAL:3d:3::0 0.000 -10000.000
+DEAL:3d:3::1 0.000 -10000.000
+DEAL:3d:3::2 0.000 -10000.000
+DEAL:3d:3::3 0.000 -10000.000
+DEAL:3d:3::4 0.000 -10000.000
+DEAL:3d:3::5 0.000 -10000.000
+DEAL:3d:3::6 0.000 -10000.000
+DEAL:3d:3::7 0.000 -10000.000
+DEAL:3d:3::8 0.000 -10000.000
+DEAL:3d:3::9 0.000 -10000.000
+DEAL:3d:3::10 0.000 -10000.000
+DEAL:3d:3::11 0.000 -10000.000
+DEAL:3d:3::12 0.000 -10000.000
+DEAL:3d:3::13 0.000 -10000.000
+DEAL:3d:3::16 0.000 -10000.000
+DEAL:3d:3::17 0.000 -10000.000
+DEAL:3d:3::18 0.000 -10000.000
+DEAL:3d:3::19 0.000 -10000.000
+DEAL:3d:3::20 0.000 -10000.000
+DEAL:3d:3::22 0.000 -10000.000
+DEAL:3d:3::24 0.000 -10000.000
+DEAL:3d:3::25 0.000 -10000.000
+DEAL:3d:3::26 0.000 -10000.000
+DEAL:3d:3::28 0.000 -10000.000
+DEAL:3d:3::30 0.000 -10000.000
+DEAL:3d:3::32 0.000 -10000.000
+DEAL:3d:3::35 0.000 -10000.000
+DEAL:3d:3::36 0.000 -10000.000
+DEAL:3d:3::37 0.000 -10000.000
+DEAL:3d:3::38 0.000 -10000.000
+DEAL:3d:3::39 0.000 -10000.000
+DEAL:3d:3::40 0.000 -10000.000
+DEAL:3d:3::41 0.000 -10000.000
+DEAL:3d:3::42 0.000 -10000.000
+DEAL:3d:3::43 0.000 -10000.000
+DEAL:3d:3::44 0.000 -10000.000
+DEAL:3d:3::45 0.000 -10000.000
+DEAL:3d:3::46 0.000 -10000.000
+DEAL:3d:3::47 0.000 -10000.000
+DEAL:3d:3::49 0.000 -10000.000
+DEAL:3d:3::50 0.000 -10000.000
+DEAL:3d:3::51 0.000 -10000.000
+DEAL:3d:3::52 0.000 -10000.000
+DEAL:3d:3::54 0.000 -10000.000
+DEAL:3d:3::57 0.000 -10000.000
+DEAL:3d:3::58 0.000 -10000.000
+DEAL:3d:3::59 0.000 -10000.000
+DEAL:3d:3::60 0.000 -10000.000
+DEAL:3d:3::61 0.000 -10000.000
+DEAL:3d:3::62 0.000 -10000.000
+DEAL:3d:3::63 0.000 -10000.000
+DEAL:3d:3::64 0.000 -10000.000
+DEAL:3d:3::65 0.000 -10000.000
+DEAL:3d:3::66 0.000 -10000.000
+DEAL:3d:3::67 0.000 -10000.000
+DEAL:3d:3::68 0.000 -10000.000
+DEAL:3d:3::70 0.000 -10000.000
+DEAL:3d:3::71 0.000 -10000.000
+DEAL:3d:3::72 0.000 -10000.000
+DEAL:3d:3::73 0.000 -10000.000
+DEAL:3d:3::75 0.000 -10000.000
+DEAL:3d:3::76 0.000 -10000.000
+DEAL:3d:3::79 0.000 -10000.000
+DEAL:3d:3::80 0.000 -10000.000
+DEAL:3d:3::81 0.000 -10000.000
+DEAL:3d:3::82 0.000 -10000.000
+DEAL:3d:3::83 0.000 -10000.000
+DEAL:3d:3::84 0.000 -10000.000
+DEAL:3d:3::85 0.000 -10000.000
+DEAL:3d:3::86 0.000 -10000.000
+DEAL:3d:3::87 0.000 -10000.000
+DEAL:3d:3::88 0.000 -10000.000
+DEAL:3d:3::89 0.000 -10000.000
+DEAL:3d:3::90 0.000 -10000.000
+DEAL:3d:3::93 0.000 -10000.000
+DEAL:3d:3::94 0.000 -10000.000
+DEAL:3d:3::95 0.000 -10000.000
+DEAL:3d:3::96 0.000 -10000.000
+DEAL:3d:3::97 0.000 -10000.000
+DEAL:3d:3::98 0.000 -10000.000
+DEAL:3d:3::99 0.000 -10000.000
+DEAL:3d:3::100 0.000 -10000.000
+DEAL:3d:3::101 0.000 -10000.000
+DEAL:3d:3::102 0.000 -10000.000
+DEAL:3d:3::103 0.000 -10000.000
+DEAL:3d:3::104 0.000 -10000.000
+DEAL:3d:3::105 0.000 -10000.000
+DEAL:3d:3::106 0.000 -10000.000
+DEAL:3d:3::107 0.000 -10000.000
+DEAL:3d:3::109 0.000 -10000.000
+DEAL:3d:3::111 0.000 -10000.000
+DEAL:3d:3::113 0.000 -10000.000
+DEAL:3d:3::115 0.000 -10000.000
+DEAL:3d:3::116 0.000 -10000.000
+DEAL:3d:3::117 0.000 -10000.000
+DEAL:3d:3::118 0.000 -10000.000
+DEAL:3d:3::119 0.000 -10000.000
+DEAL:3d:3::120 0.000 -10000.000
+DEAL:3d:3::121 0.000 -10000.000
+DEAL:3d:3::122 0.000 -10000.000
+DEAL:3d:3::123 0.000 -10000.000
+DEAL:3d:3::124 0.000 -10000.000
+DEAL:3d:3::125 0.000 -10000.000
+DEAL:3d:3::127 0.000 -10000.000
+DEAL:3d:3::129 0.000 -10000.000
+DEAL:3d:3::130 0.000 -10000.000
+DEAL:3d:3::131 0.000 -10000.000
+DEAL:3d:3::132 0.000 -10000.000
+DEAL:3d:3::133 0.000 -10000.000
+DEAL:3d:3::134 0.000 -10000.000
+DEAL:3d:3::135 0.000 -10000.000
+DEAL:3d:3::136 0.000 -10000.000
+DEAL:3d:3::137 0.000 -10000.000
+DEAL:3d:3::138 0.000 -10000.000
+DEAL:3d:3::140 0.000 -10000.000
+DEAL:3d:3::141 0.000 -10000.000
+DEAL:3d:3::143 0.000 -10000.000
+DEAL:3d:3::144 0.000 -10000.000
+DEAL:3d:3::145 0.000 -10000.000
+DEAL:3d:3::146 0.000 -10000.000
+DEAL:3d:3::147 0.000 -10000.000
+DEAL:3d:3::148 0.000 -10000.000
+DEAL:3d:3::149 0.000 -10000.000
+DEAL:3d:3::150 0.000 -10000.000
+DEAL:3d:3::Projected boundary values
+DEAL:3d:3::0 0.000 -10000.000
+DEAL:3d:3::1 0.000 -10000.000
+DEAL:3d:3::2 0.000 -10000.000
+DEAL:3d:3::3 0.000 -10000.000
+DEAL:3d:3::4 0.000 -10000.000
+DEAL:3d:3::5 0.000 -10000.000
+DEAL:3d:3::6 0.000 -10000.000
+DEAL:3d:3::7 0.000 -10000.000
+DEAL:3d:3::8 0.000 -10000.000
+DEAL:3d:3::9 0.000 -10000.000
+DEAL:3d:3::10 0.000 -10000.000
+DEAL:3d:3::11 0.000 -10000.000
+DEAL:3d:3::12 0.000 -10000.000
+DEAL:3d:3::13 0.000 -10000.000
+DEAL:3d:3::16 0.000 -10000.000
+DEAL:3d:3::17 0.000 -10000.000
+DEAL:3d:3::18 0.000 -10000.000
+DEAL:3d:3::19 0.000 -10000.000
+DEAL:3d:3::20 0.000 -10000.000
+DEAL:3d:3::22 0.000 -10000.000
+DEAL:3d:3::24 0.000 -10000.000
+DEAL:3d:3::25 0.000 -10000.000
+DEAL:3d:3::26 0.000 -10000.000
+DEAL:3d:3::28 0.000 -10000.000
+DEAL:3d:3::30 0.000 -10000.000
+DEAL:3d:3::32 0.000 -10000.000
+DEAL:3d:3::35 0.000 -10000.000
+DEAL:3d:3::36 0.000 -10000.000
+DEAL:3d:3::37 0.000 -10000.000
+DEAL:3d:3::38 0.000 -10000.000
+DEAL:3d:3::39 0.000 -10000.000
+DEAL:3d:3::40 0.000 -10000.000
+DEAL:3d:3::41 0.000 -10000.000
+DEAL:3d:3::42 0.000 -10000.000
+DEAL:3d:3::43 0.000 -10000.000
+DEAL:3d:3::44 0.000 -10000.000
+DEAL:3d:3::45 0.000 -10000.000
+DEAL:3d:3::46 0.000 -10000.000
+DEAL:3d:3::47 0.000 -10000.000
+DEAL:3d:3::49 0.000 -10000.000
+DEAL:3d:3::50 0.000 -10000.000
+DEAL:3d:3::51 0.000 -10000.000
+DEAL:3d:3::52 0.000 -10000.000
+DEAL:3d:3::54 0.000 -10000.000
+DEAL:3d:3::57 0.000 -10000.000
+DEAL:3d:3::58 0.000 -10000.000
+DEAL:3d:3::59 0.000 -10000.000
+DEAL:3d:3::60 0.000 -10000.000
+DEAL:3d:3::61 0.000 -10000.000
+DEAL:3d:3::62 0.000 -10000.000
+DEAL:3d:3::63 0.000 -10000.000
+DEAL:3d:3::64 0.000 -10000.000
+DEAL:3d:3::65 0.000 -10000.000
+DEAL:3d:3::66 0.000 -10000.000
+DEAL:3d:3::67 0.000 -10000.000
+DEAL:3d:3::68 0.000 -10000.000
+DEAL:3d:3::70 0.000 -10000.000
+DEAL:3d:3::71 0.000 -10000.000
+DEAL:3d:3::72 0.000 -10000.000
+DEAL:3d:3::73 0.000 -10000.000
+DEAL:3d:3::75 0.000 -10000.000
+DEAL:3d:3::76 0.000 -10000.000
+DEAL:3d:3::79 0.000 -10000.000
+DEAL:3d:3::80 0.000 -10000.000
+DEAL:3d:3::81 0.000 -10000.000
+DEAL:3d:3::82 0.000 -10000.000
+DEAL:3d:3::83 0.000 -10000.000
+DEAL:3d:3::84 0.000 -10000.000
+DEAL:3d:3::85 0.000 -10000.000
+DEAL:3d:3::86 0.000 -10000.000
+DEAL:3d:3::87 0.000 -10000.000
+DEAL:3d:3::88 0.000 -10000.000
+DEAL:3d:3::89 0.000 -10000.000
+DEAL:3d:3::90 0.000 -10000.000
+DEAL:3d:3::93 0.000 -10000.000
+DEAL:3d:3::94 0.000 -10000.000
+DEAL:3d:3::95 0.000 -10000.000
+DEAL:3d:3::96 0.000 -10000.000
+DEAL:3d:3::97 0.000 -10000.000
+DEAL:3d:3::98 0.000 -10000.000
+DEAL:3d:3::99 0.000 -10000.000
+DEAL:3d:3::100 0.000 -10000.000
+DEAL:3d:3::101 0.000 -10000.000
+DEAL:3d:3::102 0.000 -10000.000
+DEAL:3d:3::103 0.000 -10000.000
+DEAL:3d:3::104 0.000 -10000.000
+DEAL:3d:3::105 0.000 -10000.000
+DEAL:3d:3::106 0.000 -10000.000
+DEAL:3d:3::107 0.000 -10000.000
+DEAL:3d:3::109 0.000 -10000.000
+DEAL:3d:3::111 0.000 -10000.000
+DEAL:3d:3::113 0.000 -10000.000
+DEAL:3d:3::115 0.000 -10000.000
+DEAL:3d:3::116 0.000 -10000.000
+DEAL:3d:3::117 0.000 -10000.000
+DEAL:3d:3::118 0.000 -10000.000
+DEAL:3d:3::119 0.000 -10000.000
+DEAL:3d:3::120 0.000 -10000.000
+DEAL:3d:3::121 0.000 -10000.000
+DEAL:3d:3::122 0.000 -10000.000
+DEAL:3d:3::123 0.000 -10000.000
+DEAL:3d:3::124 0.000 -10000.000
+DEAL:3d:3::125 0.000 -10000.000
+DEAL:3d:3::127 0.000 -10000.000
+DEAL:3d:3::129 0.000 -10000.000
+DEAL:3d:3::130 0.000 -10000.000
+DEAL:3d:3::131 0.000 -10000.000
+DEAL:3d:3::132 0.000 -10000.000
+DEAL:3d:3::133 0.000 -10000.000
+DEAL:3d:3::134 0.000 -10000.000
+DEAL:3d:3::135 0.000 -10000.000
+DEAL:3d:3::136 0.000 -10000.000
+DEAL:3d:3::137 0.000 -10000.000
+DEAL:3d:3::138 0.000 -10000.000
+DEAL:3d:3::140 0.000 -10000.000
+DEAL:3d:3::141 0.000 -10000.000
+DEAL:3d:3::143 0.000 -10000.000
+DEAL:3d:3::144 0.000 -10000.000
+DEAL:3d:3::145 0.000 -10000.000
+DEAL:3d:3::146 0.000 -10000.000
+DEAL:3d:3::147 0.000 -10000.000
+DEAL:3d:3::148 0.000 -10000.000
+DEAL:3d:3::149 0.000 -10000.000
+DEAL:3d:3::150 0.000 -10000.000
+DEAL:3d:3::Interpolated boundary values
+DEAL:3d:3::0 0.000 -10000.000
+DEAL:3d:3::1 0.000 -10000.000
+DEAL:3d:3::2 0.000 -10000.000
+DEAL:3d:3::3 0.000 -10000.000
+DEAL:3d:3::4 0.000 -10000.000
+DEAL:3d:3::5 0.000 -10000.000
+DEAL:3d:3::6 0.000 -10000.000
+DEAL:3d:3::8 0.000 -10000.000
+DEAL:3d:3::9 0.000 -10000.000
+DEAL:3d:3::10 0.000 -10000.000
+DEAL:3d:3::11 0.000 -10000.000
+DEAL:3d:3::12 0.000 -10000.000
+DEAL:3d:3::13 0.000 -10000.000
+DEAL:3d:3::14 0.000 -10000.000
+DEAL:3d:3::15 0.000 -10000.000
+DEAL:3d:3::16 0.000 -10000.000
+DEAL:3d:3::17 0.000 -10000.000
+DEAL:3d:3::18 0.000 -10000.000
+DEAL:3d:3::19 0.000 -10000.000
+DEAL:3d:3::20 0.000 -10000.000
+DEAL:3d:3::21 0.000 -10000.000
+DEAL:3d:3::22 0.000 -10000.000
+DEAL:3d:3::23 0.000 -10000.000
+DEAL:3d:3::24 0.000 -10000.000
+DEAL:3d:3::25 0.000 -10000.000
+DEAL:3d:3::26 0.000 -10000.000
+DEAL:3d:3::Projected boundary values
+DEAL:3d:3::0 0.000 -10000.000
+DEAL:3d:3::1 0.000 -10000.000
+DEAL:3d:3::2 0.000 -10000.000
+DEAL:3d:3::3 0.000 -10000.000
+DEAL:3d:3::4 0.000 -10000.000
+DEAL:3d:3::5 0.000 -10000.000
+DEAL:3d:3::6 0.000 -10000.000
+DEAL:3d:3::8 0.000 -10000.000
+DEAL:3d:3::9 0.000 -10000.000
+DEAL:3d:3::10 0.000 -10000.000
+DEAL:3d:3::11 0.000 -10000.000
+DEAL:3d:3::12 0.000 -10000.000
+DEAL:3d:3::13 0.000 -10000.000
+DEAL:3d:3::14 0.000 -10000.000
+DEAL:3d:3::15 0.000 -10000.000
+DEAL:3d:3::16 0.000 -10000.000
+DEAL:3d:3::17 0.000 -10000.000
+DEAL:3d:3::18 0.000 -10000.000
+DEAL:3d:3::19 0.000 -10000.000
+DEAL:3d:3::20 0.000 -10000.000
+DEAL:3d:3::21 0.000 -10000.000
+DEAL:3d:3::22 0.000 -10000.000
+DEAL:3d:3::23 0.000 -10000.000
+DEAL:3d:3::24 0.000 -10000.000
+DEAL:3d:3::25 0.000 -10000.000
+DEAL:3d:3::26 0.000 -10000.000