[cig-commits] r19207 - in seismo/2D/SPECFEM2D/trunk/EXAMPLES: Abel_Brest Gmsh_example_MPI Gmsh_example_serial INDUSTRIAL_FORMAT M2_UPPA Tape2007 Tape2007_kernel Tromp2005 Tromp2005_kernel acoustic_poroelastic canyon fluid_solid init_plane noise_uniform
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Wed Nov 16 04:46:58 PST 2011
Author: dkomati1
Date: 2011-11-16 04:46:57 -0800 (Wed, 16 Nov 2011)
New Revision: 19207
Modified:
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/Par_file
seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/Par_file_Tape2007_onerec
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/Par_file_Tromp2005
seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/Par_file_acoustic_poroelastic
seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon
seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/Par_file_fluid_solid
seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave
seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_for
seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_kernel
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_1
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_2
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_3
Log:
Paul Cristini converted all the EXAMPLES/Par_files to the new format with more input parameters
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = Test for Abel Balanche UMR 6538 Brest
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# 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 @@
# 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 @@
nreceiverlines = 2 # number of receiver lines
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 line (in the ocean)
nrec = 11 # number of receivers
@@ -67,8 +70,13 @@
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)
@@ -91,6 +99,15 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job, and file that contains interface data
+# title of job, and file that contains interface data r19201
title = Shallow water 30 Hz Source middle 4 receiver lines
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# parameters concerning partitioning
nproc = 4 # number of processes
partitioning_method = 3 # ascending order = 1, Scotch = 3
+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)
@@ -25,6 +26,7 @@
# time step parameters
nt = 2400 # total number of time steps
deltat = 5.d-7 # 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]
@@ -40,6 +42,7 @@
nreceiverlines = 1 # number of receiver lines
anglerec = 0.d0 # angle to rotate components at receivers
rec_normal_to_surface = .false. # replace anglerec with the normal to the surface (external mesh and curve file needed, anglerec above ignored)
+SU_FORMAT = .false. # output seismograms in Seismic Unix format (adjoint traces will be read in the same format)
# first receiver line
nrec = 11 # number of receivers
@@ -60,8 +63,13 @@
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)
@@ -87,6 +95,15 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job, and file that contains interface data
+# title of job, and file that contains interface data r19201
title = Shallow water 30 Hz Source middle 4 receiver lines
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # ascending order = 1, Scotch = 3
+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)
@@ -25,6 +26,7 @@
# time step parameters
nt = 2400 # total number of time steps
deltat = 5.d-7 # 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]
@@ -40,6 +42,7 @@
nreceiverlines = 1 # number of receiver lines
anglerec = 0.d0 # angle to rotate components at receivers
rec_normal_to_surface = .false. # replace anglerec with the normal to the surface (external mesh and curve file needed, anglerec above ignored)
+SU_FORMAT = .false. # output seismograms in Seismic Unix format (adjoint traces will be read in the same format)
# first receiver line
nrec = 11 # number of receivers
@@ -60,8 +63,13 @@
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)
@@ -87,6 +95,15 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/Par_file
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/Par_file 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/Par_file 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = Test Adjoint Pure Acoustic
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# 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 @@
# time step parameters
nt = 10000 # total number of time steps
deltat = 5.d-4 # 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 @@
nreceiverlines = 1 # number of receiver lines
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 line (repeat these 6 lines and adjust nreceiverlines accordingly)
nrec = 351 # number of receivers
@@ -59,8 +62,13 @@
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)
@@ -85,6 +93,15 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = Test for M2 UPPA
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# 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 @@
# time step parameters
nt = 1600 # 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 @@
nreceiverlines = 1 # number of receiver lines
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 line (repeat these 6 lines and adjust nreceiverlines accordingly)
nrec = 11 # number of receivers
@@ -59,8 +62,13 @@
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)
@@ -85,6 +93,15 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = Tape-Liu-Tromp (GJI 2007)
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# 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 = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
@@ -24,6 +25,7 @@
# time step parameters
nt = 4800 # total number of time steps
deltat = 6.0d-2 # 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 @@
nreceiverlines = 132 # number of receiver lines
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 line
nrec = 1 # number of receivers
@@ -845,8 +848,13 @@
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)
@@ -868,6 +876,15 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = Tape-Liu-Tromp (GJI 2007)
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# 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 = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
@@ -24,6 +25,7 @@
# time step parameters
nt = 3000 # total number of time steps
deltat = 6.0d-2 # 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 @@
nreceiverlines = 1 # number of receiver lines
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 line
nrec = 1 # number of receivers
@@ -60,8 +63,13 @@
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 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/Par_file_Tape2007_onerec
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/Par_file_Tape2007_onerec 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/Par_file_Tape2007_onerec 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = Tape-Liu-Tromp (GJI 2007)
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# 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 = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
@@ -24,6 +25,7 @@
# time step parameters
nt = 3000 # total number of time steps
deltat = 6.0d-2 # 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 @@
nreceiverlines = 1 # number of receiver lines
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 line
nrec = 1 # number of receivers
@@ -60,8 +63,13 @@
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 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = Tromp-Tape-Liu (GJI 2005)
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# 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 = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
@@ -24,6 +25,7 @@
# time step parameters
nt = 3000 # total number of time steps
deltat = 2.0d-2 # 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 @@
nreceiverlines = 1 # number of receiver lines
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 line
nrec = 1 # number of receivers
@@ -59,8 +62,13 @@
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 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = Tromp-Tape-Liu (GJI 2005)
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# 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 = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
@@ -24,6 +25,7 @@
# time step parameters
nt = 3000 # total number of time steps
deltat = 2.0d-4 # 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 @@
nreceiverlines = 1 # number of receiver lines
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 line
nrec = 1 # number of receivers
@@ -59,8 +62,13 @@
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 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/Par_file_Tromp2005
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/Par_file_Tromp2005 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/Par_file_Tromp2005 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = Tromp-Tape-Liu (GJI 2005)
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# 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 = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
@@ -24,6 +25,7 @@
# time step parameters
nt = 3000 # total number of time steps
deltat = 2.0d-2 # 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 @@
nreceiverlines = 1 # number of receiver lines
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 line
nrec = 1 # number of receivers
@@ -59,8 +62,13 @@
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 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/Par_file_acoustic_poroelastic
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/Par_file_acoustic_poroelastic 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/Par_file_acoustic_poroelastic 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = Test for 2 layers: acoustic/poroelastic
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# 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 @@
# time step parameters
nt = 5000 # total number of time steps
deltat = 3.d-4 # 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 @@
nreceiverlines = 2 # number of receiver lines
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 line (in the acoustic layer)
nrec = 1 # number of receivers
@@ -67,8 +70,13 @@
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)
@@ -91,6 +99,15 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = Canyon Mexico Paco P-SV
# forward or adjoint simulation
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/Par_file_fluid_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/Par_file_fluid_solid 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/Par_file_fluid_solid 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = External mesh with fluid and solid
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# 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 = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
@@ -24,6 +25,7 @@
# time step parameters
nt = 18000 # total number of time steps
deltat = 0.25d-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 @@
nreceiverlines = 1 # number of receiver lines
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 line
nrec = 12 # number of receivers
@@ -59,8 +62,13 @@
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)
@@ -87,6 +95,15 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = Slave Craton
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# 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 = 4 # number of control nodes per element (4 or 9)
initialfield = .true. # use a plane wave as source or not
@@ -24,6 +25,7 @@
# time step parameters
nt = 6000 # total number of time steps
deltat = 0.01 # 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 @@
nreceiverlines = 1 # number of receiver lines
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 line
nrec = 30
@@ -59,8 +62,13 @@
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 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_for
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_for 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_for 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = Slave Craton
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# 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 = 4 # number of control nodes per element (4 or 9)
initialfield = .true. # use a plane wave as source or not
@@ -24,6 +25,7 @@
# time step parameters
nt = 6000 # total number of time steps
deltat = 0.01 # 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 @@
nreceiverlines = 1 # number of receiver lines
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 line
nrec = 30
@@ -59,8 +62,13 @@
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 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_kernel
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_kernel 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_kernel 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = Slave Craton
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# 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 = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
@@ -24,6 +25,7 @@
# time step parameters
nt = 6000 # total number of time steps
deltat = 0.01 # 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 @@
nreceiverlines = 1 # number of receiver lines
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 line
nrec = 30
@@ -59,8 +62,13 @@
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 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_1
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_1 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_1 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = Noise_2D
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# 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 @@
# time step parameters
nt = 3000 # total number of time steps
deltat = 4.d-2 # 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 @@
nreceiverlines = 1 # number of receiver lines
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 line (repeat these 6 lines and adjust nreceiverlines accordingly)
nrec = 1 # number of receivers
@@ -59,8 +62,13 @@
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 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_2
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_2 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_2 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = Noise_2D
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# 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 @@
# time step parameters
nt = 3000 # total number of time steps
deltat = 4.d-2 # 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 @@
nreceiverlines = 1 # number of receiver lines
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 line (repeat these 6 lines and adjust nreceiverlines accordingly)
nrec = 1 # number of receivers
@@ -59,8 +62,13 @@
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 @@
# 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
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_3
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_3 2011-11-16 11:57:37 UTC (rev 19206)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_3 2011-11-16 12:46:57 UTC (rev 19207)
@@ -1,4 +1,4 @@
-# title of job
+# title of job r19201
title = Noise_2D
# forward or adjoint simulation
@@ -9,6 +9,7 @@
# 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 @@
# time step parameters
nt = 3000 # total number of time steps
deltat = 4.d-2 # 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 @@
nreceiverlines = 1 # number of receiver lines
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 line (repeat these 6 lines and adjust nreceiverlines accordingly)
nrec = 1 # number of receivers
@@ -59,8 +62,13 @@
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 @@
# 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