[cig-commits] [commit] release_1.2: Clarifications on Instaseis databases (d75b01c)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/axisem

On branch  : release_1.2
Link       : https://github.com/geodynamics/axisem/compare/ec5ac1f7d2b600f32d2c7d315ba9d7161b51fca3...d75b01c37b7729b58e4e6b0fc36db2eefb92d0e2

>---------------------------------------------------------------

commit d75b01c37b7729b58e4e6b0fc36db2eefb92d0e2
Author: Simon Stähler <staehler at geophysik.uni-muenchen.de>
Date:   Wed Jan 14 17:55:45 2015 +0100

    Clarifications on Instaseis databases


>---------------------------------------------------------------

d75b01c37b7729b58e4e6b0fc36db2eefb92d0e2
 MANUAL/manual_axisem1.2.pdf | Bin 1927600 -> 1927838 bytes
 MANUAL/manual_axisem1.2.tex |  15 ++++++++++-----
 2 files changed, 10 insertions(+), 5 deletions(-)

diff --git a/MANUAL/manual_axisem1.2.pdf b/MANUAL/manual_axisem1.2.pdf
index 7538821..22007d5 100644
Binary files a/MANUAL/manual_axisem1.2.pdf and b/MANUAL/manual_axisem1.2.pdf differ
diff --git a/MANUAL/manual_axisem1.2.tex b/MANUAL/manual_axisem1.2.tex
index 23721d8..b7c5453 100644
--- a/MANUAL/manual_axisem1.2.tex
+++ b/MANUAL/manual_axisem1.2.tex
@@ -304,12 +304,11 @@ element approach.\\
 
 AxiSEM allows to output larger datasets, especially wavefields in the NetCDF format. \\
 Unfortunately, the installation of the NetCDF libraries may still be cumbersome; below we
-list some limited support for installing NetCDF..
+list some limited support for installing NetCDF.
 
 \begin{description}
     \item[HPC systems:] The libraries should be provided by the system. Use the
-          recommended settings. AxiSEM has worked with NetCDF versions down to 4.1.2, but 
-          unpredictable errors still occur.
+          recommended settings. AxiSEM has worked reliably with NetCDF versions starting from 4.2. 
           If you want to use the new Parallel NetCDF output, choose a package named netcdf-mpi 
           or similar.
     \item[Ubuntu 12.10 and newer, Debian Wheezy and newer]: The code is working with the NetCDF libraries delivered
@@ -868,7 +867,7 @@ By default, the sampling rate is set to the time step length of the simulation.
 strongly recommend to leave it as such to avoid aliasing. The resampling can better be
 done with \textit{ObsPy} or another tool that supports filtering.
 
-\subsection{Output wavefields in netcdf format needed for \textit{Instaseis}}
+\subsection{\textit{Instaseis}-ready wavefield output in netcdf format}
 
 \begin{itemize}
     \item Compile with netcdf (see section~\ref{sssec:netcdf} for hints on the
@@ -902,13 +901,19 @@ KERNEL_RMIN         5671.
 KERNEL_RMAX         7000.
 \end{verbatim}
 
+
+    \item For high frequencies and propagation lengths over 100 wavelengths, the 
+    symplectic time scheme is recommended to keep the dispersion errors small 
+    (compare benchmarks in van Driel \& Nissen-Meyer 2014a,b):\\
+    In \verb|SOLVER/inparam_advanced| set \verb|TIME_SCHEME| to \verb|symplec4|. 
+
     \item in \verb|SOLVER/inparam_basic| set \verb|SIMULATION_TYPE| to \verb|force| to
     simulate both the vertical and the horizontal database. To generate just one of the
     two, set it to \verb|single| and select either \verb|thetaforce| or \verb|vertforce|
     in \\\verb|SOLVER/inparam_source|
 
     \item once the AxiSEM runs are finished, go to the rundirectory and run
-    \verb|./fieltransform.sh| (in case of single runs just run \verb|./fieltransform|
+    \verb|./fieldtransform.sh| (in case of single runs just run \verb|./fieldtransform|
     directly). This rechunks the wavefields in the netcdf file to optimize them for
     reading as time series instead of snapshots.