From b3bf99bdb199232ff3d31d67be611c6b63daf0cd Mon Sep 17 00:00:00 2001 From: Daniel Peter Date: Mon, 13 May 2024 20:28:40 +0200 Subject: [PATCH 1/9] updates workflow script (macOS) --- .github/workflows/CI.yml | 6 ++++++ 1 file changed, 6 insertions(+) diff --git a/.github/workflows/CI.yml b/.github/workflows/CI.yml index e31803dfd..71b97ad5e 100644 --- a/.github/workflows/CI.yml +++ b/.github/workflows/CI.yml @@ -107,6 +107,12 @@ jobs: echo "PATH=/usr/local/opt/gnu-sed/libexec/gnubin:$PATH" >> $GITHUB_ENV ## OpenMP echo "OMP_NUM_THREADS=2" >> $GITHUB_ENV + ## MPI + echo "MPI environment:" + echo "mpif90 on path: $(which mpif90)" + echo + # add OpenMPI path to have ./configure detect mpi.h + echo "MPI_INC=/opt/homebrew/include" >> $GITHUB_ENV ## avoids MPI issue with number of slots echo "OMPI_MCA_rmaps_base_oversubscribe=1" >> $GITHUB_ENV echo "OMPI_MCA_rmaps_base_inherit=1" >> $GITHUB_ENV From 90a67c4df0ce4e78075b4bca6fe1cf91badd7d01 Mon Sep 17 00:00:00 2001 From: Daniel Peter Date: Tue, 14 May 2024 10:33:49 +0200 Subject: [PATCH 2/9] updates script comments (to avoid a SyntaxWarning: invalid escape sequence '\ ') --- utils/scripts/createTomographyFile.py | 29 ++++++++++++++++----------- 1 file changed, 17 insertions(+), 12 deletions(-) diff --git a/utils/scripts/createTomographyFile.py b/utils/scripts/createTomographyFile.py index d5a897482..3455512cf 100755 --- a/utils/scripts/createTomographyFile.py +++ b/utils/scripts/createTomographyFile.py @@ -7,28 +7,31 @@ vmin o----------------> velocity (m/s) - | | / - | | / - | | / gradient : -d1 - ZMIN |___|/ - | \ + | | // + | | // + | | // gradient : -d1 + ZMIN |___|// + | \\ | - - | \ gradient : d2 + | \\ gradient : d2 | - - | \ + | \\ Depth (m) v How to use external velocity model in specfem (03/20/2017) : -To set an arbitrary velocity model you have to use the option TOMOGRAPHY_FILE that I have implemented and that you have to set like that in the Par_file when you define the velocity model: +To set an arbitrary velocity model you have to use the option TOMOGRAPHY_FILE that I have implemented, +and you have to set it like that in the Par_file when you define the velocity model: 2 -1 9999 9999 9999 9999 9999 9999 9999 0 0 0 0 0 0 # external tomography file TOMOGRAPHY_FILE = ./profile.xyz -The program understand that the velocity model number 2 has to be read in the file profile.xyz as a regular grid it will then deal with the interpolation. You can set what you want in place of the 9999, it does not matter. -The file profile.xyz has to be written under the following format: +The program understands that the velocity model number 2 has to be read in the file profile.xyz as a regular grid; +it will then deal with the interpolation. You can set what you want in place of the 9999, it doesn't matter. + +The file profile.xyz has to be written using the following format: ! The xyz file TOMOGRAPHY_FILE that describe the tomography should be located in the TOMOGRAPHY_PATH ! directory, set in the Par_file. The format of the file, as read from define_external_model_from_xyz_file.f90 looks like : @@ -70,14 +73,16 @@ ! coordinate (from ORIGIN_X to END_X with step of SPACING_X) for each given z (from ORIGIN_Z ! to END_Z, with step of SPACING_Z). -This is a working script writing such file for a 1D profile (velocity depending just on depth). Read carefully all the comments but in particular the lines containing "!!WARNING!!" that I have added for +This is a working script, writing such a file for a 1D profile (velocity depending just on depth). +Read carefully all the comments, in particular the lines containing "!!WARNING!!" that I have added for people who need range dependent models. run that script ./createTomographyFile.py and look at the file profile.xyz created. In the case of range dependent models vp, vs and rho will be 2D arrays, you have to modify that. -For the moment just one model can be read from an external file but it is not very difficult to implement that. Likewise viscoelastic parameters can not vary with position for now. +For the moment just one model can be read from an external file, but it is not very difficult to implement that. +Likewise, viscoelastic parameters can not vary with position for now. @author: alexis bottero (alexis dot bottero at gmail dot com) """ From e15fc6183130e33f6e4b67e3209d42eb8745c483 Mon Sep 17 00:00:00 2001 From: Daniel Peter Date: Tue, 14 May 2024 11:23:29 +0200 Subject: [PATCH 3/9] updates workflow (sed path for macOS) --- .github/workflows/CI.yml | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/.github/workflows/CI.yml b/.github/workflows/CI.yml index 71b97ad5e..0d90e621c 100644 --- a/.github/workflows/CI.yml +++ b/.github/workflows/CI.yml @@ -104,10 +104,11 @@ jobs: echo "" ## avoids sed -i '' issue on MacOS, using gnu sed to have the same sed command lines as in linux brew install gnu-sed - echo "PATH=/usr/local/opt/gnu-sed/libexec/gnubin:$PATH" >> $GITHUB_ENV + echo "PATH=/opt/homebrew/opt/gnu-sed/libexec/gnubin:$PATH" >> $GITHUB_ENV ## OpenMP echo "OMP_NUM_THREADS=2" >> $GITHUB_ENV ## MPI + echo echo "MPI environment:" echo "mpif90 on path: $(which mpif90)" echo From 504bb07f652f9067ee1a37a09bc4fbed3a6160d9 Mon Sep 17 00:00:00 2001 From: Daniel Peter Date: Wed, 5 Jun 2024 10:43:20 +0200 Subject: [PATCH 4/9] updates github action checkout version --- .github/workflows/CI.yml | 34 +++++++++++++++++----------------- 1 file changed, 17 insertions(+), 17 deletions(-) diff --git a/.github/workflows/CI.yml b/.github/workflows/CI.yml index 0d90e621c..2be076522 100644 --- a/.github/workflows/CI.yml +++ b/.github/workflows/CI.yml @@ -11,7 +11,7 @@ jobs: runs-on: ubuntu-latest steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Check for changes id: diff @@ -86,7 +86,7 @@ jobs: needs: changesCheck steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install dependencies run: | @@ -146,7 +146,7 @@ jobs: os: [ubuntu-latest,ubuntu-20.04] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -165,7 +165,7 @@ jobs: needs: changesCheck steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Cache Intel oneapi packages id: cache-intel-oneapi @@ -301,7 +301,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -319,7 +319,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -343,7 +343,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -368,7 +368,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -393,7 +393,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -417,7 +417,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -441,7 +441,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -466,7 +466,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -490,7 +490,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -514,7 +514,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -538,7 +538,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -563,7 +563,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -587,7 +587,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh From ca6d7537a8cba3164954e3c3d100604af6f83517 Mon Sep 17 00:00:00 2001 From: Daniel Peter Date: Wed, 5 Jun 2024 10:59:55 +0200 Subject: [PATCH 5/9] updates Par_file comments --- DATA/Par_file | 4 +- EXAMPLES/BENCHIE_water_PVC/DATA/Par_file | 5 +- .../ACOUSTIC/DATA/Par_file | 5 +- .../ACOUSTIC_ELASTIC/DATA/Par_file | 5 +- .../ELASTIC/DATA/Par_file | 5 +- EXAMPLES/Gmsh_example_CPML_MPI/DATA/Par_file | 2 +- .../Gmsh_example_Stacey_MPI/DATA/Par_file | 5 +- .../DATA/Par_file.serial | 5 +- EXAMPLES/Industrial_Format_SEP/DATA/Par_file | 5 +- .../REF_SEIS/Par_file.step_1 | 5 +- .../REF_SEIS/Par_file.step_2 | 5 +- .../LuoYang_fluid_solid_kernel/DATA/Par_file | 5 +- EXAMPLES/Marmousi2/DATA/Par_file | 5 +- .../Marmousi_mesh_of_the_model/DATA/Par_file | 2 +- .../DATA/Par_file.mesh-improved | 2 +- .../REF_SEIS/Par_file | 2 +- EXAMPLES/Rayleigh_wave_no_crack/DATA/Par_file | 5 +- .../DATA/Par_file | 5 +- .../Rayleigh_wave_with_crack/DATA/Par_file | 5 +- .../DATA/Par_file_Tape2007_132rec_checker | 5 +- .../Tape2007/DATA/Par_file_Tape2007_onerec | 5 +- EXAMPLES/Tape2007/REF_SEIS/Par_file | 5 +- EXAMPLES/Tape2007_kernel/DATA/Par_file | 5 +- EXAMPLES/Tromp2005/DATA/Par_file | 5 +- .../Tromp2005/DATA/Par_file_Tromp2005_s100 | 5 +- EXAMPLES/Tromp2005_kernel/DATA/Par_file | 5 +- .../anisotropic_isotropic_model/DATA/Par_file | 2 +- .../Par_file_homogeneous_with_attenuation | 2 +- .../anisotropic_zinc_crystal/DATA/Par_file | 5 +- .../DATA/Par_file | 5 +- .../DATA/Par_file | 5 +- ...ation_2D_source_in_point_inside_an_element | 2 +- ...2D_source_in_point_shared_by_four_elements | 2 +- .../DATA/Par_file | 5 +- .../REF_SEIS/Par_file | 5 +- EXAMPLES/canyon/DATA/Par_file | 2 +- ...e_corner_between_several_spectral_elements | 5 +- ...nuation_2D_inside_a_given_spectral_element | 5 +- .../REF_SEIS/Par_file | 5 +- ...e_corner_between_several_spectral_elements | 5 +- .../REF_SEIS/Par_file | 5 +- ...e_corner_between_several_spectral_elements | 5 +- ...nuation_2D_inside_a_given_spectral_element | 5 +- .../REF_SEIS/Par_file | 5 +- EXAMPLES/electromagnetic_GPR_TM/DATA/Par_file | 2 +- .../fluid_solid_external_mesh/DATA/Par_file | 2 +- .../DATA/Par_file | 2 +- .../DATA/Par_file | 2 +- EXAMPLES/global_Earth_ak135f/DATA/Par_file | 5 +- .../DATA/Par_file | 2 +- .../DATA/Par_file | 4 +- EXAMPLES/initial_mode_LDDRK/DATA/Par_file | 4 +- .../DATA/Par_file | 2 +- .../moving_sources_acoustic/DATA/Par_file | 4 +- EXAMPLES/noise_uniform/DATA/Par_file | 2 +- EXAMPLES/noise_uniform/REF_SEIS/Par_file | 4 +- .../paper_axisymmetry_example/DATA/Par_file | 2 +- .../Pekeris_waveguide/DATA/Par_file | 2 +- .../with_free_bottom_surface/DATA/Par_file | 2 +- .../with_rigid_bottom_surface/DATA/Par_file | 5 +- EXAMPLES/poroelastic_acoustic/DATA/Par_file | 5 +- .../DATA/Par_file | 5 +- .../CPML_homogeneous/DATA/Par_file | 2 +- .../CPML_normal_fluid_solid/DATA/Par_file | 2 +- .../CPML_normal_solid_only/DATA/Par_file | 2 +- .../Stacey_homogeneous/DATA/Par_file | 2 +- .../Stacey_normal_fluid_solid/DATA/Par_file | 2 +- .../DATA/Par_file | 5 +- .../DATA/Par_file_adj | 5 +- .../DATA/Par_file_fwd | 5 +- .../DATA/Par_file_kernel | 5 +- .../DATA/Par_file_save_forward | 5 +- .../semi_infinite_homogeneous/DATA/Par_file | 5 +- .../DATA/Par_file | 5 +- .../REF_KERNEL/Par_file | 2 +- .../REF_SEIS/Par_file | 5 +- EXAMPLES/thermocline/DATA/Par_file | 2 +- .../tomographic_ocean_model/DATA/Par_file | 2 +- .../tomographic_ocean_model/REF_SEIS/Par_file | 2 +- .../Par_file_fluid_solid | 5 +- .../Par_file_fluid_solid | 5 +- .../DATA/Par_file | 5 +- ...file_to_use_with_this_timing_routine_HOWTO | 62 ++++++++++--------- 83 files changed, 217 insertions(+), 166 deletions(-) diff --git a/DATA/Par_file b/DATA/Par_file index e322ac11d..e9716e3a8 100644 --- a/DATA/Par_file +++ b/DATA/Par_file @@ -30,7 +30,7 @@ DT = 1.1d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/BENCHIE_water_PVC/DATA/Par_file b/EXAMPLES/BENCHIE_water_PVC/DATA/Par_file index 3de437d64..3b1715bca 100644 --- a/EXAMPLES/BENCHIE_water_PVC/DATA/Par_file +++ b/EXAMPLES/BENCHIE_water_PVC/DATA/Par_file @@ -30,7 +30,7 @@ DT = 4.5d-8 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file b/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file index 645bd3305..6f2a9abd6 100644 --- a/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file +++ b/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file @@ -30,7 +30,7 @@ DT = 0.00035 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -272,6 +272,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC_ELASTIC/DATA/Par_file b/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC_ELASTIC/DATA/Par_file index b4de65fe9..b5481f9fa 100644 --- a/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC_ELASTIC/DATA/Par_file +++ b/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC_ELASTIC/DATA/Par_file @@ -30,7 +30,7 @@ DT = 0.00025 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -272,6 +272,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file b/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file index 0a1be5a8f..206ed7060 100644 --- a/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file +++ b/EXAMPLES/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file @@ -30,7 +30,7 @@ DT = 0.00035 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -272,6 +272,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Gmsh_example_CPML_MPI/DATA/Par_file b/EXAMPLES/Gmsh_example_CPML_MPI/DATA/Par_file index 5ee802b61..a4f38b451 100644 --- a/EXAMPLES/Gmsh_example_CPML_MPI/DATA/Par_file +++ b/EXAMPLES/Gmsh_example_CPML_MPI/DATA/Par_file @@ -30,7 +30,7 @@ DT = 0.002 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/Gmsh_example_Stacey_MPI/DATA/Par_file b/EXAMPLES/Gmsh_example_Stacey_MPI/DATA/Par_file index 70325daa6..241c7e0ea 100644 --- a/EXAMPLES/Gmsh_example_Stacey_MPI/DATA/Par_file +++ b/EXAMPLES/Gmsh_example_Stacey_MPI/DATA/Par_file @@ -30,7 +30,7 @@ DT = 5.d-7 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Gmsh_example_Stacey_MPI/DATA/Par_file.serial b/EXAMPLES/Gmsh_example_Stacey_MPI/DATA/Par_file.serial index 5810e7217..f0ad6eefb 100644 --- a/EXAMPLES/Gmsh_example_Stacey_MPI/DATA/Par_file.serial +++ b/EXAMPLES/Gmsh_example_Stacey_MPI/DATA/Par_file.serial @@ -30,7 +30,7 @@ DT = 5.d-7 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Industrial_Format_SEP/DATA/Par_file b/EXAMPLES/Industrial_Format_SEP/DATA/Par_file index 8fd0d342e..ddf686235 100644 --- a/EXAMPLES/Industrial_Format_SEP/DATA/Par_file +++ b/EXAMPLES/Industrial_Format_SEP/DATA/Par_file @@ -30,7 +30,7 @@ DT = 5d-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Industrial_Format_SEP/REF_SEIS/Par_file.step_1 b/EXAMPLES/Industrial_Format_SEP/REF_SEIS/Par_file.step_1 index 990b09eca..9c209b117 100644 --- a/EXAMPLES/Industrial_Format_SEP/REF_SEIS/Par_file.step_1 +++ b/EXAMPLES/Industrial_Format_SEP/REF_SEIS/Par_file.step_1 @@ -30,7 +30,7 @@ DT = 5d-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Industrial_Format_SEP/REF_SEIS/Par_file.step_2 b/EXAMPLES/Industrial_Format_SEP/REF_SEIS/Par_file.step_2 index bbe1b89f7..d4e1c258c 100644 --- a/EXAMPLES/Industrial_Format_SEP/REF_SEIS/Par_file.step_2 +++ b/EXAMPLES/Industrial_Format_SEP/REF_SEIS/Par_file.step_2 @@ -30,7 +30,7 @@ DT = 5d-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/LuoYang_fluid_solid_kernel/DATA/Par_file b/EXAMPLES/LuoYang_fluid_solid_kernel/DATA/Par_file index 273d5f79e..a9a311fd4 100644 --- a/EXAMPLES/LuoYang_fluid_solid_kernel/DATA/Par_file +++ b/EXAMPLES/LuoYang_fluid_solid_kernel/DATA/Par_file @@ -30,7 +30,7 @@ DT = 1.1d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 4 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Marmousi2/DATA/Par_file b/EXAMPLES/Marmousi2/DATA/Par_file index 9a41c9e8b..b0b27e258 100644 --- a/EXAMPLES/Marmousi2/DATA/Par_file +++ b/EXAMPLES/Marmousi2/DATA/Par_file @@ -30,7 +30,7 @@ DT = 1.d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # marmousi: model_number 1 rho Vp Vs k 0 QKappa Qmu 0 0 0 0 0 0 (similar to elastic, except that the fourth parameter k accounts for the VELOCITY GRADIENT in Marmousi2) diff --git a/EXAMPLES/Marmousi_mesh_of_the_model/DATA/Par_file b/EXAMPLES/Marmousi_mesh_of_the_model/DATA/Par_file index 37be2f307..0b18dcc8b 100644 --- a/EXAMPLES/Marmousi_mesh_of_the_model/DATA/Par_file +++ b/EXAMPLES/Marmousi_mesh_of_the_model/DATA/Par_file @@ -30,7 +30,7 @@ DT = 5.d-6 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/Marmousi_mesh_of_the_model/DATA/Par_file.mesh-improved b/EXAMPLES/Marmousi_mesh_of_the_model/DATA/Par_file.mesh-improved index 29075f783..50e4985e6 100644 --- a/EXAMPLES/Marmousi_mesh_of_the_model/DATA/Par_file.mesh-improved +++ b/EXAMPLES/Marmousi_mesh_of_the_model/DATA/Par_file.mesh-improved @@ -30,7 +30,7 @@ DT = 5.d-6 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/Marmousi_mesh_of_the_model/REF_SEIS/Par_file b/EXAMPLES/Marmousi_mesh_of_the_model/REF_SEIS/Par_file index 37be2f307..0b18dcc8b 100644 --- a/EXAMPLES/Marmousi_mesh_of_the_model/REF_SEIS/Par_file +++ b/EXAMPLES/Marmousi_mesh_of_the_model/REF_SEIS/Par_file @@ -30,7 +30,7 @@ DT = 5.d-6 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/Rayleigh_wave_no_crack/DATA/Par_file b/EXAMPLES/Rayleigh_wave_no_crack/DATA/Par_file index cde32b61a..490bff77b 100644 --- a/EXAMPLES/Rayleigh_wave_no_crack/DATA/Par_file +++ b/EXAMPLES/Rayleigh_wave_no_crack/DATA/Par_file @@ -30,7 +30,7 @@ DT = 3.90625e-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Rayleigh_wave_no_crack_simpler_way_to_generate_the_wave/DATA/Par_file b/EXAMPLES/Rayleigh_wave_no_crack_simpler_way_to_generate_the_wave/DATA/Par_file index 43a8136f7..ba0a7fd2d 100644 --- a/EXAMPLES/Rayleigh_wave_no_crack_simpler_way_to_generate_the_wave/DATA/Par_file +++ b/EXAMPLES/Rayleigh_wave_no_crack_simpler_way_to_generate_the_wave/DATA/Par_file @@ -30,7 +30,7 @@ DT = 0.85d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Rayleigh_wave_with_crack/DATA/Par_file b/EXAMPLES/Rayleigh_wave_with_crack/DATA/Par_file index edfe6d5ac..12bd791a7 100644 --- a/EXAMPLES/Rayleigh_wave_with_crack/DATA/Par_file +++ b/EXAMPLES/Rayleigh_wave_with_crack/DATA/Par_file @@ -30,7 +30,7 @@ DT = 3.90625e-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 3 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Tape2007/DATA/Par_file_Tape2007_132rec_checker b/EXAMPLES/Tape2007/DATA/Par_file_Tape2007_132rec_checker index a0ca45c0d..6f1598767 100644 --- a/EXAMPLES/Tape2007/DATA/Par_file_Tape2007_132rec_checker +++ b/EXAMPLES/Tape2007/DATA/Par_file_Tape2007_132rec_checker @@ -30,7 +30,7 @@ DT = 6.0d-2 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .false. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Tape2007/DATA/Par_file_Tape2007_onerec b/EXAMPLES/Tape2007/DATA/Par_file_Tape2007_onerec index 112ff4c37..70e9ae385 100644 --- a/EXAMPLES/Tape2007/DATA/Par_file_Tape2007_onerec +++ b/EXAMPLES/Tape2007/DATA/Par_file_Tape2007_onerec @@ -30,7 +30,7 @@ DT = 6.0d-2 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .false. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Tape2007/REF_SEIS/Par_file b/EXAMPLES/Tape2007/REF_SEIS/Par_file index a0ca45c0d..6f1598767 100644 --- a/EXAMPLES/Tape2007/REF_SEIS/Par_file +++ b/EXAMPLES/Tape2007/REF_SEIS/Par_file @@ -30,7 +30,7 @@ DT = 6.0d-2 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .false. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Tape2007_kernel/DATA/Par_file b/EXAMPLES/Tape2007_kernel/DATA/Par_file index 5e5c5420a..bf99e4cfa 100644 --- a/EXAMPLES/Tape2007_kernel/DATA/Par_file +++ b/EXAMPLES/Tape2007_kernel/DATA/Par_file @@ -30,7 +30,7 @@ DT = 6.0d-2 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .false. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Tromp2005/DATA/Par_file b/EXAMPLES/Tromp2005/DATA/Par_file index eb2537872..5414ad56d 100644 --- a/EXAMPLES/Tromp2005/DATA/Par_file +++ b/EXAMPLES/Tromp2005/DATA/Par_file @@ -30,7 +30,7 @@ DT = 2.0d-2 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .false. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Tromp2005/DATA/Par_file_Tromp2005_s100 b/EXAMPLES/Tromp2005/DATA/Par_file_Tromp2005_s100 index 1d7f0127c..e8c7651bf 100644 --- a/EXAMPLES/Tromp2005/DATA/Par_file_Tromp2005_s100 +++ b/EXAMPLES/Tromp2005/DATA/Par_file_Tromp2005_s100 @@ -30,7 +30,7 @@ DT = 2.0d-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Tromp2005_kernel/DATA/Par_file b/EXAMPLES/Tromp2005_kernel/DATA/Par_file index a26cb7a1f..98008e0f6 100644 --- a/EXAMPLES/Tromp2005_kernel/DATA/Par_file +++ b/EXAMPLES/Tromp2005_kernel/DATA/Par_file @@ -30,7 +30,7 @@ DT = 2.0d-2 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/anisotropic_isotropic_model/DATA/Par_file b/EXAMPLES/anisotropic_isotropic_model/DATA/Par_file index 72ab0c30d..745d89f1e 100644 --- a/EXAMPLES/anisotropic_isotropic_model/DATA/Par_file +++ b/EXAMPLES/anisotropic_isotropic_model/DATA/Par_file @@ -30,7 +30,7 @@ DT = 0.8e-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/anisotropic_isotropic_model/DATA/Par_file_homogeneous_with_attenuation b/EXAMPLES/anisotropic_isotropic_model/DATA/Par_file_homogeneous_with_attenuation index 4bb3c32a5..c72021169 100644 --- a/EXAMPLES/anisotropic_isotropic_model/DATA/Par_file_homogeneous_with_attenuation +++ b/EXAMPLES/anisotropic_isotropic_model/DATA/Par_file_homogeneous_with_attenuation @@ -30,7 +30,7 @@ DT = 0.8e-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/anisotropic_zinc_crystal/DATA/Par_file b/EXAMPLES/anisotropic_zinc_crystal/DATA/Par_file index 70529b1f2..44fc87a9b 100644 --- a/EXAMPLES/anisotropic_zinc_crystal/DATA/Par_file +++ b/EXAMPLES/anisotropic_zinc_crystal/DATA/Par_file @@ -30,7 +30,7 @@ DT = 55.e-9 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/attenuation/viscoacoustic_attenuation_off_versus_analytical/DATA/Par_file b/EXAMPLES/attenuation/viscoacoustic_attenuation_off_versus_analytical/DATA/Par_file index 0bb1994d4..59ff92cda 100644 --- a/EXAMPLES/attenuation/viscoacoustic_attenuation_off_versus_analytical/DATA/Par_file +++ b/EXAMPLES/attenuation/viscoacoustic_attenuation_off_versus_analytical/DATA/Par_file @@ -30,7 +30,7 @@ DT = 1.3d-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/attenuation/viscoacoustic_attenuation_on_versus_analytical/DATA/Par_file b/EXAMPLES/attenuation/viscoacoustic_attenuation_on_versus_analytical/DATA/Par_file index 5b6183071..cdabb87b5 100644 --- a/EXAMPLES/attenuation/viscoacoustic_attenuation_on_versus_analytical/DATA/Par_file +++ b/EXAMPLES/attenuation/viscoacoustic_attenuation_on_versus_analytical/DATA/Par_file @@ -30,7 +30,7 @@ DT = 1.3d-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/attenuation/viscoelastic/DATA/Par_file_attenuation_2D_source_in_point_inside_an_element b/EXAMPLES/attenuation/viscoelastic/DATA/Par_file_attenuation_2D_source_in_point_inside_an_element index d56faa5c4..318e77b01 100644 --- a/EXAMPLES/attenuation/viscoelastic/DATA/Par_file_attenuation_2D_source_in_point_inside_an_element +++ b/EXAMPLES/attenuation/viscoelastic/DATA/Par_file_attenuation_2D_source_in_point_inside_an_element @@ -30,7 +30,7 @@ DT = 4e-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/attenuation/viscoelastic/DATA/Par_file_attenuation_2D_source_in_point_shared_by_four_elements b/EXAMPLES/attenuation/viscoelastic/DATA/Par_file_attenuation_2D_source_in_point_shared_by_four_elements index 5f3e9a4f8..bb1c53c95 100644 --- a/EXAMPLES/attenuation/viscoelastic/DATA/Par_file_attenuation_2D_source_in_point_shared_by_four_elements +++ b/EXAMPLES/attenuation/viscoelastic/DATA/Par_file_attenuation_2D_source_in_point_shared_by_four_elements @@ -30,7 +30,7 @@ DT = 4e-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/axisymmetric_case_AXISYM_option/DATA/Par_file b/EXAMPLES/axisymmetric_case_AXISYM_option/DATA/Par_file index cf4274bad..5324f1b37 100644 --- a/EXAMPLES/axisymmetric_case_AXISYM_option/DATA/Par_file +++ b/EXAMPLES/axisymmetric_case_AXISYM_option/DATA/Par_file @@ -30,7 +30,7 @@ DT = 0.60d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -296,6 +296,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/axisymmetric_case_AXISYM_option/REF_SEIS/Par_file b/EXAMPLES/axisymmetric_case_AXISYM_option/REF_SEIS/Par_file index cf4274bad..5324f1b37 100644 --- a/EXAMPLES/axisymmetric_case_AXISYM_option/REF_SEIS/Par_file +++ b/EXAMPLES/axisymmetric_case_AXISYM_option/REF_SEIS/Par_file @@ -30,7 +30,7 @@ DT = 0.60d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -296,6 +296,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/canyon/DATA/Par_file b/EXAMPLES/canyon/DATA/Par_file index 2f9dda128..7c532ba24 100644 --- a/EXAMPLES/canyon/DATA/Par_file +++ b/EXAMPLES/canyon/DATA/Par_file @@ -30,7 +30,7 @@ DT = 2.e-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_elastic/DATA/Par_file_no_attenuation_2D_at_the_corner_between_several_spectral_elements b/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_elastic/DATA/Par_file_no_attenuation_2D_at_the_corner_between_several_spectral_elements index f7b9ffabf..b89987a1e 100644 --- a/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_elastic/DATA/Par_file_no_attenuation_2D_at_the_corner_between_several_spectral_elements +++ b/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_elastic/DATA/Par_file_no_attenuation_2D_at_the_corner_between_several_spectral_elements @@ -30,7 +30,7 @@ DT = 4e-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_elastic/DATA/Par_file_no_attenuation_2D_inside_a_given_spectral_element b/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_elastic/DATA/Par_file_no_attenuation_2D_inside_a_given_spectral_element index 121a74eaf..9e51b92d4 100644 --- a/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_elastic/DATA/Par_file_no_attenuation_2D_inside_a_given_spectral_element +++ b/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_elastic/DATA/Par_file_no_attenuation_2D_inside_a_given_spectral_element @@ -30,7 +30,7 @@ DT = 4e-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_elastic/REF_SEIS/Par_file b/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_elastic/REF_SEIS/Par_file index f7b9ffabf..b89987a1e 100644 --- a/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_elastic/REF_SEIS/Par_file +++ b/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_elastic/REF_SEIS/Par_file @@ -30,7 +30,7 @@ DT = 4e-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_viscoelastic/DATA/Par_file_no_attenuation_2D_at_the_corner_between_several_spectral_elements b/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_viscoelastic/DATA/Par_file_no_attenuation_2D_at_the_corner_between_several_spectral_elements index f47d65119..e16e76a11 100644 --- a/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_viscoelastic/DATA/Par_file_no_attenuation_2D_at_the_corner_between_several_spectral_elements +++ b/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_viscoelastic/DATA/Par_file_no_attenuation_2D_at_the_corner_between_several_spectral_elements @@ -30,7 +30,7 @@ DT = 4e-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_viscoelastic/REF_SEIS/Par_file b/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_viscoelastic/REF_SEIS/Par_file index f47d65119..e16e76a11 100644 --- a/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_viscoelastic/REF_SEIS/Par_file +++ b/EXAMPLES/check_absolute_amplitude_of_force_source_seismograms_viscoelastic/REF_SEIS/Par_file @@ -30,7 +30,7 @@ DT = 4e-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/check_absolute_amplitude_of_pressure_source_seismograms_acoustic/DATA/Par_file_no_attenuation_2D_at_the_corner_between_several_spectral_elements b/EXAMPLES/check_absolute_amplitude_of_pressure_source_seismograms_acoustic/DATA/Par_file_no_attenuation_2D_at_the_corner_between_several_spectral_elements index 66d2b13d9..8d32d65f8 100644 --- a/EXAMPLES/check_absolute_amplitude_of_pressure_source_seismograms_acoustic/DATA/Par_file_no_attenuation_2D_at_the_corner_between_several_spectral_elements +++ b/EXAMPLES/check_absolute_amplitude_of_pressure_source_seismograms_acoustic/DATA/Par_file_no_attenuation_2D_at_the_corner_between_several_spectral_elements @@ -30,7 +30,7 @@ DT = 4e-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/check_absolute_amplitude_of_pressure_source_seismograms_acoustic/DATA/Par_file_no_attenuation_2D_inside_a_given_spectral_element b/EXAMPLES/check_absolute_amplitude_of_pressure_source_seismograms_acoustic/DATA/Par_file_no_attenuation_2D_inside_a_given_spectral_element index c201ad286..0ab68b068 100644 --- a/EXAMPLES/check_absolute_amplitude_of_pressure_source_seismograms_acoustic/DATA/Par_file_no_attenuation_2D_inside_a_given_spectral_element +++ b/EXAMPLES/check_absolute_amplitude_of_pressure_source_seismograms_acoustic/DATA/Par_file_no_attenuation_2D_inside_a_given_spectral_element @@ -30,7 +30,7 @@ DT = 4e-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/check_absolute_amplitude_of_pressure_source_seismograms_acoustic/REF_SEIS/Par_file b/EXAMPLES/check_absolute_amplitude_of_pressure_source_seismograms_acoustic/REF_SEIS/Par_file index 66d2b13d9..8d32d65f8 100644 --- a/EXAMPLES/check_absolute_amplitude_of_pressure_source_seismograms_acoustic/REF_SEIS/Par_file +++ b/EXAMPLES/check_absolute_amplitude_of_pressure_source_seismograms_acoustic/REF_SEIS/Par_file @@ -30,7 +30,7 @@ DT = 4e-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/electromagnetic_GPR_TM/DATA/Par_file b/EXAMPLES/electromagnetic_GPR_TM/DATA/Par_file index e8bb40a0b..310e70d08 100644 --- a/EXAMPLES/electromagnetic_GPR_TM/DATA/Par_file +++ b/EXAMPLES/electromagnetic_GPR_TM/DATA/Par_file @@ -30,7 +30,7 @@ DT = 4d-12 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .false. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/fluid_solid/fluid_solid_external_mesh/DATA/Par_file b/EXAMPLES/fluid_solid/fluid_solid_external_mesh/DATA/Par_file index f7d90509c..b91f98ec3 100644 --- a/EXAMPLES/fluid_solid/fluid_solid_external_mesh/DATA/Par_file +++ b/EXAMPLES/fluid_solid/fluid_solid_external_mesh/DATA/Par_file @@ -30,7 +30,7 @@ DT = 0.25d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/DATA/Par_file b/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/DATA/Par_file index dd3cb7e6e..9078539e5 100644 --- a/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/DATA/Par_file +++ b/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/DATA/Par_file @@ -30,7 +30,7 @@ DT = 0.85d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/DATA/Par_file b/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/DATA/Par_file index 3eb709250..707992d14 100644 --- a/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/DATA/Par_file +++ b/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/DATA/Par_file @@ -30,7 +30,7 @@ DT = 0.7d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/global_Earth_ak135f/DATA/Par_file b/EXAMPLES/global_Earth_ak135f/DATA/Par_file index 59c8ce87f..c66082d0f 100644 --- a/EXAMPLES/global_Earth_ak135f/DATA/Par_file +++ b/EXAMPLES/global_Earth_ak135f/DATA/Par_file @@ -30,7 +30,7 @@ DT = 0.25 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 4 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/infinite_homogeneous_moment_tensor_vertical_dip_slip/DATA/Par_file b/EXAMPLES/infinite_homogeneous_moment_tensor_vertical_dip_slip/DATA/Par_file index a9ef1c29f..865e4aa72 100644 --- a/EXAMPLES/infinite_homogeneous_moment_tensor_vertical_dip_slip/DATA/Par_file +++ b/EXAMPLES/infinite_homogeneous_moment_tensor_vertical_dip_slip/DATA/Par_file @@ -30,7 +30,7 @@ DT = 0.8e-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/infinite_homogeneous_plane_wave/DATA/Par_file b/EXAMPLES/infinite_homogeneous_plane_wave/DATA/Par_file index 9046c7f2e..7202daa8b 100644 --- a/EXAMPLES/infinite_homogeneous_plane_wave/DATA/Par_file +++ b/EXAMPLES/infinite_homogeneous_plane_wave/DATA/Par_file @@ -30,7 +30,7 @@ DT = 20.d-8 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/initial_mode_LDDRK/DATA/Par_file b/EXAMPLES/initial_mode_LDDRK/DATA/Par_file index f1bf3ff1e..28cdaf4fb 100644 --- a/EXAMPLES/initial_mode_LDDRK/DATA/Par_file +++ b/EXAMPLES/initial_mode_LDDRK/DATA/Par_file @@ -30,7 +30,7 @@ DT = 0.4e-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 2 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/initial_plane_wave_with_free_surface/DATA/Par_file b/EXAMPLES/initial_plane_wave_with_free_surface/DATA/Par_file index b9cfd9b12..fd10df566 100644 --- a/EXAMPLES/initial_plane_wave_with_free_surface/DATA/Par_file +++ b/EXAMPLES/initial_plane_wave_with_free_surface/DATA/Par_file @@ -30,7 +30,7 @@ DT = 0.01 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/moving_sources_acoustic/DATA/Par_file b/EXAMPLES/moving_sources_acoustic/DATA/Par_file index dab2bcac5..98f968c5e 100644 --- a/EXAMPLES/moving_sources_acoustic/DATA/Par_file +++ b/EXAMPLES/moving_sources_acoustic/DATA/Par_file @@ -30,7 +30,7 @@ DT = 2.5d-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/noise_uniform/DATA/Par_file b/EXAMPLES/noise_uniform/DATA/Par_file index fd7727895..2c7855a0e 100644 --- a/EXAMPLES/noise_uniform/DATA/Par_file +++ b/EXAMPLES/noise_uniform/DATA/Par_file @@ -30,7 +30,7 @@ DT = 5.d-2 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .false. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/noise_uniform/REF_SEIS/Par_file b/EXAMPLES/noise_uniform/REF_SEIS/Par_file index fd7727895..c36596f94 100644 --- a/EXAMPLES/noise_uniform/REF_SEIS/Par_file +++ b/EXAMPLES/noise_uniform/REF_SEIS/Par_file @@ -30,7 +30,7 @@ DT = 5.d-2 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .false. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/paper_axisymmetry_example/DATA/Par_file b/EXAMPLES/paper_axisymmetry_example/DATA/Par_file index a2f8a1156..688f5017e 100644 --- a/EXAMPLES/paper_axisymmetry_example/DATA/Par_file +++ b/EXAMPLES/paper_axisymmetry_example/DATA/Par_file @@ -30,7 +30,7 @@ DT = 1.220703125d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/Pekeris_waveguide/DATA/Par_file b/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/Pekeris_waveguide/DATA/Par_file index 54e5e4c26..085a15750 100644 --- a/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/Pekeris_waveguide/DATA/Par_file +++ b/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/Pekeris_waveguide/DATA/Par_file @@ -30,7 +30,7 @@ DT = 1.6d-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/with_free_bottom_surface/DATA/Par_file b/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/with_free_bottom_surface/DATA/Par_file index 987b8e11d..7b072b760 100644 --- a/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/with_free_bottom_surface/DATA/Par_file +++ b/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/with_free_bottom_surface/DATA/Par_file @@ -30,7 +30,7 @@ DT = 2.0d-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/with_rigid_bottom_surface/DATA/Par_file b/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/with_rigid_bottom_surface/DATA/Par_file index 874cfadbe..49d00f9cc 100644 --- a/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/with_rigid_bottom_surface/DATA/Par_file +++ b/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/with_rigid_bottom_surface/DATA/Par_file @@ -30,7 +30,7 @@ DT = 2.0d-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/poroelastic_acoustic/DATA/Par_file b/EXAMPLES/poroelastic_acoustic/DATA/Par_file index d050d5217..2b77bf527 100644 --- a/EXAMPLES/poroelastic_acoustic/DATA/Par_file +++ b/EXAMPLES/poroelastic_acoustic/DATA/Par_file @@ -30,7 +30,7 @@ DT = 4.2d-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -288,6 +288,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/poroelastic_semi_infinite_homogeneous/DATA/Par_file b/EXAMPLES/poroelastic_semi_infinite_homogeneous/DATA/Par_file index 1897341f3..2bc1fd2fd 100644 --- a/EXAMPLES/poroelastic_semi_infinite_homogeneous/DATA/Par_file +++ b/EXAMPLES/poroelastic_semi_infinite_homogeneous/DATA/Par_file @@ -30,7 +30,7 @@ DT = 1.d-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_homogeneous/DATA/Par_file b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_homogeneous/DATA/Par_file index b9944c23d..fc3df43af 100644 --- a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_homogeneous/DATA/Par_file +++ b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_homogeneous/DATA/Par_file @@ -30,7 +30,7 @@ DT = 4.25e-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_normal_fluid_solid/DATA/Par_file b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_normal_fluid_solid/DATA/Par_file index 109ca69dd..e1eee148e 100644 --- a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_normal_fluid_solid/DATA/Par_file +++ b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_normal_fluid_solid/DATA/Par_file @@ -30,7 +30,7 @@ DT = 4.25e-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_normal_solid_only/DATA/Par_file b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_normal_solid_only/DATA/Par_file index d42c49c6b..706dfbf36 100644 --- a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_normal_solid_only/DATA/Par_file +++ b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_normal_solid_only/DATA/Par_file @@ -30,7 +30,7 @@ DT = 4.25e-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/Stacey_homogeneous/DATA/Par_file b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/Stacey_homogeneous/DATA/Par_file index 83f4bc422..ba135589b 100644 --- a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/Stacey_homogeneous/DATA/Par_file +++ b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/Stacey_homogeneous/DATA/Par_file @@ -30,7 +30,7 @@ DT = 4.25e-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/Stacey_normal_fluid_solid/DATA/Par_file b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/Stacey_normal_fluid_solid/DATA/Par_file index 94b5c7166..ec167848d 100644 --- a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/Stacey_normal_fluid_solid/DATA/Par_file +++ b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/Stacey_normal_fluid_solid/DATA/Par_file @@ -30,7 +30,7 @@ DT = 4.25e-4 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file b/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file index 21ed46ea1..d334a9f28 100644 --- a/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file +++ b/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file @@ -30,7 +30,7 @@ DT = 0.25d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file_adj b/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file_adj index c1a4729a9..c46602d39 100644 --- a/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file_adj +++ b/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file_adj @@ -30,7 +30,7 @@ DT = 0.25d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file_fwd b/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file_fwd index 296453aa8..38aaa1f7b 100644 --- a/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file_fwd +++ b/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file_fwd @@ -30,7 +30,7 @@ DT = 0.25d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file_kernel b/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file_kernel index cf4e0eefa..14ae4731e 100644 --- a/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file_kernel +++ b/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file_kernel @@ -30,7 +30,7 @@ DT = 0.25d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file_save_forward b/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file_save_forward index 8ff579f9d..d43196aaa 100644 --- a/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file_save_forward +++ b/EXAMPLES/salt_dome_Vadim/specfem2d_run_input_files_and_scripts/DATA/Par_file_save_forward @@ -30,7 +30,7 @@ DT = 0.25d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/semi_infinite_homogeneous/DATA/Par_file b/EXAMPLES/semi_infinite_homogeneous/DATA/Par_file index 5620837c0..7cf086989 100644 --- a/EXAMPLES/semi_infinite_homogeneous/DATA/Par_file +++ b/EXAMPLES/semi_infinite_homogeneous/DATA/Par_file @@ -30,7 +30,7 @@ DT = 0.8e-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/simple_topography_and_also_a_simple_fluid_layer/DATA/Par_file b/EXAMPLES/simple_topography_and_also_a_simple_fluid_layer/DATA/Par_file index 9f9f329a9..94c7ddfe3 100644 --- a/EXAMPLES/simple_topography_and_also_a_simple_fluid_layer/DATA/Par_file +++ b/EXAMPLES/simple_topography_and_also_a_simple_fluid_layer/DATA/Par_file @@ -30,7 +30,7 @@ DT = 1.1d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -288,6 +288,7 @@ nbmodels = 4 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/simple_topography_and_also_a_simple_fluid_layer/REF_KERNEL/Par_file b/EXAMPLES/simple_topography_and_also_a_simple_fluid_layer/REF_KERNEL/Par_file index c18e824aa..1c496fa6d 100644 --- a/EXAMPLES/simple_topography_and_also_a_simple_fluid_layer/REF_KERNEL/Par_file +++ b/EXAMPLES/simple_topography_and_also_a_simple_fluid_layer/REF_KERNEL/Par_file @@ -30,7 +30,7 @@ DT = 1.1d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/simple_topography_and_also_a_simple_fluid_layer/REF_SEIS/Par_file b/EXAMPLES/simple_topography_and_also_a_simple_fluid_layer/REF_SEIS/Par_file index 9f9f329a9..94c7ddfe3 100644 --- a/EXAMPLES/simple_topography_and_also_a_simple_fluid_layer/REF_SEIS/Par_file +++ b/EXAMPLES/simple_topography_and_also_a_simple_fluid_layer/REF_SEIS/Par_file @@ -30,7 +30,7 @@ DT = 1.1d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -288,6 +288,7 @@ nbmodels = 4 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/thermocline/DATA/Par_file b/EXAMPLES/thermocline/DATA/Par_file index b9a215e5b..de1ca68b5 100644 --- a/EXAMPLES/thermocline/DATA/Par_file +++ b/EXAMPLES/thermocline/DATA/Par_file @@ -30,7 +30,7 @@ DT = 1.d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/tomographic_ocean_model/DATA/Par_file b/EXAMPLES/tomographic_ocean_model/DATA/Par_file index 6cbd57fbd..eae6ad36f 100644 --- a/EXAMPLES/tomographic_ocean_model/DATA/Par_file +++ b/EXAMPLES/tomographic_ocean_model/DATA/Par_file @@ -30,7 +30,7 @@ DT = 1.5d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/tomographic_ocean_model/REF_SEIS/Par_file b/EXAMPLES/tomographic_ocean_model/REF_SEIS/Par_file index 6cbd57fbd..eae6ad36f 100644 --- a/EXAMPLES/tomographic_ocean_model/REF_SEIS/Par_file +++ b/EXAMPLES/tomographic_ocean_model/REF_SEIS/Par_file @@ -30,7 +30,7 @@ DT = 1.5d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation diff --git a/EXAMPLES/unresolved_bugs/bugPMLs_1_extra_phase_coming_back_from_the_left_PML/Par_file_fluid_solid b/EXAMPLES/unresolved_bugs/bugPMLs_1_extra_phase_coming_back_from_the_left_PML/Par_file_fluid_solid index 9671ec924..9d480ddc6 100644 --- a/EXAMPLES/unresolved_bugs/bugPMLs_1_extra_phase_coming_back_from_the_left_PML/Par_file_fluid_solid +++ b/EXAMPLES/unresolved_bugs/bugPMLs_1_extra_phase_coming_back_from_the_left_PML/Par_file_fluid_solid @@ -30,7 +30,7 @@ DT = 0.875d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/unresolved_bugs/bugPMLs_1_extra_phase_coming_back_from_the_left_PML/the_same_bug_appears_from_the_right_PML/Par_file_fluid_solid b/EXAMPLES/unresolved_bugs/bugPMLs_1_extra_phase_coming_back_from_the_left_PML/the_same_bug_appears_from_the_right_PML/Par_file_fluid_solid index 9671ec924..9d480ddc6 100644 --- a/EXAMPLES/unresolved_bugs/bugPMLs_1_extra_phase_coming_back_from_the_left_PML/the_same_bug_appears_from_the_right_PML/Par_file_fluid_solid +++ b/EXAMPLES/unresolved_bugs/bugPMLs_1_extra_phase_coming_back_from_the_left_PML/the_same_bug_appears_from_the_right_PML/Par_file_fluid_solid @@ -30,7 +30,7 @@ DT = 0.875d-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/unresolved_bugs/bugPMLs_2_weak_instability_in_long_times/DATA/Par_file b/EXAMPLES/unresolved_bugs/bugPMLs_2_weak_instability_in_long_times/DATA/Par_file index 536bb4b6e..e9e09ed7c 100644 --- a/EXAMPLES/unresolved_bugs/bugPMLs_2_weak_instability_in_long_times/DATA/Par_file +++ b/EXAMPLES/unresolved_bugs/bugPMLs_2_weak_instability_in_long_times/DATA/Par_file @@ -30,7 +30,7 @@ DT = 1.6e-3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta time_stepping_scheme = 1 -# set the type of calculation (P-SV or SH/membrane waves) +# set the type of calculation (P-SV (EM TE) or SH/membrane waves (EM TM)) P_SV = .true. # axisymmetric (2.5D) or Cartesian planar (2D) simulation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 3 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/utils/CPML/Par_file_to_use_with_this_timing_routine_HOWTO b/utils/CPML/Par_file_to_use_with_this_timing_routine_HOWTO index 76e5201e0..4d2af207b 100644 --- a/utils/CPML/Par_file_to_use_with_this_timing_routine_HOWTO +++ b/utils/CPML/Par_file_to_use_with_this_timing_routine_HOWTO @@ -1,8 +1,8 @@ -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # -# simulation input parameters +# Simulation input parameters # -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # title of job title = Test of SPECFEM2D with curved interfaces @@ -17,7 +17,7 @@ NOISE_TOMOGRAPHY = 0 SAVE_FORWARD = .false. # parameters concerning partitioning -NPROC = 1 +NPROC = 1 # number of processes # time step parameters # total number of time steps @@ -72,20 +72,20 @@ SAVE_MODEL = default # (used for tomography materials with negative material ids and/or MODEL==tomo settings) TOMOGRAPHY_FILE = ./DATA/tomo_file.xyz -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # # Attenuation # -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # attenuation parameters ATTENUATION_VISCOELASTIC = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model ATTENUATION_VISCOACOUSTIC = .false. # turn attenuation (viscoacousticity) on or off for non-poroelastic fluid parts of the model -# for viscoelastic attenuation +# for viscoelastic or viscoacoustic attenuation N_SLS = 3 # number of standard linear solids for attenuation (3 is usually the minimum) -ATTENUATION_f0_REFERENCE = 5.196 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file -READ_VELOCITIES_AT_f0 = .false. # shift velocities to account for physical dispersion (see user manual for more information) +ATTENUATION_f0_REFERENCE = 5.196 # in case of attenuation, reference frequency in Hz at which the velocity values in the velocity model are given (unused otherwise); relevant only if source is a Dirac or a Heaviside, otherwise it is automatically set to f0 the dominant frequency of the source in the DATA/SOURCE file +READ_VELOCITIES_AT_f0 = .false. # read seismic velocities at ATTENUATION_f0_REFERENCE instead of at infinite frequency (see user manual for more information) USE_SOLVOPT = .false. # use more precise but much more expensive way of determining the Q factor relaxation times, as in https://doi.org/10.1093/gji/ggw024 # for poroelastic attenuation @@ -96,10 +96,10 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file -# to undo attenuation for sensitivity kernel calculations or forward runs with SAVE_FORWARD -# use the flag below. It performs undoing of attenuation in an exact way for sensitivity kernel calculations +# to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD +# use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations # but requires disk space for temporary storage, and uses a significant amount of memory used as buffers for temporary storage. # When that option is on the second parameter indicates how often the code dumps restart files to disk (if in doubt, use something between 100 and 1000). UNDO_ATTENUATION_AND_OR_PML = .false. @@ -109,11 +109,11 @@ NT_DUMP_ATTENUATION = 500 # Outperforms conventional mode using a value of NTSTEP_BETWEEN_COMPUTE_KERNELS high enough. NO_BACKWARD_RECONSTRUCTION = .false. -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # # Sources # -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # source parameters NSOURCES = 1 # number of sources (source information is then read from the DATA/SOURCE file) @@ -142,14 +142,14 @@ noise_source_time_function_type = 4 # This is only useful for GPU version (for now) write_moving_sources_database = .false. -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # # Receivers # -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # receiver set parameters for recording stations (i.e. recording points) -# record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential 7=electric field +# seismotype : record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential 7=electric field seismotype = 1 # several values can be chosen. For example : 1,2,4 # interval in time steps for writing of seismograms @@ -193,11 +193,11 @@ zfin = 2200. # last receiver z in meters (ig record_at_surface_same_vertical = .true. # receivers inside the medium or at the surface (z values are ignored if this is set to true, they are replaced with the topography height) -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # -# Adjoint kernel outputs +# adjoint kernel outputs # -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # save sensitivity kernels in ASCII format (much bigger files, but compatible with current GMT scripts) or in binary format save_ASCII_kernels = .true. @@ -208,11 +208,11 @@ NTSTEP_BETWEEN_COMPUTE_KERNELS = 1 # outputs approximate Hessian for preconditioning APPROXIMATE_HESS_KL = .false. -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # # Boundary conditions # -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # Perfectly Matched Layer (PML) boundaries # absorbing boundary active or not @@ -237,11 +237,11 @@ STACEY_ABSORBING_CONDITIONS = .false. ADD_PERIODIC_CONDITIONS = .false. PERIODIC_HORIZ_DIST = 4000.d0 -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # # MESHING - Velocity and density models # -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # use an external mesh created by an external meshing tool or use the internal mesher read_external_mesh = .false. @@ -270,6 +270,7 @@ tangential_detection_curve_file = dummy # file # PARAMETERS FOR INTERNAL MESHING # #----------------------------------------------------------- + # material properties # number of model materials nbmodels = 1 @@ -279,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, @@ -309,11 +311,11 @@ nbregions = 1 # then set below the different # format of each line: nxmin nxmax nzmin nzmax material_number 1 80 1 60 1 -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # # Display parameters # -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # interval at which we output time step info and max of norm of displacement # (every how many time steps we display information about the simulation. costly, do not use a very small value) @@ -334,11 +336,11 @@ NTSTEP_BETWEEN_OUTPUT_ENERGY = 10 # the script utils/visualisation/plotIntegratedEnergyFile.py to watch. It is refreshed at the same time than the seismograms COMPUTE_INTEGRATED_ENERGY_FIELD = .false. -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # # Movies/images/snaphots visualizations # -#----------------------------------------------------------------------------- +#----------------------------------------------------------- # every how many time steps we draw JPEG or PostScript pictures of the simulation # and/or we dump results of the simulation as ASCII or binary files (costly, do not use a very small value) @@ -391,9 +393,9 @@ use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision # When that option is on, of course the number of processor cores used to start # the code in the batch system must be a multiple of NUMBER_OF_SIMULTANEOUS_RUNS, # all the individual runs must use the same number of processor cores, -# which as usual is proc in the Par_file, +# which as usual is NPROC in the Par_file, # and thus the total number of processor cores to request from the batch system -# should be NUMBER_OF_SIMULTANEOUS_RUNS * proc. +# should be NUMBER_OF_SIMULTANEOUS_RUNS * NPROC. # All the runs to perform must be placed in directories called run0001, run0002, run0003 and so on # (with exactly four digits). # From 0e640754237650c135acd783e7eaa1681b888eec Mon Sep 17 00:00:00 2001 From: Daniel Peter Date: Wed, 5 Jun 2024 11:00:34 +0200 Subject: [PATCH 6/9] code cleaning (whitespace) --- src/shared/read_material_table.f90 | 26 +++++------ src/shared/read_parameter_file.F90 | 4 +- src/shared/shared_par.F90 | 2 +- src/specfem2D/check_grid.F90 | 4 +- .../compute_add_sources_electromagnetic.f90 | 6 +-- .../compute_attenuation_electromagnetic.f90 | 2 +- .../compute_forces_electromagnetic.F90 | 8 ++-- ...forces_electromagnetic_calling_routine.F90 | 4 +- .../compute_stacey_electromagnetic.f90 | 46 +++++++++---------- src/specfem2D/invert_mass_matrix.F90 | 28 +++++------ src/specfem2D/iterate_time.F90 | 2 +- src/specfem2D/prepare_optimized_arrays.F90 | 2 +- src/specfem2D/prepare_wavefields.f90 | 4 +- src/specfem2D/read_materials.f90 | 18 ++++---- src/specfem2D/read_mesh_databases.F90 | 2 +- src/specfem2D/save_model_files.f90 | 2 +- src/specfem2D/setup_mesh.F90 | 18 ++++---- src/specfem2D/specfem2D_par.f90 | 16 +++---- 18 files changed, 97 insertions(+), 97 deletions(-) diff --git a/src/shared/read_material_table.f90 b/src/shared/read_material_table.f90 index 78a2d4114..8ed8f7f9e 100644 --- a/src/shared/read_material_table.f90 +++ b/src/shared/read_material_table.f90 @@ -45,7 +45,7 @@ subroutine read_material_table() rho_s_read,rho_f_read, & phi_read,tortuosity_read, & permxx_read,permxz_read,permzz_read,kappa_s_read,kappa_f_read,kappa_fr_read, & - eta_f_read,mu_fr_read,mu0_read,e0_read,e11_read,e33_read,sig11_read,sig33_read,& + eta_f_read,mu_fr_read,mu0_read,e0_read,e11_read,e33_read,sig11_read,sig33_read, & Qe11_read,Qe33_read,Qs11_read,Qs33_read, & compaction_grad @@ -167,7 +167,7 @@ subroutine read_material_table() call stop_the_code('Invalid poroelastic material: non-positive value of Qmu') else if (icodemat(i) == ELECTROMAGNETIC_MATERIAL) then - ! electromagnetic material + ! electromagnetic material mu0_read(i) = val0read e0_read(i) = val1read e11_read(i) = val2read @@ -248,7 +248,7 @@ subroutine initialize_material_properties() eta_f_read,mu_fr_read ! electromagnetic - use shared_parameters, only: mu0_read,e0_read,e11_read,e33_read,sig11_read,sig33_read,& + use shared_parameters, only: mu0_read,e0_read,e11_read,e33_read,sig11_read,sig33_read, & Qe11_read,Qe33_read,Qs11_read,Qs33_read implicit none @@ -289,15 +289,15 @@ subroutine initialize_material_properties() eta_f_read(nbmodels), & mu_fr_read(nbmodels)) - allocate( mu0_read(nbmodels),& - e0_read(nbmodels),& - e11_read(nbmodels),& - e33_read(nbmodels),& - sig11_read(nbmodels),& - sig33_read(nbmodels),& - Qe11_read(nbmodels),& - Qe33_read(nbmodels),& - Qs11_read(nbmodels),& + allocate( mu0_read(nbmodels), & + e0_read(nbmodels), & + e11_read(nbmodels), & + e33_read(nbmodels), & + sig11_read(nbmodels), & + sig33_read(nbmodels), & + Qe11_read(nbmodels), & + Qe33_read(nbmodels), & + Qs11_read(nbmodels), & Qs33_read(nbmodels)) allocate(compaction_grad(nbmodels)) @@ -379,7 +379,7 @@ subroutine print_materials_info() eta_f_read,mu_fr_read ! electromagnetic - use shared_parameters, only: mu0_read,e0_read,e11_read,e33_read,sig11_read,sig33_read,& + use shared_parameters, only: mu0_read,e0_read,e11_read,e33_read,sig11_read,sig33_read, & Qe11_read,Qe33_read,Qs11_read,Qs33_read implicit none diff --git a/src/shared/read_parameter_file.F90 b/src/shared/read_parameter_file.F90 index e592c066a..3ededff71 100644 --- a/src/shared/read_parameter_file.F90 +++ b/src/shared/read_parameter_file.F90 @@ -564,14 +564,14 @@ subroutine read_parameter_file_only() write(*,'(a)') 'ATTENUATION_PERMITTIVITY = .false.' write(*,*) endif - + call read_value_logical_p(ATTENUATION_CONDUCTIVITY, 'ATTENUATION_CONDUCTIVITY') if (err_occurred() /= 0) then some_parameters_missing_from_Par_file = .true. write(*,'(a)') 'ATTENUATION_CONDUCTIVITY = .false.' write(*,*) endif - + call read_value_double_precision_p(f0_electromagnetic, 'f0_electromagnetic') if (err_occurred() /= 0) then some_parameters_missing_from_Par_file = .true. diff --git a/src/shared/shared_par.F90 b/src/shared/shared_par.F90 index 5deee0578..508eb6591 100644 --- a/src/shared/shared_par.F90 +++ b/src/shared/shared_par.F90 @@ -416,7 +416,7 @@ module shared_parameters permzz_read,kappa_s_read,kappa_f_read,kappa_fr_read,eta_f_read,mu_fr_read ! electromagnetic (EM) - double precision, dimension(:),allocatable :: mu0_read,e0_read,e11_read,e33_read,sig11_read,sig33_read,& + double precision, dimension(:),allocatable :: mu0_read,e0_read,e11_read,e33_read,sig11_read,sig33_read, & Qe11_read,Qe33_read,Qs11_read,Qs33_read ! compaction gradient diff --git a/src/specfem2D/check_grid.F90 b/src/specfem2D/check_grid.F90 index 8695f6b79..982c103fe 100644 --- a/src/specfem2D/check_grid.F90 +++ b/src/specfem2D/check_grid.F90 @@ -218,7 +218,7 @@ subroutine check_grid() cpIIloc = 0.d0 endif cEloc = 0.d0 - + else cpIloc = 0.d0 cpIIloc = 0.d0 @@ -241,7 +241,7 @@ subroutine check_grid() ! ignore fluid regions with Vs = 0 if (csloc > TINYVAL) vsmin = min(vsmin,csloc) vsmax = max(vsmax,csloc) - + if (cEloc > TINYVAL) vEmin = min(vEmin,cEloc) vEmax = max(vEmax,cEloc) diff --git a/src/specfem2D/compute_add_sources_electromagnetic.f90 b/src/specfem2D/compute_add_sources_electromagnetic.f90 index 2812b3e62..f869e4032 100644 --- a/src/specfem2D/compute_add_sources_electromagnetic.f90 +++ b/src/specfem2D/compute_add_sources_electromagnetic.f90 @@ -38,7 +38,7 @@ subroutine compute_add_sources_electromagnetic(accel_electromagnetic,it,i_stage) use constants, only: CUSTOM_REAL,NGLLX,NGLLZ,NDIM,HALF,myrank use specfem_par, only: P_SV,ispec_is_electromagnetic,nglob_electromagnetic, & - NSOURCES,source_time_function,& + NSOURCES,source_time_function, & islice_selected_source,ispec_selected_source,sourcearrays, & ibool implicit none @@ -58,7 +58,7 @@ subroutine compute_add_sources_electromagnetic(accel_electromagnetic,it,i_stage) ! element containing source ispec = ispec_selected_source(i_source) - + ! source element is electromagnetic if (ispec_is_electromagnetic(ispec)) then ! source time function @@ -83,7 +83,7 @@ subroutine compute_add_sources_electromagnetic(accel_electromagnetic,it,i_stage) ! 2D: y-component uses array(1,..) accel_electromagnetic(1,iglob) = accel_electromagnetic(1,iglob) - sourcearrays(1,i,j,i_source) * stf_used enddo - enddo + enddo endif endif ! source element is electromagnetic diff --git a/src/specfem2D/compute_attenuation_electromagnetic.f90 b/src/specfem2D/compute_attenuation_electromagnetic.f90 index 97dde245d..c6c527749 100644 --- a/src/specfem2D/compute_attenuation_electromagnetic.f90 +++ b/src/specfem2D/compute_attenuation_electromagnetic.f90 @@ -66,7 +66,7 @@ subroutine compute_attenuation_electromagnetic() permlxx = spermittivitystore(1,i,j,ispec) permlzz = spermittivitystore(2,i,j,ispec) - permx_loc(i,j) = veloc_electromagnetic(1,iglob) * permlxx + permx_loc(i,j) = veloc_electromagnetic(1,iglob) * permlxx permz_loc(i,j) = veloc_electromagnetic(2,iglob) * permlzz ! time stepping diff --git a/src/specfem2D/compute_forces_electromagnetic.F90 b/src/specfem2D/compute_forces_electromagnetic.F90 index d0627cba5..9782a8610 100644 --- a/src/specfem2D/compute_forces_electromagnetic.F90 +++ b/src/specfem2D/compute_forces_electromagnetic.F90 @@ -31,14 +31,14 @@ !======================================================================== - subroutine compute_forces_electromagnetic(accel_electromagnetic,displ_electromagnetic,& + subroutine compute_forces_electromagnetic(accel_electromagnetic,displ_electromagnetic, & rx_permattenuation,rz_permattenuation,iphase) ! compute forces for the electromagnetic elements use constants, only: CUSTOM_REAL,NGLLX,NGLLZ,NGLJ,NDIM, & ONE,HALF,PI,TINYVAL,FOUR_THIRDS - use specfem_par, only: nglob,& + use specfem_par, only: nglob, & ibool,kmato,ispec_is_electromagnetic, & xix,xiz,gammax,gammaz,jacobian, & hprime_xx,hprimewgll_xx,hprime_zz,hprimewgll_zz,wxgll,wzgll, & @@ -177,7 +177,7 @@ subroutine compute_forces_electromagnetic(accel_electromagnetic,displ_electromag !--- if external medium, get electromagnetic parameters of current grid point inv_mag_permeability = inv_magpermeabilitystore(i,j,ispec) - if(P_SV) then + if (P_SV) then sigma_xx = inv_mag_permeability*(dux_dxl(i,j)+duz_dzl(i,j)) ! sigma_zz = inv_mag_permeability*(dux_dxl(i,j)+duz_dzl(i,j)) ! sigma_xz = inv_mag_permeability*(duz_dxl(i,j) - dux_dzl(i,j)) ! d_t \tilde{H}_xz = d_t H_y = mu^-1(dxEz-dzEx) @@ -255,7 +255,7 @@ subroutine compute_forces_electromagnetic(accel_electromagnetic,displ_electromag accel_electromagnetic(1,iglob) = accel_electromagnetic(1,iglob) - (wzgll(j) * tempx1l + wxgll(i) * tempx2l) & - wxgll(i) * wzgll(j) * jacobian(i,j,ispec) * rx_permattenuation(i,j,ispec) accel_electromagnetic(2,iglob) = 0.d0 - endif + endif else if (P_SV) then accel_electromagnetic(1,iglob) = accel_electromagnetic(1,iglob) - (wzgll(j) * tempx1l + wxgll(i) * tempx2l) diff --git a/src/specfem2D/compute_forces_electromagnetic_calling_routine.F90 b/src/specfem2D/compute_forces_electromagnetic_calling_routine.F90 index 99353b2fa..0df6588b8 100644 --- a/src/specfem2D/compute_forces_electromagnetic_calling_routine.F90 +++ b/src/specfem2D/compute_forces_electromagnetic_calling_routine.F90 @@ -55,7 +55,7 @@ subroutine compute_forces_electromagnetic_main() ! main solver for the electromagnetic elements ! electromagnetic term - call compute_forces_electromagnetic(accel_electromagnetic,displ_electromagnetic,& + call compute_forces_electromagnetic(accel_electromagnetic,displ_electromagnetic, & rx_permattenuation,rz_permattenuation,iphase) ! PML_BOUNDARY_CONDITIONS,e1,e11,e13,iphase) @@ -73,7 +73,7 @@ subroutine compute_forces_electromagnetic_main() ! ! PML boundary ! if (PML_BOUNDARY_CONDITIONS) then -! call pml_boundary_electromagnetic(accel_electromagnetic,veloc_electromagnetic,displ_electromagnetic,& +! call pml_boundary_electromagnetic(accel_electromagnetic,veloc_electromagnetic,displ_electromagnetic, & ! displ_electromagnetic_old) ! endif diff --git a/src/specfem2D/compute_stacey_electromagnetic.f90 b/src/specfem2D/compute_stacey_electromagnetic.f90 index 51a92ddc3..e3e14ef45 100644 --- a/src/specfem2D/compute_stacey_electromagnetic.f90 +++ b/src/specfem2D/compute_stacey_electromagnetic.f90 @@ -50,7 +50,7 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag b_absorb_electromagnetic_left,b_absorb_electromagnetic_right, & b_absorb_electromagnetic_bottom,b_absorb_electromagnetic_top, & ib_left,ib_right,ib_bottom,ib_top, & - anyabs,STACEY_ABSORBING_CONDITIONS,Qe11_electromagnetic,Qe33_electromagnetic,& + anyabs,STACEY_ABSORBING_CONDITIONS,Qe11_electromagnetic,Qe33_electromagnetic, & Qs11_electromagnetic,Qs33_electromagnetic,P_SV implicit none @@ -101,14 +101,14 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag permlxx = spermittivitystore(1,i,j,ispec) !e11 permlzz = spermittivitystore(2,i,j,ispec) !e33 condlxx = sconductivitystore(1,i,j,ispec) !sig11 - condlzz = sconductivitystore(2,i,j,ispec) !sig33 + condlzz = sconductivitystore(2,i,j,ispec) !sig33 two_inv_magpermeability = 2.d0 * inv_magpermeabilitystore(i,j,ispec) !2mu0^-1 - + call get_electromagnetic_velocities(cpxsquare,cpzsquare,econdl,eperml,ATTENUATION_CONDUCTIVITY, & - ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)),& + ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)), & Qs11_electromagnetic(kmato(ispec)),Qs33_electromagnetic(kmato(ispec)), & permlxx,permlzz,condlxx,condlzz,two_inv_magpermeability) - + cpxl = sqrt(cpxsquare) cpzl = sqrt(cpzsquare) @@ -117,7 +117,7 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag if (P_SV) then ! P_SV case - vx = veloc_electromagnetic(1,iglob) + vx = veloc_electromagnetic(1,iglob) vz = veloc_electromagnetic(2,iglob) vn = nx*vx + nz*vz @@ -128,7 +128,7 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag ! with velocity's tangential component: veloc_tangential = v - vn * normal ! total traction ! T = T_normal + T_tangential - tx = rho_cpx*vn*nx + tx = rho_cpx*vn*nx tz = rho_cpz*vn*nz else @@ -178,14 +178,14 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag permlxx = spermittivitystore(1,i,j,ispec) !e11 permlzz = spermittivitystore(2,i,j,ispec) !e33 condlxx = sconductivitystore(1,i,j,ispec) !sig11 - condlzz = sconductivitystore(2,i,j,ispec) !sig33 + condlzz = sconductivitystore(2,i,j,ispec) !sig33 two_inv_magpermeability = 2.d0 * inv_magpermeabilitystore(i,j,ispec) !2mu0^-1 - + call get_electromagnetic_velocities(cpxsquare,cpzsquare,econdl,eperml,ATTENUATION_CONDUCTIVITY, & - ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)),& + ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)), & Qs11_electromagnetic(kmato(ispec)),Qs33_electromagnetic(kmato(ispec)), & permlxx,permlzz,condlxx,condlzz,two_inv_magpermeability) - + cpxl = sqrt(cpxsquare) cpzl = sqrt(cpzsquare) @@ -194,7 +194,7 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag if (P_SV) then ! P_SV case - vx = veloc_electromagnetic(1,iglob) + vx = veloc_electromagnetic(1,iglob) vz = veloc_electromagnetic(2,iglob) vn = nx*vx + nz*vz tx = rho_cpx*vn*nx @@ -215,7 +215,7 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag ! SH case accel_electromagnetic(1,iglob) = accel_electromagnetic(1,iglob) - ty*weight endif - + if (SAVE_FORWARD .and. SIMULATION_TYPE == 1) then if (P_SV) then ! P_SV case @@ -245,14 +245,14 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag permlxx = spermittivitystore(1,i,j,ispec) !e11 permlzz = spermittivitystore(2,i,j,ispec) !e33 condlxx = sconductivitystore(1,i,j,ispec) !sig11 - condlzz = sconductivitystore(2,i,j,ispec) !sig33 + condlzz = sconductivitystore(2,i,j,ispec) !sig33 two_inv_magpermeability = 2.d0 * inv_magpermeabilitystore(i,j,ispec) !2mu0^-1 - + call get_electromagnetic_velocities(cpxsquare,cpzsquare,econdl,eperml,ATTENUATION_CONDUCTIVITY, & - ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)),& + ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)), & Qs11_electromagnetic(kmato(ispec)),Qs33_electromagnetic(kmato(ispec)), & permlxx,permlzz,condlxx,condlzz,two_inv_magpermeability) - + cpxl = sqrt(cpxsquare) cpzl = sqrt(cpzsquare) @@ -261,7 +261,7 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag if (P_SV) then ! P_SV case - vx = veloc_electromagnetic(1,iglob) + vx = veloc_electromagnetic(1,iglob) vz = veloc_electromagnetic(2,iglob) vn = nx*vx+nz*vz tx = rho_cpx*vn*nx @@ -326,14 +326,14 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag permlxx = spermittivitystore(1,i,j,ispec) !e11 permlzz = spermittivitystore(2,i,j,ispec) !e33 condlxx = sconductivitystore(1,i,j,ispec) !sig11 - condlzz = sconductivitystore(2,i,j,ispec) !sig33 + condlzz = sconductivitystore(2,i,j,ispec) !sig33 two_inv_magpermeability = 2.d0 * inv_magpermeabilitystore(i,j,ispec) !2mu0^-1 - + call get_electromagnetic_velocities(cpxsquare,cpzsquare,econdl,eperml,ATTENUATION_CONDUCTIVITY, & - ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)),& + ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)), & Qs11_electromagnetic(kmato(ispec)),Qs33_electromagnetic(kmato(ispec)), & permlxx,permlzz,condlxx,condlzz,two_inv_magpermeability) - + cpxl = sqrt(cpxsquare) cpzl = sqrt(cpzsquare) @@ -342,7 +342,7 @@ subroutine compute_stacey_electromagnetic(accel_electromagnetic,veloc_electromag if (P_SV) then ! P_SV case - vx = veloc_electromagnetic(1,iglob) + vx = veloc_electromagnetic(1,iglob) vz = veloc_electromagnetic(2,iglob) vn = nx*vx+nz*vz tx = rho_cpx*vn*nx diff --git a/src/specfem2D/invert_mass_matrix.F90 b/src/specfem2D/invert_mass_matrix.F90 index db7686c94..00557b5aa 100644 --- a/src/specfem2D/invert_mass_matrix.F90 +++ b/src/specfem2D/invert_mass_matrix.F90 @@ -43,7 +43,7 @@ subroutine invert_mass_matrix_init() rmass_inverse_elastic, & rmass_inverse_acoustic,rmass_inverse_e1,ATTENUATION_VISCOACOUSTIC,phi_nu1,N_SLS, & time_stepping_scheme,rmass_s_inverse_poroelastic,rmass_w_inverse_poroelastic, & - rmass_inverse_electromagnetic, & + rmass_inverse_electromagnetic, & nspec,ibool,wxgll,wzgll,jacobian, & ispec_is_elastic,ispec_is_acoustic,ispec_is_poroelastic,ispec_is_electromagnetic, & phistore,tortstore,rhoarraystore, & @@ -457,19 +457,19 @@ subroutine invert_mass_matrix_init() endif else if (ispec_is_electromagnetic(ispec)) then - ! electromagnetic medium + ! electromagnetic medium !!! PML NOT WORKING YET !!! this_element_has_PML = .false. if (PML_BOUNDARY_CONDITIONS) then if (ispec_is_PML(ispec)) stop 'PML not implemented yet for electromagnetic case' endif - + permlxx = spermittivitystore(1,i,j,ispec) !e11 permlzz = spermittivitystore(2,i,j,ispec) !e33 condlxx = sconductivitystore(1,i,j,ispec) !sig11 - condlzz = sconductivitystore(2,i,j,ispec) !sig33 - + condlzz = sconductivitystore(2,i,j,ispec) !sig33 + ! permittivity (Zener model) taud_x = (sqrt(Qe11_electromagnetic(kmato(ispec))**2+1.d0) +1.d0)/ & (2.d0*pi*f0_electromagnetic*Qe11_electromagnetic(kmato(ispec))) @@ -886,11 +886,11 @@ subroutine invert_mass_matrix_init_Stacey() permlxx = spermittivitystore(1,i,j,ispec) !e11 permlzz = spermittivitystore(2,i,j,ispec) !e33 condlxx = sconductivitystore(1,i,j,ispec) !sig11 - condlzz = sconductivitystore(2,i,j,ispec) !sig33 + condlzz = sconductivitystore(2,i,j,ispec) !sig33 two_inv_magpermeability = 2.d0 * inv_magpermeabilitystore(i,j,ispec) !2mu0^-1 call get_electromagnetic_velocities(cpxsquare,cpzsquare,econdl,eperml,ATTENUATION_CONDUCTIVITY, & - ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)),& + ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)), & Qs11_electromagnetic(kmato(ispec)),Qs33_electromagnetic(kmato(ispec)), & permlxx,permlzz,condlxx,condlzz,two_inv_magpermeability) @@ -924,11 +924,11 @@ subroutine invert_mass_matrix_init_Stacey() permlxx = spermittivitystore(1,i,j,ispec) !e11 permlzz = spermittivitystore(2,i,j,ispec) !e33 condlxx = sconductivitystore(1,i,j,ispec) !sig11 - condlzz = sconductivitystore(2,i,j,ispec) !sig33 + condlzz = sconductivitystore(2,i,j,ispec) !sig33 two_inv_magpermeability = 2.d0 * inv_magpermeabilitystore(i,j,ispec) !2mu0^-1 call get_electromagnetic_velocities(cpxsquare,cpzsquare,econdl,eperml,ATTENUATION_CONDUCTIVITY, & - ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)),& + ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)), & Qs11_electromagnetic(kmato(ispec)),Qs33_electromagnetic(kmato(ispec)), & permlxx,permlzz,condlxx,condlzz,two_inv_magpermeability) @@ -962,11 +962,11 @@ subroutine invert_mass_matrix_init_Stacey() permlxx = spermittivitystore(1,i,j,ispec) !e11 permlzz = spermittivitystore(2,i,j,ispec) !e33 condlxx = sconductivitystore(1,i,j,ispec) !sig11 - condlzz = sconductivitystore(2,i,j,ispec) !sig33 + condlzz = sconductivitystore(2,i,j,ispec) !sig33 two_inv_magpermeability = 2.d0 * inv_magpermeabilitystore(i,j,ispec) !2mu0^-1 call get_electromagnetic_velocities(cpxsquare,cpzsquare,econdl,eperml,ATTENUATION_CONDUCTIVITY, & - ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)),& + ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)), & Qs11_electromagnetic(kmato(ispec)),Qs33_electromagnetic(kmato(ispec)), & permlxx,permlzz,condlxx,condlzz,two_inv_magpermeability) @@ -1007,11 +1007,11 @@ subroutine invert_mass_matrix_init_Stacey() permlxx = spermittivitystore(1,i,j,ispec) !e11 permlzz = spermittivitystore(2,i,j,ispec) !e33 condlxx = sconductivitystore(1,i,j,ispec) !sig11 - condlzz = sconductivitystore(2,i,j,ispec) !sig33 + condlzz = sconductivitystore(2,i,j,ispec) !sig33 two_inv_magpermeability = 2.d0 * inv_magpermeabilitystore(i,j,ispec) !2mu0^-1 call get_electromagnetic_velocities(cpxsquare,cpzsquare,econdl,eperml,ATTENUATION_CONDUCTIVITY, & - ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)),& + ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)), & Qs11_electromagnetic(kmato(ispec)),Qs33_electromagnetic(kmato(ispec)), & permlxx,permlzz,condlxx,condlzz,two_inv_magpermeability) @@ -1041,7 +1041,7 @@ subroutine invert_mass_matrix_init_Stacey() rmass_inverse_electromagnetic(2,iglob) = rmass_inverse_electromagnetic(2,iglob) + tz*weight enddo endif ! end of top absorbing boundary - + enddo endif ! any_electromagnetic diff --git a/src/specfem2D/iterate_time.F90 b/src/specfem2D/iterate_time.F90 index 77484777c..79e267794 100644 --- a/src/specfem2D/iterate_time.F90 +++ b/src/specfem2D/iterate_time.F90 @@ -302,7 +302,7 @@ subroutine iterate_time() if (any_electromagnetic) call exit_MPI(myrank,'electromagnetic not implemented in GPU MODE yet') endif endif - + case (3) ! kernel simulations ! forward (adjoint) wavefields use switched update ordering diff --git a/src/specfem2D/prepare_optimized_arrays.F90 b/src/specfem2D/prepare_optimized_arrays.F90 index 928105c05..611948107 100644 --- a/src/specfem2D/prepare_optimized_arrays.F90 +++ b/src/specfem2D/prepare_optimized_arrays.F90 @@ -149,7 +149,7 @@ subroutine prepare_fused_array() deriv_mapping(4,INDEX_IJ,ispec) = gammaz(INDEX_IJ,ispec) deriv_mapping(5,INDEX_IJ,ispec) = jacobian(INDEX_IJ,ispec) - + if (any_acoustic) then ! for acoustic elements ! density model diff --git a/src/specfem2D/prepare_wavefields.f90 b/src/specfem2D/prepare_wavefields.f90 index 9aa13d842..5f895ff32 100644 --- a/src/specfem2D/prepare_wavefields.f90 +++ b/src/specfem2D/prepare_wavefields.f90 @@ -464,7 +464,7 @@ subroutine prepare_wavefields() modeAmplitude(:) = 0.0d0 ! -! electromagnetic material +! electromagnetic material ! if (ELECTROMAGNETIC_SIMULATION) then if (myrank == 0) then @@ -480,7 +480,7 @@ subroutine prepare_wavefields() ! dummy allocate unused arrays with fictitious size nglob_electromagnetic = 1 endif - + allocate(displ_electromagnetic(NDIM,nglob_electromagnetic), & veloc_electromagnetic(NDIM,nglob_electromagnetic), & accel_electromagnetic(NDIM,nglob_electromagnetic),stat=ier) diff --git a/src/specfem2D/read_materials.f90 b/src/specfem2D/read_materials.f90 index fe8d5ca20..1dfc30a2f 100644 --- a/src/specfem2D/read_materials.f90 +++ b/src/specfem2D/read_materials.f90 @@ -42,9 +42,9 @@ subroutine read_materials(f0) use specfem_par, only: AXISYM,density,porosity,tortuosity,anisotropycoef,permeability,poroelastcoef, & numat,myrank,QKappa_attenuationcoef,Qmu_attenuationcoef, & USE_CONDUCTIVE_DIFFUSION,ATTENUATION_CONDUCTIVITY, & - ATTENUATION_PERMITTIVITY,f0_electromagnetic,& - Qe11_electromagnetic,Qe33_electromagnetic,Qs11_electromagnetic,Qs33_electromagnetic,& - spermittivity,sconductivity,inv_magpermeability,& + ATTENUATION_PERMITTIVITY,f0_electromagnetic, & + Qe11_electromagnetic,Qe33_electromagnetic,Qs11_electromagnetic,Qs33_electromagnetic, & + spermittivity,sconductivity,inv_magpermeability, & freq0_poroelastic,Q0_poroelastic,ATTENUATION_PORO_FLUID_PART, & use_external_velocity_model,tomo_material,myrank @@ -96,7 +96,7 @@ subroutine read_materials(f0) Qs11_electromagnetic(:) = 9999. Qs33_electromagnetic(:) = 9999. - USE_CONDUCTIVE_DIFFUSION = .false. ! for conduction diffusion in electromagnetic material. + USE_CONDUCTIVE_DIFFUSION = .false. ! for conduction diffusion in electromagnetic material. ! Index of the material that will be defined by an external tomo file if needed (TOMOGRAPHY_FILE) tomo_material = 0 @@ -320,7 +320,7 @@ subroutine read_materials(f0) ! Poisson's ratio must be between -1 and +1/2 if (poisson_s < -1.d0 .or. poisson_s > 0.5d0) call stop_the_code('Poisson''s ratio for the solid phase out of range') - + else if (indic == ELECTROMAGNETIC_MATERIAL ) then ! model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 spermittivity(1,n) = val1 * val2 !e11 @@ -428,11 +428,11 @@ subroutine read_materials(f0) permlxx = spermittivity(1,n) !e11 permlzz = spermittivity(2,n) !e33 condlxx = sconductivity(1,n) !sig11 - condlzz = sconductivity(2,n) !sig33 + condlzz = sconductivity(2,n) !sig33 two_inv_magpermeability = 2.d0 * inv_magpermeability(n) !2mu0^-1 call get_electromagnetic_velocities(cpxsquare,cpzsquare,econdl,eperml,ATTENUATION_CONDUCTIVITY, & - ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(n),Qe33_electromagnetic(n),& + ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(n),Qe33_electromagnetic(n), & Qs11_electromagnetic(n),Qs33_electromagnetic(n), & permlxx,permlzz,condlxx,condlzz,two_inv_magpermeability) @@ -652,8 +652,8 @@ subroutine read_materials(f0) '-- Electromagnetic material --',/5x, & '-------------------------------',/5x, & 'Material set number. . . . . . . . (jmat) =',i6,/5x, & - 'magnetic permeability (free space). . . . . (mu0=4xpix10-7 H.m^-1) =',1pe15.8,/5x, & - 'permittivity (free space). . . . . . . (e0=8.85x10-12 F.m^-1) =',1pe15.8,/5x, & + 'magnetic permeability (free space). . . . . (mu0 = 4xpix10-7 H.m^-1) =',1pe15.8,/5x, & + 'permittivity (free space). . . . . . . (e0 = 8.85x10-12 F.m^-1) =',1pe15.8,/5x, & 'permittivity xx component. . . (e11(e0)) =',1pe15.8,/5x, & 'permittivity zz component. . . (e33(e0)) =',1pe15.8,/5x, & 'conductivity xx component. . . (sig11) =',1pe15.8,/5x, & diff --git a/src/specfem2D/read_mesh_databases.F90 b/src/specfem2D/read_mesh_databases.F90 index f15ec3893..98628e98d 100644 --- a/src/specfem2D/read_mesh_databases.F90 +++ b/src/specfem2D/read_mesh_databases.F90 @@ -241,7 +241,7 @@ subroutine read_mesh_for_init() read(IIN) local_dble,local_dble ! Q0_poroelastic,freq0_poroelastic - read(IIN) local_l1,local_l2,local_dble ! ATTENUATION_PERMITTIVITY,ATTENUATION_CONDUCTIVITY,f0_electromagnetic + read(IIN) local_l1,local_l2,local_dble ! ATTENUATION_PERMITTIVITY,ATTENUATION_CONDUCTIVITY,f0_electromagnetic read(IIN) local_l ! AXISYM if (local_l .neqv. AXISYM) then diff --git a/src/specfem2D/save_model_files.f90 b/src/specfem2D/save_model_files.f90 index a59e9cebd..1b61146e3 100644 --- a/src/specfem2D/save_model_files.f90 +++ b/src/specfem2D/save_model_files.f90 @@ -272,7 +272,7 @@ subroutine save_model_files() write(172) inv_magpermeability_save close(172) if (myrank == 0) write(IMAIN,*) ' written to file(s): ',trim(filename) - endif + endif else call stop_the_code('Save Model not implemented for external and tomo') endif !Type of model diff --git a/src/specfem2D/setup_mesh.F90 b/src/specfem2D/setup_mesh.F90 index 48ca05992..cf7ce6f2b 100644 --- a/src/specfem2D/setup_mesh.F90 +++ b/src/specfem2D/setup_mesh.F90 @@ -544,7 +544,7 @@ subroutine setup_mesh_material_properties() c33ext,c35ext,c55ext real(kind=CUSTOM_REAL), dimension(:,:,:,:), allocatable :: spermittivityext,sconductivityext real(kind=CUSTOM_REAL), dimension(:,:,:), allocatable :: inv_magpermeabilityext - + ! for shifting of velocities if needed in the case of viscoelasticity double precision :: vp,vs,rhol,mul,lambdal,kappal,qmul,qkappal double precision :: phi,tort,kappa_s,kappa_f,kappa_fr,mu_s,mu_fr @@ -865,32 +865,32 @@ subroutine setup_mesh_material_properties() permlxx = spermittivityext(1,i,j,ispec) !e11 permlzz = spermittivityext(2,i,j,ispec) !e33 condlxx = sconductivityext(1,i,j,ispec) !sig11 - condlzz = sconductivityext(2,i,j,ispec) !sig33 + condlzz = sconductivityext(2,i,j,ispec) !sig33 two_inv_magpermeability = 2.d0 * inv_magpermeabilityext(i,j,ispec) !2mu0^-1 else ! internal mesh permlxx = spermittivity(1,kmato(ispec)) !e11 permlzz = spermittivity(2,kmato(ispec)) !e33 condlxx = sconductivity(1,kmato(ispec)) !sig11 - condlzz = sconductivity(2,kmato(ispec)) !sig33 + condlzz = sconductivity(2,kmato(ispec)) !sig33 two_inv_magpermeability = 2.d0 * inv_magpermeability(kmato(ispec)) !2mu0^-1 endif call get_electromagnetic_velocities(cpxsquare,cpzsquare,econdl,eperml,ATTENUATION_CONDUCTIVITY, & - ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)),& + ATTENUATION_PERMITTIVITY,f0_electromagnetic,Qe11_electromagnetic(kmato(ispec)),Qe33_electromagnetic(kmato(ispec)), & Qs11_electromagnetic(kmato(ispec)),Qs33_electromagnetic(kmato(ispec)), & permlxx,permlzz,condlxx,condlzz,two_inv_magpermeability) - ! stores - vEstore(i,j,ispec)=max(sqrt(cpxsquare),sqrt(cpzsquare)) + ! stores + vEstore(i,j,ispec)=max(sqrt(cpxsquare),sqrt(cpzsquare)) spermittivitystore(1,i,j,ispec)=permlxx spermittivitystore(2,i,j,ispec)=permlzz sconductivitystore(1,i,j,ispec)=condlxx sconductivitystore(2,i,j,ispec)=condlzz inv_magpermeabilitystore(i,j,ispec)=two_inv_magpermeability/2.d0 - vEmin_glob = min(vEmin_glob,max(sqrt(cpxsquare),sqrt(cpzsquare))) - vEmax_glob = max(vEmax_glob,max(sqrt(cpxsquare),sqrt(cpzsquare))) + vEmin_glob = min(vEmin_glob,max(sqrt(cpxsquare),sqrt(cpzsquare))) + vEmax_glob = max(vEmax_glob,max(sqrt(cpxsquare),sqrt(cpzsquare))) endif enddo enddo @@ -985,7 +985,7 @@ subroutine setup_mesh_material_properties() deallocate(spermittivityext,sconductivityext,inv_magpermeabilityext) ! stats - ! poroelasticity + ! poroelasticity ! check if any poroelastic elements in domains call any_all_l(any_elastic,has_elasticity) if (has_elasticity) then diff --git a/src/specfem2D/specfem2D_par.f90 b/src/specfem2D/specfem2D_par.f90 index 566b9221f..0631357a2 100644 --- a/src/specfem2D/specfem2D_par.f90 +++ b/src/specfem2D/specfem2D_par.f90 @@ -108,8 +108,8 @@ module specfem_par ! electromagnetic real(kind=CUSTOM_REAL), dimension(:,:,:,:), allocatable :: spermittivitystore,sconductivitystore - real(kind=CUSTOM_REAL), dimension(:,:,:), allocatable :: inv_magpermeabilitystore,vEstore - + real(kind=CUSTOM_REAL), dimension(:,:,:), allocatable :: inv_magpermeabilitystore,vEstore + ! resolution double precision :: mesh_T_min @@ -614,20 +614,20 @@ module specfem_par ! global flag for electromagnetic simulations logical :: ELECTROMAGNETIC_SIMULATION logical :: USE_CONDUCTIVE_DIFFUSION - + logical, dimension(:), allocatable :: ispec_is_electromagnetic - + ! inverse mass matrices real(kind=CUSTOM_REAL), dimension(:,:), allocatable :: rmass_inverse_electromagnetic - ! material properties of the electromagnetic medium + ! material properties of the electromagnetic medium real(kind=CUSTOM_REAL), dimension(:,:), allocatable :: & accel_electromagnetic,veloc_electromagnetic,displ_electromagnetic - + double precision, dimension(:,:), allocatable :: spermittivity,sconductivity double precision, dimension(:), allocatable :: inv_magpermeability - ! for permittivity & conductivity attenuation + ! for permittivity & conductivity attenuation double precision, dimension(:), allocatable :: Qe11_electromagnetic,Qe33_electromagnetic double precision, dimension(:), allocatable :: Qs11_electromagnetic,Qs33_electromagnetic double precision, dimension(:,:,:,:), allocatable :: tau_e,tau_d,tau_s @@ -703,7 +703,7 @@ module specfem_par ! poroelastic fluid real(kind=CUSTOM_REAL), dimension(:,:,:,:), allocatable :: & b_absorb_poro_w_left,b_absorb_poro_w_right,b_absorb_poro_w_bottom,b_absorb_poro_w_top - ! electromagnetic + ! electromagnetic real(kind=CUSTOM_REAL), dimension(:,:,:,:), allocatable :: & b_absorb_electromagnetic_left,b_absorb_electromagnetic_right,b_absorb_electromagnetic_bottom,b_absorb_electromagnetic_top From aa370ca575dd395e1a7afb4403fd8c95c85055ba Mon Sep 17 00:00:00 2001 From: Daniel Peter Date: Wed, 5 Jun 2024 12:52:30 +0200 Subject: [PATCH 7/9] more Par_file updates --- EXAMPLES/Gmsh_example_CPML_MPI/DATA/Par_file | 3 ++- .../Marmousi_mesh_of_the_model/DATA/Par_file | 3 ++- .../DATA/Par_file.mesh-improved | 3 ++- .../REF_SEIS/Par_file | 3 ++- .../anisotropic_isotropic_model/DATA/Par_file | 3 ++- .../Par_file_homogeneous_with_attenuation | 3 ++- ...ation_2D_source_in_point_inside_an_element | 2 +- ...2D_source_in_point_shared_by_four_elements | 2 +- EXAMPLES/canyon/DATA/Par_file | 2 +- EXAMPLES/electromagnetic_GPR_TM/DATA/Par_file | 2 +- .../fluid_solid_external_mesh/DATA/Par_file | 2 +- .../DATA/Par_file | 2 +- .../DATA/Par_file | 2 +- .../DATA/Par_file | 2 +- .../DATA/Par_file | 2 +- EXAMPLES/noise_uniform/DATA/Par_file | 21 +++++++++++++++++-- .../paper_axisymmetry_example/DATA/Par_file | 3 ++- .../Pekeris_waveguide/DATA/Par_file | 2 +- .../with_free_bottom_surface/DATA/Par_file | 2 +- .../CPML_homogeneous/DATA/Par_file | 2 +- .../CPML_normal_fluid_solid/DATA/Par_file | 2 +- .../CPML_normal_solid_only/DATA/Par_file | 2 +- .../Stacey_homogeneous/DATA/Par_file | 2 +- .../Stacey_normal_fluid_solid/DATA/Par_file | 2 +- .../REF_KERNEL/Par_file | 3 ++- EXAMPLES/thermocline/DATA/Par_file | 2 +- .../tomographic_ocean_model/DATA/Par_file | 3 ++- .../tomographic_ocean_model/REF_SEIS/Par_file | 3 ++- 28 files changed, 56 insertions(+), 29 deletions(-) diff --git a/EXAMPLES/Gmsh_example_CPML_MPI/DATA/Par_file b/EXAMPLES/Gmsh_example_CPML_MPI/DATA/Par_file index a4f38b451..83766b342 100644 --- a/EXAMPLES/Gmsh_example_CPML_MPI/DATA/Par_file +++ b/EXAMPLES/Gmsh_example_CPML_MPI/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -288,6 +288,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Marmousi_mesh_of_the_model/DATA/Par_file b/EXAMPLES/Marmousi_mesh_of_the_model/DATA/Par_file index 0b18dcc8b..52585e4b4 100644 --- a/EXAMPLES/Marmousi_mesh_of_the_model/DATA/Par_file +++ b/EXAMPLES/Marmousi_mesh_of_the_model/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Marmousi_mesh_of_the_model/DATA/Par_file.mesh-improved b/EXAMPLES/Marmousi_mesh_of_the_model/DATA/Par_file.mesh-improved index 50e4985e6..c059046ac 100644 --- a/EXAMPLES/Marmousi_mesh_of_the_model/DATA/Par_file.mesh-improved +++ b/EXAMPLES/Marmousi_mesh_of_the_model/DATA/Par_file.mesh-improved @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/Marmousi_mesh_of_the_model/REF_SEIS/Par_file b/EXAMPLES/Marmousi_mesh_of_the_model/REF_SEIS/Par_file index 0b18dcc8b..52585e4b4 100644 --- a/EXAMPLES/Marmousi_mesh_of_the_model/REF_SEIS/Par_file +++ b/EXAMPLES/Marmousi_mesh_of_the_model/REF_SEIS/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 1 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/anisotropic_isotropic_model/DATA/Par_file b/EXAMPLES/anisotropic_isotropic_model/DATA/Par_file index 745d89f1e..50a04c776 100644 --- a/EXAMPLES/anisotropic_isotropic_model/DATA/Par_file +++ b/EXAMPLES/anisotropic_isotropic_model/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 3 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 0 0 # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 0 0 # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/anisotropic_isotropic_model/DATA/Par_file_homogeneous_with_attenuation b/EXAMPLES/anisotropic_isotropic_model/DATA/Par_file_homogeneous_with_attenuation index c72021169..3ea5fae6c 100644 --- a/EXAMPLES/anisotropic_isotropic_model/DATA/Par_file_homogeneous_with_attenuation +++ b/EXAMPLES/anisotropic_isotropic_model/DATA/Par_file_homogeneous_with_attenuation @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 3 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 0 0 # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 0 0 # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/attenuation/viscoelastic/DATA/Par_file_attenuation_2D_source_in_point_inside_an_element b/EXAMPLES/attenuation/viscoelastic/DATA/Par_file_attenuation_2D_source_in_point_inside_an_element index 318e77b01..dc40e65b5 100644 --- a/EXAMPLES/attenuation/viscoelastic/DATA/Par_file_attenuation_2D_source_in_point_inside_an_element +++ b/EXAMPLES/attenuation/viscoelastic/DATA/Par_file_attenuation_2D_source_in_point_inside_an_element @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/attenuation/viscoelastic/DATA/Par_file_attenuation_2D_source_in_point_shared_by_four_elements b/EXAMPLES/attenuation/viscoelastic/DATA/Par_file_attenuation_2D_source_in_point_shared_by_four_elements index bb1c53c95..fd6e32de2 100644 --- a/EXAMPLES/attenuation/viscoelastic/DATA/Par_file_attenuation_2D_source_in_point_shared_by_four_elements +++ b/EXAMPLES/attenuation/viscoelastic/DATA/Par_file_attenuation_2D_source_in_point_shared_by_four_elements @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/canyon/DATA/Par_file b/EXAMPLES/canyon/DATA/Par_file index 7c532ba24..684536278 100644 --- a/EXAMPLES/canyon/DATA/Par_file +++ b/EXAMPLES/canyon/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/electromagnetic_GPR_TM/DATA/Par_file b/EXAMPLES/electromagnetic_GPR_TM/DATA/Par_file index 310e70d08..1a64ece9e 100644 --- a/EXAMPLES/electromagnetic_GPR_TM/DATA/Par_file +++ b/EXAMPLES/electromagnetic_GPR_TM/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/fluid_solid/fluid_solid_external_mesh/DATA/Par_file b/EXAMPLES/fluid_solid/fluid_solid_external_mesh/DATA/Par_file index b91f98ec3..46f196658 100644 --- a/EXAMPLES/fluid_solid/fluid_solid_external_mesh/DATA/Par_file +++ b/EXAMPLES/fluid_solid/fluid_solid_external_mesh/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/DATA/Par_file b/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/DATA/Par_file index 9078539e5..3e652abff 100644 --- a/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/DATA/Par_file +++ b/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/DATA/Par_file b/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/DATA/Par_file index 707992d14..ff072e10a 100644 --- a/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/DATA/Par_file +++ b/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/infinite_homogeneous_moment_tensor_vertical_dip_slip/DATA/Par_file b/EXAMPLES/infinite_homogeneous_moment_tensor_vertical_dip_slip/DATA/Par_file index 865e4aa72..09e83118c 100644 --- a/EXAMPLES/infinite_homogeneous_moment_tensor_vertical_dip_slip/DATA/Par_file +++ b/EXAMPLES/infinite_homogeneous_moment_tensor_vertical_dip_slip/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/initial_plane_wave_with_free_surface/DATA/Par_file b/EXAMPLES/initial_plane_wave_with_free_surface/DATA/Par_file index fd10df566..d36ff3190 100644 --- a/EXAMPLES/initial_plane_wave_with_free_surface/DATA/Par_file +++ b/EXAMPLES/initial_plane_wave_with_free_surface/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/noise_uniform/DATA/Par_file b/EXAMPLES/noise_uniform/DATA/Par_file index 2c7855a0e..2b270cdfa 100644 --- a/EXAMPLES/noise_uniform/DATA/Par_file +++ b/EXAMPLES/noise_uniform/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -274,7 +274,24 @@ tangential_detection_curve_file = dummy # file # material properties # number of model materials nbmodels = 1 -# define models as +# available material types (see user manual for more information) +# acoustic: model_number 1 rho Vp 0 0 0 QKappa 9999 0 0 0 0 0 0 (for QKappa use 9999 to ignore it) +# elastic: model_number 1 rho Vp Vs 0 0 QKappa Qmu 0 0 0 0 0 0 (for QKappa and Qmu use 9999 to ignore them) +# anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu +# anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu +# poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 +# tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 +# +# note: When viscoelasticity or viscoacousticity is turned on, +# the Vp and Vs values that are read here are the UNRELAXED ones i.e. the values at infinite frequency +# unless the READ_VELOCITIES_AT_f0 parameter above is set to true, in which case they are the values at frequency f0. +# +# Please also note that Qmu is always equal to Qs, but Qkappa is in general not equal to Qp. +# To convert one to the other see doc/Qkappa_Qmu_versus_Qp_Qs_relationship_in_2D_plane_strain.pdf and +# utils/attenuation/conversion_from_Qkappa_Qmu_to_Qp_Qs_from_Dahlen_Tromp_959_960.f90. +# +# model formats: # I: (model_number 1 rho Vp Vs 0 0 QKappa Qmu 0 0 0 0 0 0) or # II: (model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 0 0) or # III: (model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu). diff --git a/EXAMPLES/paper_axisymmetry_example/DATA/Par_file b/EXAMPLES/paper_axisymmetry_example/DATA/Par_file index 688f5017e..393cbf417 100644 --- a/EXAMPLES/paper_axisymmetry_example/DATA/Par_file +++ b/EXAMPLES/paper_axisymmetry_example/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -280,6 +280,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/Pekeris_waveguide/DATA/Par_file b/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/Pekeris_waveguide/DATA/Par_file index 085a15750..e1399fd50 100644 --- a/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/Pekeris_waveguide/DATA/Par_file +++ b/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/Pekeris_waveguide/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/with_free_bottom_surface/DATA/Par_file b/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/with_free_bottom_surface/DATA/Par_file index 7b072b760..bf2825927 100644 --- a/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/with_free_bottom_surface/DATA/Par_file +++ b/EXAMPLES/perfect_idealized_oceanic_guide_in_underwater_acoustics/with_free_bottom_surface/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_homogeneous/DATA/Par_file b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_homogeneous/DATA/Par_file index fc3df43af..9adc8855b 100644 --- a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_homogeneous/DATA/Par_file +++ b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_homogeneous/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_normal_fluid_solid/DATA/Par_file b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_normal_fluid_solid/DATA/Par_file index e1eee148e..24be14a90 100644 --- a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_normal_fluid_solid/DATA/Par_file +++ b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_normal_fluid_solid/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_normal_solid_only/DATA/Par_file b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_normal_solid_only/DATA/Par_file index 706dfbf36..7999bdb31 100644 --- a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_normal_solid_only/DATA/Par_file +++ b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/CPML_normal_solid_only/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/Stacey_homogeneous/DATA/Par_file b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/Stacey_homogeneous/DATA/Par_file index ba135589b..379b8f8c3 100644 --- a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/Stacey_homogeneous/DATA/Par_file +++ b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/Stacey_homogeneous/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/Stacey_normal_fluid_solid/DATA/Par_file b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/Stacey_normal_fluid_solid/DATA/Par_file index ec167848d..6636e3c4b 100644 --- a/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/Stacey_normal_fluid_solid/DATA/Par_file +++ b/EXAMPLES/salt_dome_CUBIT_mesh_from_Ronan_Madec/Stacey_normal_fluid_solid/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/simple_topography_and_also_a_simple_fluid_layer/REF_KERNEL/Par_file b/EXAMPLES/simple_topography_and_also_a_simple_fluid_layer/REF_KERNEL/Par_file index 1c496fa6d..a823dad90 100644 --- a/EXAMPLES/simple_topography_and_also_a_simple_fluid_layer/REF_KERNEL/Par_file +++ b/EXAMPLES/simple_topography_and_also_a_simple_fluid_layer/REF_KERNEL/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -288,6 +288,7 @@ nbmodels = 4 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/thermocline/DATA/Par_file b/EXAMPLES/thermocline/DATA/Par_file index de1ca68b5..a99607255 100644 --- a/EXAMPLES/thermocline/DATA/Par_file +++ b/EXAMPLES/thermocline/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations diff --git a/EXAMPLES/tomographic_ocean_model/DATA/Par_file b/EXAMPLES/tomographic_ocean_model/DATA/Par_file index eae6ad36f..e248857c0 100644 --- a/EXAMPLES/tomographic_ocean_model/DATA/Par_file +++ b/EXAMPLES/tomographic_ocean_model/DATA/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -288,6 +288,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, diff --git a/EXAMPLES/tomographic_ocean_model/REF_SEIS/Par_file b/EXAMPLES/tomographic_ocean_model/REF_SEIS/Par_file index eae6ad36f..e248857c0 100644 --- a/EXAMPLES/tomographic_ocean_model/REF_SEIS/Par_file +++ b/EXAMPLES/tomographic_ocean_model/REF_SEIS/Par_file @@ -96,7 +96,7 @@ freq0_poroelastic = 10 # frequency for viscous attenua # for electromagnetic attenuation ATTENUATION_PERMITTIVITY = .false. # turn permittivity (Zener model) attenuation on or off for electromagnetic material ATTENUATION_CONDUCTIVITY = .false. # turn conductivity (Kelvin-Voigt model) attenuation on or off for electromagnetic material -f0_electromagnetic = 1d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file +f0_electromagnetic = 1.d9 # (Hz) relevant only if source is a Dirac or a Heaviside, otherwise it is f0 the dominant frequency of the source in the DATA/SOURCE file # to undo attenuation and/or PMLs for sensitivity kernel calculations or forward runs with SAVE_FORWARD # use the flag below. It performs undoing of attenuation and/or of PMLs in an exact way for sensitivity kernel calculations @@ -288,6 +288,7 @@ nbmodels = 2 # anisotropic: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 QKappa Qmu # anisotropic in AXISYM: model_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 c22 QKappa Qmu # poroelastic: model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu +# electromagnetic: model_number 4 mu0 e0 e11(e0) e33(e0) sig11 sig33 Qe11 Qe33 Qs11 Qs33 0 0 0 # tomo: model_number -1 0 0 A 0 0 0 0 0 0 0 0 0 0 # # note: When viscoelasticity or viscoacousticity is turned on, From 71525aa37fbe7e97f69b051826d23a204afb805e Mon Sep 17 00:00:00 2001 From: Daniel Peter Date: Wed, 5 Jun 2024 16:51:17 +0200 Subject: [PATCH 8/9] updates array allocations for EM --- .../compute_attenuation_electromagnetic.f90 | 10 +-- src/specfem2D/prepare_attenuation.f90 | 75 +++++++++++-------- src/specfem2D/setup_sources_receivers.F90 | 2 +- src/specfem2D/specfem2D_par.f90 | 4 +- 4 files changed, 52 insertions(+), 39 deletions(-) diff --git a/src/specfem2D/compute_attenuation_electromagnetic.f90 b/src/specfem2D/compute_attenuation_electromagnetic.f90 index c6c527749..67ab881ed 100644 --- a/src/specfem2D/compute_attenuation_electromagnetic.f90 +++ b/src/specfem2D/compute_attenuation_electromagnetic.f90 @@ -40,7 +40,7 @@ subroutine compute_attenuation_electromagnetic() use specfem_par, only: nspec,ispec_is_electromagnetic,ibool, & veloc_electromagnetic,time_stepping_scheme,time_stepping_scheme, & permx,permz,spermittivitystore,rx_permattenuation,rz_permattenuation, & - alphavalem,betavalem,gammavalem,tau_e,tau_d,ATTENUATION_PERMITTIVITY + alphaval_em,betaval_em,gammaval_em,tau_e,tau_d,ATTENUATION_PERMITTIVITY implicit none @@ -78,16 +78,16 @@ subroutine compute_attenuation_electromagnetic() (tau_d(i,j,ispec,1)*tau_d(i,j,ispec,1))*permx(i,j,ispec) Snp1 = - (1.d0 - tau_e(i,j,ispec,1)/tau_d(i,j,ispec,1))/ & (tau_d(i,j,ispec,1)*tau_d(i,j,ispec,1))*permx_loc(i,j) - rx_permattenuation(i,j,ispec) = alphavalem(i,j,ispec,1) * rx_permattenuation(i,j,ispec) & - + betavalem(i,j,ispec,1) * Sn + gammavalem(i,j,ispec,1) * Snp1 + rx_permattenuation(i,j,ispec) = alphaval_em(i,j,ispec,1) * rx_permattenuation(i,j,ispec) & + + betaval_em(i,j,ispec,1) * Sn + gammaval_em(i,j,ispec,1) * Snp1 ! evolution rz_permattenuation Sn = - (1.d0 - tau_e(i,j,ispec,2)/tau_d(i,j,ispec,2))/ & (tau_d(i,j,ispec,2)*tau_d(i,j,ispec,2))*permz(i,j,ispec) Snp1 = - (1.d0 - tau_e(i,j,ispec,2)/tau_d(i,j,ispec,2))/ & (tau_d(i,j,ispec,2)*tau_d(i,j,ispec,2))*permz_loc(i,j) - rz_permattenuation(i,j,ispec) = alphavalem(i,j,ispec,2) * rz_permattenuation(i,j,ispec) & - + betavalem(i,j,ispec,2) * Sn + gammavalem(i,j,ispec,2) * Snp1 + rz_permattenuation(i,j,ispec) = alphaval_em(i,j,ispec,2) * rz_permattenuation(i,j,ispec) & + + betaval_em(i,j,ispec,2) * Sn + gammaval_em(i,j,ispec,2) * Snp1 case (2) ! LDDRK diff --git a/src/specfem2D/prepare_attenuation.f90 b/src/specfem2D/prepare_attenuation.f90 index 37dd1d657..fafca8414 100644 --- a/src/specfem2D/prepare_attenuation.f90 +++ b/src/specfem2D/prepare_attenuation.f90 @@ -53,6 +53,9 @@ subroutine prepare_attenuation() phi_nu1_sent,phi_nu2_sent real(kind=CUSTOM_REAL), dimension(N_SLS) :: phinu,tauinvnu,temp,coef + ! EM attenuation + double precision :: tauinv(2) + ! attenuation ! user output call synchronize_all() @@ -502,16 +505,16 @@ subroutine prepare_attenuation() ! allocate memory variables for permittivity attenuation (electromagnetic media) if (ATTENUATION_PERMITTIVITY) then - allocate(rx_permattenuation(NGLLX,NGLLZ,nspec)) - allocate(rz_permattenuation(NGLLX,NGLLZ,nspec)) - allocate(permx(NGLLX,NGLLZ,nspec)) - allocate(permz(NGLLX,NGLLZ,nspec)) - allocate(tau_e(NGLLX,NGLLZ,nspec,2)) - allocate(tau_d(NGLLX,NGLLZ,nspec,2)) - allocate(tauinv(2)) - allocate(alphavalem(NGLLX,NGLLZ,nspec,2)) - allocate(betavalem(NGLLX,NGLLZ,nspec,2)) - allocate(gammavalem(NGLLX,NGLLZ,nspec,2)) + allocate(rx_permattenuation(NGLLX,NGLLZ,nspec), & + rz_permattenuation(NGLLX,NGLLZ,nspec), & + permx(NGLLX,NGLLZ,nspec), & + permz(NGLLX,NGLLZ,nspec), & + tau_e(NGLLX,NGLLZ,nspec,2), & + tau_d(NGLLX,NGLLZ,nspec,2), & + alphaval_em(NGLLX,NGLLZ,nspec,2), & + betaval_em(NGLLX,NGLLZ,nspec,2), & + gammaval_em(NGLLX,NGLLZ,nspec,2),stat=ier) + if (ier /= 0) stop 'Error allocating rx_permattenuation,.. arrays' ! initialize memory variables for attenuation rx_permattenuation(:,:,:) = 0.d0 @@ -519,11 +522,11 @@ subroutine prepare_attenuation() permx(:,:,:) = 0.d0 permz(:,:,:) = 0.d0 - ! if source is not a Dirac or Heavyside then f0_attenuation is f0 of the first - ! source - if (.not. (time_function_type(1) == 4 .or. time_function_type(1) == 5)) then - f0_electromagnetic = f0_source(1) - endif + ! if source is not a Dirac or Heavyside then f0_attenuation is f0 of the first + ! source + if (.not. (time_function_type(1) == 4 .or. time_function_type(1) == 5)) then + f0_electromagnetic = f0_source(1) + endif ! precompute Runge Kutta coefficients if viscous attenuation ! viscous attenuation is implemented following the memory variable @@ -534,25 +537,37 @@ subroutine prepare_attenuation() do j = 1,NGLLZ do i = 1,NGLLX - tau_e(i,j,ispec,1) = (sqrt(Qe11_electromagnetic(kmato(ispec))**2+1.d0) +1.d0)/ & - (2.d0*pi*f0_electromagnetic*Qe11_electromagnetic(kmato(ispec))) - tau_e(i,j,ispec,2) = (sqrt(Qe33_electromagnetic(kmato(ispec))**2+1.d0) +1.d0)/ & - (2.d0*pi*f0_electromagnetic*Qe33_electromagnetic(kmato(ispec))) - tau_d(i,j,ispec,1) = (sqrt(Qe11_electromagnetic(kmato(ispec))**2+1.d0) -1.d0)/ & - (2.d0*pi*f0_electromagnetic*Qe11_electromagnetic(kmato(ispec))) - tau_d(i,j,ispec,2) = (sqrt(Qe33_electromagnetic(kmato(ispec))**2+1.d0) -1.d0)/ & - (2.d0*pi*f0_electromagnetic*Qe33_electromagnetic(kmato(ispec))) - - tauinv(:) = - 1.d0 / tau_d(i,j,ispec,:) - alphavalem(i,j,ispec,:) = 1.d0 + deltat*tauinv(:) + deltat**2*tauinv(:)**2 / 2.d0 & - + deltat**3*tauinv(:)**3 / 6.d0 + deltat**4*tauinv(:)**4 / 24.d0 - betavalem(i,j,ispec,:) = deltat / 2.d0 + deltat**2*tauinv(:) / 3.d0 + deltat**3*tauinv(:)**2 / 8.d0 & - + deltat**4*tauinv(:)**3 / 24.d0 - gammavalem(i,j,ispec,:) = deltat / 2.d0 + deltat**2*tauinv(:) / 6.d0 + deltat**3*tauinv(:)**2 / 24.d0 + tau_e(i,j,ispec,1) = (sqrt(Qe11_electromagnetic(kmato(ispec))**2+1.d0) +1.d0)/ & + (2.d0*pi*f0_electromagnetic*Qe11_electromagnetic(kmato(ispec))) + tau_e(i,j,ispec,2) = (sqrt(Qe33_electromagnetic(kmato(ispec))**2+1.d0) +1.d0)/ & + (2.d0*pi*f0_electromagnetic*Qe33_electromagnetic(kmato(ispec))) + tau_d(i,j,ispec,1) = (sqrt(Qe11_electromagnetic(kmato(ispec))**2+1.d0) -1.d0)/ & + (2.d0*pi*f0_electromagnetic*Qe11_electromagnetic(kmato(ispec))) + tau_d(i,j,ispec,2) = (sqrt(Qe33_electromagnetic(kmato(ispec))**2+1.d0) -1.d0)/ & + (2.d0*pi*f0_electromagnetic*Qe33_electromagnetic(kmato(ispec))) + + tauinv(:) = - 1.d0 / tau_d(i,j,ispec,:) + alphaval_em(i,j,ispec,:) = 1.d0 + deltat*tauinv(:) + deltat**2*tauinv(:)**2 / 2.d0 & + + deltat**3*tauinv(:)**3 / 6.d0 + deltat**4*tauinv(:)**4 / 24.d0 + betaval_em(i,j,ispec,:) = deltat / 2.d0 + deltat**2*tauinv(:) / 3.d0 + deltat**3*tauinv(:)**2 / 8.d0 & + + deltat**4*tauinv(:)**3 / 24.d0 + gammaval_em(i,j,ispec,:) = deltat / 2.d0 + deltat**2*tauinv(:) / 6.d0 + deltat**3*tauinv(:)**2 / 24.d0 enddo enddo enddo + else + ! dummy allocation + ! (rx_permattenuation,rz_permattenuation allocation needed for subroutine call) + allocate(rx_permattenuation(1,1,1), & + rz_permattenuation(1,1,1), & + permx(1,1,1), & + permz(1,1,1), & + tau_e(1,1,1,1), & + tau_d(1,1,1,1), & + alphaval_em(1,1,1,1), & + betaval_em(1,1,1,1), & + gammaval_em(1,1,1,1)) endif ! synchronizes all processes diff --git a/src/specfem2D/setup_sources_receivers.F90 b/src/specfem2D/setup_sources_receivers.F90 index 6cc9747f4..b4f083d22 100644 --- a/src/specfem2D/setup_sources_receivers.F90 +++ b/src/specfem2D/setup_sources_receivers.F90 @@ -519,7 +519,7 @@ subroutine setup_receivers_check_acoustic() 'on the free surface because pressure is zero there') else write(IMAIN,*) '**********************************************************************' - write(IMAIN,*) '*** Warning: the acoustic receivers located on the free surface ***' + write(IMAIN,*) '*** Warning: acoustic receivers located on the free surface ***' write(IMAIN,*) '*** Warning: will record 0 !! ***' write(IMAIN,*) '**********************************************************************' write(IMAIN,*) diff --git a/src/specfem2D/specfem2D_par.f90 b/src/specfem2D/specfem2D_par.f90 index 0631357a2..5b3f86f6c 100644 --- a/src/specfem2D/specfem2D_par.f90 +++ b/src/specfem2D/specfem2D_par.f90 @@ -631,9 +631,7 @@ module specfem_par double precision, dimension(:), allocatable :: Qe11_electromagnetic,Qe33_electromagnetic double precision, dimension(:), allocatable :: Qs11_electromagnetic,Qs33_electromagnetic double precision, dimension(:,:,:,:), allocatable :: tau_e,tau_d,tau_s - double precision, dimension(:,:,:,:), allocatable :: alphavalem,betavalem,gammavalem - double precision, dimension(:), allocatable :: tauinv - + double precision, dimension(:,:,:,:), allocatable :: alphaval_em,betaval_em,gammaval_em double precision, dimension(:,:,:), allocatable :: rx_permattenuation,rz_permattenuation,permx,permz !for backward simulation in adjoint inversion From e524e89e0f4d5d0ef3a46a0079e5fbf47e5fdb96 Mon Sep 17 00:00:00 2001 From: Daniel Peter Date: Wed, 5 Jun 2024 17:16:13 +0200 Subject: [PATCH 9/9] adds github action tests (for Gmsh & EM) --- .github/workflows/CI.yml | 48 ++++++++++++++++++++++++++++++++++++++++ 1 file changed, 48 insertions(+) diff --git a/.github/workflows/CI.yml b/.github/workflows/CI.yml index 2be076522..07c9405d5 100644 --- a/.github/workflows/CI.yml +++ b/.github/workflows/CI.yml @@ -605,4 +605,52 @@ jobs: run: ./.github/scripts/run_tests.sh shell: bash + linuxTest_13: + name: Test run example 13 - Gmsh w/ Stacey + runs-on: ubuntu-latest + needs: [linuxCheck] + + steps: + - uses: actions/checkout@v4 + + - name: Install packages + run: ./.github/scripts/run_install.sh + shell: bash + + - name: Run build + env: + TESTFLAGS: --with-mpi + run: ./.github/scripts/run_build.sh + shell: bash + + - name: Run test + env: + TESTDIR: EXAMPLES/Gmsh_example_Stacey_MPI + run: ./.github/scripts/run_tests.sh + shell: bash + + linuxTest_14: + name: Test run example 14 - EM + runs-on: ubuntu-latest + needs: [linuxCheck] + + steps: + - uses: actions/checkout@v4 + + - name: Install packages + run: ./.github/scripts/run_install.sh + shell: bash + + - name: Run build + env: + TESTFLAGS: --enable-debug + run: ./.github/scripts/run_build.sh + shell: bash + + - name: Run test + env: + TESTDIR: EXAMPLES/electromagnetic_GPR_TM + DEBUG: true + run: ./.github/scripts/run_tests.sh + shell: bash