[cig-commits] [commit] devel: doc: Spelling and grammar fixes. (700dd71)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Thu Feb 6 06:09:59 PST 2014
Repository : ssh://geoshell/specfem2d
On branch : devel
Link : https://github.com/geodynamics/specfem2d/compare/24dca31bf087ee51d2d1db912640d58ade2df07c...a3880316a04510a6c88ed36ceee2fab3c0169f6e
>---------------------------------------------------------------
commit 700dd71314d927d1f40349898ecf2e050d244ee4
Author: Elliott Sales de Andrade <esalesde at physics.utoronto.ca>
Date: Mon Jan 27 21:30:43 2014 -0500
doc: Spelling and grammar fixes.
>---------------------------------------------------------------
700dd71314d927d1f40349898ecf2e050d244ee4
doc/USER_MANUAL/manual_SPECFEM2D.pdf | Bin 3181256 -> 3181271 bytes
doc/USER_MANUAL/manual_SPECFEM2D.tex | 40 +++++++++++++++++------------------
2 files changed, 20 insertions(+), 20 deletions(-)
diff --git a/doc/USER_MANUAL/manual_SPECFEM2D.pdf b/doc/USER_MANUAL/manual_SPECFEM2D.pdf
index 616cded..50d0f3b 100644
Binary files a/doc/USER_MANUAL/manual_SPECFEM2D.pdf and b/doc/USER_MANUAL/manual_SPECFEM2D.pdf differ
diff --git a/doc/USER_MANUAL/manual_SPECFEM2D.tex b/doc/USER_MANUAL/manual_SPECFEM2D.tex
index b1cfb1b..fc36764 100644
--- a/doc/USER_MANUAL/manual_SPECFEM2D.tex
+++ b/doc/USER_MANUAL/manual_SPECFEM2D.tex
@@ -179,7 +179,7 @@ both in serial and in parallel. See SPECFEM2D
\urlwithparentheses{http://www.geodynamics.org/cig/software/packages/seismo/specfem2d}
for the source code.
-The SEM is a continuous Galerkin technique \citep{TrKoLi08,PeKoLuMaLeCaLeMaLiBlNiBaTr11}, which can easily be made discontinous \citep{BeMaPa94,Ch00,KoWoHu02,ChCaVi03,LaWaBe05,Kop06,WiStBuGh10,AcKo11}; it is then close to a particular case of the discontinuous Galerkin technique \citep{ReHi73,LeRa74,Arn82,JoPi86,BoMaHe91,FaRi99,HuHuRa99,CoKaSh00,GiHeWa02,RiWh03,MoRi05,GrScSc06,AiMoMu06,BeLaPi06,DuKa06,DeSeWh08,PuAmKa09,WiStBuGh10,DeSe10,EtChViGl10}, with optimized efficiency because of its tensorized basis functions \citep{WiStBuGh10,AcKo11}.
+The SEM is a continuous Galerkin technique \citep{TrKoLi08,PeKoLuMaLeCaLeMaLiBlNiBaTr11}, which can easily be made discontinuous \citep{BeMaPa94,Ch00,KoWoHu02,ChCaVi03,LaWaBe05,Kop06,WiStBuGh10,AcKo11}; it is then close to a particular case of the discontinuous Galerkin technique \citep{ReHi73,LeRa74,Arn82,JoPi86,BoMaHe91,FaRi99,HuHuRa99,CoKaSh00,GiHeWa02,RiWh03,MoRi05,GrScSc06,AiMoMu06,BeLaPi06,DuKa06,DeSeWh08,PuAmKa09,WiStBuGh10,DeSe10,EtChViGl10}, with optimized efficiency because of its tensorized basis functions \citep{WiStBuGh10,AcKo11}.
In particular, it can accurately handle very distorted mesh elements \citep{OlSe11}.
It has very good accuracy and convergence properties \citep{MaPa89,SePr94,DeFiMu02,Coh02,DeSe07,SeOl08,AiWa09,AiWa10,MeStTh12}.
@@ -204,7 +204,7 @@ and has been successfully adapted to address problems in seismic wave propagatio
Early seismic wave propagation applications of the SEM, utilizing Legendre basis functions and a
perfectly diagonal mass matrix, include \cite{CoJoTo93}, \cite{Kom97},
\cite{FaMaPaQu97}, \cite{CaGa97}, \cite{KoVi98} and \cite{KoTr99},
-whereas applications involving Chebyshev basis functions and a nondiagonal mass matrix
+whereas applications involving Chebyshev basis functions and a non-diagonal mass matrix
include \cite{SePr94}, \cite{PrCaSe94} and \cite{SePrPr95}.
All SPECFEM2D software is written in Fortran2003 with full portability
@@ -216,7 +216,7 @@ parallel programming based upon the Message Passing Interface (MPI)
The next release of the code will include support for GPU graphics card acceleration \citep{KoMiEr09,KoErGoMi10,MiKo10,Kom11}.
The code uses the plane strain convention when the standard P-SV equation case is used, i.e.,
-the off-plane strain $\epsilon_{zz}$ is zero by definition of the plane strain convention but thus the off-plane stress $\sigma_{zz}$ is not equal to zero,
+the off-plane strain $\epsilon_{zz}$ is zero by definition of the plane strain convention but the off-plane stress $\sigma_{zz}$ is not equal to zero,
one has $\sigma_{zz} = \lambda (\epsilon_{xx} + \epsilon_{yy})$.
%------------------------------------------------------------------------------------------------%
@@ -327,7 +327,7 @@ You will need to replace it with \texttt{clock\_get\_time} if you want to use \t
The SPECFEM2D software package relies on the SCOTCH library to partition meshes.
The SCOTCH library \citep{PeRo96}
-provides efficent static mapping, graph and mesh partitioning routines. SCOTCH is a free software package developed by
+provides efficient static mapping, graph and mesh partitioning routines. SCOTCH is a free software package developed by
Fran\c{c}ois Pellegrini et al. from LaBRI and Inria in Bordeaux, France, downloadable from the web page \url{https://gforge.inria.fr/projects/scotch/}.
In case no SCOTCH libraries can be found on the system, the configuration will bundle the version provided with the source code for compilation.
The path to an existing SCOTCH installation can to be set explicitly with the option \texttt{-{}-with-scotch-dir}.
@@ -406,7 +406,7 @@ but you can find useful information in the manuals of the 3D versions, since man
\urlwithparentheses{http://geodynamics.org/wsvn/cig/seismo/3D} in subdirectories \texttt{USER\_MANUAL/}.
To create acoustic (fluid) regions, just set the S wave speed to zero and the code will see that these elements are fluid and switch to the right equations there automatically, and automatically match them with the solid regions
-\item if you are using an external mesher (like GID or CUBIT / Trelis), you should set \texttt{read\_external\_mesh} to \texttt{.true.}:
+\item if you are using an external mesher (like GiD or CUBIT / Trelis), you should set \texttt{read\_external\_mesh} to \texttt{.true.}:
\begin{description}
\item[\texttt{mesh\_file}] is the file describing the mesh : first line is the number of elements, then a list of 4 nodes (quadrilaterals only) forming each elements on each line.
@@ -425,7 +425,7 @@ If one of your elements has more than one edge along a given absorbing contour
putting the first edge on the first line and the second edge on the second line.
Do not list the same element with the same absorbing edge twice or more, otherwise absorption will not be correct because the edge integral
will be improperly subtracted several times.
-If one of your elements has a single point along the absording contour rather than a full edge, do NOT list it
+If one of your elements has a single point along the absorbing contour rather than a full edge, do NOT list it
(it would have no weight in the contour integral anyway because it would consist of a single point).
If you use 9-node elements, list only the first and last points of the edge and not the intermediate point
located around the middle of the edge; the right 9-node curvature will be restored automatically by the code.
@@ -482,7 +482,7 @@ for generating meshes which can be processed by SPECFEM2D. Only two
modules of Gmsh are of interest for the SPECFEM2D users : the geometry
and the mesh modules. An example is given in directory \texttt{EXAMPLES/Gmsh\_example}
which illustrates the generation of an external mesh using these two
-modules. The model, which is considered, consists of a homogeneous
+modules. The model that is considered consists of a homogeneous
square containing two circles filled with a different material.
The geometry is generated by loading file \texttt{SqrCirc.geo} into
@@ -694,7 +694,7 @@ To run the solver, type:
./bin/xspecfem2D
\end{verbatim}
%
-to run the main solver (use \texttt{mpirun} or equivalent if you compiled with parallel support). This will output the seismograms and snapshots of the wave fronts at different time steps in directory \texttt{OUTPUT\_FILES/}. To visualize them, type "\texttt{gs OUTPUT\_FILES/vect*.ps}" to see the Postscript files (in which the wave field is represented with small arrows, fluid/solid matching interfaces with a thick pink line, and absorbing edges with a thick green line) and "\texttt{gimp OUTPUT\_FILES/image*.gif}" to see the color snapshot showing a pixelized image of one of the two components of the wave field (or pressure, depending on what you have selected for the output in \texttt{DATA/Par\_file}).
+to run the main solver (use \texttt{mpirun} or equivalent if you compiled with parallel support). This will output the seismograms and snapshots of the wave fronts at different time steps in directory \texttt{OUTPUT\_FILES/}. To visualize them, type "\texttt{gs OUTPUT\_FILES/vect*.ps}" to see the Postscript files (in which the wave field is represented with small arrows, fluid/solid matching interfaces with a thick pink line, and absorbing edges with a thick green line) and "\texttt{gimp OUTPUT\_FILES/image*.gif}" to see the colour snapshot showing a pixelized image of one of the two components of the wave field (or pressure, depending on what you have selected for the output in \texttt{DATA/Par\_file}).
%%
\begin{figure}[htbp]
@@ -856,7 +856,7 @@ To run a P-SV waves calculation propagating in the x-z plane,
set \texttt{p\_sv = .true.} in the \texttt{Par\_file}.
\item[SH:]
-To run a SH (membrane) waves calculation traveling in the x-z plane with a
+To run a SH (membrane) waves calculation travelling in the x-z plane with a
y-component of motion, set \texttt{p\_sv = .false.}
\end{description}
@@ -940,8 +940,8 @@ see \cite{MoTr08}.
(3) adjoint method and kernels calculation
\item In section \textbf{"\# source parameters"}:\\
-The code now support multi sources.
-\texttt{NSOURCE} is the number of source.
+The code now support multiple sources.
+\texttt{NSOURCE} is the number of sources.
Parameters of the sources are displayed in the file \texttt{SOURCE}, which must be
in the directory \texttt{DATA/}. The components of a moment tensor source must be given in N.m,
not in dyne.cm as in the \texttt{DATA/CMTSOLUTION} source file of the 3D version of the code.
@@ -969,7 +969,7 @@ There are three possible types of models:\\
\texttt{II}: (model\_number 2 rho c11 c13 c15 c33 c35 c55 c12 c23 c25 0 0 0) or\\
\texttt{III}: (model\_number 3 rhos rhof phi c kxx kxz kzz Ks Kf Kfr etaf mufr Qmu). \\
-For istropic elastic/acoustic material use \texttt{I} and set Vs to zero to make a given model acoustic, for anisotropic elastic use \texttt{II},
+For isotropic elastic/acoustic material use \texttt{I} and set Vs to zero to make a given model acoustic, for anisotropic elastic use \texttt{II},
and for isotropic poroelastic material use \texttt{III}. The mesh can contain acoustic, elastic, and poroelastic models simultaneously.
For anisotropic elastic media the last three parameters, c12 c23 c25, are used only when the user asks the code to compute pressure for display
@@ -1004,7 +1004,7 @@ permeability.
%------------------------------------------------------------------------------------------------%
\section{How to set plane waves as initial conditions}
-To simulate propagation of incoming plane waves in the simulation domain, initial conditions based on analytical formulae of plane waves in homogenous model need to be set. No additional body or boundary forces are required. To set up this senario:
+To simulate propagation of incoming plane waves in the simulation domain, initial conditions based on analytical formulae of plane waves in homogeneous model need to be set. No additional body or boundary forces are required. To set up this scenario:
%
\begin{description}
\item{\verb+Par_file+:}
@@ -1027,8 +1027,8 @@ To simulate propagation of incoming plane waves in the simulation domain, initia
Three different explicit conditionally-stable time schemes can be used for elastic, acoustic (fluid) or coupled elastic/acoustic media:
the Newmark method, the low-dissipation and low-dispersion fourth-order six-stage Runge-Kutta method (LDDRK4-6) presented in \cite{BeBoBa06},
and the classical fourth-order four-stage Runge-Kutta (RK4) method.
-Currently the last two methods are not implemented for poro-elastic media.
-According to \cite{DeSe10} and \cite{BeBoBa06}, with different degrees $N=NGLLX-1$ of the GLL basis funtions the CFL bounds are given in the following tables.
+Currently the last two methods are not implemented for poroelastic media.
+According to \cite{DeSe10} and \cite{BeBoBa06}, with different degrees $N=NGLLX-1$ of the GLL basis functions the CFL bounds are given in the following tables.
Note that by default the SPECFEM solver uses $NGLLX = 5$ and thus a degree $N = 4$, which is thus the value you should use
in most cases in the following tables.
You can directly compare these values with the value given in sentence 'Max stability for P wave velocity' in file
@@ -1121,8 +1121,8 @@ Degree $N$ & Newmark & LDDRK4-6 & RK4 \\ [0.5ex]
\item \texttt{SAVE\_FORWARD = .true.}
\item \texttt{seismotype = 1} (we need to save the displacement fields to later on derive the
adjoint source. Note: if the user forgets it, the program corrects it when reading the proper
-\texttt{SIMULATION\_TYPE} and \texttt{SAVE\_FORWARD} combination and a warning message appears in the ouput
-file)
+\texttt{SIMULATION\_TYPE} and \texttt{SAVE\_FORWARD} combination and a warning
+message appears in the output file)
\end{itemize}
Important output files (for example, for the elastic case, P-SV waves):
@@ -1215,14 +1215,14 @@ poroelastic case.
%------------------------------------------------------------------------------------------------%
-The SPECFEM2D package provides compatibilities in industrial (oil and gas industry) types of simulations.
+The SPECFEM2D package provides compatibility with industrial (oil and gas industry) types of simulations.
These features include importing Seismic Unix (SU) format wavespeed models into SPECFEM2D,
-outputing seismograms also in SU format with a few key parameters defined in the trace headers
+output of seismograms in SU format with a few key parameters defined in the trace headers
and reading adjoint sources in SU format etc.
There is one example given in EXAMPLES/INDUSTRIAL\_FORMAT, which you can follow.
We also changed the relationship between adjoint potential and adjoint displacement in fluid region
-(the relationship between forward potential and forward displacement remains the same as previsouly defined).
+(the relationship between forward potential and forward displacement remains the same as previously defined).
The new definition is critical when there are adjoint sources (in other words, receivers) in the acoustic domain,
and is the direct consequence of the optimization problem.
%
More information about the CIG-COMMITS
mailing list