par.inp

Spectral Element Method for wave propagation and rupture dynamics.

# Parameter file for SEM2DPACK 2.0 
# Anisotropic test from figure 7 of Komatitsch, Barnes and Tromp (Geophysics, 2000)

#----- Some general parameters ----------------
&ECHO title = 'Anisotropic test', verbose='1111'
    , ItInfo = 200, ItSnapshots = 200 /
&GENERAL iexec=0, ngll=6, fmax=400.d3 , ndof=2 /

#----- Build the mesh ---------------------------

# a box:
&MESH_DEF  method = 'CARTESIAN' /
&MESH_CART xlim=-.325d0,.325d0 ,zlim=-.325d0,.325d0 , nelem=130,130 /

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

# define domains, limits (ex,ez) given as element indices
# ex=1 : left
# ez=1 : bottom
&MESH_CART_DOMAIN tag=1, ex= 1,65, ez=1,130 /
&MESH_CART_DOMAIN tag=2, ex= 66,130, ez=1,130 /

# define material properties:

&MATERIAL tag=1, mode='ANISO' /
7100d0 16.5e10 5.0e10 6.2e10 3.96e10
# rho  c11 c13 c33 c44 

&MATERIAL tag=2, mode='ISOTR' /
7100d0 4820.729884542206d0 2361.665680982894d0
# rho vP vS

#----- Boundary conditions ---------------------

# if none given: free boundary

# for periodic boundaries:
#&BC_DEF  tags =1,3 ,kind = 'PERIOD' /
#&BC_DEF  tags =2,4 ,kind = 'PERIOD' /

# for absorbing boundaries:
#&BC_DEF  tag = 1 , kind = 'ABSORB' /
#&BC_ABSORB side = 'D' , stacey=F/

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

#&BC_DEF  tag = 3 , kind = 'ABSORB' /  
#&BC_ABSORB side = 'U' , stacey=F/

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

#---- Time scheme settings ----------------------
&TIME  TotalTime=1.d-4, dt= 5d-8 /

#---- Sources ----------------------------------
&SRC_DEF TimeFunction= 'RICKER' ,coord= -0.02d0,0.d0 
mechanism= 'FORCE' /
&STF_RICKER  f0= 170.d3, onset = 6.d-6, ampli = 1.d13 /
&SRC_FORCE angle = 180d0/

#----- Receivers ---------------------------------
&REC_LINE field='D', number = 49 , first = -12.d-2,-8.d-2, last = 12.d-2,-8.d-2, isamp=1 / 

#--------- Plots settings ----------------------
&PLOTS fields='D' /
&PLOTS_POSTSCRIPT  vectors =F, color=T, interpol= F / 
#,ScaleField = 1.d-2 /