[cig-commits] r8482 - seismo/2D/SPECFEM2D/trunk/DATA

[walter at geodynamics.org]

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