[cig-commits] [commit] devel, master: replaced the OCEANS flag with APPROXIMATE_OCEAN_LOAD to avoid confusion (1364e08)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jun 18 16:50:58 PDT 2014
Repository : https://github.com/geodynamics/specfem3d
On branches: devel,master
Link : https://github.com/geodynamics/specfem3d/compare/6026e367984905ab133865f62fa6293b343759b9...47f703851338234f96397e7da9fbff63d8178b8a
>---------------------------------------------------------------
commit 1364e08c38e4ecfe7099af157811b8bd392f9eb4
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date: Wed Mar 13 00:23:42 2013 +0000
replaced the OCEANS flag with APPROXIMATE_OCEAN_LOAD to avoid confusion
>---------------------------------------------------------------
1364e08c38e4ecfe7099af157811b8bd392f9eb4
BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file | 6 ++--
BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file | 6 ++--
CPML/HOMO8_FREE_SURFACE/DATA/Par_file | 6 ++--
CPML/HOMO8_NOFREESURFACE/DATA/Par_file | 6 ++--
Mount_StHelens/DATA/Par_file | 6 ++--
homogeneous_halfspace_HEX27/DATA/Par_file | 6 ++--
homogeneous_halfspace_HEX8/DATA/Par_file | 6 ++--
homogeneous_poroelastic/DATA/Par_file | 6 ++--
layered_halfspace/DATA/Par_file | 6 ++--
meshfem3D_examples/many_interfaces/DATA/Par_file | 6 ++--
meshfem3D_examples/simple_model/DATA/Par_file | 6 ++--
meshfem3D_examples/socal1D/DATA/Par_file | 6 ++--
.../socal1D/example_utm/Par_file_utm | 2 +-
noise_tomography/DATA/Par_file_step1 | 2 +-
noise_tomography/DATA/Par_file_step2 | 2 +-
noise_tomography/DATA/Par_file_step3 | 2 +-
splay_faults/DATA/Par_file | 6 ++--
splay_faults/DATA/Par_file_faults.barrier | 34 +++++++++++-----------
splay_faults/DATA/Par_file_faults.homogeneous | 34 +++++++++++-----------
tomographic_model/DATA/Par_file | 6 ++--
tpv102/DATA/Par_file | 6 ++--
tpv102/DATA/Par_file_faults | 2 +-
tpv103/DATA/Par_file | 6 ++--
tpv103/DATA/Par_file_faults | 2 +-
tpv15/DATA/Par_file | 6 ++--
tpv15/DATA/Par_file_faults | 14 ++++-----
tpv16/DATA/Par_file | 6 ++--
tpv16/DATA/Par_file_faults | 2 +-
tpv5/DATA/Par_file | 6 ++--
tpv5/DATA/Par_file_faults | 4 +--
waterlayered_halfspace/DATA/Par_file | 6 ++--
31 files changed, 110 insertions(+), 110 deletions(-)
diff --git a/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file b/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file
index 444956a..509c3ef 100644
--- a/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file
+++ b/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.001
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .true.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file b/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file
index 0c0a648..5fcb60b 100644
--- a/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file
+++ b/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.0006
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .true.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/CPML/HOMO8_FREE_SURFACE/DATA/Par_file b/CPML/HOMO8_FREE_SURFACE/DATA/Par_file
index abf71c0..ac6128c 100644
--- a/CPML/HOMO8_FREE_SURFACE/DATA/Par_file
+++ b/CPML/HOMO8_FREE_SURFACE/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.05
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .false.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .true.
+PML_CONDITIONS = .true.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/CPML/HOMO8_NOFREESURFACE/DATA/Par_file b/CPML/HOMO8_NOFREESURFACE/DATA/Par_file
index c311c67..0f5de5a 100644
--- a/CPML/HOMO8_NOFREESURFACE/DATA/Par_file
+++ b/CPML/HOMO8_NOFREESURFACE/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.05
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .false.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .true.
+PML_CONDITIONS = .true.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .true.
diff --git a/Mount_StHelens/DATA/Par_file b/Mount_StHelens/DATA/Par_file
index 7be25fa..57b2d69 100644
--- a/Mount_StHelens/DATA/Par_file
+++ b/Mount_StHelens/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.005d0
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .true.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/homogeneous_halfspace_HEX27/DATA/Par_file b/homogeneous_halfspace_HEX27/DATA/Par_file
index a2ff73e..4997d08 100644
--- a/homogeneous_halfspace_HEX27/DATA/Par_file
+++ b/homogeneous_halfspace_HEX27/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.05
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 27
+NGNOD = 27
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 27
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .false.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/homogeneous_halfspace_HEX8/DATA/Par_file b/homogeneous_halfspace_HEX8/DATA/Par_file
index 943b0cd..3221264 100644
--- a/homogeneous_halfspace_HEX8/DATA/Par_file
+++ b/homogeneous_halfspace_HEX8/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.05
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .false.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/homogeneous_poroelastic/DATA/Par_file b/homogeneous_poroelastic/DATA/Par_file
index 22b663c..df02935 100644
--- a/homogeneous_poroelastic/DATA/Par_file
+++ b/homogeneous_poroelastic/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 5.d-5 #0.5d-3
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .false.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/layered_halfspace/DATA/Par_file b/layered_halfspace/DATA/Par_file
index 2a882f5..35bf60f 100644
--- a/layered_halfspace/DATA/Par_file
+++ b/layered_halfspace/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.0075
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .true.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/meshfem3D_examples/many_interfaces/DATA/Par_file b/meshfem3D_examples/many_interfaces/DATA/Par_file
index f2cf55b..0b0365e 100644
--- a/meshfem3D_examples/many_interfaces/DATA/Par_file
+++ b/meshfem3D_examples/many_interfaces/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.03
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .true.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/meshfem3D_examples/simple_model/DATA/Par_file b/meshfem3D_examples/simple_model/DATA/Par_file
index e14db13..959ef09 100644
--- a/meshfem3D_examples/simple_model/DATA/Par_file
+++ b/meshfem3D_examples/simple_model/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.03
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .false.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/meshfem3D_examples/socal1D/DATA/Par_file b/meshfem3D_examples/socal1D/DATA/Par_file
index 581fe54..01ffb0f 100644
--- a/meshfem3D_examples/socal1D/DATA/Par_file
+++ b/meshfem3D_examples/socal1D/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.012
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .true.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/meshfem3D_examples/socal1D/example_utm/Par_file_utm b/meshfem3D_examples/socal1D/example_utm/Par_file_utm
index f3b88a4..022193d 100644
--- a/meshfem3D_examples/socal1D/example_utm/Par_file_utm
+++ b/meshfem3D_examples/socal1D/example_utm/Par_file_utm
@@ -28,7 +28,7 @@ DT = 0.012
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
USE_OLSEN_ATTENUATION = .false.
diff --git a/noise_tomography/DATA/Par_file_step1 b/noise_tomography/DATA/Par_file_step1
index 673018c..8c1a1d6 100644
--- a/noise_tomography/DATA/Par_file_step1
+++ b/noise_tomography/DATA/Par_file_step1
@@ -29,7 +29,7 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
USE_OLSEN_ATTENUATION = .false.
diff --git a/noise_tomography/DATA/Par_file_step2 b/noise_tomography/DATA/Par_file_step2
index 30f66a8..cc4e063 100644
--- a/noise_tomography/DATA/Par_file_step2
+++ b/noise_tomography/DATA/Par_file_step2
@@ -29,7 +29,7 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
USE_OLSEN_ATTENUATION = .false.
diff --git a/noise_tomography/DATA/Par_file_step3 b/noise_tomography/DATA/Par_file_step3
index f5e5df8..62851f5 100644
--- a/noise_tomography/DATA/Par_file_step3
+++ b/noise_tomography/DATA/Par_file_step3
@@ -29,7 +29,7 @@ NTSTEP_BETWEEN_READ_ADJSRC = 0
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
USE_OLSEN_ATTENUATION = .false.
diff --git a/splay_faults/DATA/Par_file b/splay_faults/DATA/Par_file
index 6e3117a..a7562a1 100644
--- a/splay_faults/DATA/Par_file
+++ b/splay_faults/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.00325d0
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .true.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/splay_faults/DATA/Par_file_faults.barrier b/splay_faults/DATA/Par_file_faults.barrier
index 4e1e0fd..7641477 100644
--- a/splay_faults/DATA/Par_file_faults.barrier
+++ b/splay_faults/DATA/Par_file_faults.barrier
@@ -2,7 +2,7 @@
0.00134 ! #tag #eta(damping) FAULT A
0.00134 ! #tag #eta(damping) FAULT BC
1 ! 1 = dyn 2=kin
-100 ! NTOUT : Number of time steps
+100 ! NTOUT : Number of time steps
1000 ! NTSNAP: time interation of snapshots
-1e0 ! V_HEALING (-1 : Healing off)
1e-3 ! V_RUPT
@@ -10,28 +10,28 @@
&BEGIN_FAULT /
&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-4.78e6,n1=0,n2=3,n3=3 /
-&DIST2D shape='rectangle', val = -2.16e6, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = -2.16e6, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
&DIST2D shape='rectangle', val = -2.88e6, xc = 166521.863e0, yc =148.0e3, zc= -7962.09773e0,lx=23794.68e0, ly=24.0e3, lz=3132.628e0 /
-&DIST2D shape='rectangle-taper', val = -2.16e6, valh =-0.216e6,xc = 226407.564e0, yc =120.0e3, zc= -4.50e3,lx=22787.2623e0, ly=80.0e3, lz=3.0e3 /
-&DIST2D shape='rectangle', val =-0.478e6, xc = 113.5e3, yc =120.0e3, zc= -15000.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle', val = -4.78e6, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle-taper', val = -4.78e6, valh=-0.478e6, xc = 226407.564e0, yc =120.0e3, zc= -4.50e3, lx=22787.2623e0, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle-taper', val = -2.16e6, valh =-0.216e6,xc = 226407.564e0, yc =120.0e3, zc= -4.50e3,lx=22787.2623e0, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val =-0.478e6, xc = 113.5e3, yc =120.0e3, zc= -15000.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = -4.78e6, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle-taper', val = -4.78e6, valh=-0.478e6, xc = 226407.564e0, yc =120.0e3, zc= -4.50e3, lx=22787.2623e0, ly=80.0e3, lz=3.0e3 /
&SWF mus=1000.0e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1,ndc=0 /
-&DIST2D shape='rectangle', val = 0.6e0, xc = 113.5e3, yc =135.0e3, zc= -15.0e3, lx=227.0e3, ly=50.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle', val = 0.2e0, xc = 113.5e3, yc =135.0e3, zc= -15.0e3, lx=227.0e3, ly=50.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = 0.6e0, xc = 113.5e3, yc =135.0e3, zc= -15.0e3, lx=227.0e3, ly=50.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = 0.2e0, xc = 113.5e3, yc =135.0e3, zc= -15.0e3, lx=227.0e3, ly=50.0e3, lz=30000.0e0 /
&BEGIN_FAULT /
&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-100.0e6,n1=0,n2=3,n3=4 /
-&DIST2D shape='rectangle', val = -5.69e6, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
-&DIST2D shape='rectangle', val = -3.54e6, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
-&DIST2D shape='rectangle-taper', val = -5.69e6, valh=-0.569e6, xc = 82205.7714e0, yc =120.0e3, zc= -4.5e3, lx=5196.15242e0, ly=80.0e3, lz=3.0e3 /
-&DIST2D shape='rectangle', val =-1.238e6, xc = 74411.5427e0, yc =120.0e3, zc= -9.0e3, lx=31176.9145e0, ly=80.0e3, lz=18.0e3 /
-&DIST2D shape='rectangle', val =-12.38e6, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
-&DIST2D shape='rectangle-taper', val =-12.38e6, valh=-1.238e6, xc = 82205.7714e0, yc =120.0e3, zc= -4.5e3, lx=5196.15242e0, ly=80.0e3, lz=3.0e3 /
-&DIST2D shape='rectangle', val = -9.04e6, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
+&DIST2D shape='rectangle', val = -5.69e6, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
+&DIST2D shape='rectangle', val = -3.54e6, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
+&DIST2D shape='rectangle-taper', val = -5.69e6, valh=-0.569e6, xc = 82205.7714e0, yc =120.0e3, zc= -4.5e3, lx=5196.15242e0, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val =-1.238e6, xc = 74411.5427e0, yc =120.0e3, zc= -9.0e3, lx=31176.9145e0, ly=80.0e3, lz=18.0e3 /
+&DIST2D shape='rectangle', val =-12.38e6, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
+&DIST2D shape='rectangle-taper', val =-12.38e6, valh=-1.238e6, xc = 82205.7714e0, yc =120.0e3, zc= -4.5e3, lx=5196.15242e0, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val = -9.04e6, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
&SWF mus=1000.e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1,ndc=0 /
-&DIST2D shape='rectangle', val = 0.6e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
-&DIST2D shape='rectangle', val = 0.2e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
+&DIST2D shape='rectangle', val = 0.6e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
+&DIST2D shape='rectangle', val = 0.2e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
diff --git a/splay_faults/DATA/Par_file_faults.homogeneous b/splay_faults/DATA/Par_file_faults.homogeneous
index bb1678c..94df165 100644
--- a/splay_faults/DATA/Par_file_faults.homogeneous
+++ b/splay_faults/DATA/Par_file_faults.homogeneous
@@ -2,7 +2,7 @@
0.00134 ! #tag #eta(damping) FAULT A
0.00134 ! #tag #eta(damping) FAULT BC
1 ! 1 = dyn 2=kin
-100 ! NTOUT : Number of time steps
+100 ! NTOUT : Number of time steps
1000 ! NTSNAP: time interation of snapshots
-1e0 ! V_HEALING (-1 : Healing off)
1e-3 ! V_RUPT
@@ -10,28 +10,28 @@
&BEGIN_FAULT /
&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-4.78e6,n1=0,n2=3,n3=3 /
-&DIST2D shape='rectangle', val = -2.16e6, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = -2.16e6, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
&DIST2D shape='rectangle', val = -2.88e6, xc = 166521.863e0, yc =148.0e3, zc= -7962.09773e0,lx=23794.68e0, ly=24.0e3, lz=3132.628e0 /
-&DIST2D shape='rectangle-taper', val = -2.16e6, valh =-0.216e6,xc = 226407.564e0, yc =120.0e3, zc= -4.50e3,lx=22787.2623e0, ly=80.0e3, lz=3.0e3 /
-&DIST2D shape='rectangle', val =-0.478e6, xc = 113.5e3, yc =120.0e3, zc= -15000.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle', val = -4.78e6, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle-taper', val = -4.78e6, valh=-0.478e6, xc = 226407.564e0, yc =120.0e3, zc= -4.50e3, lx=22787.2623e0, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle-taper', val = -2.16e6, valh =-0.216e6,xc = 226407.564e0, yc =120.0e3, zc= -4.50e3,lx=22787.2623e0, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val =-0.478e6, xc = 113.5e3, yc =120.0e3, zc= -15000.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = -4.78e6, xc = 113.5e3, yc =120.0e3, zc= -19500.0e0, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle-taper', val = -4.78e6, valh=-0.478e6, xc = 226407.564e0, yc =120.0e3, zc= -4.50e3, lx=22787.2623e0, ly=80.0e3, lz=3.0e3 /
&SWF mus=1000.0e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1,ndc=0 /
-&DIST2D shape='rectangle', val = 0.6e0, xc = 113.5e3, yc =120.0e3, zc= -15.0e3, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
-&DIST2D shape='rectangle', val = 0.2e0, xc = 113.5e3, yc =120.0e3, zc= -15.0e3, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = 0.6e0, xc = 113.5e3, yc =120.0e3, zc= -15.0e3, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
+&DIST2D shape='rectangle', val = 0.2e0, xc = 113.5e3, yc =120.0e3, zc= -15.0e3, lx=227.0e3, ly=80.0e3, lz=30000.0e0 /
&BEGIN_FAULT /
&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-100.0e6,n1=0,n2=3,n3=4 /
-&DIST2D shape='rectangle', val = -5.69e6, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
-&DIST2D shape='rectangle', val = -3.54e6, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
-&DIST2D shape='rectangle-taper', val = -5.69e6, valh=-0.569e6, xc = 82205.7714e0, yc =120.0e3, zc= -4.5e3, lx=5196.15242e0, ly=80.0e3, lz=3.0e3 /
-&DIST2D shape='rectangle', val =-1.238e6, xc = 74411.5427e0, yc =120.0e3, zc= -9.0e3, lx=31176.9145e0, ly=80.0e3, lz=18.0e3 /
-&DIST2D shape='rectangle', val =-12.38e6, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
-&DIST2D shape='rectangle-taper', val =-12.38e6, valh=-1.238e6, xc = 82205.7714e0, yc =120.0e3, zc= -4.5e3, lx=5196.15242e0, ly=80.0e3, lz=3.0e3 /
-&DIST2D shape='rectangle', val = -9.04e6, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
+&DIST2D shape='rectangle', val = -5.69e6, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
+&DIST2D shape='rectangle', val = -3.54e6, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
+&DIST2D shape='rectangle-taper', val = -5.69e6, valh=-0.569e6, xc = 82205.7714e0, yc =120.0e3, zc= -4.5e3, lx=5196.15242e0, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val =-1.238e6, xc = 74411.5427e0, yc =120.0e3, zc= -9.0e3, lx=31176.9145e0, ly=80.0e3, lz=18.0e3 /
+&DIST2D shape='rectangle', val =-12.38e6, xc = 70514.4284e0, yc =120.0e3, zc=-11.25e3, lx=23382.6859e0, ly=80.0e3, lz=13.5e3 /
+&DIST2D shape='rectangle-taper', val =-12.38e6, valh=-1.238e6, xc = 82205.7714e0, yc =120.0e3, zc= -4.5e3, lx=5196.15242e0, ly=80.0e3, lz=3.0e3 /
+&DIST2D shape='rectangle', val = -9.04e6, xc = 43894.8822e0, yc =120.0e3, zc=-22.0e3, lx=29856.4065e0, ly=80.0e3, lz= 8.0e3 /
&SWF mus=1000.e0,mud=1000.0e0,dc=0.5e0,nmus=1,nmud=1,ndc=0 /
-&DIST2D shape='rectangle', val = 0.6e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
-&DIST2D shape='rectangle', val = 0.2e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
+&DIST2D shape='rectangle', val = 0.6e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
+&DIST2D shape='rectangle', val = 0.2e0, xc = 60.1915e3, yc =120.0e3, zc=-13.0e3, lx=59616.98e0, ly=80.0e3, lz=26.0e3 /
diff --git a/tomographic_model/DATA/Par_file b/tomographic_model/DATA/Par_file
index 2d27ec6..f7e19a6 100644
--- a/tomographic_model/DATA/Par_file
+++ b/tomographic_model/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.01
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .true.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/tpv102/DATA/Par_file b/tpv102/DATA/Par_file
index 7fe6985..093325b 100644
--- a/tpv102/DATA/Par_file
+++ b/tpv102/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.0025d0
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .true.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/tpv102/DATA/Par_file_faults b/tpv102/DATA/Par_file_faults
index 04f54b8..6eb78ab 100644
--- a/tpv102/DATA/Par_file_faults
+++ b/tpv102/DATA/Par_file_faults
@@ -1,7 +1,7 @@
1
0.00051 ! #tag #eta(damping) FAULT A
1 ! 1 = dyn 2=kin
-6000 ! NTOUT : Number of time steps
+6000 ! NTOUT : Number of time steps
200 ! NTSNAP: time interation of snapshots
-1e0 ! V_HEALING (-1 : Healing off)
1e-3 ! V_RUPT
diff --git a/tpv103/DATA/Par_file b/tpv103/DATA/Par_file
index 7fe6985..093325b 100644
--- a/tpv103/DATA/Par_file
+++ b/tpv103/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.0025d0
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .true.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/tpv103/DATA/Par_file_faults b/tpv103/DATA/Par_file_faults
index f50b5ef..a270133 100644
--- a/tpv103/DATA/Par_file_faults
+++ b/tpv103/DATA/Par_file_faults
@@ -1,7 +1,7 @@
1
0.00134 ! #tag #eta(damping) FAULT A
1 ! 1 = dyn 2=kin
-1200 ! NTOUT : Number of time steps
+1200 ! NTOUT : Number of time steps
200 ! NTSNAP: time interation of snapshots
-1e0 ! V_HEALING (-1 : Healing off)
1e-3 ! V_RUPT
diff --git a/tpv15/DATA/Par_file b/tpv15/DATA/Par_file
index 8cb3837..fb15a43 100644
--- a/tpv15/DATA/Par_file
+++ b/tpv15/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.001d0
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .true.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/tpv15/DATA/Par_file_faults b/tpv15/DATA/Par_file_faults
index de5adc7..1f2daaf 100644
--- a/tpv15/DATA/Par_file_faults
+++ b/tpv15/DATA/Par_file_faults
@@ -2,7 +2,7 @@
0.00067 ! #tag #eta(damping) FAULT Main fault
0.00067 ! #tag #eta(damping) FAULT Branch
1 ! 1 = dyn 2=kin
-100 ! NTOUT : Number of time steps
+100 ! NTOUT : Number of time steps
1000 ! NTSNAP: time interation of snapshots
-1e0 ! V_HEALING (-1 : Healing off)
1e-3 ! V_RUPT
@@ -10,18 +10,18 @@
&BEGIN_FAULT /
&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-120.0e6,n1=2,n2=0,n3=0 /
-&DIST2D shape='rectangle', val = -70.0e6, xc = -2.0e3, yc =0.0e3, zc= -7.5e3, lx=28.0e3, ly=100.0e0, lz=15.0e3 /
+&DIST2D shape='rectangle', val = -70.0e6, xc = -2.0e3, yc =0.0e3, zc= -7.5e3, lx=28.0e3, ly=100.0e0, lz=15.0e3 /
&DIST2D shape='rectangle', val = -81.6e6, xc = -8.0e3, yc =0.0e3, zc= -7.5e3, lx=3.1e3, ly=100.0e0, lz=3.0e3 /
&SWF mus=1000.0e0,mud=1000.0e0,dc=0.4e0,nmus=1,nmud=1,ndc=0 /
-&DIST2D shape='rectangle', val = 0.677e0, xc = -2.0e3, yc =0.0e3, zc= -7.5e3, lx=28.0e3, ly=100.0e0, lz=15.0e3 /
-&DIST2D shape='rectangle', val = 0.525e0, xc = -2.0e3, yc =0.0e3, zc= -7.5e3, lx=28.0e3, ly=100.0e0, lz=15.0e3 /
+&DIST2D shape='rectangle', val = 0.677e0, xc = -2.0e3, yc =0.0e3, zc= -7.5e3, lx=28.0e3, ly=100.0e0, lz=15.0e3 /
+&DIST2D shape='rectangle', val = 0.525e0, xc = -2.0e3, yc =0.0e3, zc= -7.5e3, lx=28.0e3, ly=100.0e0, lz=15.0e3 /
&BEGIN_FAULT /
&INIT_STRESS S1=0.0e6,S2=0.0e0,S3=-120.0e6,n1=1,n2=0,n3=0 /
-&DIST2D shape='rectangle', val = -78.0e6, xc = 5152.85115e0, yc =-3025.0e0, zc=-7.5e3, lx=10478.9074e0, ly=6050.0e0, lz=15.0e3 /
+&DIST2D shape='rectangle', val = -78.0e6, xc = 5152.85115e0, yc =-3025.0e0, zc=-7.5e3, lx=10478.9074e0, ly=6050.0e0, lz=15.0e3 /
&SWF mus=1000.e0,mud=1000.0e0,dc=0.4e0,nmus=1,nmud=1,ndc=0 /
-&DIST2D shape='rectangle', val = 0.677e0, xc = 5152.85115e0, yc =-3025.0e0, zc=-7.5e3, lx=10478.9074e0, ly=6050.0e0, lz=15.0e3 /
-&DIST2D shape='rectangle', val = 0.525e0, xc = 5152.85115e0, yc =-3025.0e0, zc=-7.5e3, lx=10478.9074e0, ly=6050.0e0, lz=15.0e3 /
+&DIST2D shape='rectangle', val = 0.677e0, xc = 5152.85115e0, yc =-3025.0e0, zc=-7.5e3, lx=10478.9074e0, ly=6050.0e0, lz=15.0e3 /
+&DIST2D shape='rectangle', val = 0.525e0, xc = 5152.85115e0, yc =-3025.0e0, zc=-7.5e3, lx=10478.9074e0, ly=6050.0e0, lz=15.0e3 /
diff --git a/tpv16/DATA/Par_file b/tpv16/DATA/Par_file
index b67f6cb..7e36939 100644
--- a/tpv16/DATA/Par_file
+++ b/tpv16/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.0025d0
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .true.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/tpv16/DATA/Par_file_faults b/tpv16/DATA/Par_file_faults
index 27e66ea..5d2f457 100644
--- a/tpv16/DATA/Par_file_faults
+++ b/tpv16/DATA/Par_file_faults
@@ -1,7 +1,7 @@
1
0.00134 ! #tag #eta(damping) FAULT A
1 ! 1 = dyn 2=kin
-1000 ! NTOUT : Number of time steps
+1000 ! NTOUT : Number of time steps
200 ! NTSNAP: time interation of snapshots
-1e0 ! V_HEALING (-1 : Healing off)
1e-3 ! V_RUPT
diff --git a/tpv5/DATA/Par_file b/tpv5/DATA/Par_file
index 7722794..aeb8fec 100644
--- a/tpv5/DATA/Par_file
+++ b/tpv5/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.002d0
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .true.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
diff --git a/tpv5/DATA/Par_file_faults b/tpv5/DATA/Par_file_faults
index 1ea86d3..16646ca 100644
--- a/tpv5/DATA/Par_file_faults
+++ b/tpv5/DATA/Par_file_faults
@@ -1,7 +1,7 @@
1
0.00134 ! #tag #eta(damping)
1 ! 1 = dyn 2=kin
-1000 ! NTOUT : Number of time steps
+1000 ! NTOUT : Number of time steps
100 ! NTSNAP: time interation of snapshots
-1.0e0 ! V_HEALING (-1 : Healing off)
1e-3 ! V_RUPT
@@ -11,4 +11,4 @@
&DIST2D shape='square', val = 81.6e6, xc = 0e0, yc =0e0, zc= -7500.0e0, l=3000.0e0 /
&DIST2D shape='square', val = 62.0e6, xc = 7500.0e0, yc =0e0, zc= -7500.0e0, l=3000.0e0 /
&SWF mus=10000.0e0, nmus=1, mud=0.525e0,dc=0.4e0/
-&DIST2D shape='rectangle', val = 0.677e0, xc = 0e0, yc =0e0, zc= -7500.0e0, lx=30000.0e0, ly=1.0e0, lz=15000.0e0 /
+&DIST2D shape='rectangle', val = 0.677e0, xc = 0e0, yc =0e0, zc= -7500.0e0, lx=30000.0e0, ly=1.0e0, lz=15000.0e0 /
diff --git a/waterlayered_halfspace/DATA/Par_file b/waterlayered_halfspace/DATA/Par_file
index 1531f40..d95bf5b 100644
--- a/waterlayered_halfspace/DATA/Par_file
+++ b/waterlayered_halfspace/DATA/Par_file
@@ -23,7 +23,7 @@ DT = 0.005
# If you use our internal mesher, the only option is 8-node bricks (27-node elements are not supported)
# CUBIT does not support HEX27 elements either (it can generate them, but they are flat, i.e. identical to HEX8).
# To generate HEX27 elements with curvature properly taken into account, you can use Gmsh http://geuz.org/gmsh/
-NGNOD = 8
+NGNOD = 8
# models:
# available options are:
@@ -35,7 +35,7 @@ NGNOD = 8
MODEL = default
# parameters describing the model
-OCEANS = .false.
+APPROXIMATE_OCEAN_LOAD = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
ANISOTROPY = .false.
@@ -55,7 +55,7 @@ ABSORBING_CONDITIONS = .true.
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
# C-PML boundary conditions for a regional simulation
-PML_CONDITIONS = .false.
+PML_CONDITIONS = .false.
# C-PML top surface
PML_INSTEAD_OF_FREE_SURFACE = .false.
More information about the CIG-COMMITS
mailing list