[cig-commits] [commit] devel,master: committed new version of Paul Cristini's Python script to convert old Par_files to the new format; and used it to convert old examples that had not been converted before (93d2f61)

cig_noreply at geodynamics.org cig_noreply at geodynamics.org
Wed Jun 18 15:20:51 PDT 2014


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

On branches: devel,master
Link       : https://github.com/geodynamics/specfem2d/compare/fc67e6fd7ad890705b2b72b4b3c509accb22249e...e9ca46c40131588d89d7b0883250bc6584ce6b4c

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

commit 93d2f6166c7d7fe9daea30d8f6b4ccbfe0981fb5
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Thu Jan 12 23:31:23 2012 +0000

    committed new version of Paul Cristini's Python script to convert old Par_files to the new format; and used it to convert old examples that had not been converted before


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

93d2f6166c7d7fe9daea30d8f6b4ccbfe0981fb5
 noise_layered/model_0/Par_file_fair | 19 ++++++++++++++++++-
 noise_layered/model_0/Par_file_good | 19 ++++++++++++++++++-
 noise_layered/model_1/Par_file_best | 19 ++++++++++++++++++-
 noise_layered/model_1/Par_file_fair | 19 ++++++++++++++++++-
 noise_layered/model_1/Par_file_good | 19 ++++++++++++++++++-
 noise_layered/model_2/Par_file_fair | 19 ++++++++++++++++++-
 noise_layered/model_2/Par_file_good | 19 ++++++++++++++++++-
 7 files changed, 126 insertions(+), 7 deletions(-)

diff --git a/noise_layered/model_0/Par_file_fair b/noise_layered/model_0/Par_file_fair
index 6eba8af..03c40f8 100644
--- a/noise_layered/model_0/Par_file_fair
+++ b/noise_layered/model_0/Par_file_fair
@@ -9,6 +9,7 @@ SAVE_FORWARD                    = .false.
 # parameters concerning partitioning
 nproc                           = 1              # number of processes
 partitioning_method             = 3              # SCOTCH = 3, ascending order (very bad idea) = 1
+PERFORM_CUTHILL_MCKEE           = .true.         # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
 
 ngnod                           = 9              # number of control nodes per element (4 or 9)
 initialfield                    = .false.        # use a plane wave as source or not
@@ -24,6 +25,7 @@ p_sv                            = .false.        # set the type of calculation (
 # time step parameters
 nt                              = 3000           # total number of time steps
 deltat                          = 1.d-3          # duration of a time step
+USER_T0                         = 0.0d0          # use this t0 as earliest starting time rather than the automatically calculated one
 
 # source parameters
 NSOURCES                        = 1              # number of sources [source info read in CMTSOLUTION file]
@@ -39,6 +41,7 @@ generate_STATIONS               = .true.         # creates a STATION file in ./D
 nreceiversets                   = 1              # number of receiver sets
 anglerec                        = 0.d0           # angle to rotate components at receivers
 rec_normal_to_surface           = .false.        # base anglerec normal to surface (external mesh and curve file needed)
+SU_FORMAT                       = .false.        # output seismograms in Seismic Unix format (adjoint traces will be read in the same format)
 
 # first receiver set (repeat these 6 lines and adjust nreceiversets  accordingly)
 nrec                            = 29
@@ -59,8 +62,13 @@ modelvect                       = .false.        # display velocity model on vec
 boundvect                       = .true.         # display boundary conditions on plots
 interpol                        = .true.         # interpolation of the display or not
 pointsdisp                      = 6              # points for interpolation of display (set to 1 for lower-left corner only)
-subsamp                         = 1              # subsampling of color snapshots
+subsamp_postscript              = 1              # subsampling of color snapshots
+factor_subsample_image          = 1              # factor to subsample color images output by the code (useful for very large models)
+POWER_DISPLAY_COLOR             = 0.30d0         # non linear display to enhance small amplitudes in color images
+DRAW_WATER_CONSTANT_BLUE_IN_JPG = .true.         # display acoustic layers as constant blue in JPEG images, because they likely correspond to water
 sizemax_arrows                  = 1.d0           # maximum size of arrows on vector plots in cm
+US_LETTER                       = .false.        # US letter paper or European A4
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false.        # use snapshot number in the file name of JPEG color snapshots instead of the time step
 gnuplot                         = .false.        # generate a GNUPLOT file for the grid
 output_grid                     = .false.        # save the grid in a text file or not
 output_energy                   = .false.        # compute and output acoustic and elastic energy (slows down the code significantly)
