ValueError: Requested fields not available for output

PyLith
1

Dear pylith team, I am working based on step03 in subduction-2d. I am trying to simulate several seismic cycles, but encountered the following error:

ValueError: Requested fields not available for output.
Data provider: ‘faultcohesivekin’
Field type: ‘vertex’
Data type: ‘info’
Available fields: ‘normal_dir’ ‘final_slip’ ‘slip_time’ ‘final_slip_eq1’ ‘slip_time_eq1’ ‘final_slip_eq2’ ‘slip_time_eq2’ ‘final_slip_creep’ ‘slip_time_creep’ ‘strike_dir’
Fields not available: ‘final_slip_rupture’ ‘slip_time_rupture’
application called MPI_Abort(MPI_COMM_WORLD, -1) - process 0
/home/crl/Desktop/pylith/pylith-2.2.2-linux-x86_64/bin/nemesis: mpirun: exit 255

My .cfg file:

Set earthquake sources to an array consisting of creep and 3 ruptures.

eq_srcs = [creep,eq1,eq2]
eq_srcs.creep.origin_time = 0.0year
eq_srcs.eq1.origin_time = 150.0year
eq_srcs.eq2.origin_time = 300.0*year

Coseismic eq1 ---------

[pylithapp.timedependent.interfaces.fault_locked.eq_srcs.eq1.slip_function]
slip.iohandler.filename = fault_slip_coseismic.spatialdb
slip.query_type = linear
slip.label = Final_slip_eq1

Slip time (relative to origin time) is uniform, so use UniformDB for convenience

slip_time = spatialdata.spatialdb.UniformDB
slip_time.label = Slip_time_eq1
slip_time.values = [slip-time]
slip_time.data = [0.0*year]

Coseismic eq2 ---------

[pylithapp.timedependent.interfaces.fault_locked.eq_srcs.eq2.slip_function]
slip.iohandler.filename = fault_slip_coseismic.spatialdb
slip.query_type = linear
slip.label = Final_slip_eq2

Slip time (relative to origin time) is uniform, so use UniformDB for convenience

slip_time = spatialdata.spatialdb.UniformDB
slip_time.label = Slip_time_eq2
slip_time.values = [slip-time]
slip_time.data = [0.0*year]

Creep ---------

[pylithapp.timedependent.interfaces.fault_locked.eq_srcs.creep]
slip_function = pylith.faults.ConstRateSlipFn

The slip time and final slip are defined in spatial databases.

[pylithapp.timedependent.interfaces.fault_locked.eq_srcs.creep.slip_function]
slip_rate.iohandler.filename = fault_locked.spatialdb
slip_rate.query_type = linear
slip_rate.label = Final slip

Slip time is uniform, so use UniformDB for convenience

slip_time = spatialdata.spatialdb.UniformDB
slip_time.label = Slip time
slip_time.values = [slip-time]
slip_time.data = [0.0*year]

----------------------------------------------------------------------

output

----------------------------------------------------------------------

Domain

[pylithapp.problem.formulation.output.domain]
writer.filename = output/step03.h5

Ground surface

[pylithapp.problem.formulation.output.subdomain]
writer.filename = output/step03-groundsurf.h5

Faults

#[pylithapp.problem.interfaces.fault_locked.output]
writer = pylith.meshio.DataWriterHDF5
writer.filename = output/step03-fault-locked.h5

Materials

[pylithapp.timedependent.materials.concrust.output]
writer.filename = output/step03-concrust.h5

[pylithapp.timedependent.materials.conmantle.output]
writer.filename = output/step03-conmantle.h5

[pylithapp.timedependent.materials.oceancrust.output]
writer.filename = output/step03-oceancrust.h5

[pylithapp.timedependent.materials.oceanmantle.output]
writer.filename = output/step03-oceanmantle.h5

End of file

Thank you very much for any help.

2

Available fields: ‘normal_dir’ ‘final_slip’ ‘slip_time’ ‘final_slip_eq1’ ‘slip_time_eq1’ ‘final_slip_eq2’ ‘slip_time_eq2’ ‘final_slip_creep’ ‘slip_time_creep’ ‘strike_dir’
Fields not available: ‘final_slip_rupture’ ‘slip_time_rupture’

In PyLith v2 wIth multiple ruptures, the default values for the info_fields will not work, so you should specify what values you want in the info files. I would put all of the ones available (as given in the error message).

[pylithapp.problem.interfaces.fault_locked.output]
vertex_info_fields = [ normal_dir, final_slip, slip_time, final_slip_eq1, slip_time_eq1, final_slip_eq2, slip_time_eq2 ]
3

Thank you for your valuable reply. The question has been resolved.
Now I noticed that the groundsurf.h5 was not generated. What causes this problem?
Thanks!

4

I have solved this problem. Thanks again.