[cig-commits] r8540 - seismo/2D/SPECFEM2D/trunk/DATA
walter at geodynamics.org
walter at geodynamics.org
Fri Dec 7 15:54:49 PST 2007
Author: walter
Date: 2007-12-07 15:54:48 -0800 (Fri, 07 Dec 2007)
New Revision: 8540
Added:
seismo/2D/SPECFEM2D/trunk/DATA/Par_file_unstruct
Log:
added Par_file_unstruct for unstructured example.
Added: seismo/2D/SPECFEM2D/trunk/DATA/Par_file_unstruct
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file_unstruct 2007-06-26 01:02:08 UTC (rev 8539)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file_unstruct 2007-12-07 23:54:48 UTC (rev 8540)
@@ -0,0 +1,98 @@
+
+# title of job, and file that contains interface data
+title = Test for M2 UPPA
+interfacesfile = interfaces_M2_UPPA_curved.dat
+
+# data concerning mesh, when generated using third-party app (more info in README)
+read_external_mesh = .true.
+mesh_file = ./DATA/unstructured_fluide_solide_test/mesh # file containing the mesh
+nodes_coords_file = ./DATA/unstructured_fluide_solide_test/nodes_coords # file containing the nodes coordinates
+materials_file = ./DATA/unstructured_fluide_solide_test/mat # file containing the material number for each element
+free_surface_file = ./DATA/unstructured_fluide_solide_test/surface_free # file containing the free surface
+absorbing_surface_file = ./DATA/unstructured_fluide_solide_test/surface_abs # file containing the absorbing surface
+
+# parameters concerning partitionning
+nproc = 8 # number of processes
+partionning_method = 2 # ascending order = 1, Metis = 2, Scotch = 3
+partitionning_strategy = 01110 #b{strat=m{asc=g,low=g,rat=1.0,type=h,vert=4}} #01110 # options concerning partitionning strategy.
+
+# 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 = 4 # number of control nodes per element (4 or 9)
+initialfield = .false. # use a plane wave as source or not
+assign_external_model = .false. # define external earth model or not
+TURN_ANISOTROPY_ON = .false. # turn anisotropy on or off for solid medium
+TURN_ATTENUATION_ON = .false. # turn attenuation on or off for solid medium
+
+# absorbing boundaries parameters
+absorbing_conditions = .true. # absorbing boundary active or not
+absorbbottom = .true.
+absorbright = .true.
+absorbtop = .false.
+absorbleft = .true.
+
+# time step parameters
+nt = 18000 # total number of time steps
+deltat = 0.25d-3 # duration of a time step
+
+# source parameters
+source_surf = .false. # source inside the medium or at the surface
+xs = 0. # source location x in meters
+zs = -100. # source location z in meters
+source_type = 1 # elastic force or acoustic pressure = 1 or moment tensor = 2
+time_function_type = 1 # Ricker = 1, first derivative = 2, Gaussian = 3, Dirac = 4, Heaviside = 5
+f0 = 5.d0 # dominant source frequency (Hz) if not Dirac or Heaviside
+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 for seismograms
+seismotype = 1 # record 1=displ 2=veloc 3=accel 4=pressure
+generate_STATIONS = .true. # creates a STATION file in ./DATA
+nreceiverlines = 1 # number of receiver lines
+anglerec = 0.d0 # angle to rotate components at receivers
+
+# first receiver line
+nrec = 100 # number of receivers
+xdeb = 0. # first receiver x in meters
+zdeb = -30. # first receiver z in meters
+xfin = 13000. # last receiver x in meters (ignored if onlyone receiver)
+zfin = -30. # last receiver z in meters (ignored if onlyone receiver)
+enreg_surf = .false. # receivers inside the medium or at the surface
+
+# display parameters
+NTSTEP_BETWEEN_OUTPUT_INFO = 2000 # display frequency in time steps
+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
+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
+
+# velocity and density models
+nbmodels = 4 # nb of different models
+# define models as (model_number,1,rho,vp,vs,0,0) or (model_number,2,rho,c11,c13,c33,c44)
+# set vs to zero to make a given model acoustic
+# the mesh can contain both acoustic and elastic models simultaneously
+1 1 1000.d0 1500.d0 0.0d0 0 0
+2 1 2500.d0 5000.d0 2500.0d0 0 0
+3 1 2200.d0 2500.d0 1443.375d0 0 0
+4 1 2200.d0 2200.d0 1343.375d0 0 0
+# define the different regions of the model in the (nx,nz) spectral element mesh
+nbregions = 5 # nb of regions and model number for each
+1 80 1 20 1
+1 80 21 40 4
+1 80 41 60 3
+60 80 21 40 4
+30 40 50 60 2
\ No newline at end of file
More information about the cig-commits
mailing list