[cig-commits] r18555 - short/3D/PyLith/trunk/examples/3d/hex8

[willic3 at geodynamics.org]

Author: willic3
Date: 2011-06-08 14:57:26 -0700 (Wed, 08 Jun 2011)
New Revision: 18555

Modified:
   short/3D/PyLith/trunk/examples/3d/hex8/step06.cfg
   short/3D/PyLith/trunk/examples/3d/hex8/step07.cfg
   short/3D/PyLith/trunk/examples/3d/hex8/step08.cfg
   short/3D/PyLith/trunk/examples/3d/hex8/step09.cfg
Log:
Altered examples to use HDF5 output.



Modified: short/3D/PyLith/trunk/examples/3d/hex8/step06.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/step06.cfg	2011-06-08 21:40:09 UTC (rev 18554)
+++ short/3D/PyLith/trunk/examples/3d/hex8/step06.cfg	2011-06-08 21:57:26 UTC (rev 18555)
@@ -168,54 +168,53 @@
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
+# Give basename for output of solution over domain.
 [pylithapp.problem.formulation.output.domain]
 # We specify that output occurs in terms of a given time frequency, and
-# ask for output every 50 years. The time stamps of the output files are
-# in years (rather than the default of seconds), and we give a format for
-# the time stamp.
+# ask for output every 50 years.
 output_freq = time_step
 time_step = 50.0*year
-writer.filename = output/step06.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
 
-# Give basename for VTK domain output of solution over ground surface.
+# We are using HDF5 output so we must change the default writer.
+writer = pylith.meshio.DataWriterHDF5Mesh
+writer.filename = output/step06.h5
+
+# Give basename for output of solution over ground surface.
 [pylithapp.problem.formulation.output.subdomain]
 # Name of nodeset for ground surface.
 label = face_zpos
 # We keep the default output frequency behavior (skip every n steps), and
 # ask to skip 0 steps between output, so that we get output every time step.
+# We again switch the writer to produce HDF5 output.
+# Note that we specifically ask for a submesh writer.
 skip = 0
-writer.filename = output/step06-groundsurf.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+writer = pylith.meshio.DataWriterHDF5SubMesh
+writer.filename = output/step06-groundsurf.h5
 
-# Give basename for VTK fault rupture output.
+# Give basename for fault rupture output.
 [pylithapp.problem.interfaces.fault.output]
 # We keep the default output frequency behavior (skip every n steps), and
 # ask to skip 0 steps between output, so that we get output every time step.
+# We again switch the writer to produce HDF5 output.
+# Note that we specifically ask for a subsubmesh writer.
 skip = 0
-writer.filename = output/step06-fault.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+writer = pylith.meshio.DataWriterHDF5SubSubMesh
+writer.filename = output/step06-fault.h5
 
-# Give basename for VTK output of upper_crust state variables.
+# Give basename for output of upper_crust state variables.
 [pylithapp.timedependent.materials.upper_crust.output]
 # Average values over quadrature points.
 cell_filter = pylith.meshio.CellFilterAvgMesh
 output_freq = time_step
 time_step = 50.0*year
-writer.filename = output/step06-upper_crust.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+writer = pylith.meshio.DataWriterHDF5Mesh
+writer.filename = output/step06-upper_crust.h5
 
-# Give basename for VTK output of lower_crust state variables.
+# Give basename for output of lower_crust state variables.
 [pylithapp.timedependent.materials.lower_crust.output]
 # Average values over quadrature points.
 cell_filter = pylith.meshio.CellFilterAvgMesh
 output_freq = time_step
 time_step = 50.0*year
-writer.filename = output/step06-lower_crust.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+writer = pylith.meshio.DataWriterHDF5Mesh
+writer.filename = output/step06-lower_crust.h5

Modified: short/3D/PyLith/trunk/examples/3d/hex8/step07.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/step07.cfg	2011-06-08 21:40:09 UTC (rev 18554)
+++ short/3D/PyLith/trunk/examples/3d/hex8/step07.cfg	2011-06-08 21:57:26 UTC (rev 18555)
@@ -182,21 +182,20 @@
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
+# Give basename for output of solution over domain.
 [pylithapp.problem.formulation.output.domain]
 # We specify that output occurs in terms of a given time frequency, and
-# ask for output every 50 years. The time stamps of the output files are
-# in years (rather than the default of seconds), and we give a format for
-# the time stamp.
+# ask for output every 50 years.
 # We also request velocity output in addition to displacements.
 vertex_data_fields = [displacement,velocity]
 output_freq = time_step
 time_step = 50.0*year
-writer.filename = output/step07.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
 
-# Give basename for VTK domain output of solution over ground surface.
+# We are using HDF5 output so we must change the default writer.
+writer = pylith.meshio.DataWriterHDF5Mesh
+writer.filename = output/step07.h5
+
+# Give basename for output of solution over ground surface.
 [pylithapp.problem.formulation.output.subdomain]
 # Name of nodeset for ground surface.
 label = face_zpos
@@ -205,35 +204,41 @@
 # We keep the default output frequency behavior (skip every n steps), and
 # ask to skip 0 steps between output, so that we get output every time step.
 skip = 0
-writer.filename = output/step07-groundsurf.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+# We again switch the writer to produce HDF5 output.
+# Note that we specifically ask for a submesh writer.
+writer = pylith.meshio.DataWriterHDF5SubMesh
+writer.filename = output/step07-groundsurf.h5
 
-# Give basename for VTK fault rupture output.
+# Give basename for fault rupture output.
 [pylithapp.problem.interfaces.fault.output]
 # We keep the default output frequency behavior (skip every n steps), and
 # ask to skip 0 steps between output, so that we get output every time step.
 skip = 0
-writer.filename = output/step07-fault.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+# We again switch the writer to produce HDF5 output.
+# Note that we specifically ask for a subsubmesh writer.
+writer = pylith.meshio.DataWriterHDF5SubSubMesh
+writer.filename = output/step07-fault.h5
 
-# Give basename for VTK output of upper_crust state variables.
+# Give basename for output of upper_crust state variables.
 [pylithapp.timedependent.materials.upper_crust.output]
 # Average values over quadrature points.
 cell_filter = pylith.meshio.CellFilterAvgMesh
+# We specify that output occurs in terms of a given time frequency, and
+# ask for output every 50 years.
 output_freq = time_step
 time_step = 50.0*year
-writer.filename = output/step07-upper_crust.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+# We again switch the writer to produce HDF5 output.
+writer = pylith.meshio.DataWriterHDF5Mesh
+writer.filename = output/step07-upper_crust.h5
 
-# Give basename for VTK output of lower_crust state variables.
+# Give basename for output of lower_crust state variables.
 [pylithapp.timedependent.materials.lower_crust.output]
 # Average values over quadrature points.
 cell_filter = pylith.meshio.CellFilterAvgMesh
+# We specify that output occurs in terms of a given time frequency, and
+# ask for output every 50 years.
 output_freq = time_step
 time_step = 50.0*year
-writer.filename = output/step07-lower_crust.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+# We again switch the writer to produce HDF5 output.
+writer = pylith.meshio.DataWriterHDF5Mesh
+writer.filename = output/step07-lower_crust.h5

Modified: short/3D/PyLith/trunk/examples/3d/hex8/step08.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/step08.cfg	2011-06-08 21:40:09 UTC (rev 18554)
+++ short/3D/PyLith/trunk/examples/3d/hex8/step08.cfg	2011-06-08 21:57:26 UTC (rev 18555)
@@ -200,54 +200,58 @@
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
+# Give basename for output of solution over domain.
 [pylithapp.problem.formulation.output.domain]
 # We specify that output occurs in terms of a given time frequency, and
-# ask for output every 50 years. The time stamps of the output files are
-# in years (rather than the default of seconds), and we give a format for
-# the time stamp.
+# ask for output every 50 years.
 output_freq = time_step
 time_step = 50.0*year
-writer.filename = output/step08.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+# We are using HDF5 output so we must change the default writer.
+writer = pylith.meshio.DataWriterHDF5Mesh
+writer.filename = output/step08.h5
 
-# Give basename for VTK domain output of solution over ground surface.
+# Give basename for output of solution over ground surface.
 [pylithapp.problem.formulation.output.subdomain]
 # Name of nodeset for ground surface.
 label = face_zpos
 # We keep the default output frequency behavior (skip every n steps), and
 # ask to skip 0 steps between output, so that we get output every time step.
 skip = 0
-writer.filename = output/step08-groundsurf.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+# We again switch the writer to produce HDF5 output.
+# Note that we specifically ask for a submesh writer.
+writer = pylith.meshio.DataWriterHDF5SubMesh
+writer.filename = output/step08-groundsurf.h5
 
-# Give basename for VTK fault rupture output.
+# Give basename for fault rupture output.
 [pylithapp.problem.interfaces.fault.output]
 # We keep the default output frequency behavior (skip every n steps), and
 # ask to skip 0 steps between output, so that we get output every time step.
 skip = 0
-writer.filename = output/step08-fault.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+# We again switch the writer to produce HDF5 output.
+# Note that we specifically ask for a subsubmesh writer.
+writer = pylith.meshio.DataWriterHDF5SubSubMesh
+writer.filename = output/step08-fault.h5
 
-# Give basename for VTK output of upper_crust state variables.
+# Give basename for output of upper_crust state variables.
 [pylithapp.timedependent.materials.upper_crust.output]
 # Average values over quadrature points.
 cell_filter = pylith.meshio.CellFilterAvgMesh
