Thu Nov 6 08:02:01 PST 2014

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

On branches: devel,master
Link       : https://github.com/geodynamics/specfem3d_globe/compare/bc58e579b3b0838a0968725a076f5904845437ca...be63f20cbb6f462104e949894dbe205d2398cd7f

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

commit 3fcd5f07861110f0293943ff040cea621114215c
Author: Bernhard Schubert <mail at bernhard-schuberth.de>
Date:   Wed Jul 4 11:11:01 2007 +0000

    - added possibility of writing SAC binary format in write_seismograms.f90
      this involved adding a new function in write_c_binary.c
    - added flags OUTPUT_SEISMOS_ASCII_TEXT, OUTPUT_SEISMOS_SAC_ALPHANUM,
      and OUTPUT_SEISMOS_SAC_BINARY to Par_file, which are used to decide
      on the output format
    - added flag ROTATE_SEISMOGRAMS_RT which is used to decide whether the
      horizontal components of the seismograms are rotated to radial and
      transverse directions
    - added file get_backazimuth.f90 which computes the backazimuth from
      source and receiver coordinates needed to rotate the horizontal
      components
    - added get_backazimuth.f90 as object and target in Makefile.in
    - changed read_compute_parameters.f90 according to the new flags in Par_file
    - updated all calls to read_compute_parameters.f90 and the respective
      MPI_BAST call in specfem3D.f90 (not in meshfem3D.f90, the new flags are
      not used there)
    - Added information on new output flags to the manual


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

3fcd5f07861110f0293943ff040cea621114215c
 Par_file     | 20 ++++++++++++++------
 Par_file.cfg |  8 ++++++++
 2 files changed, 22 insertions(+), 6 deletions(-)

diff --git a/Par_file b/Par_file
index 1a4c14c..1e5f687 100644
--- a/Par_file
+++ b/Par_file
@@ -15,12 +15,12 @@ GAMMA_ROTATION_AZIMUTH        = 20.d0
 
 # number of elements at the surface along the two sides of the first chunk
 # (must be 16 * multiple of NPROC below)
-NEX_XI                          = 160 
-NEX_ETA                         = 160
+NEX_XI                          = 64  
+NEX_ETA                         = 64 
 
 # number of MPI processors along the two sides of the first chunk
-NPROC_XI                        = 5
-NPROC_ETA                       = 5
+NPROC_XI                        = 1
+NPROC_ETA                       = 1
 
 # 1D models with real structure:
 # 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135
@@ -44,7 +44,7 @@ ATTENUATION                     = .false.
 ABSORBING_CONDITIONS            = .false.
 
 # record length in minutes
-RECORD_LENGTH_IN_MINUTES        = 40.0d0
+RECORD_LENGTH_IN_MINUTES        = 1.0d0
 
 # save AVS or OpenDX movies
 MOVIE_SURFACE                   = .false.
@@ -60,7 +60,7 @@ NUMBER_OF_RUNS                  = 1
 NUMBER_OF_THIS_RUN              = 1
 
 # path to store the local database files on each node
-LOCAL_PATH                      = /scratch/komatits/DATABASES_MPI
+LOCAL_PATH                      = /import/tethys-data/bernhard/SEM/TEST_SPECFEM3D
 
 # interval at which we output time step info and max of norm of displacement
 NTSTEP_BETWEEN_OUTPUT_INFO      = 100
@@ -69,6 +69,14 @@ NTSTEP_BETWEEN_OUTPUT_INFO      = 100
 NTSTEP_BETWEEN_OUTPUT_SEISMOS   = 5000000
 NTSTEP_BETWEEN_READ_ADJSRC      = 1000
 
+# output format for the seismograms (one can use either or all of the three formats)
+OUTPUT_SEISMOS_ASCII_TEXT Â      = .true.
+OUTPUT_SEISMOS_SAC_ALPHANUM     = .false.
+OUTPUT_SEISMOS_SAC_BINARY Â      = .false.
+
+# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
+ROTATE_SEISMOGRAMS_RT           = .false.
+
 # flag to impose receivers at the surface or allow them to be buried
 RECEIVERS_CAN_BE_BURIED         = .false.
 
diff --git a/Par_file.cfg b/Par_file.cfg
index 0eac3d8..80762f6 100644
--- a/Par_file.cfg
+++ b/Par_file.cfg
@@ -86,6 +86,14 @@ ntstep-between-output-info = 200
 # interval in time steps for temporary writing of seismograms
 ntstep-between-output-seismos = 5000000
 
+# output format for the seismograms (one can use either or all of the three formats)
+output-seismos-ascii-text       = .true.
+output-seismos-sac-alphanum     = .false.
+output-seismos-sac-binary       = .false.
+
+# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
+rotate-seimograms-rt           = .false.
+
 # flag to impose receivers at the surface or allow them to be buried
 receivers-can-be-buried = False
 

