[cig-commits] r20325 - in seismo/2D/SPECFEM2D/trunk/EXAMPLES: . Abel_Brest Gmsh_example_MPI Gmsh_example_serial INDUSTRIAL_FORMAT M2_UPPA Rayleigh_wave_no_crack Rayleigh_wave_with_crack Tape2007 Tape2007_kernel Tromp2005 Tromp2005_kernel acoustic_poroelastic attenuation canyon fluid_solid/fluid_solid_external_mesh fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom init_plane noise_layered/model_0 noise_layered/model_1 noise_layered/model_2 noise_uniform semi_infinite_homo
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Wed Jun 6 16:28:42 PDT 2012
Author: dkomati1
Date: 2012-06-06 16:28:41 -0700 (Wed, 06 Jun 2012)
New Revision: 20325
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/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
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/attenuation/Par_file_attenuation_2D
seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon
seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid
seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/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_layered/model_0/Par_file_fair
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_good
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_best
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_fair
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_good
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_fair
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_good
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
seismo/2D/SPECFEM2D/trunk/EXAMPLES/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py
seismo/2D/SPECFEM2D/trunk/EXAMPLES/semi_infinite_homo/Par_file_elastic_2D
Log:
used Paul Cristini's script to process all the Par_files to convert them to the new format
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 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = Test for Abel Balanche UMR 6538 Brest
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 4 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 4 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
nreceiversets = 2 # 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 (in the ocean)
nrec = 11 # number of receivers
@@ -61,27 +66,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 300 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 300 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 2 # nb of different models
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 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job, and file that contains interface data r19201
+# title of job
title = Shallow water 30 Hz Source middle 4 receiver sets
# forward or adjoint simulation
@@ -38,12 +38,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 4 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 4 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
nreceiversets = 1 # number of receiver sets
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 set
nrec = 11 # number of receivers
@@ -54,27 +59,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 300 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 300 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 10 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .false. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 4 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 2 # nb of different models
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 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job, and file that contains interface data r19201
+# title of job
title = Shallow water 30 Hz Source middle 4 receiver sets
# forward or adjoint simulation
@@ -38,12 +38,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 4 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 4 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
nreceiversets = 1 # number of receiver sets
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 set
nrec = 11 # number of receivers
@@ -54,27 +59,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 300 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 300 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 10 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .false. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 4 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 2 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/Par_file
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/Par_file 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/Par_file 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = Test Adjoint Pure Acoustic
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 6 # record 1=displ 2=veloc 3=accel 4=pressure 6=potential
+seismotype = 6 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential 6=potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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 = 351 # number of receivers
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 100 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 100 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .false. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .true. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 4 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = Test for M2 UPPA
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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 = 11 # number of receivers
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .true. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 100 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 100 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 4 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D 2012-06-06 23:28:41 UTC (rev 20325)
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
nrec = 5 # number of receivers
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 200 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 200 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_postscript = 1 # subsampling of background velocity model in PostScript 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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 1 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D 2012-06-06 23:28:41 UTC (rev 20325)
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
nrec = 5 # number of receivers
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 200 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 200 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_postscript = 1 # subsampling of background velocity model in PostScript 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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 3 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = Tape-Liu-Tromp (GJI 2007)
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
nreceiversets = 132 # 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
nrec = 1 # number of receivers
@@ -839,27 +844,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 400 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 400 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .false. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .true. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 1 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = Tape-Liu-Tromp (GJI 2007)
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
nrec = 1 # number of receivers
@@ -54,27 +59,40 @@
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 400 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 400 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .false. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .true. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 1 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/Par_file_Tape2007_onerec
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/Par_file_Tape2007_onerec 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/Par_file_Tape2007_onerec 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = Tape-Liu-Tromp (GJI 2007)
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
nrec = 1 # number of receivers
@@ -54,27 +59,40 @@
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 400 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 400 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .false. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 1 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = Tromp-Tape-Liu (GJI 2005)
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
nrec = 1 # number of receivers
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 400 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 400 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .true. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 1 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = Tromp-Tape-Liu (GJI 2005)
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
nrec = 1 # number of receivers
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 400 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 400 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 1 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/Par_file_Tromp2005
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/Par_file_Tromp2005 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/Par_file_Tromp2005 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = Tromp-Tape-Liu (GJI 2005)
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
nrec = 1 # number of receivers
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 400 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 400 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 1 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/Par_file_acoustic_poroelastic
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/Par_file_acoustic_poroelastic 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/Par_file_acoustic_poroelastic 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = Test for 2 layers: acoustic/poroelastic
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 4 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 4 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
nreceiversets = 2 # 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 (in the acoustic layer)
nrec = 1 # number of receivers
@@ -61,27 +66,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 300 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 300 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .false. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .false. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 2 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/attenuation/Par_file_attenuation_2D
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/attenuation/Par_file_attenuation_2D 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/attenuation/Par_file_attenuation_2D 2012-06-06 23:28:41 UTC (rev 20325)
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
nrec = 1 # number of receivers
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 100 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 100 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 1 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = Canyon Mexico Paco P-SV
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
nrec = 5 # number of receivers
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 250 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 250 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .false. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .false. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_postscript = 1 # subsampling of background velocity model in PostScript 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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 1 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = External mesh with fluid and solid
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
nrec = 12 # number of receivers
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 2000 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 2000 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 2 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid 2012-06-06 23:28:41 UTC (rev 20325)
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 2 # record 1=displ 2=veloc 3=accel 4=pressure 6=potential
+seismotype = 2 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential 6=potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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 = 110 # number of receivers
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 200 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 200 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results
+imagetype_JPEG = 6 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output JPEG color image of the results
-imagetype = 2 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_postscript = 1 # subsampling of background velocity model in PostScript 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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 2 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid 2012-06-06 23:28:41 UTC (rev 20325)
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 2 # record 1=displ 2=veloc 3=accel 4=pressure 6=potential
+seismotype = 2 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential 6=potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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 = 110 # number of receivers
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 200 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 200 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results
+imagetype_JPEG = 6 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output JPEG color image of the results
-imagetype = 2 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_postscript = 1 # subsampling of background velocity model in PostScript 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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 2 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = Slave Craton
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# eceiver line parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
nrec = 30
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 200 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 200 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 5. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 5. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 2 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_for
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_for 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_for 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = Slave Craton
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# eceiver line parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
nrec = 30
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 200 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 200 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 5. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 5. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 2 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_kernel
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_kernel 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_kernel 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = Slave Craton
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# eceiver line parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
nrec = 30
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 200 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 200 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 5. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 5. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 2 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_fair
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_fair 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_fair 2012-06-06 23:28:41 UTC (rev 20325)
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .true. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 250 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 250 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .false. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .false. # output Postscript snapshot of the results
-output_color_image = .false. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 1 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_good
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_good 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_good 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-s# title of job
+# title of job
4itle = Passive Imaging - 0 Interfaces
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .true. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 250 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 250 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .false. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .false. # output Postscript snapshot of the results
-output_color_image = .false. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 1 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_best
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_best 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_best 2012-06-06 23:28:41 UTC (rev 20325)
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .true. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 250 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 250 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .false. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .false. # output Postscript snapshot of the results
-output_color_image = .false. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 2 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_fair
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_fair 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_fair 2012-06-06 23:28:41 UTC (rev 20325)
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .true. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 250 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 250 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .false. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .false. # output Postscript snapshot of the results
-output_color_image = .false. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 2 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_good
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_good 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_good 2012-06-06 23:28:41 UTC (rev 20325)
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .true. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 250 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 250 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .false. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .false. # output Postscript snapshot of the results
-output_color_image = .false. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 2 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_fair
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_fair 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_fair 2012-06-06 23:28:41 UTC (rev 20325)
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .true. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 250 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 250 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .false. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .false. # output Postscript snapshot of the results
-output_color_image = .false. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 3 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_good
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_good 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_good 2012-06-06 23:28:41 UTC (rev 20325)
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .true. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 250 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 250 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .false. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .false. # output Postscript snapshot of the results
-output_color_image = .false. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 3 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_1
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_1 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_1 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = Noise_2D
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .false. # creates a STATION file in ./DATA
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 = 1 # number of receivers
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 250 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 250 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .false. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .false. # output Postscript snapshot of the results
-output_color_image = .false. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 1 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_2
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_2 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_2 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = Noise_2D
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .false. # creates a STATION file in ./DATA
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 = 1 # number of receivers
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 250 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 250 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .false. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .false. # output Postscript snapshot of the results
-output_color_image = .false. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 1 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_3
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_3 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_3 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,4 +1,4 @@
-# title of job r19201
+# title of job
title = Noise_2D
# forward or adjoint simulation
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .false. # creates a STATION file in ./DATA
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 = 1 # number of receivers
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 250 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 250 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .false. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .false. # output Postscript snapshot of the results
-output_color_image = .false. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 1 # nb of different models
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py 2012-06-06 23:28:41 UTC (rev 20325)
@@ -1,7 +1,7 @@
# -*- coding: utf-8 -*-
"""
Created on Tue Nov 15 09:00:00 2011
-Updated on Fri Jan 13 2012
+Updated on Fri Jun 5 2012
Processing of Par_file to update them to new format
@@ -12,8 +12,40 @@
"""
import os, string, sys
from os import listdir, walk
+from string import capitalize
from shutil import copy, move
from os.path import exists
+
+class InfoSEM:
+ def __init__(self,titre='Simulation',f0=10.,source_type=1):
+ self.title = titre
+ self.f0 = f0
+ self.source_type = source_type
+ def __str__(self):
+ return self.title+'\n'+'f0= '+str(self.f0)+' Hz'
+#
+def OuvreParFile(fic):
+ SEM=InfoSEM()
+ f = file(fic,'r')
+ lignes= f.readlines()
+ f.close()
+ # Lists of variables
+ if exists(fic):
+ # Numerical variables
+ VariableNum=['imagetype']
+ else:
+ print 'No Par_file found !'
+ return
+ #
+ for var in VariableNum:
+ print var
+ for ligne in lignes:
+ lsplit=string.split(ligne)
+ if lsplit!= []:
+ if lsplit[0]==var:
+ exec 'SEM.'+var+'='+string.replace(string.split(''.join(ligne))[2],'d','e')
+ break
+ return SEM
#------------------------------------------------------------------------------
def LoadLig(Fichier):
f = open(Fichier,'r')
@@ -25,7 +57,6 @@
for file in os.listdir(currdir):
path=os.path.join(currdir, file)
if not os.path.isdir(path):
- #print path
Fichiers.append(path)
else:
mylister(path)
@@ -33,28 +64,39 @@
def ProcessParfile_r19201(fic):
# Open the file and get all lines from Par_file
ligs= LoadLig(fic)
- # Test pour voir si le traitement a déjà été fait
+ # Test if already processed
for lig in ligs:
if 'ADD_PERIODIC_CONDITIONS' in lig:
print '----> '+fic+' already processed to r19201'
return
# New additions to the Par_file
- a1='PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)\n'
- a2='USER_T0 = 0.0d0 # use this t0 as earliest starting time rather than the automatically calculated one\n'
- a3='SU_FORMAT = .false. # output seismograms in Seismic Unix format (adjoint traces will be read in the same format)\n'
- a4='factor_subsample_image = 1 # factor to subsample color images output by the code (useful for very large models)\n'+ \
- 'POWER_DISPLAY_COLOR = 0.30d0 # non linear display to enhance small amplitudes in color images\n'+ \
- 'DRAW_WATER_CONSTANT_BLUE_IN_JPG = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water\n'
- a5='US_LETTER = .false. # US letter paper or European A4\n'+ \
- 'USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step\n'
- a6='\n# for horizontal periodic conditions: detect common points between left and right edges\n'+ \
+ a1='PERFORM_CUTHILL_MCKEE = .false. # perform ' + \
+ 'inverse Cuthill-McKee (1969) optimization/permutation for mesh ' + \
+ 'numbering (can be very costly and not very useful)\n'
+ a2='USER_T0 = 0.0d0 # use this t0' + \
+ ' as earliest starting time rather than the automatically calculated one\n'
+ a3='SU_FORMAT = .false. # output ' + \
+ 'seismograms in Seismic Unix format (adjoint traces will be read' + \
+ ' in the same format)\n'
+ a4='factor_subsample_image = 1 # factor to' + \
+ ' subsample color images output by the code (useful for very large models)\n'+ \
+ 'POWER_DISPLAY_COLOR = 0.30d0 # non linear ' + \
+ 'display to enhance small amplitudes in color images\n'+ \
+ 'DRAW_WATER_CONSTANT_BLUE_IN_JPG = .true. # display acoustic' + \
+ ' layers as constant blue in JPEG images, because they likely correspond to water\n'
+ a5='US_LETTER = .false. # US letter ' + \
+ 'paper or European A4\n'+ \
+ 'USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot ' + \
+ 'number in the file name of JPEG color snapshots instead of the time step\n'
+ a6='\n# for horizontal periodic conditions: detect common points' + \
+ ' between left and right edges\n'+ \
'ADD_PERIODIC_CONDITIONS = .false.\n\n'+ \
'# horizontal periodicity distance for periodic conditions\n'+ \
'PERIODIC_horiz_dist = 0.3597d0\n\n'+ \
'# grid point detection tolerance for periodic conditions\n'+ \
'PERIODIC_DETECT_TOL = 3.3334d-6\n'
#--------------------------------------------------------------------------
- # Ajout des parametres supplementaires
+ # Add new parameters
#
for ilg, lig in enumerate(ligs):
if lig.startswith('partitioning'):
@@ -87,13 +129,13 @@
def ProcessParfile_r19340(fic):
# Open the file and get all lines from Par_file
ligs= LoadLig(fic)
- # Teste si le traitement a déjà été fait
+ # Test if already processed
for lig in ligs:
if 'nreceiversets' in lig:
print '----> '+fic+' already processed to r19340'
return
#
- # Ajout des parametres supplementaires
+ # Add new parameters
#
for ilg, lig in enumerate(ligs):
if lig.startswith('nreceiverlines'):
@@ -111,13 +153,13 @@
def ProcessParfile_r19346(fic):
# Open the file and get all lines from Par_file
ligs= LoadLig(fic)
- # Teste si le traitement a déjà été fait
+ # Test if already processed
for lig in ligs:
if 'ATTENUATION_PORO_FLUID_PART' in lig:
print '----> '+fic+' already processed to r19346'
return
#--------------------------------------------------------------------------
- # Ajout des parametres supplementaires
+ # Add new parameters
#
for ilg, lig in enumerate(ligs):
if lig.startswith('TURN_ATTENUATION_ON'):
@@ -136,13 +178,13 @@
print 'xxxxx------> '+fic+' processed to r19346'
return
#------------------------------------------------------------------------------
-def ProcessParfile_r19xxx(fic):
+def ProcessParfile_r19521(fic):
# Open the file and get all lines from Par_file
ligs= LoadLig(fic)
- # Teste si le traitement a déjà été fait
+ # Test if already processed
for lig in ligs:
if 'time_stepping_scheme' in lig:
- print '----> '+fic+' already processed to r19xxx'
+ print '----> '+fic+' already processed to r19521'
return
#
a1='time_stepping_scheme = 1 # 1 = Newmark (2nd order), \
@@ -150,23 +192,201 @@
3 = classical 4th-order 4-stage Runge-Kutta\n'
#--------------------------------------------------------------------------
- # Ajout des parametres supplementaires
+ # Add new parameters
#
for ilg, lig in enumerate(ligs):
if lig.startswith('USER_T0'):
ligs.insert(ilg+1,a1)
#
- move(fic,fic+'.before_update_to_r19xxx')
+ move(fic,fic+'.before_update_to_r19521')
#
fm = open(fic,'w')
fm.writelines(ligs)
fm.close()
#
- print 'xxxxx------> '+fic+' processed to r19xxx'
- return
+ print 'xxxxx------> '+fic+' processed to r19521'
+ return
#------------------------------------------------------------------------------
+def ProcessParfile_r19804(fic):
+ # Open the file and get all lines from Par_file
+ ligs= LoadLig(fic)
+
+ for ilg, lig in enumerate(ligs):
+ if lig.startswith('PERFORM_CUTHILL_MCKEE'):
+ ligs[ilg]=ligs[ilg].replace('.true.','.false.')
+
+ # Test if already processed
+ for lig in ligs:
+ if 'ADD_SPRING_TO_STACEY' in lig:
+ print '----> '+fic+' already processed to r19804'
+ return
+ #
+ a1='ADD_SPRING_TO_STACEY = .true.\n'
+
+ #--------------------------------------------------------------------------
+ # Add new parameters
+ #
+ for ilg, lig in enumerate(ligs):
+ if lig.startswith('absorbing_conditions'):
+ ligs.insert(ilg+1,a1)
+ #
+ move(fic,fic+'.before_update_to_r19804')
+ #
+ fm = open(fic,'w')
+ fm.writelines(ligs)
+ fm.close()
+ #
+ print 'xxxxx------> '+fic+' processed to r19804'
+ return
+#------------------------------------------------------------------------------
+def ProcessParfile_r20307(fic):
+ # Open the file and get all lines from Par_file
+ ligs= LoadLig(fic)
+
+ # Test if already processed
+ for lig in ligs:
+ if 'NSTEP_BETWEEN_OUTPUT_SEISMOS' in lig:
+ print '----> '+fic+' already processed to r20307'
+ return
+
+ # Change the jpeg imagetype accordingly to the original values
+ if Data.imagetype==1: imagetype_new = 3
+ if Data.imagetype==2: imagetype_new = 6
+ if Data.imagetype==3: imagetype_new = 9
+ if Data.imagetype==4: imagetype_new = 10
+
+ # Change first line
+ ligs[0]='# title of job\n'
+
+ # New parameters
+ a0='NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many' + \
+ ' time steps we save the seismograms (costly, do not use a very small' + \
+ ' value; if you use a very large value that is larger than the total ' + \
+ 'number of time steps of the run, the seismograms will automatically ' + \
+ 'be saved once at the end of the run anyway)\n' + \
+ 'save_ASCII_seismograms = .true. # save seismograms ' + \
+ 'in ASCII format or not\n' + \
+ 'save_binary_seismograms_single = .true. # save seismograms ' + \
+ 'in single precision binary format or not (can be used jointly with' + \
+ ' ASCII above to save both)\n' + \
+ 'save_binary_seismograms_double = .false. # save seismograms' + \
+ ' in double precision binary format or not (can be used jointly with' + \
+ ' both flags above to save all)\n'
+ #
+ a1='subsamp_seismos = 1 # subsampling of ' + \
+ 'the seismograms to create smaller files (but less accurately ' + \
+ 'sampled in time)\n'
+ #
+ a2='NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how ' + \
+ 'many time steps we draw JPEG or PostScript pictures of the ' + \
+ 'simulation (costly, do not use a very small value)\n'
+ #
+ a3='imagetype_JPEG = '+str(imagetype_new)+' ' + \
+ ' # display ' + \
+ '1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz ' + \
+ '6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure\n'
+ #
+ a4='DRAW_SOURCES_AND_RECEIVERS = .true. # display ' + \
+ 'sources as orange crosses and receivers as green squares in' + \
+ ' JPEG images or not\n'
+ #
+ a5='imagetype_postscript = 1 ' + \
+ ' # display 1=displ vector 2=veloc vector 3=accel vector;' + \
+ ' small arrows are displayed for the vectors\n'
+ #
+ a6='NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how ' + \
+ 'many time steps we dump results of the simulation as ASCII or ' + \
+ 'binary files (costly, do not use a very small value)\n' + \
+ 'output_wavefield_dumps = .false. # output wave' + \
+ ' field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time' + \
+ ' steps (creates very big files)\n' + \
+ 'imagetype_wavefield_dumps = 1 # display 1=displ' + \
+ ' vector 2=veloc vector 3=accel vector 4=pressure\n' + \
+ 'use_binary_for_wavefield_dumps = .false. # use ASCII' + \
+ ' or single-precision binary format for the wave field dumps\n####\n'
+ a7='output_grid_ASCII = .false. # dump the' + \
+ ' grid in an ASCII text file consisting of a set of X,Y,Z points or not\n'
+ #
+ for ilg, lig in enumerate(ligs):
+ if lig.startswith('seismotype'):
+ ligs[ilg] = ligs[ilg].replace('# record 1=displ 2=veloc' + \
+ ' 3=accel 4=pressure','# record 1=displ 2=veloc 3=accel ' + \
+ '4=pressure 5=curl of displ 6=the fluid potential')
+ ligs.insert(ilg+1,a0)
+ ligs.insert(ilg+2,a1)
+
+ if lig.startswith('SU_FORMAT'):
+ ligs[ilg] = ligs[ilg].replace('# output seismograms in Seismic' + \
+ ' Unix format','# output single precision binary seismograms' + \
+ ' in Seismic Unix format')
+ ligs.insert(ilg-5,ligs.pop(ilg))
+
+ if lig.startswith('NTSTEP_BETWEEN_OUTPUT_INFO'):
+ ligs[ilg] = ligs[ilg].replace('NTSTEP','NSTEP')
+ ligs[ilg] = ligs[ilg].replace('# display frequency in time ' + \
+ 'steps','# every how many time steps we display information' + \
+ ' about the simulation (costly, do not use a very small value)')
+ ligs.insert(ilg+1,a2)
+ ligs.insert(ilg+5,ligs.pop(ilg+2)) #move output_postscript-snapshot
+ ligs.insert(ilg+4,ligs.pop(ilg+2)) #move output_color_image
+ ligs[ilg+4] = ligs[ilg+4].replace('# output color image of' + \
+ ' the results','# output JPEG color image of the results' +\
+ ' every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not')
+ ligs.insert(ilg+4,'\n#### for JPEG color images ####\n')
+ ligs.pop(ilg+2)
+ ligs.insert(ilg+5,a3)
+ ligs.insert(ilg+6,ligs.pop(ilg+13)) # move factor_subsample_image
+ ligs.insert(ilg+7,ligs.pop(ilg+14)) # move POWER_DISPLAY_COLOR
+ ligs.insert(ilg+8,ligs.pop(ilg+15)) # move DRAW_WATER...
+ ligs[ilg+8] = ligs[ilg+8].replace('DRAW_WATER_CONSTANT' + \
+ '_BLUE_IN_JPG =','DRAW_WATER_IN_BLUE =')
+ ligs.insert(ilg+8,a4)
+ ligs.insert(ilg+10,ligs.pop(ilg+19)) # move USE_SNAPSHOT-NUMBER_IN
+ ligs.insert(ilg+11,'\n#### for PostScript snapshots ####\n')
+ ligs.insert(ilg+13,a5)
+
+ if lig.startswith('US_LETTER'):
+ ligs.insert(ilg+1,'\n#### for wavefield dumps ####\n')
+ ligs.insert(ilg+2,a6)
+
+ if lig.startswith('gnuplot ='):
+ ligs.pop(ilg)
+
+ for ilg, lig in enumerate(ligs):
+ if lig.startswith('cutsnaps ='):
+ ligs[ilg] = ligs[ilg].replace('# minimum amplitude in % for' + \
+ ' snapshots','# minimum amplitude kept in % for the JPEG and' + \
+ ' PostScript snapshots; amplitudes below that are muted')
+
+ if lig.startswith('DRAW_WATER_IN_BLUE ='):
+ ligs[ilg] = ligs[ilg].replace('# display acoustic layers as' + \
+ ' constant blue in JPEG images, because they likely correspond' + \
+ ' to water','# display acoustic layers as constant blue in JPEG'+ \
+ ' images, because they likely correspond to water in the case' + \
+ ' of ocean acoustics or in the case of offshore oil industry' + \
+ ' experiments (if off, display them as greyscale, as for ' + \
+ 'elastic or poroelastic elements, for instance for acoustic-' + \
+ 'only oil industry models of solid media)')
+
+ if lig.startswith('output_grid '):
+ ligs[ilg] = ligs[ilg].replace('output_grid ','output_grid_Gnuplot')
+ ligs.insert(ilg+1,a7)
+
+ if lig.startswith('output_wavefield_snapshot ='):
+ ligs.pop(ilg)
+
+ #
+ move(fic,fic+'.before_update_to_r20307')
+ #
+ fm = open(fic,'w')
+ fm.writelines(ligs)
+ fm.close()
+ #
+ print 'xxxxx------> '+fic+' processed to r20307'
+ return
+#------------------------------------------------------------------------------
if __name__=='__main__':
- ## Liste de tous les fichiers à partir du répertoire courant
+ ## List of all files of current directory
Fichiers=[]
mylister('.')
#
@@ -177,13 +397,16 @@
if not( ('.svn' in repert) or ('unused' in repert) or \
'.before_update_to_' in ficname):
if ficname.startswith('Par_file'):
- print 'Analysis of file : '+fic
if not (ficname.endswith('~')):
+ print 'Analysis of file : '+fic
Ct_Par_file+=1
+ Data = OuvreParFile(fic)
ProcessParfile_r19201(fic)
ProcessParfile_r19340(fic)
ProcessParfile_r19346(fic)
- ProcessParfile_r19xxx(fic)
+ ProcessParfile_r19521(fic)
+ ProcessParfile_r19804(fic)
+ ProcessParfile_r20307(fic)
print '~'*80
#
print 'Number of Par_file analysed : ', Ct_Par_file
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/semi_infinite_homo/Par_file_elastic_2D
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/semi_infinite_homo/Par_file_elastic_2D 2012-06-06 23:23:17 UTC (rev 20324)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/semi_infinite_homo/Par_file_elastic_2D 2012-06-06 23:28:41 UTC (rev 20325)
@@ -37,12 +37,17 @@
f0_attenuation = 5.196152422706633 # (Hz) relevant only if source is a Dirac or a Heaviside, else it is f0
# receiver set parameters for seismograms
-seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure 5=curl of displ 6=the fluid potential
+NSTEP_BETWEEN_OUTPUT_SEISMOS = 5000000 # every how many time steps we save the seismograms (costly, do not use a very small value; if you use a very large value that is larger than the total number of time steps of the run, the seismograms will automatically be saved once at the end of the run anyway)
+save_ASCII_seismograms = .true. # save seismograms in ASCII format or not
+save_binary_seismograms_single = .true. # save seismograms in single precision binary format or not (can be used jointly with ASCII above to save both)
+save_binary_seismograms_double = .false. # save seismograms in double precision binary format or not (can be used jointly with both flags above to save all)
+SU_FORMAT = .false. # output single precision binary seismograms in Seismic Unix format (adjoint traces will be read in the same format)
+subsamp_seismos = 1 # subsampling of the seismograms to create smaller files (but less accurately sampled in time)
generate_STATIONS = .true. # creates a STATION file in ./DATA
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
nrec = 1 # number of receivers
@@ -53,27 +58,40 @@
enreg_surf_same_vertical = .false. # receivers inside the medium or at the surface
# display parameters
-NTSTEP_BETWEEN_OUTPUT_INFO = 100 # display frequency in time steps
+NSTEP_BETWEEN_OUTPUT_INFO = 100 # every how many time steps we display information about the simulation (costly, do not use a very small value)
+NSTEP_BETWEEN_OUTPUT_IMAGES = 100 # every how many time steps we draw JPEG or PostScript pictures of the simulation (costly, do not use a very small value)
+cutsnaps = 1. # minimum amplitude kept in % for the JPEG and PostScript snapshots; amplitudes below that are muted
+
+#### for JPEG color images ####
+output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps or not
+imagetype_JPEG = 3 # display 1=displ_Ux 2=displ_Uz 3=displ_norm 4=veloc_Vx 5=veloc_Vz 6=veloc_norm 7=accel_Ax 8=accel_Az 9=accel_norm 10=pressure
+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_SOURCES_AND_RECEIVERS = .true. # display sources as orange crosses and receivers as green squares in JPEG images or not
+DRAW_WATER_IN_BLUE = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water in the case of ocean acoustics or in the case of offshore oil industry experiments (if off, display them as greyscale, as for elastic or poroelastic elements, for instance for acoustic-only oil industry models of solid media)
+USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step
+
+#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results
-output_color_image = .true. # output color image of the results
-imagetype = 1 # display 1=displ 2=veloc 3=accel 4=pressure
-cutsnaps = 1. # minimum amplitude in % for snapshots
+imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector; small arrows are displayed for the vectors
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
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_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
+
+#### for wavefield dumps ####
+NSTEP_BETWEEN_OUTPUT_WAVE_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII or binary files (costly, do not use a very small value)
+output_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+imagetype_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
+use_binary_for_wavefield_dumps = .false. # use ASCII or single-precision binary format for the wave field dumps
+####
+output_grid_Gnuplot = .false. # save the grid in a text file or not
+output_grid_ASCII = .false. # dump the grid in an ASCII text file consisting of a set of X,Y,Z points or not
output_energy = .false. # compute and output acoustic and elastic energy (slows down the code significantly)
-output_wavefield_snapshot = .false. # output Ux,Uy,Uz text file for each output time (big files)
# velocity and density models
nbmodels = 1 # nb of different models
More information about the CIG-COMMITS
mailing list