Skip to content

updates mass matrix usage; updates Par_files #1871

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
danielpeter merged 4 commits into from
May 14, 2026
Merged

updates mass matrix usage; updates Par_files #1871

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
2bacee7
updates Par_files (adding IS_WAVEFIELD_DISCONTINUITY & ROTATION)
danielpeter May 12, 2026
49389d7
updates Par_file ordering & reading of ROTATION parameter
danielpeter May 12, 2026
95e1028
updates xmeshfem3D stability criterion (to be consistent with check_m…
danielpeter May 14, 2026
ea26cf6
updates mass matrix storage and usages
danielpeter May 14, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
3 changes: 2 additions & 1 deletion DATA/Par_file
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,9 +75,10 @@ SEP_MODEL_DIRECTORY = DATA/my_SEP_model/
# parameters describing the model
APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.

# in case of attenuation, reference frequency in Hz at which the velocity values in the velocity model are given (unused otherwise)
ATTENUATION_f0_REFERENCE = 18.d0
Expand Down
10 changes: 9 additions & 1 deletion EXAMPLES/applications/CPML_examples/CPML_realistic_mesh_with_topography/DATA/Par_file
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,9 +75,10 @@ SEP_MODEL_DIRECTORY = ./DATA/my_SEP_model/
# parameters describing the model
APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.

# in case of attenuation, reference frequency in Hz at which the velocity values in the velocity model are given (unused otherwise)
ATTENUATION_f0_REFERENCE = 0.33333d0
Expand Down Expand Up @@ -324,6 +325,13 @@ TRACTION_PATH = ./DATA/AxiSEM_tractions/3/
FKMODEL_FILE = FKmodel
RECIPROCITY_AND_KH_INTEGRAL = .false. # does not work yet

#-----------------------------------------------------------
##
## Prescribed wavefield discontinuity on an interface
##
##-----------------------------------------------------------
IS_WAVEFIELD_DISCONTINUITY = .false.

#-----------------------------------------------------------
#
# Run modes
Expand Down
10 changes: 9 additions & 1 deletion ...ons/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_5sides/DATA/Par_file
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,9 +75,10 @@ SEP_MODEL_DIRECTORY = ./DATA/my_SEP_model/
# parameters describing the model
APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.

# in case of attenuation, reference frequency in Hz at which the velocity values in the velocity model are given (unused otherwise)
ATTENUATION_f0_REFERENCE = 0.33333d0
Expand Down Expand Up @@ -324,6 +325,13 @@ TRACTION_PATH = ./DATA/AxiSEM_tractions/3/
FKMODEL_FILE = FKmodel
RECIPROCITY_AND_KH_INTEGRAL = .false. # does not work yet

#-----------------------------------------------------------
##
## Prescribed wavefield discontinuity on an interface
##
##-----------------------------------------------------------
IS_WAVEFIELD_DISCONTINUITY = .false.

#-----------------------------------------------------------
#
# Run modes
Expand Down
10 changes: 9 additions & 1 deletion ...ons/CPML_examples/homogeneous_halfspace_HEX8_acoustic_absorbing_CPML_6sides/DATA/Par_file
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,9 +75,10 @@ SEP_MODEL_DIRECTORY = ./DATA/my_SEP_model/
# parameters describing the model
APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.

# in case of attenuation, reference frequency in Hz at which the velocity values in the velocity model are given (unused otherwise)
ATTENUATION_f0_REFERENCE = 0.33333d0
Expand Down Expand Up @@ -324,6 +325,13 @@ TRACTION_PATH = ./DATA/AxiSEM_tractions/3/
FKMODEL_FILE = FKmodel
RECIPROCITY_AND_KH_INTEGRAL = .false. # does not work yet

#-----------------------------------------------------------
##
## Prescribed wavefield discontinuity on an interface
##
##-----------------------------------------------------------
IS_WAVEFIELD_DISCONTINUITY = .false.

#-----------------------------------------------------------
#
# Run modes
Expand Down
10 changes: 9 additions & 1 deletion ..._examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_5sides/DATA/Par_file
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,9 +75,10 @@ SEP_MODEL_DIRECTORY = ./DATA/my_SEP_model/
# parameters describing the model
APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.

# in case of attenuation, reference frequency in Hz at which the velocity values in the velocity model are given (unused otherwise)
ATTENUATION_f0_REFERENCE = 0.33333d0
Expand Down Expand Up @@ -324,6 +325,13 @@ TRACTION_PATH = ./DATA/AxiSEM_tractions/3/
FKMODEL_FILE = FKmodel
RECIPROCITY_AND_KH_INTEGRAL = .false. # does not work yet

#-----------------------------------------------------------
##
## Prescribed wavefield discontinuity on an interface
##
##-----------------------------------------------------------
IS_WAVEFIELD_DISCONTINUITY = .false.

#-----------------------------------------------------------
#
# Run modes
Expand Down
10 changes: 9 additions & 1 deletion ..._examples/homogeneous_halfspace_HEX8_acoustic_elastic_absorbing_CPML_6sides/DATA/Par_file
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,9 +75,10 @@ SEP_MODEL_DIRECTORY = ./DATA/my_SEP_model/
# parameters describing the model
APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.

# in case of attenuation, reference frequency in Hz at which the velocity values in the velocity model are given (unused otherwise)
ATTENUATION_f0_REFERENCE = 0.33333d0
Expand Down Expand Up @@ -324,6 +325,13 @@ TRACTION_PATH = ./DATA/AxiSEM_tractions/3/
FKMODEL_FILE = FKmodel
RECIPROCITY_AND_KH_INTEGRAL = .false. # does not work yet

#-----------------------------------------------------------
##
## Prescribed wavefield discontinuity on an interface
##
##-----------------------------------------------------------
IS_WAVEFIELD_DISCONTINUITY = .false.

#-----------------------------------------------------------
#
# Run modes
Expand Down
10 changes: 9 additions & 1 deletion ...ions/CPML_examples/homogeneous_halfspace_HEX8_elastic_absorbing_CPML_5sides/DATA/Par_file
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,9 +75,10 @@ SEP_MODEL_DIRECTORY = ./DATA/my_SEP_model/
# parameters describing the model
APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.

# in case of attenuation, reference frequency in Hz at which the velocity values in the velocity model are given (unused otherwise)
ATTENUATION_f0_REFERENCE = 0.33333d0
Expand Down Expand Up @@ -324,6 +325,13 @@ TRACTION_PATH = ./DATA/AxiSEM_tractions/3/
FKMODEL_FILE = FKmodel
RECIPROCITY_AND_KH_INTEGRAL = .false. # does not work yet

#-----------------------------------------------------------
##
## Prescribed wavefield discontinuity on an interface
##
##-----------------------------------------------------------
IS_WAVEFIELD_DISCONTINUITY = .false.

#-----------------------------------------------------------
#
# Run modes
Expand Down
10 changes: 9 additions & 1 deletion ...ions/CPML_examples/homogeneous_halfspace_HEX8_elastic_absorbing_CPML_6sides/DATA/Par_file
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,9 +75,10 @@ SEP_MODEL_DIRECTORY = ./DATA/my_SEP_model/
# parameters describing the model
APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.

# in case of attenuation, reference frequency in Hz at which the velocity values in the velocity model are given (unused otherwise)
ATTENUATION_f0_REFERENCE = 0.33333d0
Expand Down Expand Up @@ -324,6 +325,13 @@ TRACTION_PATH = ./DATA/AxiSEM_tractions/3/
FKMODEL_FILE = FKmodel
RECIPROCITY_AND_KH_INTEGRAL = .false. # does not work yet

#-----------------------------------------------------------
##
## Prescribed wavefield discontinuity on an interface
##
##-----------------------------------------------------------
IS_WAVEFIELD_DISCONTINUITY = .false.

#-----------------------------------------------------------
#
# Run modes
Expand Down
10 changes: 9 additions & 1 deletion EXAMPLES/applications/Gmsh_simple_box_hex27/DATA/Par_file
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,9 +75,10 @@ SEP_MODEL_DIRECTORY = ./DATA/my_SEP_model/
# parameters describing the model
APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .true.
ANISOTROPY = .false.
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .true.

# in case of attenuation, reference frequency in Hz at which the velocity values in the velocity model are given (unused otherwise)
ATTENUATION_f0_REFERENCE = 1.0d0
Expand Down Expand Up @@ -324,6 +325,13 @@ TRACTION_PATH = ./DATA/AxiSEM_tractions/3/
FKMODEL_FILE = FKmodel
RECIPROCITY_AND_KH_INTEGRAL = .false. # does not work yet

#-----------------------------------------------------------
##
## Prescribed wavefield discontinuity on an interface
##
##-----------------------------------------------------------
IS_WAVEFIELD_DISCONTINUITY = .false.

#-----------------------------------------------------------
#
# Run modes
Expand Down
10 changes: 9 additions & 1 deletion EXAMPLES/applications/Gmsh_simple_lddrk/DATA/Par_file
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,9 +75,10 @@ SEP_MODEL_DIRECTORY = ./DATA/my_SEP_model/
# parameters describing the model
APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .true.
ANISOTROPY = .false.
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .true.

# in case of attenuation, reference frequency in Hz at which the velocity values in the velocity model are given (unused otherwise)
ATTENUATION_f0_REFERENCE = 0.33333d0
Expand Down Expand Up @@ -324,6 +325,13 @@ TRACTION_PATH = ./DATA/AxiSEM_tractions/3/
FKMODEL_FILE = FKmodel
RECIPROCITY_AND_KH_INTEGRAL = .false. # does not work yet

#-----------------------------------------------------------
##
## Prescribed wavefield discontinuity on an interface
##
##-----------------------------------------------------------
IS_WAVEFIELD_DISCONTINUITY = .false.

#-----------------------------------------------------------
#
# Run modes
Expand Down
10 changes: 9 additions & 1 deletion EXAMPLES/applications/LTS_homogeneous_halfspace_HEX8/DATA/Par_file
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,9 +75,10 @@ SEP_MODEL_DIRECTORY = DATA/my_SEP_model/
# parameters describing the model
APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.

# in case of attenuation, reference frequency in Hz at which the velocity values in the velocity model are given (unused otherwise)
ATTENUATION_f0_REFERENCE = 18.d0
Expand Down Expand Up @@ -324,6 +325,13 @@ TRACTION_PATH = DATA/AxiSEM_tractions/3/
FKMODEL_FILE = FKmodel
RECIPROCITY_AND_KH_INTEGRAL = .false. # does not work yet

#-----------------------------------------------------------
##
## Prescribed wavefield discontinuity on an interface
##
##-----------------------------------------------------------
IS_WAVEFIELD_DISCONTINUITY = .false.

#-----------------------------------------------------------
#
# Run modes
Expand Down
10 changes: 9 additions & 1 deletion EXAMPLES/applications/cmt_force_hybrid/DATA/Par_file
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,9 +75,10 @@ SEP_MODEL_DIRECTORY = ./DATA/my_SEP_model/
# parameters describing the model
APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.

# in case of attenuation, reference frequency in Hz at which the velocity values in the velocity model are given (unused otherwise)
ATTENUATION_f0_REFERENCE = 1.d0
Expand Down Expand Up @@ -330,6 +331,13 @@ TRACTION_PATH = ./DATA/AxiSEM_tractions/3/
FKMODEL_FILE = FKmodel
RECIPROCITY_AND_KH_INTEGRAL = .false. # does not work yet

#-----------------------------------------------------------
##
## Prescribed wavefield discontinuity on an interface
##
##-----------------------------------------------------------
IS_WAVEFIELD_DISCONTINUITY = .false.

#-----------------------------------------------------------
#
# Run modes
Expand Down
10 changes: 9 additions & 1 deletion EXAMPLES/applications/coffee_mug_with_fluid_inside/DATA/Par_file
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,9 +75,10 @@ SEP_MODEL_DIRECTORY = DATA/my_SEP_model/
# parameters describing the model
APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.

# in case of attenuation, reference frequency in Hz at which the velocity values in the velocity model are given (unused otherwise)
ATTENUATION_f0_REFERENCE = 18.d0
Expand Down Expand Up @@ -324,6 +325,13 @@ TRACTION_PATH = DATA/AxiSEM_tractions/3/
FKMODEL_FILE = FKmodel
RECIPROCITY_AND_KH_INTEGRAL = .false. # does not work yet

#-----------------------------------------------------------
##
## Prescribed wavefield discontinuity on an interface
##
##-----------------------------------------------------------
IS_WAVEFIELD_DISCONTINUITY = .false.

#-----------------------------------------------------------
#
# Run modes
Expand Down
10 changes: 9 additions & 1 deletion EXAMPLES/applications/decompose_mesh_MPI/DATA/Par_file
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,9 +75,10 @@ SEP_MODEL_DIRECTORY = ./DATA/my_SEP_model/
# parameters describing the model
APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.

# in case of attenuation, reference frequency in Hz at which the velocity values in the velocity model are given (unused otherwise)
ATTENUATION_f0_REFERENCE = 18.d0
Expand Down Expand Up @@ -324,6 +325,13 @@ TRACTION_PATH = ./DATA/AxiSEM_tractions/3/
FKMODEL_FILE = FKmodel
RECIPROCITY_AND_KH_INTEGRAL = .false. # does not work yet

#-----------------------------------------------------------
##
## Prescribed wavefield discontinuity on an interface
##
##-----------------------------------------------------------
IS_WAVEFIELD_DISCONTINUITY = .false.

#-----------------------------------------------------------
#
# Run modes
Expand Down
10 changes: 9 additions & 1 deletion EXAMPLES/applications/decompose_mesh_MPI_with_faults/DATA/Par_file
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -75,9 +75,10 @@ SEP_MODEL_DIRECTORY = ./DATA/my_SEP_model/
# parameters describing the model
APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
ROTATION = .false.
ATTENUATION = .false.

# in case of attenuation, reference frequency in Hz at which the velocity values in the velocity model are given (unused otherwise)
ATTENUATION_f0_REFERENCE = 0.33333d0
Expand Down Expand Up @@ -324,6 +325,13 @@ TRACTION_PATH = ./DATA/AxiSEM_tractions/3/
FKMODEL_FILE = FKmodel
RECIPROCITY_AND_KH_INTEGRAL = .false. # does not work yet

#-----------------------------------------------------------
##
## Prescribed wavefield discontinuity on an interface
##
##-----------------------------------------------------------
IS_WAVEFIELD_DISCONTINUITY = .false.

#-----------------------------------------------------------
#
# Run modes
Expand Down
Loading
Loading