+# We specify that output occurs in terms of a given time frequency, and
+# ask for output every 50 years.
 output_freq = time_step
 time_step = 50.0*year
-writer.filename = output/step08-upper_crust.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+# We again switch the writer to produce HDF5 output.
+writer = pylith.meshio.DataWriterHDF5Mesh
+writer.filename = output/step08-upper_crust.h5
 
-# Give basename for VTK output of lower_crust state variables.
+# Give basename for output of lower_crust state variables.
 [pylithapp.timedependent.materials.lower_crust.output]
 # Average values over quadrature points.
 cell_filter = pylith.meshio.CellFilterAvgMesh
+# We specify that output occurs in terms of a given time frequency, and
+# ask for output every 50 years.
 output_freq = time_step
 time_step = 50.0*year
-writer.filename = output/step08-lower_crust.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+# We again switch the writer to produce HDF5 output.
+writer = pylith.meshio.DataWriterHDF5Mesh
+writer.filename = output/step08-lower_crust.h5

Modified: short/3D/PyLith/trunk/examples/3d/hex8/step09.cfg
===================================================================
--- short/3D/PyLith/trunk/examples/3d/hex8/step09.cfg	2011-06-08 21:40:09 UTC (rev 18554)
+++ short/3D/PyLith/trunk/examples/3d/hex8/step09.cfg	2011-06-08 21:57:26 UTC (rev 18555)
@@ -200,54 +200,58 @@
 # ----------------------------------------------------------------------
 # output
 # ----------------------------------------------------------------------
-# Give basename for VTK domain output of solution over domain.
+# Give basename for output of solution over domain.
 [pylithapp.problem.formulation.output.domain]
 # We specify that output occurs in terms of a given time frequency, and
-# ask for output every 50 years. The time stamps of the output files are
-# in years (rather than the default of seconds), and we give a format for
-# the time stamp.
+# ask for output every 50 years.
 output_freq = time_step
 time_step = 50.0*year
-writer.filename = output/step09.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+# We are using HDF5 output so we must change the default writer.
+writer = pylith.meshio.DataWriterHDF5Mesh
+writer.filename = output/step09.h5
 
-# Give basename for VTK domain output of solution over ground surface.
+# Give basename for output of solution over ground surface.
 [pylithapp.problem.formulation.output.subdomain]
 # Name of nodeset for ground surface.
 label = face_zpos
 # We keep the default output frequency behavior (skip every n steps), and
 # ask to skip 0 steps between output, so that we get output every time step.
 skip = 0
-writer.filename = output/step09-groundsurf.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+# We again switch the writer to produce HDF5 output.
+# Note that we specifically ask for a submesh writer.
+writer = pylith.meshio.DataWriterHDF5SubMesh
+writer.filename = output/step09-groundsurf.h5
 
-# Give basename for VTK fault rupture output.
+# Give basename for fault rupture output.
 [pylithapp.problem.interfaces.fault.output]
 # We keep the default output frequency behavior (skip every n steps), and
 # ask to skip 0 steps between output, so that we get output every time step.
 skip = 0
-writer.filename = output/step09-fault.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+# We again switch the writer to produce HDF5 output.
+# Note that we specifically ask for a subsubmesh writer.
+writer = pylith.meshio.DataWriterHDF5SubSubMesh
+writer.filename = output/step09-fault.h5
 
-# Give basename for VTK output of upper_crust state variables.
+# Give basename for output of upper_crust state variables.
 [pylithapp.timedependent.materials.upper_crust.output]
 # Average values over quadrature points.
 cell_filter = pylith.meshio.CellFilterAvgMesh
+# We specify that output occurs in terms of a given time frequency, and
+# ask for output every 50 years.
 output_freq = time_step
 time_step = 50.0*year
-writer.filename = output/step09-upper_crust.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+# We again switch the writer to produce HDF5 output.
+writer = pylith.meshio.DataWriterHDF5Mesh
+writer.filename = output/step09-upper_crust.h5
 
-# Give basename for VTK output of lower_crust state variables.
+# Give basename for output of lower_crust state variables.
 [pylithapp.timedependent.materials.lower_crust.output]
 # Average values over quadrature points.
 cell_filter = pylith.meshio.CellFilterAvgMesh
+# We specify that output occurs in terms of a given time frequency, and
+# ask for output every 50 years.
 output_freq = time_step
 time_step = 50.0*year
-writer.filename = output/step09-lower_crust.vtk
-writer.time_format = %04.0f
-writer.time_constant = 1.0*year
+# We again switch the writer to produce HDF5 output.
+writer = pylith.meshio.DataWriterHDF5Mesh
+writer.filename = output/step09-lower_crust.h5