[CIG-SHORT] running TPV102 benchmark with version 1.9
Brad Aagaard
baagaard at usgs.gov
Wed Aug 28 17:05:15 PDT 2013
Yuval,
The benchmarks repository is setup for the PyLith development version
(v2.0). To run the benchmark with PyLith 1.9, use the attached
pylithapp.cfg file, which uses data writer settings available in v1.9
(things get simpler in v2.0).
Regards,
Brad
On 08/28/2013 04:40 PM, Yuval Tal wrote:
> Hello,
>
> I downloaded the TPV102 benchmark from the website, and when I'm trying to run it I get the following massage:
> -- pyre.inventory(error)
> -- meshimporter.partitioner.data_writer <- 'pylith.meshio.DataWriterHDF5Ext'
> -- could not bind facility 'data_writer': no factory named 'data_writer' in 'pylith.meshio.DataWriterHDF5Ext'
> etc.
>
> Is it possible to run TPV102 benchmark with version 1.9?
> Do I have to change some parameters to run it with version 1.9?
>
> Thanks,
> Yuval
>
>
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>
-------------- next part --------------
[pylithapp]
#initialize_only = True
# ----------------------------------------------------------------------
# journal
# ----------------------------------------------------------------------
[pylithapp.journal.info]
timedependent = 1
explicitlumped = 1
petsc = 1
meshiocubit = 1
elasticityexplicit = 1
elasticityexplicittet4 = 1
meshimporter = 1
faultcohesivedyn = 1
[pylithapp.journal.debug]
#pylithapp = 1
#problem = 1
#explicit = 1
# ----------------------------------------------------------------------
# mesh_generator
# ----------------------------------------------------------------------
[pylithapp.mesh_generator]
reader = pylith.meshio.MeshIOCubit
reorder_mesh = True
distributor.partitioner = parmetis
distributor.write_partition = True
distributor.data_writer = pylith.meshio.DataWriterHDF5ExtMesh
distributor.data_writer.filename = output/partition.h5
[pylithapp.mesh_generator.reader]
coordsys.space_dim = 3
# ----------------------------------------------------------------------
# problem
# ----------------------------------------------------------------------
[pylithapp.timedependent]
dimension = 3
elastic_prestep = False
formulation = pylith.problems.Explicit
formulation.norm_viscosity = 0.3
normalizer = spatialdata.units.NondimElasticDynamic
normalizer.shear_wave_speed = 3333*m/s
[pylithapp.timedependent.formulation.time_step]
total_time = 12.001*s
dt = 0.0125*s
# ----------------------------------------------------------------------
# materials
# ----------------------------------------------------------------------
[pylithapp.timedependent]
materials = [elastic]
materials.elastic = pylith.materials.ElasticIsotropic3D
[pylithapp.timedependent.materials.elastic]
label = Elastic material
id = 1
db_properties.label = Elastic properties
db_properties.iohandler.filename = matprops.spatialdb
# ----------------------------------------------------------------------
# boundary conditions
# ----------------------------------------------------------------------
[pylithapp.timedependent]
bc = [x_neg,x_pos,y_neg,y_pos,z_neg]
bc.x_neg = pylith.bc.AbsorbingDampers
bc.x_pos = pylith.bc.AbsorbingDampers
bc.y_neg = pylith.bc.AbsorbingDampers
bc.y_pos = pylith.bc.AbsorbingDampers
bc.z_neg = pylith.bc.AbsorbingDampers
[pylithapp.timedependent.bc.x_pos]
label = face_xpos
db.label = Absorbing BC +x
db.iohandler.filename = matprops.spatialdb
[pylithapp.timedependent.bc.x_neg]
label = face_xneg
db.label = Absorbing BC -x
db.iohandler.filename = matprops.spatialdb
[pylithapp.timedependent.bc.y_pos]
label = face_ypos
db.label = Absorbing BC +y
db.iohandler.filename = matprops.spatialdb
[pylithapp.timedependent.bc.y_neg]
label = face_yneg
db.label = Absorbing BC -y
db.iohandler.filename = matprops.spatialdb
[pylithapp.timedependent.bc.z_neg]
label = face_zneg
up_dir = [0, 1, 0]
db.label = Absorbing BC -z
db.iohandler.filename = matprops.spatialdb
# ----------------------------------------------------------------------
# faults
# ----------------------------------------------------------------------
[pylithapp.timedependent]
interfaces = [fault]
interfaces.fault = pylith.faults.FaultCohesiveDyn
[pylithapp.timedependent.interfaces.fault]
id = 10
label = fault
friction = pylith.friction.RateStateAgeing
friction.label = Rate-state w/ageing
friction.db_properties = spatialdata.spatialdb.SimpleGridDB
friction.db_properties.label = Rate-state parameters
friction.db_properties.filename = friction.spatialdb
friction.db_properties.query_type = linear
friction.db_initial_state = spatialdata.spatialdb.SimpleGridDB
friction.db_initial_state.label = Initial value for state variable
friction.db_initial_state.filename = friction.spatialdb
friction.db_initial_state.query_type = linear
traction_perturbation = pylith.faults.TractPerturbation
[pylithapp.timedependent.interfaces.fault.traction_perturbation]
db_initial = spatialdata.spatialdb.UniformDB
db_initial.label = Initial fault tractions
db_initial.values = [traction-shear-leftlateral, traction-shear-updip, traction-normal]
db_initial.data = [-75.0*MPa, 0.0*MPa, -120.0*MPa]
db_change = spatialdata.spatialdb.SimpleGridDB
db_change.label = Traction perturbation
db_change.filename = traction_change.spatialdb
db_change.query_type = linear
th_change = spatialdata.spatialdb.TimeHistory
th_change.label = Traction perturbation time history
th_change.filename = traction_change.timedb
# ----------------------------------------------------------------------
# output
# ----------------------------------------------------------------------
[pylithapp.timedependent.formulation]
output = [subdomain]
output.subdomain = pylith.meshio.OutputSolnSubset
[pylithapp.problem.formulation.output.subdomain]
label = face_zpos
vertex_data_fields=[displacement,velocity]
output_freq = time_step
time_step = 0.04999*s
writer = pylith.meshio.DataWriterHDF5ExtSubMesh
[pylithapp.timedependent.interfaces.fault.output]
vertex_info_fields = [strike_dir,normal_dir,traction_initial_value,traction_change_in_value,traction_change_start_time,reference_friction_coefficient,reference_slip_rate,characteristic_slip_distance,constitutive_parameter_a,constitutive_parameter_b,cohesion]
vertex_data_fields = [slip,slip_rate,traction,state_variable]
output_freq = time_step
time_step = 0.04999*s
writer = pylith.meshio.DataWriterHDF5ExtSubSubMesh
[pylithapp.timedependent.materials.elastic.output]
cell_info_fields = []
cell_data_fields = []
cell_filter = pylith.meshio.CellFilterAvgMesh
writer = pylith.meshio.DataWriterHDF5ExtMesh
# ----------------------------------------------------------------------
# PETSc
# ----------------------------------------------------------------------
[pylithapp.petsc]
log_summary = true
More information about the CIG-SHORT
mailing list