[cig-commits] r8482 - seismo/2D/SPECFEM2D/trunk/DATA
walter at geodynamics.org
walter at geodynamics.org
Fri Dec 7 15:49:44 PST 2007
Author: walter
Date: 2007-12-07 15:49:43 -0800 (Fri, 07 Dec 2007)
New Revision: 8482
Added:
seismo/2D/SPECFEM2D/trunk/DATA/Par_file_benchmark_Hertzberg
Log:
added Par_file for Jacques Hertzberg's benchmark at University of Pau, France
Added: seismo/2D/SPECFEM2D/trunk/DATA/Par_file_benchmark_Hertzberg
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file_benchmark_Hertzberg 2006-01-18 15:25:13 UTC (rev 8481)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file_benchmark_Hertzberg 2007-12-07 23:49:43 UTC (rev 8482)
@@ -0,0 +1,92 @@
+# ----------------------------------------------------------------
+#
+# 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
+#
+# 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
+#
+# absorbing boundaries parameters
+#
+absorbhaut = .true. ! 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
+#
+# 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
+#
+# receiver line parameters
+#
+sismostype = 4 ! record 1=displ 2=veloc 3=accel
+nreceiverlines = 1
+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
+
+#
+# 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
+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
+#
+# 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
More information about the cig-commits
mailing list