@@ -82,6 +90,15 @@ read_external_mesh              = .false.
 # absorbing boundary active or not
 absorbing_conditions            = .true.
 
+# for horizontal periodic conditions: detect common points between left and right edges
+ADD_PERIODIC_CONDITIONS         = .false.
+
+# horizontal periodicity distance for periodic conditions
+PERIODIC_horiz_dist             = 0.3597d0
+
+# grid point detection tolerance for periodic conditions
+PERIODIC_DETECT_TOL             = 3.3334d-6
+
 #-----------------------------------------------------------------------------
 # PARAMETERS FOR EXTERNAL MESHING
 
diff --git a/noise_layered/model_0/Par_file_good b/noise_layered/model_0/Par_file_good
index 8743a9c..4c3edd2 100644
--- a/noise_layered/model_0/Par_file_good
+++ b/noise_layered/model_0/Par_file_good
@@ -9,6 +9,7 @@ SAVE_FORWARD                    = .false.
 # parameters concerning partitioning
 nproc                           = 1              # number of processes
 partitioning_method             = 3              # SCOTCH = 3, ascending order (very bad idea) = 1
+PERFORM_CUTHILL_MCKEE           = .true.         # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
 
 ngnod                           = 9              # number of control nodes per element (4 or 9)
 initialfield                    = .false.        # use a plane wave as source or not
@@ -24,6 +25,7 @@ p_sv                            = .false.        # set the type of calculation (
 # time step parameters
 nt                              = 3000           # total number of time steps
 deltat                          = 1.d-3          # duration of a time step
+USER_T0                         = 0.0d0          # use this t0 as earliest starting time rather than the automatically calculated one
 
 # source parameters
 NSOURCES                        = 1              # number of sources [source info read in CMTSOLUTION file]
@@ -39,6 +41,7 @@ generate_STATIONS               = .true.         # creates a STATION file in ./D
 nreceiversets                   = 1              # number of receiver sets
 anglerec                        = 0.d0           # angle to rotate components at receivers
 rec_normal_to_surface           = .false.        # base anglerec normal to surface (external mesh and curve file needed)
+SU_FORMAT                       = .false.        # output seismograms in Seismic Unix format (adjoint traces will be read in the same format)
 
 # first receiver set (repeat these 6 lines and adjust nreceiversets  accordingly)
 nrec                            = 29
@@ -59,8 +62,13 @@ modelvect                       = .false.        # display velocity model on vec
 boundvect                       = .true.         # display boundary conditions on plots
 interpol                        = .true.         # interpolation of the display or not
 pointsdisp                      = 6              # points for interpolation of display (set to 1 for lower-left corner only)
-subsamp                         = 1              # subsampling of color snapshots
+subsamp_postscript              = 1              # subsampling of color snapshots
+factor_subsample_image          = 1              # factor to subsample color images output by the code (useful for very large models)
+POWER_DISPLAY_COLOR             = 0.30d0         # non linear display to enhance small amplitudes in color images
+DRAW_WATER_CONSTANT_BLUE_IN_JPG = .true.         # display acoustic layers as constant blue in JPEG images, because they likely correspond to water
 sizemax_arrows                  = 1.d0           # maximum size of arrows on vector plots in cm
+US_LETTER                       = .false.        # US letter paper or European A4
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false.        # use snapshot number in the file name of JPEG color snapshots instead of the time step
 gnuplot                         = .false.        # generate a GNUPLOT file for the grid
 output_grid                     = .false.        # save the grid in a text file or not
 output_energy                   = .false.        # compute and output acoustic and elastic energy (slows down the code significantly)
@@ -82,6 +90,15 @@ read_external_mesh              = .false.
 # absorbing boundary active or not
 absorbing_conditions            = .true.
 
+# for horizontal periodic conditions: detect common points between left and right edges
+ADD_PERIODIC_CONDITIONS         = .false.
+
+# horizontal periodicity distance for periodic conditions
+PERIODIC_horiz_dist             = 0.3597d0
+
+# grid point detection tolerance for periodic conditions
+PERIODIC_DETECT_TOL             = 3.3334d-6
+
 #-----------------------------------------------------------------------------
 # PARAMETERS FOR EXTERNAL MESHING
 
