[cig-commits] r20944 - in seismo/3D/SPECFEM3D/trunk: . doc/USER_MANUAL examples/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA examples/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA examples/Mount_StHelens/DATA examples/homogeneous_halfspace_HEX27/DATA examples/homogeneous_halfspace_HEX8/DATA examples/homogeneous_poroelastic/DATA examples/layered_halfspace/DATA examples/meshfem3D_examples/many_interfaces/DATA examples/meshfem3D_examples/simple_model/DATA examples/meshfem3D_examples/socal1D/DATA examples/tomographic_model/DATA examples/waterlayered_halfspace/DATA src/shared
joseph.charles at geodynamics.org
joseph.charles at geodynamics.org
Fri Oct 26 15:42:57 PDT 2012
Author: joseph.charles
Date: 2012-10-26 15:42:56 -0700 (Fri, 26 Oct 2012)
New Revision: 20944
Removed:
seismo/3D/SPECFEM3D/trunk/doc/USER_MANUAL/manual_SPECFEM3D_Cartesian.pdf
Modified:
seismo/3D/SPECFEM3D/trunk/doc/USER_MANUAL/manual_SPECFEM3D_Cartesian.tex
seismo/3D/SPECFEM3D/trunk/examples/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/Mount_StHelens/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX27/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/homogeneous_poroelastic/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/layered_halfspace/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/many_interfaces/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/simple_model/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/tomographic_model/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/DATA/Par_file
seismo/3D/SPECFEM3D/trunk/src/shared/read_parameter_file.f90
seismo/3D/SPECFEM3D/trunk/todo_list_please_dont_remove.txt
Log:
updates the user manual; modifies all parameter files
Deleted: seismo/3D/SPECFEM3D/trunk/doc/USER_MANUAL/manual_SPECFEM3D_Cartesian.pdf
===================================================================
(Binary files differ)
Modified: seismo/3D/SPECFEM3D/trunk/doc/USER_MANUAL/manual_SPECFEM3D_Cartesian.tex
===================================================================
--- seismo/3D/SPECFEM3D/trunk/doc/USER_MANUAL/manual_SPECFEM3D_Cartesian.tex 2012-10-26 21:01:27 UTC (rev 20943)
+++ seismo/3D/SPECFEM3D/trunk/doc/USER_MANUAL/manual_SPECFEM3D_Cartesian.tex 2012-10-26 22:42:56 UTC (rev 20944)
@@ -1059,38 +1059,42 @@
Please note that the Vp- and Vs-velocities of your model are given for a reference frequency.
To change this reference frequency, you change the value of
\texttt{ATTENUATION\_f0\_REFERENCE} in the main constants file \texttt{constants.h} found in subdirectory \texttt{src/shared/}.
+\item [{\texttt{ANISOTROPY}}] Set to \texttt{.true.} if you
+want to use an anisotropy model. Please see the file \texttt{model\_aniso.f90} in subdirectory \texttt{src/generate\_databases/} for the current implementation of anisotropic models.
\item [{\texttt{USE\_OLSEN\_ATTENUATION}}] Set to \texttt{.true.} if you
want to use the attenuation model that scaled from the S-wave speed model
using Olsen's empirical relation (see \citet{OlDaBr03}).
-\item [{\texttt{ANISOTROPY}}] Set to \texttt{.true.} if you
-want to use an anisotropy model. Please see the file \texttt{model\_aniso.f90} in subdirectory \texttt{src/generate\_databases/} for the current implementation of anisotropic models.
+\item [{\texttt{OLSEN\_ATTENUATION\_RATIO}}] Determines the Olsen's constant in Olsen's empirical relation (see \citet{OlDaBr03}).
\item [{\texttt{ABSORBING\_CONDITIONS}}] Set to \texttt{.true.} to turn
on Clayton-Enquist absorbing boundary conditions (see \citet{KoTr99}).
\item [{\texttt{ABSORB\_INSTEAD\_OF\_FREE\_SURFACE}}] Set to \texttt{.true.} to turn on
absorbing boundary conditions on the top surface which by default constitutes a free surface of the model.
+\item [{\texttt{CREATE\_SHAKEMAP}}] Set this flag to \texttt{.true.} to
+create a ShakeMap\textregistered, i.e., a peak ground velocity map
+of the maximum absolute value of the two horizontal components of the velocity vector.
\item [{\texttt{MOVIE\_SURFACE}}] Set to \texttt{.false.}, unless you want
to create a movie of seismic wave propagation on the Earth's surface.
Turning this option on generates large output files. See Section~\ref{sec:Movies}
for a discussion on the generation of movies. This feature is only relevant for the solver.
+\item [{\texttt{MOVIE\_TYPE}}] Set this flag to 1 to show the top
+ surface (tomography + oceans) only, to 2 to show all external faces
+ of the mesh (i.e. topography + vertical edges + bottom) in shakemaps and surface movies.
\item [{\texttt{MOVIE\_VOLUME}}] Set to \texttt{.false.}, unless you want
to create a movie of seismic wave propagation in the Earth's interior.
Turning this option on generates huge output files. See Section~\ref{sec:Movies}
for a discussion on the generation of movies. This feature is only relevant for the solver.
-\item [{\texttt{NTSTEP\_BETWEEN\_FRAMES}}] Determines the number of timesteps
-between movie frames. Typically you want to save a snapshot every
-100 timesteps. The smaller you make this number the more output will
-be generated! See Section~\ref{sec:Movies} for a discussion on the
-generation of movies. This feature is only relevant for the solver.
-\item [{\texttt{CREATE\_SHAKEMAP}}] Set this flag to \texttt{.true.} to
-create a ShakeMap\textregistered, i.e., a peak ground velocity map
-of the maximum absolute value of the two horizontal components of the velocity vector.
\item [{\texttt{SAVE\_DISPLACEMENT}}] Set this flag to \texttt{.true.}
if you want to save the displacement instead of velocity for the movie
frames.
\item [{\texttt{USE\_HIGHRES\_FOR\_MOVIES}}] Set this flag to \texttt{.true.}
if you want to save the values at all the NGLL grid points for the
movie frames.
-\item [{\texttt{HDUR\_MOVIE}}] determines the half duration of the source time function for the movie simulations.
+\item [{\texttt{NTSTEP\_BETWEEN\_FRAMES}}] Determines the number of timesteps
+between movie frames. Typically you want to save a snapshot every
+100 timesteps. The smaller you make this number the more output will
+be generated! See Section~\ref{sec:Movies} for a discussion on the
+generation of movies. This feature is only relevant for the solver.
+\item [{\texttt{HDUR\_MOVIE}}] Determines the half duration of the source time function for the movie simulations.
When this parameter is set to be 0, a default half duration that corresponds to the accuracy of the simulation is provided.
Otherwise, it adds this half duration to the half duration specified in the source file \texttt{CMTSOLUTION},
thus simulates longer periods to make the movie images look smoother.
@@ -1195,7 +1199,7 @@
\item the movie control parameters \texttt{MOVIE\_SURFACE}, \texttt{MOVIE\_VOLUME},
and \texttt{NTSTEPS\_BETWEEN\_FRAMES}
\item the ShakeMap\textregistered option \texttt{CREATE\_SHAKEMAP}
-\item the output information parameters \texttt{NTSTEP\_BETWEEN\_OUTPUT\_INFO}
+\item the output information parameters \texttt{MOVIE\_TYPE}, \texttt{NTSTEP\_BETWEEN\_OUTPUT\_INFO}
and \texttt{NTSTEP\_BETWEEN\_OUTPUT\_}~\\
\texttt{SEISMOS}
\item the \texttt{PRINT\_SOURCE\_TIME\_FUNCTION} flags
@@ -1942,7 +1946,7 @@
\section{\label{sec:Movies}Movies}
To make a surface or volume movie of the simulation, set parameters
-\texttt{MOVIE\_SURFACE}, \texttt{MOVIE\_VOLUME}, and \texttt{NTSTEP\_BETWEEN\_FRAMES}
+\texttt{MOVIE\_SURFACE}, \texttt{MOVIE\_VOLUME}, \texttt{MOVIE\_TYPE}, and \texttt{NTSTEP\_BETWEEN\_FRAMES}
in the \texttt{Par\_file}. Turning on the movie flags, in particular
\texttt{MOVIE\_VOLUME}, produces large output files. \texttt{MOVIE\_VOLUME}
files are saved in the \texttt{LOCAL\_PATH} directory, whereas \texttt{MOVIE\_SURFACE}
@@ -1974,7 +1978,7 @@
When running \texttt{xspecfem3D} with the \texttt{MOVIE\_SURFACE}
flag turned on, the code outputs \texttt{moviedata??????} files in
the \texttt{OUTPUT\_FILES} directory. There are several flags
-in the main parameter file \texttt{Par\_file} that control the output of these moviedata files
+in the main parameter file \texttt{Par\_file} that control the output of these movie data files
(see section \ref{cha:Main-Parameter} for details):
\texttt{\small NTSTEP\_BETWEEN\_FRAMES} to set the timesteps between frames,
\texttt{\small SAVE\_DISPLACEMENT} to save displacement instead of velocity,
Modified: seismo/3D/SPECFEM3D/trunk/examples/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file 2012-10-26 21:01:27 UTC (rev 20943)
+++ seismo/3D/SPECFEM3D/trunk/examples/BENCHMARK_CLAERBOUT_ADJOINT/ACOUSTIC/DATA/Par_file 2012-10-26 22:42:56 UTC (rev 20944)
@@ -38,24 +38,25 @@
OCEANS = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
+GRAVITY = .false.
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
+
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
-# Olsen's constant for Q_mu = constant * v_s attenuation rule
-OLSEN_ATTENUATION_RATIO = 0.05
-
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
-MOVIE_TYPE = 1
CREATE_SHAKEMAP = .false.
MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
MOVIE_VOLUME = .false.
SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
Modified: seismo/3D/SPECFEM3D/trunk/examples/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file 2012-10-26 21:01:27 UTC (rev 20943)
+++ seismo/3D/SPECFEM3D/trunk/examples/BENCHMARK_CLAERBOUT_ADJOINT/ELASTIC/DATA/Par_file 2012-10-26 22:42:56 UTC (rev 20944)
@@ -38,24 +38,25 @@
OCEANS = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
+GRAVITY = .false.
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
+
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
-# Olsen's constant for Q_mu = constant * v_s attenuation rule
-OLSEN_ATTENUATION_RATIO = 0.05
-
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
-MOVIE_TYPE = 1
CREATE_SHAKEMAP = .false.
MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
MOVIE_VOLUME = .false.
SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
Modified: seismo/3D/SPECFEM3D/trunk/examples/Mount_StHelens/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/Mount_StHelens/DATA/Par_file 2012-10-26 21:01:27 UTC (rev 20943)
+++ seismo/3D/SPECFEM3D/trunk/examples/Mount_StHelens/DATA/Par_file 2012-10-26 22:42:56 UTC (rev 20944)
@@ -38,24 +38,25 @@
OCEANS = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
+GRAVITY = .false.
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
+
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
-# Olsen's constant for Q_mu = constant * v_s attenuation rule
-OLSEN_ATTENUATION_RATIO = 0.05
-
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
-MOVIE_TYPE = 1
CREATE_SHAKEMAP = .false.
MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
MOVIE_VOLUME = .false.
SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
Modified: seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX27/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX27/DATA/Par_file 2012-10-26 21:01:27 UTC (rev 20943)
+++ seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX27/DATA/Par_file 2012-10-26 22:42:56 UTC (rev 20944)
@@ -38,25 +38,25 @@
OCEANS = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
+
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
-# Olsen's constant for Q_mu = constant * v_s attenuation rule
-OLSEN_ATTENUATION_RATIO = 0.05
-
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
-MOVIE_TYPE = 1
CREATE_SHAKEMAP = .false.
MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
MOVIE_VOLUME = .false.
SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
Modified: seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8/DATA/Par_file 2012-10-26 21:01:27 UTC (rev 20943)
+++ seismo/3D/SPECFEM3D/trunk/examples/homogeneous_halfspace_HEX8/DATA/Par_file 2012-10-26 22:42:56 UTC (rev 20944)
@@ -38,25 +38,25 @@
OCEANS = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
+
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
-# Olsen's constant for Q_mu = constant * v_s attenuation rule
-OLSEN_ATTENUATION_RATIO = 0.05
-
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
-MOVIE_TYPE = 1
CREATE_SHAKEMAP = .false.
MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
MOVIE_VOLUME = .false.
SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
Modified: seismo/3D/SPECFEM3D/trunk/examples/homogeneous_poroelastic/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/homogeneous_poroelastic/DATA/Par_file 2012-10-26 21:01:27 UTC (rev 20943)
+++ seismo/3D/SPECFEM3D/trunk/examples/homogeneous_poroelastic/DATA/Par_file 2012-10-26 22:42:56 UTC (rev 20944)
@@ -38,25 +38,25 @@
OCEANS = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
+
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
-# Olsen's constant for Q_mu = constant * v_s attenuation rule
-OLSEN_ATTENUATION_RATIO = 0.05
-
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
-MOVIE_TYPE = 1
CREATE_SHAKEMAP = .false.
MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
MOVIE_VOLUME = .false.
SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
Modified: seismo/3D/SPECFEM3D/trunk/examples/layered_halfspace/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/layered_halfspace/DATA/Par_file 2012-10-26 21:01:27 UTC (rev 20943)
+++ seismo/3D/SPECFEM3D/trunk/examples/layered_halfspace/DATA/Par_file 2012-10-26 22:42:56 UTC (rev 20944)
@@ -38,25 +38,25 @@
OCEANS = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
+
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
-# Olsen's constant for Q_mu = constant * v_s attenuation rule
-OLSEN_ATTENUATION_RATIO = 0.05
-
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
-MOVIE_TYPE = 1
CREATE_SHAKEMAP = .false.
MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
MOVIE_VOLUME = .false.
SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
Modified: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/many_interfaces/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/many_interfaces/DATA/Par_file 2012-10-26 21:01:27 UTC (rev 20943)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/many_interfaces/DATA/Par_file 2012-10-26 22:42:56 UTC (rev 20944)
@@ -38,25 +38,25 @@
OCEANS = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
+
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
-# Olsen's constant for Q_mu = constant * v_s attenuation rule
-OLSEN_ATTENUATION_RATIO = 0.05
-
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
-MOVIE_TYPE = 1
CREATE_SHAKEMAP = .false.
MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
MOVIE_VOLUME = .false.
SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
Modified: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/simple_model/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/simple_model/DATA/Par_file 2012-10-26 21:01:27 UTC (rev 20943)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/simple_model/DATA/Par_file 2012-10-26 22:42:56 UTC (rev 20944)
@@ -38,25 +38,25 @@
OCEANS = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
+
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .false.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
-# Olsen's constant for Q_mu = constant * v_s attenuation rule
-OLSEN_ATTENUATION_RATIO = 0.05
-
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
-MOVIE_TYPE = 1
CREATE_SHAKEMAP = .false.
MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
MOVIE_VOLUME = .false.
SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
Modified: seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/DATA/Par_file 2012-10-26 21:01:27 UTC (rev 20943)
+++ seismo/3D/SPECFEM3D/trunk/examples/meshfem3D_examples/socal1D/DATA/Par_file 2012-10-26 22:42:56 UTC (rev 20944)
@@ -38,25 +38,25 @@
OCEANS = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
+
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
-# Olsen's constant for Q_mu = constant * v_s attenuation rule
-OLSEN_ATTENUATION_RATIO = 0.05
-
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
-MOVIE_TYPE = 1
CREATE_SHAKEMAP = .false.
MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
MOVIE_VOLUME = .false.
SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
Modified: seismo/3D/SPECFEM3D/trunk/examples/tomographic_model/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/tomographic_model/DATA/Par_file 2012-10-26 21:01:27 UTC (rev 20943)
+++ seismo/3D/SPECFEM3D/trunk/examples/tomographic_model/DATA/Par_file 2012-10-26 22:42:56 UTC (rev 20944)
@@ -38,25 +38,25 @@
OCEANS = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
+
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
-# Olsen's constant for Q_mu = constant * v_s attenuation rule
-OLSEN_ATTENUATION_RATIO = 0.05
-
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
-MOVIE_TYPE = 1
CREATE_SHAKEMAP = .false.
MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
MOVIE_VOLUME = .false.
SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
Modified: seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/DATA/Par_file 2012-10-26 21:01:27 UTC (rev 20943)
+++ seismo/3D/SPECFEM3D/trunk/examples/waterlayered_halfspace/DATA/Par_file 2012-10-26 22:42:56 UTC (rev 20944)
@@ -38,25 +38,25 @@
OCEANS = .false.
TOPOGRAPHY = .false.
ATTENUATION = .false.
-USE_OLSEN_ATTENUATION = .false.
ANISOTROPY = .false.
GRAVITY = .false.
+# Olsen's constant for Q_mu = constant * v_s attenuation rule
+USE_OLSEN_ATTENUATION = .false.
+OLSEN_ATTENUATION_RATIO = 0.05
+
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
# absorbing top surface (defined in mesh as 'free_surface_file')
ABSORB_INSTEAD_OF_FREE_SURFACE = .false.
-# Olsen's constant for Q_mu = constant * v_s attenuation rule
-OLSEN_ATTENUATION_RATIO = 0.05
-
# save AVS or OpenDX movies
# MOVIE_TYPE = 1 to show the top surface
# MOVIE_TYPE = 2 to show all the external faces of the mesh
-MOVIE_TYPE = 1
CREATE_SHAKEMAP = .false.
MOVIE_SURFACE = .false.
+MOVIE_TYPE = 1
MOVIE_VOLUME = .false.
SAVE_DISPLACEMENT = .false.
USE_HIGHRES_FOR_MOVIES = .false.
Modified: seismo/3D/SPECFEM3D/trunk/src/shared/read_parameter_file.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/src/shared/read_parameter_file.f90 2012-10-26 21:01:27 UTC (rev 20943)
+++ seismo/3D/SPECFEM3D/trunk/src/shared/read_parameter_file.f90 2012-10-26 22:42:56 UTC (rev 20944)
@@ -114,22 +114,22 @@
if(err_occurred() /= 0) return
call read_value_logical(ATTENUATION, 'model.ATTENUATION')
if(err_occurred() /= 0) return
+ call read_value_logical(ANISOTROPY, 'model.ANISOTROPY')
+ if(err_occurred() /= 0) return
call read_value_logical(USE_OLSEN_ATTENUATION, 'model.USE_OLSEN_ATTENUATION')
if(err_occurred() /= 0) return
- call read_value_logical(ANISOTROPY, 'model.ANISOTROPY')
+ call read_value_double_precision(OLSEN_ATTENUATION_RATIO, 'model.OLSEN_ATTENUATION_RATIO')
if(err_occurred() /= 0) return
call read_value_logical(ABSORBING_CONDITIONS, 'solver.ABSORBING_CONDITIONS')
if(err_occurred() /= 0) return
call read_value_logical(ABSORB_INSTEAD_OF_FREE_SURFACE, 'model.ABSORB_INSTEAD_OF_FREE_SURFACE')
if(err_occurred() /= 0) return
- call read_value_double_precision(OLSEN_ATTENUATION_RATIO, 'model.OLSEN_ATTENUATION_RATIO')
- if(err_occurred() /= 0) return
- call read_value_integer(MOVIE_TYPE, 'solver.MOVIE_TYPE')
- if(err_occurred() /= 0) return
call read_value_logical(CREATE_SHAKEMAP, 'solver.CREATE_SHAKEMAP')
if(err_occurred() /= 0) return
call read_value_logical(MOVIE_SURFACE, 'solver.MOVIE_SURFACE')
if(err_occurred() /= 0) return
+ call read_value_integer(MOVIE_TYPE, 'solver.MOVIE_TYPE')
+ if(err_occurred() /= 0) return
call read_value_logical(MOVIE_VOLUME, 'solver.MOVIE_VOLUME')
if(err_occurred() /= 0) return
call read_value_logical(SAVE_DISPLACEMENT, 'solver.SAVE_DISPLACEMENT')
Modified: seismo/3D/SPECFEM3D/trunk/todo_list_please_dont_remove.txt
===================================================================
--- seismo/3D/SPECFEM3D/trunk/todo_list_please_dont_remove.txt 2012-10-26 21:01:27 UTC (rev 20943)
+++ seismo/3D/SPECFEM3D/trunk/todo_list_please_dont_remove.txt 2012-10-26 22:42:56 UTC (rev 20944)
@@ -8,62 +8,9 @@
CPML :
------------------------------------------------
-- suggestion 00: about adding new CPML flags and also new (or modified) movie and shakemap flags + VERCE project modifications
+- suggestion 00: about VERCE project modifications
----------------
-From Dimitri: here are two more parameters to rename and/or move from the constants.h file to the input Par_file:
-
-could you also rename "USE_RICKER_IPATI" to "USE_RICKER" or "USE_RICKER_TIME_FUNCTION" everywhere in the source code and in the input file?
-Because "Ipati" is the name of an industrial field and thus we should make the name of that parameter less specific.
-
-Let us also move that parameter from the constants.h file to the input file (Par_file), to avoid having to recompile the code
-(see my previous email about this)
-
-================
-
-As discussed yesterday during our Skype call, let us make the
-modifications below in the input Par_file of SPECFEM3D.
-NGNOD should be added and set to 8 or 27, and then when reading the
-Par_file the code should check that it is equal to one or the other, and
-otherwise exit with an error message.
-The code should then automatically set NGNOD2D to 4 if NGNOD == 8 and to 9 if NGNOD == 27.
-
-For the new movie flags, I suggest changing MOVIE_SURFACE and
-EXTERNAL_MESH_MOVIE_SURFACE, which are confusing, to MOVIE_SURFACE only,
-and then define a second (integer) parameter, called for instance
-MOVIE_TYPE, which could be equal to 1 to show the top surface
-(topography + oceans) only, to 2 to show all the external faces of the
-mesh (i.e. topography + vertical edges + bottom), and so on.
-
-Let me add this to the todo list. Let us also decide who does it
-(since new flags will be needed in the Par_file for PML, I guess it will
-be easy and better to do all the modifications at the same time).
-
-Let us also modify the code to allow for shakemaps and movies at the
-same time, rather than having to run the code twice.
-
-> add NGNOD = 27 or 8 to the input Par_file (move it from shared/constants.h.in to the Par_file, since now either value can be used depending on the mesh,
-whether it is HEX8 or HEX27)
->
->
-> In the VERCE framework they ask some modifications in order to improve
-> the flexibility of the code if installed as a module....
->
-> Move from constants.h to Par_file the flags:
-> EXTERNAL_MESH_MOVIE_SURFACE
-> EXTERNAL_MESH_CREATE_SHAKEMAP
-> OLSEN_ATTENUATION_RATIO
->
-> Allow shakemaps and surface movies to be generated during the same run
-> (EXTERNAL_MESH_MOVIE_SURFACE and EXTERNAL_MESH_CREATE_SHAKEMAP not
-> mutually exclusive)
->
-> We believe also that
->
-> Explain in the manual the difference in using
-> EXTERNAL_MESH_MOVIE_SURFACE instead of MOVIE_SURFACE, or
-> EXTERNAL_MESH_CREATE_SHAKEMAP instead of CREATE_SHAKEMAP
->
> Possibly set the name of the tomography file in the Par_file instead
> of hard-coding it in model_tomography.f90 or creating a separate
> parfile
@@ -73,7 +20,7 @@
>
> ---- multifiles tomography
> Piero told me that the run for turkey/europe mesh had this possibility
-> so probably it is matter to commint that files in the svn version
+> so probably it is matter to commit that files in the svn version
>
@@ -639,6 +586,64 @@
Par_file modifications:
------------------------------------------------
+- suggestion 00: about adding new CPML flags and also new (or modified) movie and shakemap flags + VERCE project modifications
+----------------
+
+From Dimitri: here are two more parameters to rename and/or move from the constants.h file to the input Par_file:
+
+could you also rename "USE_RICKER_IPATI" to "USE_RICKER" or "USE_RICKER_TIME_FUNCTION" everywhere in the source code and in the input file?
+Because "Ipati" is the name of an industrial field and thus we should make the name of that parameter less specific.
+
+Let us also move that parameter from the constants.h file to the input file (Par_file), to avoid having to recompile the code
+(see my previous email about this)
+
+================
+
+As discussed yesterday during our Skype call, let us make the
+modifications below in the input Par_file of SPECFEM3D.
+NGNOD should be added and set to 8 or 27, and then when reading the
+Par_file the code should check that it is equal to one or the other, and
+otherwise exit with an error message.
+The code should then automatically set NGNOD2D to 4 if NGNOD == 8 and to 9 if NGNOD == 27.
+
+For the new movie flags, I suggest changing MOVIE_SURFACE and
+EXTERNAL_MESH_MOVIE_SURFACE, which are confusing, to MOVIE_SURFACE only,
+and then define a second (integer) parameter, called for instance
+MOVIE_TYPE, which could be equal to 1 to show the top surface
+(topography + oceans) only, to 2 to show all the external faces of the
+mesh (i.e. topography + vertical edges + bottom), and so on.
+
+Let me add this to the todo list. Let us also decide who does it
+(since new flags will be needed in the Par_file for PML, I guess it will
+be easy and better to do all the modifications at the same time).
+
+Let us also modify the code to allow for shakemaps and movies at the
+same time, rather than having to run the code twice.
+
+> add NGNOD = 27 or 8 to the input Par_file (move it from shared/constants.h.in to the Par_file, since now either value can be used depending on the mesh,
+whether it is HEX8 or HEX27)
+>
+>
+> In the VERCE framework they ask some modifications in order to improve
+> the flexibility of the code if installed as a module....
+>
+> Move from constants.h to Par_file the flags:
+> EXTERNAL_MESH_MOVIE_SURFACE
+> EXTERNAL_MESH_CREATE_SHAKEMAP
+> OLSEN_ATTENUATION_RATIO
+>
+> Allow shakemaps and surface movies to be generated during the same run
+> (EXTERNAL_MESH_MOVIE_SURFACE and EXTERNAL_MESH_CREATE_SHAKEMAP not
+> mutually exclusive)
+>
+> We believe also that
+>
+> Explain in the manual the difference in using
+> EXTERNAL_MESH_MOVIE_SURFACE instead of MOVIE_SURFACE, or
+> EXTERNAL_MESH_CREATE_SHAKEMAP instead of CREATE_SHAKEMAP
+
+
+
- suggestion 05:
----------------
More information about the CIG-COMMITS
mailing list