[cig-commits] [commit] devel, master: Use Paul's script to add in Par_file about how to choose time step (2f18e6d)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 15:21:53 PDT 2014
Repository : https://github.com/geodynamics/specfem2d
On branches: devel,master
Link : https://github.com/geodynamics/specfem2d/compare/fc67e6fd7ad890705b2b72b4b3c509accb22249e...e9ca46c40131588d89d7b0883250bc6584ce6b4c
>---------------------------------------------------------------
commit 2f18e6d2bf50b762bd4da8dd6618fba4b299c8fa
Author: Xie Zhinan <xiezhinan1984 at gmail.com>
Date: Mon Oct 22 09:04:34 2012 +0000
Use Paul's script to add in Par_file about how to choose time step
>---------------------------------------------------------------
2f18e6d2bf50b762bd4da8dd6618fba4b299c8fa
.../Par_file_Abel_Balanche_bathy_source_solid | 2 +-
Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in | 2 +-
Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in | 2 +-
INDUSTRIAL_FORMAT/Par_file | 2 +-
M2_UPPA/Par_file_M2_UPPA | 2 +-
Rayleigh_wave_no_crack/Par_file_Rayleigh_2D | 2 +-
Rayleigh_wave_with_crack/Par_file_Rayleigh_2D | 2 +-
Tape2007/Par_file_Tape2007_132rec_checker | 2 +-
Tape2007/Par_file_Tape2007_onerec | 2 +-
Tape2007_kernel/Par_file_Tape2007_onerec | 2 +-
Tromp2005/Par_file_Tromp2005 | 2 +-
Tromp2005/Par_file_Tromp2005_s100 | 2 +-
Tromp2005_kernel/Par_file_Tromp2005 | 2 +-
acoustic_poroelastic/Par_file_acoustic_poroelastic | 2 +-
attenuation/Par_file_attenuation_2D | 2 +-
canyon/Par_file_canyon | 2 +-
.../fluid_solid_external_mesh/Par_file_fluid_solid | 2 +-
.../Par_file_fluid_solid | 2 +-
.../Par_file_fluid_solid | 2 +-
init_plane/Par_file_Slave | 2 +-
init_plane/Par_file_Slave_for | 2 +-
init_plane/Par_file_Slave_kernel | 2 +-
noise_layered/model_0/Par_file_fair | 2 +-
noise_layered/model_0/Par_file_good | 2 +-
noise_layered/model_1/Par_file_best | 2 +-
noise_layered/model_1/Par_file_fair | 2 +-
noise_layered/model_1/Par_file_good | 2 +-
noise_layered/model_2/Par_file_fair | 2 +-
noise_layered/model_2/Par_file_good | 2 +-
noise_uniform/Par_file_noise_1 | 2 +-
noise_uniform/Par_file_noise_2 | 2 +-
noise_uniform/Par_file_noise_3 | 2 +-
...e_their_format_when_new_parameters_are_added.py | 28 ++++++++++++++++++++++
salt_dome/Stacey_homogeneous/Par_file | 2 +-
salt_dome/Stacey_normal/Par_file | 2 +-
semi_infinite_homo/Par_file_elastic_2D | 2 +-
36 files changed, 63 insertions(+), 35 deletions(-)
diff --git a/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid b/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
index bb9a460..7ae86f9 100644
--- a/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
+++ b/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
-deltat = 1.d-3 # duration of a time step
+deltat = 1.d-3 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in b/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
index ef108e6..065c9da 100644
--- a/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
+++ b/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
@@ -25,7 +25,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 2400 # total number of time steps
-deltat = 5.d-7 # duration of a time step
+deltat = 5.d-7 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in b/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
index 7c53c91..2b8002d 100644
--- a/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
+++ b/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
@@ -25,7 +25,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 2400 # total number of time steps
-deltat = 5.d-7 # duration of a time step
+deltat = 5.d-7 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/INDUSTRIAL_FORMAT/Par_file b/INDUSTRIAL_FORMAT/Par_file
index bc5bb94..3deec23 100644
--- a/INDUSTRIAL_FORMAT/Par_file
+++ b/INDUSTRIAL_FORMAT/Par_file
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 10000 # total number of time steps
-deltat = 5.d-4 # duration of a time step
+deltat = 5.d-4 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/M2_UPPA/Par_file_M2_UPPA b/M2_UPPA/Par_file_M2_UPPA
index cf64618..fe9b5d5 100644
--- a/M2_UPPA/Par_file_M2_UPPA
+++ b/M2_UPPA/Par_file_M2_UPPA
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 1600 # total number of time steps
-deltat = 1.d-3 # duration of a time step
+deltat = 1.d-3 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D b/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
index b3251a9..7e52b96 100644
--- a/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
+++ b/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 5000 # total number of time steps
-deltat = 3.90625e-3 # duration of a time step
+deltat = 3.90625e-3 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D b/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
index 6ecd0c6..f804dec 100644
--- a/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
+++ b/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 5000 # total number of time steps
-deltat = 3.90625e-3 # duration of a time step
+deltat = 3.90625e-3 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/Tape2007/Par_file_Tape2007_132rec_checker b/Tape2007/Par_file_Tape2007_132rec_checker
index 07531d8..646fa3d 100644
--- a/Tape2007/Par_file_Tape2007_132rec_checker
+++ b/Tape2007/Par_file_Tape2007_132rec_checker
@@ -24,7 +24,7 @@ p_sv = .false. # set the type of calculation (
# time step parameters
nt = 4800 # total number of time steps
-deltat = 6.0d-2 # duration of a time step
+deltat = 6.0d-2 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
USER_T0 = 48.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
diff --git a/Tape2007/Par_file_Tape2007_onerec b/Tape2007/Par_file_Tape2007_onerec
index d94cc74..2faaeeb 100644
--- a/Tape2007/Par_file_Tape2007_onerec
+++ b/Tape2007/Par_file_Tape2007_onerec
@@ -24,7 +24,7 @@ p_sv = .false. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
-deltat = 6.0d-2 # duration of a time step
+deltat = 6.0d-2 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
USER_T0 = 48.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
diff --git a/Tape2007_kernel/Par_file_Tape2007_onerec b/Tape2007_kernel/Par_file_Tape2007_onerec
index 4979d59..3bab0c1 100644
--- a/Tape2007_kernel/Par_file_Tape2007_onerec
+++ b/Tape2007_kernel/Par_file_Tape2007_onerec
@@ -24,7 +24,7 @@ p_sv = .false. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
-deltat = 6.0d-2 # duration of a time step
+deltat = 6.0d-2 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/Tromp2005/Par_file_Tromp2005 b/Tromp2005/Par_file_Tromp2005
index adeef85..f0df889 100644
--- a/Tromp2005/Par_file_Tromp2005
+++ b/Tromp2005/Par_file_Tromp2005
@@ -24,7 +24,7 @@ p_sv = .false. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
-deltat = 2.0d-2 # duration of a time step
+deltat = 2.0d-2 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
USER_T0 = 8.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
diff --git a/Tromp2005/Par_file_Tromp2005_s100 b/Tromp2005/Par_file_Tromp2005_s100
index d906477..01c8bbd 100644
--- a/Tromp2005/Par_file_Tromp2005_s100
+++ b/Tromp2005/Par_file_Tromp2005_s100
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
-deltat = 2.0d-4 # duration of a time step
+deltat = 2.0d-4 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
USER_T0 = 8.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
diff --git a/Tromp2005_kernel/Par_file_Tromp2005 b/Tromp2005_kernel/Par_file_Tromp2005
index 617825a..0c654cd 100644
--- a/Tromp2005_kernel/Par_file_Tromp2005
+++ b/Tromp2005_kernel/Par_file_Tromp2005
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
-deltat = 2.0d-2 # duration of a time step
+deltat = 2.0d-2 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/acoustic_poroelastic/Par_file_acoustic_poroelastic b/acoustic_poroelastic/Par_file_acoustic_poroelastic
index 463b040..bd7f008 100644
--- a/acoustic_poroelastic/Par_file_acoustic_poroelastic
+++ b/acoustic_poroelastic/Par_file_acoustic_poroelastic
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 2500 # total number of time steps
-deltat = 4.d-4 # duration of a time step
+deltat = 4.d-4 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/attenuation/Par_file_attenuation_2D b/attenuation/Par_file_attenuation_2D
index f35e3ae..7451b23 100644
--- a/attenuation/Par_file_attenuation_2D
+++ b/attenuation/Par_file_attenuation_2D
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 1400 # total number of time steps
-deltat = 1e-3 # duration of a time step
+deltat = 1e-3 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/canyon/Par_file_canyon b/canyon/Par_file_canyon
index be7e6f3..d656c2b 100644
--- a/canyon/Par_file_canyon
+++ b/canyon/Par_file_canyon
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 5000 # total number of time steps
-deltat = 2.e-3 # duration of a time step
+deltat = 2.e-3 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid b/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
index 058dddd..1bda086 100644
--- a/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
+++ b/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 18000 # total number of time steps
-deltat = 0.25d-3 # duration of a time step
+deltat = 0.25d-3 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid b/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid
index f62ca48..44f4c29 100644
--- a/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid
+++ b/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 4500 # total number of time steps
-deltat = 0.85d-3 # duration of a time step
+deltat = 0.85d-3 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
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 3c4123d..1683ab0 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
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 2000 # total number of time steps
-deltat = 0.7d-3 # duration of a time step
+deltat = 0.7d-3 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/init_plane/Par_file_Slave b/init_plane/Par_file_Slave
index 8a62491..e4c9abb 100644
--- a/init_plane/Par_file_Slave
+++ b/init_plane/Par_file_Slave
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 6000 # total number of time steps
-deltat = 0.01 # duration of a time step
+deltat = 0.01 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/init_plane/Par_file_Slave_for b/init_plane/Par_file_Slave_for
index dde9f7b..6bd7da8 100644
--- a/init_plane/Par_file_Slave_for
+++ b/init_plane/Par_file_Slave_for
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 6000 # total number of time steps
-deltat = 0.01 # duration of a time step
+deltat = 0.01 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/init_plane/Par_file_Slave_kernel b/init_plane/Par_file_Slave_kernel
index 5adf56e..14ec765 100644
--- a/init_plane/Par_file_Slave_kernel
+++ b/init_plane/Par_file_Slave_kernel
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 6000 # total number of time steps
-deltat = 0.01 # duration of a time step
+deltat = 0.01 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/noise_layered/model_0/Par_file_fair b/noise_layered/model_0/Par_file_fair
index 4c8742d..11ddef0 100644
--- a/noise_layered/model_0/Par_file_fair
+++ b/noise_layered/model_0/Par_file_fair
@@ -24,7 +24,7 @@ p_sv = .false. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
-deltat = 1.d-3 # duration of a time step
+deltat = 1.d-3 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/noise_layered/model_0/Par_file_good b/noise_layered/model_0/Par_file_good
index 254df64..57afc70 100644
--- a/noise_layered/model_0/Par_file_good
+++ b/noise_layered/model_0/Par_file_good
@@ -24,7 +24,7 @@ p_sv = .false. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
-deltat = 1.d-3 # duration of a time step
+deltat = 1.d-3 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/noise_layered/model_1/Par_file_best b/noise_layered/model_1/Par_file_best
index 09b2e98..2f3fcd4 100644
--- a/noise_layered/model_1/Par_file_best
+++ b/noise_layered/model_1/Par_file_best
@@ -24,7 +24,7 @@ p_sv = .false. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
-deltat = 1.d-3 # duration of a time step
+deltat = 1.d-3 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/noise_layered/model_1/Par_file_fair b/noise_layered/model_1/Par_file_fair
index 2cc307f..65f6178 100644
--- a/noise_layered/model_1/Par_file_fair
+++ b/noise_layered/model_1/Par_file_fair
@@ -24,7 +24,7 @@ p_sv = .false. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
-deltat = 1.d-3 # duration of a time step
+deltat = 1.d-3 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/noise_layered/model_1/Par_file_good b/noise_layered/model_1/Par_file_good
index d7f63ed..fe10e8e 100644
--- a/noise_layered/model_1/Par_file_good
+++ b/noise_layered/model_1/Par_file_good
@@ -24,7 +24,7 @@ p_sv = .false. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
-deltat = 1.d-3 # duration of a time step
+deltat = 1.d-3 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/noise_layered/model_2/Par_file_fair b/noise_layered/model_2/Par_file_fair
index 2c0ae75..e96dd9c 100644
--- a/noise_layered/model_2/Par_file_fair
+++ b/noise_layered/model_2/Par_file_fair
@@ -24,7 +24,7 @@ p_sv = .false. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
-deltat = 1.d-3 # duration of a time step
+deltat = 1.d-3 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/noise_layered/model_2/Par_file_good b/noise_layered/model_2/Par_file_good
index c85fb7f..52d2f52 100644
--- a/noise_layered/model_2/Par_file_good
+++ b/noise_layered/model_2/Par_file_good
@@ -24,7 +24,7 @@ p_sv = .false. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
-deltat = 1.d-3 # duration of a time step
+deltat = 1.d-3 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/noise_uniform/Par_file_noise_1 b/noise_uniform/Par_file_noise_1
index 19a538b..8707b10 100644
--- a/noise_uniform/Par_file_noise_1
+++ b/noise_uniform/Par_file_noise_1
@@ -24,7 +24,7 @@ p_sv = .false. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
-deltat = 4.d-2 # duration of a time step
+deltat = 4.d-2 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/noise_uniform/Par_file_noise_2 b/noise_uniform/Par_file_noise_2
index d501b6d..5d52758 100644
--- a/noise_uniform/Par_file_noise_2
+++ b/noise_uniform/Par_file_noise_2
@@ -24,7 +24,7 @@ p_sv = .false. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
-deltat = 4.d-2 # duration of a time step
+deltat = 4.d-2 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/noise_uniform/Par_file_noise_3 b/noise_uniform/Par_file_noise_3
index 317c44e..c6f692f 100644
--- a/noise_uniform/Par_file_noise_3
+++ b/noise_uniform/Par_file_noise_3
@@ -24,7 +24,7 @@ p_sv = .false. # set the type of calculation (
# time step parameters
nt = 3000 # total number of time steps
-deltat = 4.d-2 # duration of a time step
+deltat = 4.d-2 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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
diff --git a/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py b/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py
index 1d0917a..2930379 100755
--- a/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py
+++ b/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py
@@ -444,6 +444,33 @@ def ProcessParfile_r20561(fic):
print 'xxxxx------> '+fic+' processed to r20561'
return
#------------------------------------------------------------------------------
+def ProcessParfile_r20859(fic):
+ # Open the file and get all lines from Par_file
+ ligs= LoadLig(fic)
+
+ # Test if already processed
+ for lig in ligs:
+ if lig.endswith('section 4.5 of the user manual)\n'):
+ print '----> '+fic+' already processed to r20859'
+ return
+ #--------------------------------------------------------------------------
+ # Add new parameters
+ #
+ for ilg, lig in enumerate(ligs):
+ if lig.startswith('deltat'):
+ ligs[ilg] = ligs[ilg].replace('# duration of a time step','#'+ \
+ ' duration of a time step (for the choice of deltat please refer'+ \
+ ' to section 4.5 of the user manual)')
+ #
+ move(fic,fic+'.before_update_to_r20859')
+ #
+ fm = open(fic,'w')
+ fm.writelines(ligs)
+ fm.close()
+ #
+ print 'xxxxx------> '+fic+' processed to r20859'
+ return
+#------------------------------------------------------------------------------
if __name__=='__main__':
## List of all files of current directory
Fichiers=[]
@@ -468,6 +495,7 @@ if __name__=='__main__':
ProcessParfile_r20307(fic)
ProcessParfile_r20359(fic)
ProcessParfile_r20561(fic)
+ ProcessParfile_r20859(fic)
print '~'*80
#
print 'Number of Par_file analysed : ', Ct_Par_file
diff --git a/salt_dome/Stacey_homogeneous/Par_file b/salt_dome/Stacey_homogeneous/Par_file
index 7b55a1f..5a07307 100644
--- a/salt_dome/Stacey_homogeneous/Par_file
+++ b/salt_dome/Stacey_homogeneous/Par_file
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 10000 # total number of time steps
-deltat = 4.25e-4 # duration of a time step
+deltat = 4.25e-4 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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 RK4 4th-order 4-stage Runge-Kutta
diff --git a/salt_dome/Stacey_normal/Par_file b/salt_dome/Stacey_normal/Par_file
index 3f45770..02228c1 100644
--- a/salt_dome/Stacey_normal/Par_file
+++ b/salt_dome/Stacey_normal/Par_file
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 10000 # total number of time steps
-deltat = 4.25e-4 # duration of a time step
+deltat = 4.25e-4 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
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 RK4 4th-order 4-stage Runge-Kutta
diff --git a/semi_infinite_homo/Par_file_elastic_2D b/semi_infinite_homo/Par_file_elastic_2D
index afebc0f..e5df40b 100644
--- a/semi_infinite_homo/Par_file_elastic_2D
+++ b/semi_infinite_homo/Par_file_elastic_2D
@@ -24,7 +24,7 @@ p_sv = .true. # set the type of calculation (
# time step parameters
nt = 15000 # total number of time steps
-deltat = 0.8e-3 # duration of a time step
+deltat = 0.8e-3 # duration of a time step (for the choice of deltat please refer to section 4.5 of the user manual)
USER_T0 = 0.15d0 # 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
More information about the CIG-COMMITS
mailing list