diff --git a/noise_layered/model_1/Par_file_best b/noise_layered/model_1/Par_file_best
index 381ec26..902d952 100644
--- a/noise_layered/model_1/Par_file_best
+++ b/noise_layered/model_1/Par_file_best
@@ -9,6 +9,7 @@ SAVE_FORWARD                    = .false.
 # parameters concerning partitioning
 nproc                           = 1              # number of processes
 partitioning_method             = 3              # SCOTCH = 3, ascending order (very bad idea) = 1
+PERFORM_CUTHILL_MCKEE           = .true.         # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
 
 ngnod                           = 9              # number of control nodes per element (4 or 9)
 initialfield                    = .false.        # use a plane wave as source or not
@@ -24,6 +25,7 @@ p_sv                            = .false.        # set the type of calculation (
 # time step parameters
 nt                              = 3000           # total number of time steps
 deltat                          = 1.d-3          # duration of a time step
+USER_T0                         = 0.0d0          # use this t0 as earliest starting time rather than the automatically calculated one
 
 # source parameters
 NSOURCES                        = 1              # number of sources [source info read in CMTSOLUTION file]
@@ -39,6 +41,7 @@ generate_STATIONS               = .true.         # creates a STATION file in ./D
 nreceiversets                   = 1              # number of receiver sets
 anglerec                        = 0.d0           # angle to rotate components at receivers
 rec_normal_to_surface           = .false.        # base anglerec normal to surface (external mesh and curve file needed)
+SU_FORMAT                       = .false.        # output seismograms in Seismic Unix format (adjoint traces will be read in the same format)
 
 # first receiver set (repeat these 6 lines and adjust nreceiversets  accordingly)
 nrec                            = 29
@@ -59,8 +62,13 @@ modelvect                       = .false.        # display velocity model on vec
 boundvect                       = .true.         # display boundary conditions on plots
 interpol                        = .true.         # interpolation of the display or not
 pointsdisp                      = 6              # points for interpolation of display (set to 1 for lower-left corner only)
-subsamp                         = 1              # subsampling of color snapshots
+subsamp_postscript              = 1              # subsampling of color snapshots
+factor_subsample_image          = 1              # factor to subsample color images output by the code (useful for very large models)
+POWER_DISPLAY_COLOR             = 0.30d0         # non linear display to enhance small amplitudes in color images
+DRAW_WATER_CONSTANT_BLUE_IN_JPG = .true.         # display acoustic layers as constant blue in JPEG images, because they likely correspond to water
 sizemax_arrows                  = 1.d0           # maximum size of arrows on vector plots in cm
+US_LETTER                       = .false.        # US letter paper or European A4
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false.        # use snapshot number in the file name of JPEG color snapshots instead of the time step
 gnuplot                         = .false.        # generate a GNUPLOT file for the grid
 output_grid                     = .false.        # save the grid in a text file or not
 output_energy                   = .false.        # compute and output acoustic and elastic energy (slows down the code significantly)
@@ -83,6 +91,15 @@ read_external_mesh              = .false.
 # absorbing boundary active or not
 absorbing_conditions            = .true.
 
+# for horizontal periodic conditions: detect common points between left and right edges
+ADD_PERIODIC_CONDITIONS         = .false.
+
+# horizontal periodicity distance for periodic conditions
+PERIODIC_horiz_dist             = 0.3597d0
+
+# grid point detection tolerance for periodic conditions
+PERIODIC_DETECT_TOL             = 3.3334d-6
+
 #-----------------------------------------------------------------------------
 # PARAMETERS FOR EXTERNAL MESHING
 
diff --git a/noise_layered/model_1/Par_file_fair b/noise_layered/model_1/Par_file_fair
index 6a4ca29..b0107db 100644
--- a/noise_layered/model_1/Par_file_fair
+++ b/noise_layered/model_1/Par_file_fair
@@ -9,6 +9,7 @@ SAVE_FORWARD                    = .false.
 # parameters concerning partitioning
 nproc                           = 1              # number of processes
 partitioning_method             = 3              # SCOTCH = 3, ascending order (very bad idea) = 1
+PERFORM_CUTHILL_MCKEE           = .true.         # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
 
 ngnod                           = 9              # number of control nodes per element (4 or 9)
 initialfield                    = .false.        # use a plane wave as source or not
@@ -24,6 +25,7 @@ p_sv                            = .false.        # set the type of calculation (
 # time step parameters
 nt                              = 3000           # total number of time steps
 deltat                          = 1.d-3          # duration of a time step
+USER_T0                         = 0.0d0          # use this t0 as earliest starting time rather than the automatically calculated one
 
 # source parameters
 NSOURCES                        = 1              # number of sources [source info read in CMTSOLUTION file]
