dynamic/kinematic faults enhancements

[mpbl]

• fault_solver_common:
  • make ordered version of subroutine assemble_MPI_vector_blocking, and use it in subroutine initialize_fault
• fault_solver_dynamic:
  • many hard-coded ad hoc features need to be set instead as user options (e.g. rate-and-state friction)
• rate-and-state friction:
  • make it a user-friendly option (currently a flag in fault_solver)
  • add documentation
• cubit interface:
  • consolidate python scripts (eliminate redundancy)
• verify consistency of fault edge nodes in kinematic solver: currently, these nodes are not split but they are included in the fault database
• simplify mesh generation with faults:
  • eliminate fault edge nodes? Currently fault edge nodes must be defined as closed and non-split in CUBIT, SESAME asigns a single ibulk, but it adds it to the fault database.
  • develop a recipe for mesh coarsening in CUBIT based on half "cubed sphere" mesh
  • test tet-to-hex strategy
• add snapshot outputs to kinematic solver (e.g. to export fault stresses for adjoint source inversion)
• prepare automated tests

[komatits]

Only Jean-Paul Ampuero's @jpampuero group at Caltech could address and implement that; it is maybe already done, we should check with them.

[jpampuero]

We have completed many of those in an older version, but my team still has to "migrate" to Git.

[komatits]

Great, thanks. This Wiki can help: https://github.com/geodynamics/specfem3d/wiki/Git-tips-for-developers

Best wishes, Dimitri.

On 10/04/2014 18:21, jpampuero wrote:

We have completed many of those in an older version, but my team still has to "migrate" to Git.

— Reply to this email directly or view it on GitHub https://github.com/geodynamics/specfem3d/issues/30#issuecomment-40106260.

Dimitri Komatitsch CNRS Research Director (DR CNRS), Laboratory of Mechanics and Acoustics, UPR 7051, Marseille, France http://komatitsch.free.fr

[komatits]

Done, now Jean-Paul's group has committed their nice improvements to that part. Thanks!