[cig-commits] r20313 - in seismo/2D/SPECFEM2D/trunk: DATA src/meshfem2D
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Mon Jun 4 17:45:18 PDT 2012
Author: dkomati1
Date: 2012-06-04 17:45:17 -0700 (Mon, 04 Jun 2012)
New Revision: 20313
Modified:
seismo/2D/SPECFEM2D/trunk/DATA/Par_file
seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90
Log:
moved cutsnaps a few lines above
Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file 2012-06-05 00:39:19 UTC (rev 20312)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file 2012-06-05 00:45:17 UTC (rev 20313)
@@ -26,10 +26,10 @@
nt = 1600 # total number of time steps
deltat = 1.d-3 # duration of a time step
USER_T0 = 0.0d0 # use this t0 as earliest starting time rather than the automatically calculated one
-time_stepping_scheme = 1 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical 4th-order 4-stage Runge-Kutta
+time_stepping_scheme = 1 # 1 = Newmark (2nd order), 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), 3 = classical RK4 4th-order 4-stage Runge-Kutta
# source parameters
-NSOURCES = 1 # number of sources [source info read in CMTSOLUTION file]
+NSOURCES = 1 # number of sources (source info read from DATA/CMTSOLUTION file)
force_normal_to_surface = .false. # angleforce normal to surface (external mesh and curve file needed)
# constants for attenuation
@@ -60,6 +60,7 @@
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 in % for JPEG and PostScript snapshots
#### for JPEG color images ####
output_color_image = .true. # output JPEG color image of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps
imagetype_JPEG = 2 # 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
@@ -71,7 +72,6 @@
#### for PostScript snapshots ####
output_postscript_snapshot = .true. # output Postscript snapshot of the results every NSTEP_BETWEEN_OUTPUT_IMAGES time steps
imagetype_postscript = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
-cutsnaps = 1. # minimum amplitude in % for JPEG and PostScript snapshots
meshvect = .true. # display mesh on PostScript plots or not
modelvect = .false. # display velocity model on PostScript plots
boundvect = .true. # display boundary conditions on PostScript plots
@@ -79,10 +79,10 @@
pointsdisp = 6 # points for interpolation of PostScript snapshots (set to 1 for lower-left corner only)
subsamp_postscript = 1 # subsampling of background velocity model in PostScript snapshots
sizemax_arrows = 1.d0 # maximum size of arrows on PostScript plots in centimeters
-US_LETTER = .false. # US letter or European A4 paper for PostScript plots
+US_LETTER = .false. # use US letter or European A4 paper for PostScript plots
####
NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS = 100 # every how many time steps we dump results of the simulation as ASCII text files (costly, do not use a very small value)
-output_TEXT_wavefield_dumps = .false. # output wave field to text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
+output_TEXT_wavefield_dumps = .false. # output wave field to a text file every NSTEP_BETWEEN_OUTPUT_TEXT_DUMPS time steps (creates very big files)
imagetype_TEXT_wavefield_dumps = 1 # display 1=displ vector 2=veloc vector 3=accel vector 4=pressure
####
output_grid_gnuplot = .false. # generate a GNUPLOT file containing the grid, and a script to plot it
Modified: seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90 2012-06-05 00:39:19 UTC (rev 20312)
+++ seismo/2D/SPECFEM2D/trunk/src/meshfem2D/read_parameter_file.F90 2012-06-05 00:45:17 UTC (rev 20313)
@@ -333,6 +333,9 @@
call read_value_integer_p(NSTEP_BETWEEN_OUTPUT_IMAGES, 'solver.NSTEP_BETWEEN_OUTPUT_IMAGES')
if(err_occurred() /= 0) stop 'error reading parameter 33c in Par_file'
+ call read_value_double_precision_p(cutsnaps, 'solver.cutsnaps')
+ if(err_occurred() /= 0) stop 'error reading parameter 37 in Par_file'
+
call read_value_logical_p(output_color_image, 'solver.output_color_image')
if(err_occurred() /= 0) stop 'error reading parameter 35 in Par_file'
@@ -360,9 +363,6 @@
call read_value_integer_p(imagetype_postscript, 'solver.imagetype_postscript')
if(err_occurred() /= 0) stop 'error reading parameter 36 in Par_file'
- call read_value_double_precision_p(cutsnaps, 'solver.cutsnaps')
- if(err_occurred() /= 0) stop 'error reading parameter 37 in Par_file'
-
call read_value_logical_p(meshvect, 'solver.meshvect')
if(err_occurred() /= 0) stop 'error reading parameter 38 in Par_file'
More information about the CIG-COMMITS
mailing list