@@ -39,6 +41,7 @@ generate_STATIONS               = .true.         # creates a STATION file in ./D
 nreceiversets                   = 1              # number of receiver sets
 anglerec                        = 0.d0           # angle to rotate components at receivers
 rec_normal_to_surface           = .false.        # base anglerec normal to surface (external mesh and curve file needed)
+SU_FORMAT                       = .false.        # output seismograms in Seismic Unix format (adjoint traces will be read in the same format)
 
 # first receiver set (repeat these 6 lines and adjust nreceiversets  accordingly)
 nrec                            = 29
@@ -59,8 +62,13 @@ modelvect                       = .false.        # display velocity model on vec
 boundvect                       = .true.         # display boundary conditions on plots
 interpol                        = .true.         # interpolation of the display or not
 pointsdisp                      = 6              # points for interpolation of display (set to 1 for lower-left corner only)
-subsamp                         = 1              # subsampling of color snapshots
+subsamp_postscript              = 1              # subsampling of color snapshots
+factor_subsample_image          = 1              # factor to subsample color images output by the code (useful for very large models)
+POWER_DISPLAY_COLOR             = 0.30d0         # non linear display to enhance small amplitudes in color images
+DRAW_WATER_CONSTANT_BLUE_IN_JPG = .true.         # display acoustic layers as constant blue in JPEG images, because they likely correspond to water
 sizemax_arrows                  = 1.d0           # maximum size of arrows on vector plots in cm
+US_LETTER                       = .false.        # US letter paper or European A4
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false.        # use snapshot number in the file name of JPEG color snapshots instead of the time step
 gnuplot                         = .false.        # generate a GNUPLOT file for the grid
 output_grid                     = .false.        # save the grid in a text file or not
 output_energy                   = .false.        # compute and output acoustic and elastic energy (slows down the code significantly)
@@ -83,6 +91,15 @@ read_external_mesh              = .false.
 # absorbing boundary active or not
 absorbing_conditions            = .true.
 
+# for horizontal periodic conditions: detect common points between left and right edges
+ADD_PERIODIC_CONDITIONS         = .false.
+
+# horizontal periodicity distance for periodic conditions
+PERIODIC_horiz_dist             = 0.3597d0
+
+# grid point detection tolerance for periodic conditions
+PERIODIC_DETECT_TOL             = 3.3334d-6
+
 #-----------------------------------------------------------------------------
 # PARAMETERS FOR EXTERNAL MESHING
 
diff --git a/noise_layered/model_1/Par_file_good b/noise_layered/model_1/Par_file_good
index 918c7a5..b5a9730 100644
--- a/noise_layered/model_1/Par_file_good
+++ b/noise_layered/model_1/Par_file_good
@@ -9,6 +9,7 @@ SAVE_FORWARD                    = .false.
 # parameters concerning partitioning
 nproc                           = 1              # number of processes
 partitioning_method             = 3              # SCOTCH = 3, ascending order (very bad idea) = 1
+PERFORM_CUTHILL_MCKEE           = .true.         # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
 
 ngnod                           = 9              # number of control nodes per element (4 or 9)
 initialfield                    = .false.        # use a plane wave as source or not
@@ -24,6 +25,7 @@ p_sv                            = .false.        # set the type of calculation (
 # time step parameters
 nt                              = 3000           # total number of time steps
 deltat                          = 1.d-3          # duration of a time step
+USER_T0                         = 0.0d0          # use this t0 as earliest starting time rather than the automatically calculated one
 
 # source parameters
 NSOURCES                        = 1              # number of sources [source info read in CMTSOLUTION file]
@@ -39,6 +41,7 @@ generate_STATIONS               = .true.         # creates a STATION file in ./D
 nreceiversets                   = 1              # number of receiver sets
 anglerec                        = 0.d0           # angle to rotate components at receivers
 rec_normal_to_surface           = .false.        # base anglerec normal to surface (external mesh and curve file needed)
+SU_FORMAT                       = .false.        # output seismograms in Seismic Unix format (adjoint traces will be read in the same format)
 
 # first receiver set (repeat these 6 lines and adjust nreceiversets  accordingly)
 nrec                            = 29
@@ -59,8 +62,13 @@ modelvect                       = .false.        # display velocity model on vec
 boundvect                       = .true.         # display boundary conditions on plots
 interpol                        = .true.         # interpolation of the display or not
 pointsdisp                      = 6              # points for interpolation of display (set to 1 for lower-left corner only)
