[cig-commits] [commit] devel, master: Modified the way MOVIE_VOLUME is done. Added several new parameters to the Par_file to allow saving movie in sub volume. Moved all movie writing to subroutines in write_movie_volume.f90 (similarly to write_seismograms). Added new parameter MOVIE_VOLUME_TYPE to allow for adding new output values/formats without modifying par_file and by simply adding new sub routines. (b5b14ee)

[cig_noreply at geodynamics.org]

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

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

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

commit b5b14eef42408a759ee95071487834bee30417a3
Author: Vala Hjörleifsdóttir <vala at gps.caltech.edu>
Date:   Sat Dec 22 02:27:39 2007 +0000

    Modified the way MOVIE_VOLUME is done.  Added several new parameters
    to the Par_file to allow saving movie in sub volume.  Moved all movie writing
    to subroutines in write_movie_volume.f90 (similarly to write_seismograms).  Added new parameter
    MOVIE_VOLUME_TYPE to allow for adding new output values/formats without
    modifying par_file and by simply adding new sub routines.
    
    Also added combine_paraview_strain_data to allow for combining the files that
    are not stored in whole slices, but rather subvolumes of slices.


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

b5b14eef42408a759ee95071487834bee30417a3
 Par_file | 24 +++++++++++++++++++++---
 1 file changed, 21 insertions(+), 3 deletions(-)

diff --git a/Par_file b/Par_file
index 21f7785..f522264 100644
--- a/Par_file
+++ b/Par_file
@@ -53,6 +53,24 @@ MOVIE_VOLUME                    = .false.
 NTSTEP_BETWEEN_FRAMES           = 100
 HDUR_MOVIE                      = 0.d0
 
+# save movie in volume.  Will save element if center of element is in prescribed volume
+# top/bottom: depth in KM, use MOVIE_TOP = -100 to make sure the surface is stored.
+# west/east: longitude, degrees East [-180/180] top/bottom: latitute, degrees North [-90/90]
+# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
+# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
+# type 4 saves the trace and deviatoric stress in the whole volume
+#MOVIE_VOLUME_COARSE saves movie only at corners of elements
+MOVIE_VOLUME_TYPE               = 2
+MOVIE_VOLUME_COARSE             = .true.
+MOVIE_TOP_KM                    = -100.0
+MOVIE_BOTTOM_KM                 = 1000.0
+MOVIE_WEST_DEG                  = -90.0
+MOVIE_EAST_DEG                  = 90.0
+MOVIE_NORTH_DEG                 = 90.0
+MOVIE_SOUTH_DEG                 = -90.0
+MOVIE_START                     = 0
+MOVIE_STOP                      = 40000
+
 # save mesh files to check the mesh
 SAVE_MESH_FILES                 = .false.
 
@@ -72,18 +90,18 @@ 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     = .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           = .true.
+ROTATE_SEISMOGRAMS_RT           = .false.
 
 # decide wheter master process writes all the seismograms or if all processes do it in parallel
 WRITE_SEISMOGRAMS_BY_MASTER     = .true.
 
 # save all seismograms in one large combined file instead of one file per seismogram
 # to avoid overloading shared non-local file systems such as GPFS for instance
-SAVE_ALL_SEISMOS_IN_ONE_FILE    = .true.
+SAVE_ALL_SEISMOS_IN_ONE_FILE    = .false.
 USE_BINARY_FOR_LARGE_FILE       = .false.
 
 # flag to impose receivers at the surface or allow them to be buried