[cig-commits] r8485 - seismo/2D/SPECFEM2D/trunk/DATA
walter at geodynamics.org
walter at geodynamics.org
Fri Dec 7 15:49:59 PST 2007
Author: walter
Date: 2007-12-07 15:49:58 -0800 (Fri, 07 Dec 2007)
New Revision: 8485
Modified:
seismo/2D/SPECFEM2D/trunk/DATA/Par_file
seismo/2D/SPECFEM2D/trunk/DATA/Par_file_benchmark_Hertzberg
Log:
changed the default Par_file to have elastic simulations instead of acoustic.
also updated the Par_file for the UPPA benchmark performed by Jacques Hertzberg.
Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file 2006-03-20 16:57:27 UTC (rev 8484)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file 2007-12-07 23:49:58 UTC (rev 8485)
@@ -1,18 +1,18 @@
# title of job, and file that contains interface data
-title = Test sinusoide acoustique
-interfacesfile = interface_sinus.dat
+title = Test pour cours M2 UPPA
+interfacesfile = interfaces_cours_M2_UPPA_curved.dat
# geometry of the model (origin lower-left corner = 0,0) and mesh description
-xmin = 0.d0 # abscissa of left side of the model
-xmax = 4000.d0 # abscissa of right side of the model
-nx = 80 # number of elements along X
-ngnod = 9 # noeuds de controle pour blocs (4 ou 9)
-initialfield = .false. # use a plane wave as source or not
-readmodel = .false. # read external earth model or not
-ELASTIC = .false. # elastic or acoustic simulation
-TURN_ANISOTROPY_ON = .false. # turn anisotropy on or off
-TURN_ATTENUATION_ON = .false. # turn attenuation on or off
+xmin = 0.d0 # abscissa of left side of the model
+xmax = 4000.d0 # abscissa of right side of the model
+nx = 80 # number of elements along X
+ngnod = 9 # noeuds de controle pour blocs (4 ou 9)
+initialfield = .false. # use a plane wave as source or not
+readmodel = .false. # read external earth model or not
+ELASTIC = .true. # elastic or acoustic simulation
+TURN_ANISOTROPY_ON = .false. # turn anisotropy on or off
+TURN_ATTENUATION_ON = .false. # turn attenuation on or off
# absorbing boundaries parameters
absorbhaut = .false. # absorbing boundary active or not
@@ -21,40 +21,40 @@
absorbdroite = .true.
# time step parameters
-nt = 1400 # nb total de pas de temps
+nt = 1600 # nb total de pas de temps
deltat = 1.d-3 # valeur du pas de temps
# source parameters
-source_surf = .false. # source dans le volume ou a la surface
+source_surf = .true. # source dans le volume ou a la surface
xs = 2000. # source location x in meters
-zs = 1200. # source location z in meters
+zs = 1500. # source location z in meters
source_type = 1 # force = 1 or moment tensor = 2
-time_function_type = 2 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
+time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
f0 = 10.0 # dominant source frequency (Hz) if not Dirac
angleforce = 0. # angle of the source (for a force only)
Mxx = 1. # Mxx component (for a moment tensor source only)
Mzz = 1. # Mzz component (for a moment tensor source only)
Mxz = 0. # Mxz component (for a moment tensor source only)
-factor = 1.d4 # amplification factor
+factor = 1.d10 # amplification factor
# receiver line parameters
-sismostype = 4 # record 1=displ 2=veloc 3=accel 4=pressure
+sismostype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
nreceiverlines = 1 # number of receiver lines
anglerec = 0.d0 # angle to rotate components at receivers
# first receiver line
-nrec = 101 # number of receivers
+nrec = 11 # number of receivers
xdeb = 300. # first receiver x in meters
-zdeb = 1700. # first receiver z in meters
+zdeb = 2200. # first receiver z in meters
xfin = 3700. # last receiver x in meters (ignored if onlyone receiver)
-zfin = 1700. # last receiver z in meters (ignored if onlyone receiver)
-enreg_surf = .false. # enregistrement volume ou surface
+zfin = 2200. # last receiver z in meters (ignored if onlyone receiver)
+enreg_surf = .true. # enregistrement volume ou surface
# display parameters
itaff = 100 # display frequency in time steps
output_postscript_snapshot = .true. # output Postscript image of the results
output_PNM_image = .true. # output PNM image of the results
-vecttype = 2 # display 1=displ 2=veloc 3=accel
+vecttype = 1 # display 1=displ 2=veloc 3=accel
cutvect = 1. # amplitude min en % pour vector plots
meshvect = .true. # display mesh on vector plots or not
modelvect = .false. # display velocity model on vector plots
@@ -67,11 +67,11 @@
outputgrid = .false. # save the grid in a text file or not
# velocity and density model (nx,nz)
-nbmodels = 2 # nb de modeles differents (0,rho,vp,vs,0,0)
-1 0 2700.d0 3000.d0 0.d0 0 0
-2 0 1800.d0 2000.d0 0.d0 0 0
-#3 0 2200.d0 2500.d0 0.d0 0 0
-nbzone = 2 # nb of zones and model number for each
-1 80 1 40 1
-1 80 31 40 2
-#1 80 41 60 3
+nbmodels = 3 # nb de modeles differents (0,rho,vp,vs,0,0)
+1 0 2700.d0 3000.d0 1732.051d0 0 0
+2 0 2500.d0 2700.d0 1558.845d0 0 0
+3 0 2200.d0 2500.d0 1443.375d0 0 0
+nbzone = 3 # nb of zones and model number for each
+1 80 1 20 1
+1 80 21 40 2
+1 80 41 60 3
Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file_benchmark_Hertzberg
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file_benchmark_Hertzberg 2006-03-20 16:57:27 UTC (rev 8484)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file_benchmark_Hertzberg 2007-12-07 23:49:58 UTC (rev 8485)
@@ -1,92 +1,77 @@
-# ----------------------------------------------------------------
-#
-# This is the parameter file
-# Put variable names first and actual value after 34th column
-#
-# ----------------------------------------------------------------
-#
+
# title of job, and file that contains interface data
-#
-title = Interface Cercle 40 60 10Hz
-interfacesfile = InterfaceCercleMod40_60.dat
-#
+title = Test pour cours M2 UPPA
+interfacesfile = interfaces_cours_M2_UPPA_curved.dat
+
# geometry of the model (origin lower-left corner = 0,0) and mesh description
-#
-xmin = 0.d0 ! abscissa of left side of the model
-xmax = 4000.d0 ! abscissa of right side of the model
-nx = 300 ! number of elements along X
-ngnod = 9 ! noeuds de controle pour blocs (4 ou 9)
-initialfield = .false. ! use a plane wave as source or not
-readmodel = .false. ! read external earth model or not
-ELASTIC = .false. ! elastic or acoustic simulation
-TURN_ANISOTROPY_ON = .false. ! turn anisotropy on or off
-TURN_ATTENUATION_ON = .false. ! turn attenuation on or off
-#
+xmin = 0.d0 # abscissa of left side of the model
+xmax = 4000.d0 # abscissa of right side of the model
+nx = 80 # number of elements along X
+ngnod = 9 # noeuds de controle pour blocs (4 ou 9)
+initialfield = .false. # use a plane wave as source or not
+readmodel = .false. # read external earth model or not
+ELASTIC = .true. # elastic or acoustic simulation
+TURN_ANISOTROPY_ON = .false. # turn anisotropy on or off
+TURN_ATTENUATION_ON = .false. # turn attenuation on or off
+
# absorbing boundaries parameters
-#
-absorbhaut = .true. ! absorbing boundary active or not
+absorbhaut = .false. # absorbing boundary active or not
absorbbas = .true.
absorbgauche = .true.
absorbdroite = .true.
-#
+
# time step parameters
-#
-nt = 15000 ! nb total de pas de temps
-dt = 1.d-4 ! valeur du pas de temps
-#
+nt = 5000 # nb total de pas de temps
+deltat = 1.d-3 # valeur du pas de temps
+
# source parameters
-#
-source_surf = .false. ! source dans le volume ou a la surface
-xs = 2000. ! source location x in meters
-zs = 2850. ! source location z in meters
-source_type = 1 ! force = 1 or explosion = 2
-time_function_type = 1 ! Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4
-f0 = 10 ! dominant source frequency (Hz) if not Dirac
-angleforce = 0. ! angle of the source (for a force only)
-Mxx = 1. ! Mxx component (for a moment tensor source only)
-Mzz = 1. ! Mzz component (for a moment tensor source only)
-Mxz = 0. ! Mxz component (for a moment tensor source only)
-factor = 1.d4 ! amplification factor
-#
+source_surf = .true. # source dans le volume ou a la surface
+xs = 2000. # source location x in meters
+zs = 1500. # source location z in meters
+source_type = 1 # force = 1 or moment tensor = 2
+time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
+f0 = 10.0 # dominant source frequency (Hz) if not Dirac
+angleforce = 0. # angle of the source (for a force only)
+Mxx = 1. # Mxx component (for a moment tensor source only)
+Mzz = 1. # Mzz component (for a moment tensor source only)
+Mxz = 0. # Mxz component (for a moment tensor source only)
+factor = 1.d10 # amplification factor
+
# receiver line parameters
-#
-sismostype = 4 ! record 1=displ 2=veloc 3=accel
-nreceiverlines = 1
-anglerec = 0.d0 ! angle to rotate components at receivers
+sismostype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+nreceiverlines = 1 # number of receiver lines
+anglerec = 0.d0 # angle to rotate components at receivers
# first receiver line
-nrec = 101 ! number of receivers
-xdeb = 300. ! first receiver x in meters
-zdeb = 2750. ! first receiver z in meters
-xfin = 3700. ! last receiver x in meters
-zfin = 2750. ! last receiver z in meters
-enreg_surf = .false. ! enregistrement volume ou surface
+nrec = 11 # number of receivers
+xdeb = 300. # first receiver x in meters
+zdeb = 2200. # first receiver z in meters
+xfin = 3700. # last receiver x in meters (ignored if onlyone receiver)
+zfin = 2200. # last receiver z in meters (ignored if onlyone receiver)
+enreg_surf = .true. # enregistrement volume ou surface
-#
# display parameters
-#
-itaff = 400 ! display frequency in time steps
-output_postscript_snapshot = .false. # output Postscript image of the results
-output_PNM_image = .false. # output PNM image of the results
-vecttype = 2 ! display 1=displ 2=veloc 3=accel
-cutvect = 1. ! amplitude min en % pour vector plots
-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
-subsamp = 2 ! subsampling of color snapshots
+itaff = 10000 # display frequency in time steps
+output_postscript_snapshot = .true. # output Postscript image of the results
+output_PNM_image = .true. # output PNM image of the results
+vecttype = 1 # display 1=displ 2=veloc 3=accel
+cutvect = 1. # amplitude min en % pour vector plots
+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 = 1 # subsampling of color snapshots
sizemax_arrows = 1.d0 # maximum size of arrows on vector plots in cm
-gnuplot = .false. ! generate a GNUPLOT file for the grid
-outputgrid = .false. ! save the grid in a text file or not
-#
+gnuplot = .false. # generate a GNUPLOT file for the grid
+outputgrid = .false. # save the grid in a text file or not
+
# velocity and density model (nx,nz)
-#
-nbmodels = 2 ! nb de modeles differents (0,rho,vp,vs,0,0)
-1 0 700.d0 3000.d0 0.d0 0 0
-2 0 1800.d0 2000.d0 0.d0 0 0
-#3 0 2200.d0 2500.d0 0.d0 0 0
-nbzone = 2 ! nb of zones and model number for each
-1 300 1 150 1
-1 300 151 323 2
-#1 80 41 60 3
+nbmodels = 3 # nb de modeles differents (0,rho,vp,vs,0,0)
+1 0 2700.d0 3000.d0 1732.051d0 0 0
+2 0 2500.d0 2700.d0 1558.845d0 0 0
+3 0 2200.d0 2500.d0 1443.375d0 0 0
+nbzone = 3 # nb of zones and model number for each
+1 80 1 20 1
+1 80 21 40 2
+1 80 41 60 3
More information about the cig-commits
mailing list