-subsamp                         = 1              # subsampling of color snapshots
+subsamp_postscript              = 1              # subsampling of color snapshots
+factor_subsample_image          = 1              # factor to subsample color images output by the code (useful for very large models)
+POWER_DISPLAY_COLOR             = 0.30d0         # non linear display to enhance small amplitudes in color images
+DRAW_WATER_CONSTANT_BLUE_IN_JPG = .true.         # display acoustic layers as constant blue in JPEG images, because they likely correspond to water
 sizemax_arrows                  = 1.d0           # maximum size of arrows on vector plots in cm
+US_LETTER                       = .false.        # US letter paper or European A4
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false.        # use snapshot number in the file name of JPEG color snapshots instead of the time step
 gnuplot                         = .false.        # generate a GNUPLOT file for the grid
 output_grid                     = .false.        # save the grid in a text file or not
 output_energy                   = .false.        # compute and output acoustic and elastic energy (slows down the code significantly)
@@ -83,6 +91,15 @@ read_external_mesh              = .false.
 # absorbing boundary active or not
 absorbing_conditions            = .true.
 
+# for horizontal periodic conditions: detect common points between left and right edges
+ADD_PERIODIC_CONDITIONS         = .false.
+
+# horizontal periodicity distance for periodic conditions
+PERIODIC_horiz_dist             = 0.3597d0
+
+# grid point detection tolerance for periodic conditions
+PERIODIC_DETECT_TOL             = 3.3334d-6
+
 #-----------------------------------------------------------------------------
 # PARAMETERS FOR EXTERNAL MESHING
 
diff --git a/noise_layered/model_2/Par_file_fair b/noise_layered/model_2/Par_file_fair
index aaeadde..d35603d 100644
--- a/noise_layered/model_2/Par_file_fair
+++ b/noise_layered/model_2/Par_file_fair
@@ -9,6 +9,7 @@ SAVE_FORWARD                    = .false.
 # parameters concerning partitioning
 nproc                           = 1              # number of processes
 partitioning_method             = 3              # SCOTCH = 3, ascending order (very bad idea) = 1
+PERFORM_CUTHILL_MCKEE           = .true.         # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
 
 ngnod                           = 9              # number of control nodes per element (4 or 9)
 initialfield                    = .false.        # use a plane wave as source or not
@@ -24,6 +25,7 @@ p_sv                            = .false.        # set the type of calculation (
 # time step parameters
 nt                              = 3000           # total number of time steps
 deltat                          = 1.d-3          # duration of a time step
+USER_T0                         = 0.0d0          # use this t0 as earliest starting time rather than the automatically calculated one
 
 # source parameters
 NSOURCES                        = 1              # number of sources [source info read in CMTSOLUTION file]
@@ -39,6 +41,7 @@ generate_STATIONS               = .true.         # creates a STATION file in ./D
 nreceiversets                   = 1              # number of receiver sets
 anglerec                        = 0.d0           # angle to rotate components at receivers
 rec_normal_to_surface           = .false.        # base anglerec normal to surface (external mesh and curve file needed)
+SU_FORMAT                       = .false.        # output seismograms in Seismic Unix format (adjoint traces will be read in the same format)
 
 # first receiver set (repeat these 6 lines and adjust nreceiversets  accordingly)
 nrec                            = 29             # number of receivers
@@ -59,8 +62,13 @@ modelvect                       = .false.        # display velocity model on vec
 boundvect                       = .true.         # display boundary conditions on plots
 interpol                        = .true.         # interpolation of the display or not
 pointsdisp                      = 6              # points for interpolation of display (set to 1 for lower-left corner only)
-subsamp                         = 1              # subsampling of color snapshots
+subsamp_postscript              = 1              # subsampling of color snapshots
+factor_subsample_image          = 1              # factor to subsample color images output by the code (useful for very large models)
+POWER_DISPLAY_COLOR             = 0.30d0         # non linear display to enhance small amplitudes in color images
+DRAW_WATER_CONSTANT_BLUE_IN_JPG = .true.         # display acoustic layers as constant blue in JPEG images, because they likely correspond to water
 sizemax_arrows                  = 1.d0           # maximum size of arrows on vector plots in cm
