[cig-commits] [commit] devel, master: RK and LDDRK time schemes not supported for adjoint inversion, mixed elastic/poroelastic forword simulations and poroelastic forword simulations with attenuation (2eb7d1f)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 15:21:00 PDT 2014
Repository : https://github.com/geodynamics/specfem2d
On branches: devel,master
Link : https://github.com/geodynamics/specfem2d/compare/fc67e6fd7ad890705b2b72b4b3c509accb22249e...e9ca46c40131588d89d7b0883250bc6584ce6b4c
>---------------------------------------------------------------
commit 2eb7d1f75f466afcc959bed29fcbf6629ec7a7ee
Author: Xie Zhinan <xiezhinan1984 at gmail.com>
Date: Tue Jan 31 12:44:59 2012 +0000
RK and LDDRK time schemes not supported for
adjoint inversion, mixed elastic/poroelastic forword simulations
and poroelastic forword simulations with attenuation
>---------------------------------------------------------------
2eb7d1f75f466afcc959bed29fcbf6629ec7a7ee
acoustic_poroelastic/Par_file_acoustic_poroelastic | 16 ++++++++--------
acoustic_poroelastic/interfaces_acoustic_poroelastic.dat | 4 ++--
attenuation/Par_file_attenuation_2D | 6 +++---
.../Par_file_fluid_solid | 6 +++---
semi_infinite_homo/Par_file_elastic_2D | 4 ++--
5 files changed, 18 insertions(+), 18 deletions(-)
diff --git a/acoustic_poroelastic/Par_file_acoustic_poroelastic b/acoustic_poroelastic/Par_file_acoustic_poroelastic
index afcdd00..9a45f50 100644
--- a/acoustic_poroelastic/Par_file_acoustic_poroelastic
+++ b/acoustic_poroelastic/Par_file_acoustic_poroelastic
@@ -19,12 +19,12 @@ READ_EXTERNAL_SEP_FILE = .false. # Read external SEP file from D
ATTENUATION_VISCOELASTIC_SOLID = .false. # turn attenuation (viscoelasticity) on or off for non-poroelastic solid parts of the model
ATTENUATION_PORO_FLUID_PART = .false. # turn viscous attenuation on or off for the fluid part of poroelastic parts of the model
Q0 = 1 # quality factor for viscous attenuation
-freq0 = 1000 # frequency for viscous attenuation
+freq0 = 1000 # frequency for viscous attenuation
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
# time step parameters
-nt = 5000 # total number of time steps
-deltat = 1.2d-4 # duration of a time step
+nt = 2500 # total number of time steps
+deltat = 4.d-4 # duration of a time step
USER_T0 = 0.0d0 # use this t0 as earliest starting time rather than the automatically calculated one
time_stepping_scheme = 2 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical 4th-order 4-stage Runge-Kutta
@@ -91,8 +91,8 @@ nbmodels = 2 # nb of different models
# III: (model_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu).
# For istropic elastic/acoustic material use I and set Vs to zero to make a given model acoustic, for anisotropic elastic use II,
# and for isotropic poroelastic material use III. The mesh can contain acoustic, elastic, and poroelastic models simultaneously.
-1 3 2650.d0 880.d0 0.3d0 1.2 1d-9 0.d0 1d-9 1.22d10 1.985d9 9.6d9 10d-4 0.0001 10000
-2 1 2500.d0 3400.d0 1963.d0 0 0 9999 9999 0 0 0 0 0 0
+1 3 2650.d0 880.d0 0.3d0 1.2 1d-9 0.d0 1d-9 1.22d10 1.985d9 9.6d9 10d-4 0.0001 9999999999
+2 1 1020.d0 1500.d0 0.d0 0 0 10.d0 10.d0 0 0 0 0 0 0
# external mesh or not
read_external_mesh = .false.
@@ -130,7 +130,7 @@ interfacesfile = ../interfaces_acoustic_poroelastic.dat
# geometry of the model (origin lower-left corner = 0,0) and mesh description
xmin = 0.d0 # abscissa of left side of the model
xmax = 4800.d0 # abscissa of right side of the model
-nx = 260 # number of elements along X
+nx = 130 # number of elements along X
# absorbing boundary parameters (see absorbing_conditions above)
absorbbottom = .true.
@@ -140,5 +140,5 @@ absorbleft = .true.
# define the different regions of the model in the (nx,nz) spectral element mesh
nbregions = 2 # nb of regions and model number for each
-1 260 1 110 1
-1 260 111 220 2
+1 130 1 55 1
+1 130 56 110 2
diff --git a/acoustic_poroelastic/interfaces_acoustic_poroelastic.dat b/acoustic_poroelastic/interfaces_acoustic_poroelastic.dat
index c42ba85..cea3ac1 100644
--- a/acoustic_poroelastic/interfaces_acoustic_poroelastic.dat
+++ b/acoustic_poroelastic/interfaces_acoustic_poroelastic.dat
@@ -29,6 +29,6 @@
#
# layer number 1 (bottom layer)
#
- 110
+ 55
#
- 110
+ 55
diff --git a/attenuation/Par_file_attenuation_2D b/attenuation/Par_file_attenuation_2D
index 23ae87b..0ad3adb 100644
--- a/attenuation/Par_file_attenuation_2D
+++ b/attenuation/Par_file_attenuation_2D
@@ -23,10 +23,10 @@ freq0 = 10 # frequency for viscous attenua
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
# time step parameters
-nt = 3000 # total number of time steps
-deltat = 5e-4 # duration of a time step
+nt = 1400 # total number of time steps
+deltat = 1e-3 # duration of a time step
USER_T0 = 0.0d0 # use this t0 as earliest starting time rather than the automatically calculated one
-time_stepping_scheme = 1 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical 4th-order 4-stage Runge-Kutta
+time_stepping_scheme = 2 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical 4th-order 4-stage Runge-Kutta
# source parameters
NSOURCES = 1 # number of sources [source info read in CMTSOLUTION file]
diff --git a/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid b/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid
index db10342..11077e6 100644
--- a/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid
+++ b/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid
@@ -23,10 +23,10 @@ freq0 = 10 # frequency for viscous attenua
p_sv = .true. # set the type of calculation (P-SV or SH/membrane waves)
# time step parameters
-nt = 5000 # total number of time steps
-deltat = 0.55d-3 # duration of a time step
+nt = 2000 # total number of time steps
+deltat = 0.7d-3 # duration of a time step
USER_T0 = 0.0d0 # use this t0 as earliest starting time rather than the automatically calculated one
-time_stepping_scheme = 1 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical 4th-order 4-stage Runge-Kutta
+time_stepping_scheme = 2 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical 4th-order 4-stage Runge-Kutta
# source parameters
NSOURCES = 1 # number of sources [source info read in CMTSOLUTION file]
diff --git a/semi_infinite_homo/Par_file_elastic_2D b/semi_infinite_homo/Par_file_elastic_2D
index e33bd55..cd837cb 100644
--- a/semi_infinite_homo/Par_file_elastic_2D
+++ b/semi_infinite_homo/Par_file_elastic_2D
@@ -24,9 +24,9 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 1500 # total number of time steps
-deltat = 1.9e-3 # duration of a time step
+deltat = 0.8e-3 # duration of a time step
USER_T0 = 0.15d0 # use this t0 as earliest starting time rather than the automatically calculated one
-time_stepping_scheme = 3 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical 4th-order 4-stage Runge-Kutta
+time_stepping_scheme = 1 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical 4th-order 4-stage Runge-Kutta
# source parameters
NSOURCES = 1 # number of sources [source info read in CMTSOLUTION file]
More information about the CIG-COMMITS
mailing list