par.inp

Spectral Element Method for wave propagation and rupture dynamics.

# Lambs' problem, same as in Mercerat et al. 2006

#----- Some general parameters ----------------
&ECHO title = "Lamb's problem", verbose='1111'
    , ItInfo = 200, ItSnapshots = 5000 /
&GENERAL iexec=1, ngll=9, fmax=25d0 , ndof=2 /

#----- Build the mesh ---------------------------
&MESH_DEF  method = 'CARTESIAN' /
&MESH_CART xlim=0.d3,4.d3 ,zlim=-2.d3,0.d3 , nelem=40,20 /

#---- Elastic material parameters --------------

# homogeneous medium: rho, cp ,cs
&MATERIAL tag=1, mode='ISOTR' /
2000.d0 3200.d0 1847.5d0

#----- Boundary conditions ---------------------
# Assign tags and kinds to boundaries: (comment it if free) 
#The default kind is FREE SURFACE

&BC_DEF  tag = 1 , kind = 'ABSORB' /
&BC_ABSORB side = 'D' , stacey=T/

&BC_DEF  tag = 2 , kind = 'ABSORB' /
&BC_ABSORB side = 'R' , stacey=T/

&BC_DEF  tag = 4 , kind = 'ABSORB' /
&BC_ABSORB side = 'L' , stacey=T/

#---- Time scheme settings ----------------------
&TIME  TotalTime=1.5d0, Dt=0.5d-3 /

#---- Sources ----------------------------------

# a point force
&SRC_DEF TimeFunction= 'RICKER' ,coord= 1500.d0,-50.d0, mechanism= 'FORCE' /
&STF_RICKER  f0= 10.d0, onset = 0.15d0, ampli = 1.d0 /
&SRC_FORCE angle = 0d0/

#----- Receivers ---------------------------------
&REC_LINE number = 2 , field='D', first = 2200.d0,0d0, last = 2700.d0,0d0 / 

#--------- Plots settings ----------------------
&PLOTS fields='V'/
&PLOTS_POSTSCRIPT  vectors =F, color=T, interpol= T /