+US_LETTER                       = .false.        # US letter paper or European A4
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false.        # use snapshot number in the file name of JPEG color snapshots instead of the time step
 gnuplot                         = .false.        # generate a GNUPLOT file for the grid
 output_grid                     = .false.        # save the grid in a text file or not
 output_energy                   = .false.        # compute and output acoustic and elastic energy (slows down the code significantly)
@@ -84,6 +92,15 @@ read_external_mesh              = .false.
 # absorbing boundary active or not
 absorbing_conditions            = .true.
 
+# for horizontal periodic conditions: detect common points between left and right edges
+ADD_PERIODIC_CONDITIONS         = .false.
+
+# horizontal periodicity distance for periodic conditions
+PERIODIC_horiz_dist             = 0.3597d0
+
+# grid point detection tolerance for periodic conditions
+PERIODIC_DETECT_TOL             = 3.3334d-6
+
 #-----------------------------------------------------------------------------
 # PARAMETERS FOR EXTERNAL MESHING
 
diff --git a/noise_layered/model_2/Par_file_good b/noise_layered/model_2/Par_file_good
index 0673fc6..5102d0b 100644
--- a/noise_layered/model_2/Par_file_good
+++ b/noise_layered/model_2/Par_file_good
@@ -9,6 +9,7 @@ SAVE_FORWARD                    = .false.
 # parameters concerning partitioning
 nproc                           = 1              # number of processes
 partitioning_method             = 3              # SCOTCH = 3, ascending order (very bad idea) = 1
+PERFORM_CUTHILL_MCKEE           = .true.         # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
 
 ngnod                           = 9              # number of control nodes per element (4 or 9)
 initialfield                    = .false.        # use a plane wave as source or not
@@ -24,6 +25,7 @@ p_sv                            = .false.        # set the type of calculation (
 # time step parameters
 nt                              = 3000           # total number of time steps
 deltat                          = 1.d-3          # duration of a time step
+USER_T0                         = 0.0d0          # use this t0 as earliest starting time rather than the automatically calculated one
 
 # source parameters
 NSOURCES                        = 1              # number of sources [source info read in CMTSOLUTION file]
@@ -39,6 +41,7 @@ generate_STATIONS               = .true.         # creates a STATION file in ./D
 nreceiversets                   = 1              # number of receiver sets
 anglerec                        = 0.d0           # angle to rotate components at receivers
 rec_normal_to_surface           = .false.        # base anglerec normal to surface (external mesh and curve file needed)
+SU_FORMAT                       = .false.        # output seismograms in Seismic Unix format (adjoint traces will be read in the same format)
 
 # first receiver set (repeat these 6 lines and adjust nreceiversets  accordingly)
 nrec                            = 29             # number of receivers
@@ -59,8 +62,13 @@ modelvect                       = .false.        # display velocity model on vec
 boundvect                       = .true.         # display boundary conditions on plots
 interpol                        = .true.         # interpolation of the display or not
 pointsdisp                      = 6              # points for interpolation of display (set to 1 for lower-left corner only)
-subsamp                         = 1              # subsampling of color snapshots
+subsamp_postscript              = 1              # subsampling of color snapshots
+factor_subsample_image          = 1              # factor to subsample color images output by the code (useful for very large models)
+POWER_DISPLAY_COLOR             = 0.30d0         # non linear display to enhance small amplitudes in color images
+DRAW_WATER_CONSTANT_BLUE_IN_JPG = .true.         # display acoustic layers as constant blue in JPEG images, because they likely correspond to water
 sizemax_arrows                  = 1.d0           # maximum size of arrows on vector plots in cm
+US_LETTER                       = .false.        # US letter paper or European A4
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false.        # use snapshot number in the file name of JPEG color snapshots instead of the time step
 gnuplot                         = .false.        # generate a GNUPLOT file for the grid
 output_grid                     = .false.        # save the grid in a text file or not
 output_energy                   = .false.        # compute and output acoustic and elastic energy (slows down the code significantly)
@@ -84,6 +92,15 @@ read_external_mesh              = .false.
 # absorbing boundary active or not
 absorbing_conditions            = .true.
 
+# for horizontal periodic conditions: detect common points between left and right edges
+ADD_PERIODIC_CONDITIONS         = .false.
+
+# horizontal periodicity distance for periodic conditions
+PERIODIC_horiz_dist             = 0.3597d0
+
+# grid point detection tolerance for periodic conditions
+PERIODIC_DETECT_TOL             = 3.3334d-6
+
 #-----------------------------------------------------------------------------
 # PARAMETERS FOR EXTERNAL MESHING
 



More information about the CIG-COMMITS mailing list