diff --git a/.github/workflows/CI.yml b/.github/workflows/CI.yml index e31803dfd..07c9405d5 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: | @@ -104,9 +104,16 @@ 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 + # 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 @@ -139,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 @@ -158,7 +165,7 @@ jobs: needs: changesCheck steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Cache Intel oneapi packages id: cache-intel-oneapi @@ -294,7 +301,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -312,7 +319,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -336,7 +343,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -361,7 +368,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -386,7 +393,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -410,7 +417,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -434,7 +441,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -459,7 +466,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -483,7 +490,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -507,7 +514,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -531,7 +538,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -556,7 +563,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -580,7 +587,7 @@ jobs: needs: [linuxCheck] steps: - - uses: actions/checkout@v3 + - uses: actions/checkout@v4 - name: Install packages run: ./.github/scripts/run_install.sh @@ -598,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 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..83766b342 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 @@ -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/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..52585e4b4 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 @@ -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 29075f783..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 @@ -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 @@ -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 37be2f307..52585e4b4 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 @@ -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/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..50a04c776 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 @@ -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 4bb3c32a5..3ea5fae6c 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 @@ -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_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..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 @@ -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 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..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 @@ -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 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..684536278 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 @@ -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/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..1a64ece9e 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 @@ -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 f7d90509c..46f196658 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 @@ -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 dd3cb7e6e..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 @@ -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 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..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 @@ -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 @@ -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/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..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 @@ -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 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..d36ff3190 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 @@ -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/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..2b270cdfa 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 @@ -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/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..393cbf417 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 @@ -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 54e5e4c26..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 @@ -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 @@ -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 987b8e11d..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 @@ -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 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..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 @@ -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 @@ -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 109ca69dd..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 @@ -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 @@ -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 d42c49c6b..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 @@ -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 @@ -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 83f4bc422..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 @@ -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 @@ -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 94b5c7166..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 @@ -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 @@ -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_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..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 @@ -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_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..a99607255 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 @@ -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 6cbd57fbd..e248857c0 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 @@ -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 6cbd57fbd..e248857c0 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 @@ -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/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/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..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 @@ -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 @@ -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/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_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/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/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 566b9221f..5b3f86f6c 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,26 +614,24 @@ 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 - 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 @@ -703,7 +701,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 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). # 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) """