Skip to content

Update Fortran 3D Mesher #1673

@lsawade

Description

@lsawade
  • Add IDOMAIN_COSSERAT = 4 constant in fortran/meshfem3d/setup/constants.h.in alongside IDOMAIN_ACOUSTIC=1, IDOMAIN_ELASTIC=2, IDOMAIN_POROELASTIC=3
  • Add Cosserat material parsing branch in fortran/meshfem3d/meshfem3D/read_value_mesh_parameters.f90 (~lines 253–476). Define a new input line format distinguishable by domain_id == 4. Read the 8 parameters: rho, kappa, mu, nu, j, lambda_c, mu_c, nu_c (same as 2D in fortran/meshfem2d/read_material_table.f90 lines 143–150)
  • Increase NUMBER_OF_MATERIAL_PROPERTIES in fortran/meshfem3d/meshfem3D/constants_meshfem3D.h from 18 if needed to accommodate the 8 Cosserat parameters + domain ID
  • Extend material writing in fortran/meshfem3d/meshfem3D/save_databases.F90 (~lines 170–250) to serialize Cosserat parameters into the matpropl(:) array with the IDOMAIN_COSSERAT identifier
  • Ensure material_properties array in fortran/meshfem3d/meshfem3D/meshfem3D_par.f90 can store the Cosserat parameters

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions