[cig-commits] commit 2303 by bangerth to /var/svn/dealii/aspect
dealii.demon at gmail.com
dealii.demon at gmail.com
Sun Feb 16 19:23:45 PST 2014
Revision 2303
Make the depth average something that can be visualized properly in a visualization program.
U trunk/aspect/include/aspect/postprocess/depth_average.h
U trunk/aspect/source/postprocess/depth_average.cc
A trunk/aspect/tests/depth_average_01/depth_average.gnuplot
D trunk/aspect/tests/depth_average_01/depthaverage.plt
A trunk/aspect/tests/depth_average_02/depth_average.gnuplot
D trunk/aspect/tests/depth_average_02/depthaverage.plt
A trunk/aspect/tests/depth_average_03/depth_average.gnuplot
D trunk/aspect/tests/depth_average_03/depthaverage.plt
A trunk/aspect/tests/depth_average_04/
A trunk/aspect/tests/depth_average_04/depth_average.vtk
D trunk/aspect/tests/depth_average_04/depthaverage.plt
A trunk/aspect/tests/depth_average_04.prm
http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=2303&peg=2303
Diff:
Modified: trunk/aspect/include/aspect/postprocess/depth_average.h
===================================================================
--- trunk/aspect/include/aspect/postprocess/depth_average.h 2014-02-15 13:32:17 UTC (rev 2302)
+++ trunk/aspect/include/aspect/postprocess/depth_average.h 2014-02-17 03:23:42 UTC (rev 2303)
@@ -121,12 +121,12 @@
unsigned int n_depth_zones;
/**
- * A structure for a single time/depth/value record.
+ * A structure for a single time step record.
*/
struct DataPoint
{
- double time, depth;
- std::vector<double> values;
+ double time;
+ std::vector<std::vector<double> > values;
template <class Archive>
void serialize (Archive &ar, const unsigned int version);
Modified: trunk/aspect/source/postprocess/depth_average.cc
===================================================================
--- trunk/aspect/source/postprocess/depth_average.cc 2014-02-15 13:32:17 UTC (rev 2302)
+++ trunk/aspect/source/postprocess/depth_average.cc 2014-02-17 03:23:42 UTC (rev 2303)
@@ -24,7 +24,11 @@
#include <aspect/simulator_access.h>
#include <aspect/global.h>
-#include <deal.II/dofs/dof_tools.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/fe/fe_dgq.h>
+#include <deal.II/dofs/dof_handler.h>
+#include <deal.II/numerics/data_out_stack.h>
#include <math.h>
@@ -37,7 +41,7 @@
template <class Archive>
void DepthAverage<dim>::DataPoint::serialize (Archive &ar, const unsigned int version)
{
- ar &time &depth &values;
+ ar &time &values;
}
@@ -69,71 +73,101 @@
return std::pair<std::string,std::string>();
const unsigned int n_statistics = 7+this->n_compositional_fields();
- std::vector<std::vector<double> > temp(n_statistics,
- std::vector<double> (n_depth_zones));
+
+ DataPoint data_point;
+ data_point.time = this->get_time();
+ data_point.values.resize(n_statistics, std::vector<double> (n_depth_zones));
+
+ // add temperature and the compositional fields that follow
+ // it immediately
{
- // add temperature and the compositional fields that follow
- // it immediately
- this->get_depth_average_temperature(temp[0]);
+ this->get_depth_average_temperature(data_point.values[0]);
for (unsigned int c=0; c<this->n_compositional_fields(); ++c)
- this->get_depth_average_composition(c, temp[1+c]);
- this->get_adiabatic_conditions().get_adiabatic_temperature_profile(temp[1+this->n_compositional_fields()]);
- this->get_depth_average_velocity_magnitude(temp[2+this->n_compositional_fields()]);
- this->get_depth_average_sinking_velocity(temp[3+this->n_compositional_fields()]);
- this->get_depth_average_Vs(temp[4+this->n_compositional_fields()]);
- this->get_depth_average_Vp(temp[5+this->n_compositional_fields()]);
- this->get_depth_average_viscosity(temp[6+this->n_compositional_fields()]);
+ this->get_depth_average_composition(c, data_point.values[1+c]);
+ this->get_adiabatic_conditions().get_adiabatic_temperature_profile(data_point.values[1+this->n_compositional_fields()]);
+ this->get_depth_average_velocity_magnitude(data_point.values[2+this->n_compositional_fields()]);
+ this->get_depth_average_sinking_velocity(data_point.values[3+this->n_compositional_fields()]);
+ this->get_depth_average_Vs(data_point.values[4+this->n_compositional_fields()]);
+ this->get_depth_average_Vp(data_point.values[5+this->n_compositional_fields()]);
+ this->get_depth_average_viscosity(data_point.values[6+this->n_compositional_fields()]);
}
+ entries.push_back (data_point);
const double max_depth = this->get_geometry_model().maximal_depth();
+ // On the root process, write out the file. do this using the DataOutStack
+ // class on a piecewise constant finite element space on
+ // a 1d mesh with the correct subdivisions
if (Utilities::MPI::this_mpi_process(this->get_mpi_communicator()) == 0)
{
- // store data from the current step
- for (unsigned int i=0; i<n_depth_zones; ++i)
- {
- DataPoint data_point;
- data_point.time = this->get_time();
- data_point.depth = max_depth*i/temp[0].size();
- data_point.values.resize(n_statistics);
- for (unsigned int j=0; j<n_statistics; ++j)
- data_point.values[j]=temp[j][i];
- entries.push_back(data_point);
- }
+ Triangulation<1> mesh;
+ GridGenerator::subdivided_hyper_cube(mesh, n_depth_zones, 0, max_depth);
- // leave a marker for the end of this time step. We'll later
- // use this to leave a blank line in the output file at this
- // position.
- DataPoint data_point;
- entries.push_back(data_point);
+ FE_DGQ<1> fe(0);
+ DoFHandler<1> dof_handler (mesh);
+ dof_handler.distribute_dofs(fe);
+ Assert (dof_handler.n_dofs() == n_depth_zones, ExcInternalError());
- // write out the file
- const std::string filename = this->get_output_directory() + "depthaverage.plt";
- std::ofstream f (filename.c_str());
- f << "# time, depth, avg T";
+ DataOutStack<1> data_out_stack;
+ std::vector<string> variables;
+ variables.push_back ("temperature");
for (unsigned int c=0; c<this->n_compositional_fields(); ++c)
- f << ", avg C_" << c+1;
- f << ", adiabatic T, velocity magnitude, avg sinking velocity, avg Vs, avg Vp, avg viscosity" << std::endl;
- f << std::scientific;
+ variables.push_back(std::string("C_") + Utilities::int_to_string(c));
+ variables.push_back ("adiabatic_temperature");
+ variables.push_back ("velocity_magnitude");
+ variables.push_back ("sinking_velocity");
+ variables.push_back ("Vs");
+ variables.push_back ("Vp");
+ variables.push_back ("viscosity");
+ Assert (variables.size() == n_statistics, ExcInternalError());
+ for (unsigned int j=0; j<n_statistics; ++j)
+ data_out_stack.declare_data_vector (variables[j],
+ DataOutStack<1>::cell_vector);
+
for (unsigned int i=0; i<entries.size(); ++i)
{
- if (entries[i].values.empty())
+ data_out_stack.new_parameter_value ((this->convert_output_to_years()
+ ?
+ entries[i].time / year_in_seconds
+ :
+ entries[i].time),
+ // declare the time step, which here is the difference
+ // between successive output times. we don't have anything
+ // for the first time step, however. we could do a zero
+ // delta, but that leads to invisible output. rather, we
+ // use an artificial value of one tenth of the first interval,
+ // if available
+ (i == 0 ?
+ (entries.size() > 1 ? (entries[1].time - entries[0].time)/10 : 0) :
+ entries[i].time - entries[i-1].time) /
+ (this->convert_output_to_years()
+ ?
+ year_in_seconds
+ :
+ 1));
+
+ data_out_stack.attach_dof_handler (dof_handler);
+
+ Vector<double> tmp(n_depth_zones);
+ for (unsigned int j=0; j<n_statistics; ++j)
{
- f << std::endl;
- continue;
+ std::copy (entries[i].values[j].begin(),
+ entries[i].values[j].end(),
+ tmp.begin());
+ data_out_stack.add_data_vector (tmp,
+ variables[j]);
}
- f << (this->convert_output_to_years()
- ?
- entries[i].time / year_in_seconds
- :
- entries[i].time)
- << " "
- << entries[i].depth;
- for (unsigned int j=0; j<n_statistics; ++j)
- f << " " << entries[i].values[j];
- f << std::endl;
+ data_out_stack.build_patches ();
+ data_out_stack.finish_parameter_value ();
}
+
+
+ const std::string filename = (this->get_output_directory() +
+ "depth_average" +
+ DataOutBase::default_suffix(output_format));
+ std::ofstream f (filename.c_str());
+ data_out_stack.write (f, output_format);
}
@@ -142,7 +176,9 @@
// return what should be printed to the screen. note that we had
// just incremented the number, so use the previous value
return std::make_pair (std::string ("Writing depth average"),
- this->get_output_directory() + "depthaverage.plt");
+ this->get_output_directory() +
+ "depth_average" +
+ DataOutBase::default_suffix(output_format));
}
@@ -173,6 +209,11 @@
"less than this default. It may also make computations slightly "
"faster. On the other hand, if you have an extremely highly "
"resolved mesh, choosing more zones might also make sense.");
+ prm.declare_entry ("Output format", "gnuplot",
+ Patterns::Selection(DataOutBase::get_output_format_names()),
+ "The format in which the output shall be produced. The "
+ "format in which the output is generated also determiens "
+ "the extension of the file into which data is written.");
}
prm.leave_subsection();
}
@@ -190,6 +231,7 @@
{
output_interval = prm.get_double ("Time between graphical output");
n_depth_zones = prm.get_integer ("Number of zones");
+ output_format = DataOutBase::parse_output_format(prm.get("Output format"));
}
prm.leave_subsection();
}
@@ -244,9 +286,9 @@
if (output_interval > 0)
{
// the current time is always in seconds, so we need to convert the output_interval to the same unit
- double output_interval_in_s = (this->convert_output_to_years() ?
- (output_interval*year_in_seconds) :
- output_interval);
+ const double output_interval_in_s = (this->convert_output_to_years() ?
+ (output_interval*year_in_seconds) :
+ output_interval);
// we need to compute the smallest integer that is bigger than current_time/my_output_interval,
// even if it is a whole number already (otherwise we output twice in a row)
@@ -265,9 +307,11 @@
ASPECT_REGISTER_POSTPROCESSOR(DepthAverage,
"depth average",
"A postprocessor that computes depth averaged "
- "quantities and writes them into a file "
- "in the output directory. A number of parameters influence "
+ "quantities and writes them into a file <depth_average.ext> "
+ "in the output directory, where the extension of the file "
+ "is determined by the output format you select. In addition "
+ "to the output format, a number of other parameters also influence "
"this postprocessor, and they can be set in the section "
- "\texttt{Postprocess/Depth average} in the input file.")
+ "\texttt{Postprocess/Depth average} in the input file. ")
}
}
Copied: trunk/aspect/tests/depth_average_01/depth_average.gnuplot (from rev 2302, trunk/aspect/tests/depth_average_01/depthaverage.plt)
===================================================================
--- trunk/aspect/tests/depth_average_01/depth_average.gnuplot (rev 0)
+++ trunk/aspect/tests/depth_average_01/depth_average.gnuplot 2014-02-17 03:23:42 UTC (rev 2303)
@@ -0,0 +1,4207 @@
+# This file was generated by the deal.II library.
+
+
+#
+# For a description of the GNUPLOT format see the GNUPLOT manual.
+#
+# <x> <y> <temperature> <adiabatic_temperature> <velocity_magnitude> <sinking_velocity> <Vs> <Vp> <viscosity>
+0 -0.0100187 0 0 0 0 0 0 0
+0.01 -0.0100187 0 0 0 0 0 0 0
+
+0 0 0 0 0 0 0 0 0
+0.01 0 0 0 0 0 0 0 0
+
+
+0.01 -0.0100187 0.0125 0 3.20592 1.03518e+06 -1 -1 1
+0.02 -0.0100187 0.0125 0 3.20592 1.03518e+06 -1 -1 1
+
+0.01 0 0.0125 0 3.20592 1.03518e+06 -1 -1 1
+0.02 0 0.0125 0 3.20592 1.03518e+06 -1 -1 1
+
+
+0.02 -0.0100187 0 0 0 0 0 0 0
+0.03 -0.0100187 0 0 0 0 0 0 0
+
+0.02 0 0 0 0 0 0 0 0
+0.03 0 0 0 0 0 0 0 0
+
+
+0.03 -0.0100187 0.0375 0 3.19899 3.06349e+06 -1 -1 1
+0.04 -0.0100187 0.0375 0 3.19899 3.06349e+06 -1 -1 1
+
+0.03 0 0.0375 0 3.19899 3.06349e+06 -1 -1 1
+0.04 0 0.0375 0 3.19899 3.06349e+06 -1 -1 1
+
+
+0.04 -0.0100187 0 0 0 0 0 0 0
+0.05 -0.0100187 0 0 0 0 0 0 0
+
+0.04 0 0 0 0 0 0 0 0
+0.05 0 0 0 0 0 0 0 0
+
+
+0.05 -0.0100187 0 0 0 0 0 0 0
+0.06 -0.0100187 0 0 0 0 0 0 0
+
+0.05 0 0 0 0 0 0 0 0
+0.06 0 0 0 0 0 0 0 0
+
+
+0.06 -0.0100187 0.0625 0 3.19608 5.03572e+06 -1 -1 1
+0.07 -0.0100187 0.0625 0 3.19608 5.03572e+06 -1 -1 1
+
+0.06 0 0.0625 0 3.19608 5.03572e+06 -1 -1 1
+0.07 0 0.0625 0 3.19608 5.03572e+06 -1 -1 1
+
+
+0.07 -0.0100187 0 0 0 0 0 0 0
+0.08 -0.0100187 0 0 0 0 0 0 0
+
+0.07 0 0 0 0 0 0 0 0
+0.08 0 0 0 0 0 0 0 0
+
+
+0.08 -0.0100187 0.0875 0 3.19524 6.95189e+06 -1 -1 1
+0.09 -0.0100187 0.0875 0 3.19524 6.95189e+06 -1 -1 1
+
+0.08 0 0.0875 0 3.19524 6.95189e+06 -1 -1 1
+0.09 0 0.0875 0 3.19524 6.95189e+06 -1 -1 1
+
+
+0.09 -0.0100187 0 0 0 0 0 0 0
+0.1 -0.0100187 0 0 0 0 0 0 0
+
+0.09 0 0 0 0 0 0 0 0
+0.1 0 0 0 0 0 0 0 0
+
+
+0.1 -0.0100187 0 0 0 0 0 0 0
+0.11 -0.0100187 0 0 0 0 0 0 0
+
+0.1 0 0 0 0 0 0 0 0
+0.11 0 0 0 0 0 0 0 0
+
+
+0.11 -0.0100187 0.1125 0 3.19519 8.81198e+06 -1 -1 1
+0.12 -0.0100187 0.1125 0 3.19519 8.81198e+06 -1 -1 1
+
+0.11 0 0.1125 0 3.19519 8.81198e+06 -1 -1 1
+0.12 0 0.1125 0 3.19519 8.81198e+06 -1 -1 1
+
+
+0.12 -0.0100187 0 0 0 0 0 0 0
+0.13 -0.0100187 0 0 0 0 0 0 0
+
+0.12 0 0 0 0 0 0 0 0
+0.13 0 0 0 0 0 0 0 0
+
+
+0.13 -0.0100187 0.1375 0 3.19533 1.0616e+07 -1 -1 1
+0.14 -0.0100187 0.1375 0 3.19533 1.0616e+07 -1 -1 1
+
+0.13 0 0.1375 0 3.19533 1.0616e+07 -1 -1 1
+0.14 0 0.1375 0 3.19533 1.0616e+07 -1 -1 1
+
+
+0.14 -0.0100187 0 0 0 0 0 0 0
+0.15 -0.0100187 0 0 0 0 0 0 0
+
+0.14 0 0 0 0 0 0 0 0
+0.15 0 0 0 0 0 0 0 0
+
+
+0.15 -0.0100187 0 0 0 0 0 0 0
+0.16 -0.0100187 0 0 0 0 0 0 0
+
+0.15 0 0 0 0 0 0 0 0
+0.16 0 0 0 0 0 0 0 0
+
+
+0.16 -0.0100187 0.1625 0 3.19571 1.23639e+07 -1 -1 1
+0.17 -0.0100187 0.1625 0 3.19571 1.23639e+07 -1 -1 1
+
+0.16 0 0.1625 0 3.19571 1.23639e+07 -1 -1 1
+0.17 0 0.1625 0 3.19571 1.23639e+07 -1 -1 1
+
+
+0.17 -0.0100187 0 0 0 0 0 0 0
+0.18 -0.0100187 0 0 0 0 0 0 0
+
+0.17 0 0 0 0 0 0 0 0
+0.18 0 0 0 0 0 0 0 0
+
+
+0.18 -0.0100187 0.1875 0 3.19704 1.40558e+07 -1 -1 1
+0.19 -0.0100187 0.1875 0 3.19704 1.40558e+07 -1 -1 1
+
+0.18 0 0.1875 0 3.19704 1.40558e+07 -1 -1 1
+0.19 0 0.1875 0 3.19704 1.40558e+07 -1 -1 1
+
+
+0.19 -0.0100187 0 0 0 0 0 0 0
+0.2 -0.0100187 0 0 0 0 0 0 0
+
+0.19 0 0 0 0 0 0 0 0
+0.2 0 0 0 0 0 0 0 0
+
+
+0.2 -0.0100187 0 0 0 0 0 0 0
+0.21 -0.0100187 0 0 0 0 0 0 0
+
+0.2 0 0 0 0 0 0 0 0
+0.21 0 0 0 0 0 0 0 0
+
+
+0.21 -0.0100187 0.2125 0 3.2007 1.56916e+07 -1 -1 1
+0.22 -0.0100187 0.2125 0 3.2007 1.56916e+07 -1 -1 1
+
+0.21 0 0.2125 0 3.2007 1.56916e+07 -1 -1 1
+0.22 0 0.2125 0 3.2007 1.56916e+07 -1 -1 1
+
+
+0.22 -0.0100187 0 0 0 0 0 0 0
+0.23 -0.0100187 0 0 0 0 0 0 0
+
+0.22 0 0 0 0 0 0 0 0
+0.23 0 0 0 0 0 0 0 0
+
+
+0.23 -0.0100187 0.2375 0 3.20874 1.72713e+07 -1 -1 1
+0.24 -0.0100187 0.2375 0 3.20874 1.72713e+07 -1 -1 1
+
+0.23 0 0.2375 0 3.20874 1.72713e+07 -1 -1 1
+0.24 0 0.2375 0 3.20874 1.72713e+07 -1 -1 1
+
+
+0.24 -0.0100187 0 0 0 0 0 0 0
+0.25 -0.0100187 0 0 0 0 0 0 0
+
+0.24 0 0 0 0 0 0 0 0
+0.25 0 0 0 0 0 0 0 0
+
+
+0.25 -0.0100187 0 0 0 0 0 0 0
+0.26 -0.0100187 0 0 0 0 0 0 0
+
+0.25 0 0 0 0 0 0 0 0
+0.26 0 0 0 0 0 0 0 0
+
+
+0.26 -0.0100187 0.2625 0 3.2087 1.87368e+07 -1 -1 1
+0.27 -0.0100187 0.2625 0 3.2087 1.87368e+07 -1 -1 1
+
+0.26 0 0.2625 0 3.2087 1.87368e+07 -1 -1 1
+0.27 0 0.2625 0 3.2087 1.87368e+07 -1 -1 1
+
+
+0.27 -0.0100187 0 0 0 0 0 0 0
+0.28 -0.0100187 0 0 0 0 0 0 0
+
+0.27 0 0 0 0 0 0 0 0
+0.28 0 0 0 0 0 0 0 0
+
+
+0.28 -0.0100187 0.2875 0 3.20041 2.00278e+07 -1 -1 1
+0.29 -0.0100187 0.2875 0 3.20041 2.00278e+07 -1 -1 1
+
+0.28 0 0.2875 0 3.20041 2.00278e+07 -1 -1 1
+0.29 0 0.2875 0 3.20041 2.00278e+07 -1 -1 1
+
+
+0.29 -0.0100187 0 0 0 0 0 0 0
+0.3 -0.0100187 0 0 0 0 0 0 0
+
+0.29 0 0 0 0 0 0 0 0
+0.3 0 0 0 0 0 0 0 0
+
+
+0.3 -0.0100187 0 0 0 0 0 0 0
+0.31 -0.0100187 0 0 0 0 0 0 0
+
+0.3 0 0 0 0 0 0 0 0
+0.31 0 0 0 0 0 0 0 0
+
+
+0.31 -0.0100187 0.3125 0 3.19639 2.11824e+07 -1 -1 1
+0.32 -0.0100187 0.3125 0 3.19639 2.11824e+07 -1 -1 1
+
+0.31 0 0.3125 0 3.19639 2.11824e+07 -1 -1 1
+0.32 0 0.3125 0 3.19639 2.11824e+07 -1 -1 1
+
+
+0.32 -0.0100187 0 0 0 0 0 0 0
+0.33 -0.0100187 0 0 0 0 0 0 0
+
+0.32 0 0 0 0 0 0 0 0
+0.33 0 0 0 0 0 0 0 0
+
+
+0.33 -0.0100187 0.3375 0 3.19474 2.22006e+07 -1 -1 1
+0.34 -0.0100187 0.3375 0 3.19474 2.22006e+07 -1 -1 1
+
+0.33 0 0.3375 0 3.19474 2.22006e+07 -1 -1 1
+0.34 0 0.3375 0 3.19474 2.22006e+07 -1 -1 1
+
+
+0.34 -0.0100187 0 0 0 0 0 0 0
+0.35 -0.0100187 0 0 0 0 0 0 0
+
+0.34 0 0 0 0 0 0 0 0
+0.35 0 0 0 0 0 0 0 0
+
+
+0.35 -0.0100187 0 0 0 0 0 0 0
+0.36 -0.0100187 0 0 0 0 0 0 0
+
+0.35 0 0 0 0 0 0 0 0
+0.36 0 0 0 0 0 0 0 0
+
+
+0.36 -0.0100187 0.3625 0 3.19419 2.30825e+07 -1 -1 1
+0.37 -0.0100187 0.3625 0 3.19419 2.30825e+07 -1 -1 1
+
+0.36 0 0.3625 0 3.19419 2.30825e+07 -1 -1 1
+0.37 0 0.3625 0 3.19419 2.30825e+07 -1 -1 1
+
+
+0.37 -0.0100187 0 0 0 0 0 0 0
+0.38 -0.0100187 0 0 0 0 0 0 0
+
+0.37 0 0 0 0 0 0 0 0
+0.38 0 0 0 0 0 0 0 0
+
+
+0.38 -0.0100187 0.3875 0 3.19415 2.3828e+07 -1 -1 1
+0.39 -0.0100187 0.3875 0 3.19415 2.3828e+07 -1 -1 1
+
+0.38 0 0.3875 0 3.19415 2.3828e+07 -1 -1 1
+0.39 0 0.3875 0 3.19415 2.3828e+07 -1 -1 1
+
+
+0.39 -0.0100187 0 0 0 0 0 0 0
+0.4 -0.0100187 0 0 0 0 0 0 0
+
+0.39 0 0 0 0 0 0 0 0
+0.4 0 0 0 0 0 0 0 0
+
+
+0.4 -0.0100187 0 0 0 0 0 0 0
+0.41 -0.0100187 0 0 0 0 0 0 0
+
+0.4 0 0 0 0 0 0 0 0
+0.41 0 0 0 0 0 0 0 0
+
+
+0.41 -0.0100187 0.4125 0 3.19465 2.44372e+07 -1 -1 1
+0.42 -0.0100187 0.4125 0 3.19465 2.44372e+07 -1 -1 1
+
+0.41 0 0.4125 0 3.19465 2.44372e+07 -1 -1 1
+0.42 0 0.4125 0 3.19465 2.44372e+07 -1 -1 1
+
+
+0.42 -0.0100187 0 0 0 0 0 0 0
+0.43 -0.0100187 0 0 0 0 0 0 0
+
+0.42 0 0 0 0 0 0 0 0
+0.43 0 0 0 0 0 0 0 0
+
+
+0.43 -0.0100187 0.4375 0 3.1964 2.491e+07 -1 -1 1
+0.44 -0.0100187 0.4375 0 3.1964 2.491e+07 -1 -1 1
+
+0.43 0 0.4375 0 3.1964 2.491e+07 -1 -1 1
+0.44 0 0.4375 0 3.1964 2.491e+07 -1 -1 1
+
+
+0.44 -0.0100187 0 0 0 0 0 0 0
+0.45 -0.0100187 0 0 0 0 0 0 0
+
+0.44 0 0 0 0 0 0 0 0
+0.45 0 0 0 0 0 0 0 0
+
+
+0.45 -0.0100187 0 0 0 0 0 0 0
+0.46 -0.0100187 0 0 0 0 0 0 0
+
+0.45 0 0 0 0 0 0 0 0
+0.46 0 0 0 0 0 0 0 0
+
+
+0.46 -0.0100187 0.4625 0 3.20075 2.52464e+07 -1 -1 1
+0.47 -0.0100187 0.4625 0 3.20075 2.52464e+07 -1 -1 1
+
+0.46 0 0.4625 0 3.20075 2.52464e+07 -1 -1 1
+0.47 0 0.4625 0 3.20075 2.52464e+07 -1 -1 1
+
+
+0.47 -0.0100187 0 0 0 0 0 0 0
+0.48 -0.0100187 0 0 0 0 0 0 0
+
+0.47 0 0 0 0 0 0 0 0
+0.48 0 0 0 0 0 0 0 0
+
+
+0.48 -0.0100187 0.4875 0 3.20971 2.54465e+07 -1 -1 1
+0.49 -0.0100187 0.4875 0 3.20971 2.54465e+07 -1 -1 1
+
+0.48 0 0.4875 0 3.20971 2.54465e+07 -1 -1 1
+0.49 0 0.4875 0 3.20971 2.54465e+07 -1 -1 1
+
+
+0.49 -0.0100187 0 0 0 0 0 0 0
+0.5 -0.0100187 0 0 0 0 0 0 0
+
+0.49 0 0 0 0 0 0 0 0
+0.5 0 0 0 0 0 0 0 0
+
+
+0.5 -0.0100187 0 0 0 0 0 0 0
+0.51 -0.0100187 0 0 0 0 0 0 0
+
+0.5 0 0 0 0 0 0 0 0
+0.51 0 0 0 0 0 0 0 0
+
+
+0.51 -0.0100187 0.5125 0 3.20963 2.54466e+07 -1 -1 1
+0.52 -0.0100187 0.5125 0 3.20963 2.54466e+07 -1 -1 1
+
+0.51 0 0.5125 0 3.20963 2.54466e+07 -1 -1 1
+0.52 0 0.5125 0 3.20963 2.54466e+07 -1 -1 1
+
+
+0.52 -0.0100187 0 0 0 0 0 0 0
+0.53 -0.0100187 0 0 0 0 0 0 0
+
+0.52 0 0 0 0 0 0 0 0
+0.53 0 0 0 0 0 0 0 0
+
+
+0.53 -0.0100187 0.5375 0 3.20057 2.52465e+07 -1 -1 1
+0.54 -0.0100187 0.5375 0 3.20057 2.52465e+07 -1 -1 1
+
+0.53 0 0.5375 0 3.20057 2.52465e+07 -1 -1 1
+0.54 0 0.5375 0 3.20057 2.52465e+07 -1 -1 1
+
+
+0.54 -0.0100187 0 0 0 0 0 0 0
+0.55 -0.0100187 0 0 0 0 0 0 0
+
+0.54 0 0 0 0 0 0 0 0
+0.55 0 0 0 0 0 0 0 0
+
+
+0.55 -0.0100187 0 0 0 0 0 0 0
+0.56 -0.0100187 0 0 0 0 0 0 0
+
+0.55 0 0 0 0 0 0 0 0
+0.56 0 0 0 0 0 0 0 0
+
+
+0.56 -0.0100187 0.5625 0 3.19617 2.491e+07 -1 -1 1
+0.57 -0.0100187 0.5625 0 3.19617 2.491e+07 -1 -1 1
+
+0.56 0 0.5625 0 3.19617 2.491e+07 -1 -1 1
+0.57 0 0.5625 0 3.19617 2.491e+07 -1 -1 1
+
+
+0.57 -0.0100187 0 0 0 0 0 0 0
+0.58 -0.0100187 0 0 0 0 0 0 0
+
+0.57 0 0 0 0 0 0 0 0
+0.58 0 0 0 0 0 0 0 0
+
+
+0.58 -0.0100187 0.5875 0 3.19445 2.44369e+07 -1 -1 1
+0.59 -0.0100187 0.5875 0 3.19445 2.44369e+07 -1 -1 1
+
+0.58 0 0.5875 0 3.19445 2.44369e+07 -1 -1 1
+0.59 0 0.5875 0 3.19445 2.44369e+07 -1 -1 1
+
+
+0.59 -0.0100187 0 0 0 0 0 0 0
+0.6 -0.0100187 0 0 0 0 0 0 0
+
+0.59 0 0 0 0 0 0 0 0
+0.6 0 0 0 0 0 0 0 0
+
+
+0.6 -0.0100187 0 0 0 0 0 0 0
+0.61 -0.0100187 0 0 0 0 0 0 0
+
+0.6 0 0 0 0 0 0 0 0
+0.61 0 0 0 0 0 0 0 0
+
+
+0.61 -0.0100187 0.6125 0 3.19405 2.38274e+07 -1 -1 1
+0.62 -0.0100187 0.6125 0 3.19405 2.38274e+07 -1 -1 1
+
+0.61 0 0.6125 0 3.19405 2.38274e+07 -1 -1 1
+0.62 0 0.6125 0 3.19405 2.38274e+07 -1 -1 1
+
+
+0.62 -0.0100187 0 0 0 0 0 0 0
+0.63 -0.0100187 0 0 0 0 0 0 0
+
+0.62 0 0 0 0 0 0 0 0
+0.63 0 0 0 0 0 0 0 0
+
+
+0.63 -0.0100187 0.6375 0 3.19427 2.30813e+07 -1 -1 1
+0.64 -0.0100187 0.6375 0 3.19427 2.30813e+07 -1 -1 1
+
+0.63 0 0.6375 0 3.19427 2.30813e+07 -1 -1 1
+0.64 0 0.6375 0 3.19427 2.30813e+07 -1 -1 1
+
+
+0.64 -0.0100187 0 0 0 0 0 0 0
+0.65 -0.0100187 0 0 0 0 0 0 0
+
+0.64 0 0 0 0 0 0 0 0
+0.65 0 0 0 0 0 0 0 0
+
+
+0.65 -0.0100187 0 0 0 0 0 0 0
+0.66 -0.0100187 0 0 0 0 0 0 0
+
+0.65 0 0 0 0 0 0 0 0
+0.66 0 0 0 0 0 0 0 0
+
+
+0.66 -0.0100187 0.6625 0 3.19508 2.21987e+07 -1 -1 1
+0.67 -0.0100187 0.6625 0 3.19508 2.21987e+07 -1 -1 1
+
+0.66 0 0.6625 0 3.19508 2.21987e+07 -1 -1 1
+0.67 0 0.6625 0 3.19508 2.21987e+07 -1 -1 1
+
+
+0.67 -0.0100187 0 0 0 0 0 0 0
+0.68 -0.0100187 0 0 0 0 0 0 0
+
+0.67 0 0 0 0 0 0 0 0
+0.68 0 0 0 0 0 0 0 0
+
+
+0.68 -0.0100187 0.6875 0 3.19707 2.11796e+07 -1 -1 1
+0.69 -0.0100187 0.6875 0 3.19707 2.11796e+07 -1 -1 1
+
+0.68 0 0.6875 0 3.19707 2.11796e+07 -1 -1 1
+0.69 0 0.6875 0 3.19707 2.11796e+07 -1 -1 1
+
+
+0.69 -0.0100187 0 0 0 0 0 0 0
+0.7 -0.0100187 0 0 0 0 0 0 0
+
+0.69 0 0 0 0 0 0 0 0
+0.7 0 0 0 0 0 0 0 0
+
+
+0.7 -0.0100187 0 0 0 0 0 0 0
+0.71 -0.0100187 0 0 0 0 0 0 0
+
+0.7 0 0 0 0 0 0 0 0
+0.71 0 0 0 0 0 0 0 0
+
+
+0.71 -0.0100187 0.7125 0 3.2015 2.00239e+07 -1 -1 1
+0.72 -0.0100187 0.7125 0 3.2015 2.00239e+07 -1 -1 1
+
+0.71 0 0.7125 0 3.2015 2.00239e+07 -1 -1 1
+0.72 0 0.7125 0 3.2015 2.00239e+07 -1 -1 1
+
+
+0.72 -0.0100187 0 0 0 0 0 0 0
+0.73 -0.0100187 0 0 0 0 0 0 0
+
+0.72 0 0 0 0 0 0 0 0
+0.73 0 0 0 0 0 0 0 0
+
+
+0.73 -0.0100187 0.7375 0 3.21028 1.87318e+07 -1 -1 1
+0.74 -0.0100187 0.7375 0 3.21028 1.87318e+07 -1 -1 1
+
+0.73 0 0.7375 0 3.21028 1.87318e+07 -1 -1 1
+0.74 0 0.7375 0 3.21028 1.87318e+07 -1 -1 1
+
+
+0.74 -0.0100187 0 0 0 0 0 0 0
+0.75 -0.0100187 0 0 0 0 0 0 0
+
+0.74 0 0 0 0 0 0 0 0
+0.75 0 0 0 0 0 0 0 0
+
+
+0.75 -0.0100187 0 0 0 0 0 0 0
+0.76 -0.0100187 0 0 0 0 0 0 0
+
+0.75 0 0 0 0 0 0 0 0
+0.76 0 0 0 0 0 0 0 0
+
+
+0.76 -0.0100187 0.7625 0 3.20997 1.72653e+07 -1 -1 1
+0.77 -0.0100187 0.7625 0 3.20997 1.72653e+07 -1 -1 1
+
+0.76 0 0.7625 0 3.20997 1.72653e+07 -1 -1 1
+0.77 0 0.7625 0 3.20997 1.72653e+07 -1 -1 1
+
+
+0.77 -0.0100187 0 0 0 0 0 0 0
+0.78 -0.0100187 0 0 0 0 0 0 0
+
+0.77 0 0 0 0 0 0 0 0
+0.78 0 0 0 0 0 0 0 0
+
+
+0.78 -0.0100187 0.7875 0 3.20074 1.5685e+07 -1 -1 1
+0.79 -0.0100187 0.7875 0 3.20074 1.5685e+07 -1 -1 1
+
+0.78 0 0.7875 0 3.20074 1.5685e+07 -1 -1 1
+0.79 0 0.7875 0 3.20074 1.5685e+07 -1 -1 1
+
+
+0.79 -0.0100187 0 0 0 0 0 0 0
+0.8 -0.0100187 0 0 0 0 0 0 0
+
+0.79 0 0 0 0 0 0 0 0
+0.8 0 0 0 0 0 0 0 0
+
+
+0.8 -0.0100187 0 0 0 0 0 0 0
+0.81 -0.0100187 0 0 0 0 0 0 0
+
+0.8 0 0 0 0 0 0 0 0
+0.81 0 0 0 0 0 0 0 0
+
+
+0.81 -0.0100187 0.8125 0 3.19598 1.40489e+07 -1 -1 1
+0.82 -0.0100187 0.8125 0 3.19598 1.40489e+07 -1 -1 1
+
+0.81 0 0.8125 0 3.19598 1.40489e+07 -1 -1 1
+0.82 0 0.8125 0 3.19598 1.40489e+07 -1 -1 1
+
+
+0.82 -0.0100187 0 0 0 0 0 0 0
+0.83 -0.0100187 0 0 0 0 0 0 0
+
+0.82 0 0 0 0 0 0 0 0
+0.83 0 0 0 0 0 0 0 0
+
+
+0.83 -0.0100187 0.8375 0 3.19365 1.2357e+07 -1 -1 1
+0.84 -0.0100187 0.8375 0 3.19365 1.2357e+07 -1 -1 1
+
+0.83 0 0.8375 0 3.19365 1.2357e+07 -1 -1 1
+0.84 0 0.8375 0 3.19365 1.2357e+07 -1 -1 1
+
+
+0.84 -0.0100187 0 0 0 0 0 0 0
+0.85 -0.0100187 0 0 0 0 0 0 0
+
+0.84 0 0 0 0 0 0 0 0
+0.85 0 0 0 0 0 0 0 0
+
+
+0.85 -0.0100187 0 0 0 0 0 0 0
+0.86 -0.0100187 0 0 0 0 0 0 0
+
+0.85 0 0 0 0 0 0 0 0
+0.86 0 0 0 0 0 0 0 0
+
+
+0.86 -0.0100187 0.8625 0 3.19238 1.06093e+07 -1 -1 1
+0.87 -0.0100187 0.8625 0 3.19238 1.06093e+07 -1 -1 1
+
+0.86 0 0.8625 0 3.19238 1.06093e+07 -1 -1 1
+0.87 0 0.8625 0 3.19238 1.06093e+07 -1 -1 1
+
+
+0.87 -0.0100187 0 0 0 0 0 0 0
+0.88 -0.0100187 0 0 0 0 0 0 0
+
+0.87 0 0 0 0 0 0 0 0
+0.88 0 0 0 0 0 0 0 0
+
+
+0.88 -0.0100187 0.8875 0 3.19146 8.80585e+06 -1 -1 1
+0.89 -0.0100187 0.8875 0 3.19146 8.80585e+06 -1 -1 1
+
+0.88 0 0.8875 0 3.19146 8.80585e+06 -1 -1 1
+0.89 0 0.8875 0 3.19146 8.80585e+06 -1 -1 1
+
+
+0.89 -0.0100187 0 0 0 0 0 0 0
+0.9 -0.0100187 0 0 0 0 0 0 0
+
+0.89 0 0 0 0 0 0 0 0
+0.9 0 0 0 0 0 0 0 0
+
+
+0.9 -0.0100187 0 0 0 0 0 0 0
+0.91 -0.0100187 0 0 0 0 0 0 0
+
+0.9 0 0 0 0 0 0 0 0
+0.91 0 0 0 0 0 0 0 0
+
+
+0.91 -0.0100187 0.9125 0 3.19083 6.94661e+06 -1 -1 1
+0.92 -0.0100187 0.9125 0 3.19083 6.94661e+06 -1 -1 1
+
+0.91 0 0.9125 0 3.19083 6.94661e+06 -1 -1 1
+0.92 0 0.9125 0 3.19083 6.94661e+06 -1 -1 1
+
+
+0.92 -0.0100187 0 0 0 0 0 0 0
+0.93 -0.0100187 0 0 0 0 0 0 0
+
+0.92 0 0 0 0 0 0 0 0
+0.93 0 0 0 0 0 0 0 0
+
+
+0.93 -0.0100187 0.9375 0 3.1911 5.03159e+06 -1 -1 1
+0.94 -0.0100187 0.9375 0 3.1911 5.03159e+06 -1 -1 1
+
+0.93 0 0.9375 0 3.1911 5.03159e+06 -1 -1 1
+0.94 0 0.9375 0 3.1911 5.03159e+06 -1 -1 1
+
+
+0.94 -0.0100187 0 0 0 0 0 0 0
+0.95 -0.0100187 0 0 0 0 0 0 0
+
+0.94 0 0 0 0 0 0 0 0
+0.95 0 0 0 0 0 0 0 0
+
+
+0.95 -0.0100187 0 0 0 0 0 0 0
+0.96 -0.0100187 0 0 0 0 0 0 0
+
+0.95 0 0 0 0 0 0 0 0
+0.96 0 0 0 0 0 0 0 0
+
+
+0.96 -0.0100187 0.9625 0 3.19353 3.06079e+06 -1 -1 1
+0.97 -0.0100187 0.9625 0 3.19353 3.06079e+06 -1 -1 1
+
+0.96 0 0.9625 0 3.19353 3.06079e+06 -1 -1 1
+0.97 0 0.9625 0 3.19353 3.06079e+06 -1 -1 1
+
+
+0.97 -0.0100187 0 0 0 0 0 0 0
+0.98 -0.0100187 0 0 0 0 0 0 0
+
+0.97 0 0 0 0 0 0 0 0
+0.98 0 0 0 0 0 0 0 0
+
+
+0.98 -0.0100187 0.9875 0 3.20005 1.03421e+06 -1 -1 1
+0.99 -0.0100187 0.9875 0 3.20005 1.03421e+06 -1 -1 1
+
+0.98 0 0.9875 0 3.20005 1.03421e+06 -1 -1 1
+0.99 0 0.9875 0 3.20005 1.03421e+06 -1 -1 1
+
+
+0.99 -0.0100187 0 0 0 0 0 0 0
+1 -0.0100187 0 0 0 0 0 0 0
+
+0.99 0 0 0 0 0 0 0 0
+1 0 0 0 0 0 0 0 0
+
+
+0 0 0 0 0 0 0 0 0
+0.01 0 0 0 0 0 0 0 0
+
+0 0.100187 0 0 0 0 0 0 0
+0.01 0.100187 0 0 0 0 0 0 0
+
+
+0.01 0 0.0514965 0 2908.1 2.81237e+07 -1 -1 1
+0.02 0 0.0514965 0 2908.1 2.81237e+07 -1 -1 1
+
+0.01 0.100187 0.0514965 0 2908.1 2.81237e+07 -1 -1 1
+0.02 0.100187 0.0514965 0 2908.1 2.81237e+07 -1 -1 1
+
+
+0.02 0 0 0 0 0 0 0 0
+0.03 0 0 0 0 0 0 0 0
+
+0.02 0.100187 0 0 0 0 0 0 0
+0.03 0.100187 0 0 0 0 0 0 0
+
+
+0.03 0 0.14652 0 2871.62 8.40585e+07 -1 -1 1
+0.04 0 0.14652 0 2871.62 8.40585e+07 -1 -1 1
+
+0.03 0.100187 0.14652 0 2871.62 8.40585e+07 -1 -1 1
+0.04 0.100187 0.14652 0 2871.62 8.40585e+07 -1 -1 1
+
+
+0.04 0 0 0 0 0 0 0 0
+0.05 0 0 0 0 0 0 0 0
+
More information about the CIG-COMMITS
mailing list