par.inp

Spectral Element Method for wave propagation and rupture dynamics.

# Parameter file for SEM2DPACK 2.2 
# Test for the Slip Weakening Fault boundary condition

#----- Some general parameters ----------------
&ECHO title = 'SCEC Test#3 for slip weakening fault', verbose='1011'
    , ItInfo = 400, ItSnapshots = 100000 /
&GENERAL iexec=1, ngll=5, fmax=3.d0 , ndof=2 /

#----- Build the mesh ---------------------------
&MESH_DEF  method = 'CARTESIAN' /
#&MESH_CART xlim=0d3,32.d3, zlim=0d3,50.d3, nelem=102,160 /
#&MESH_CART xlim=0d3,25.6d3, zlim=0d3,40.d3, nelem=102,160 /
&MESH_CART xlim=0d3,25.5d3, zlim=0d3,40.d3, nelem=102,160 /
#&MESH_CART xlim=0d3,25.5d3, zlim=0d3,40.d3, nelem=50,80/

#---- Material parameters --------------
&MATERIAL tag=1, mode='ISOTR' /
2670.d0 6000.d0 3464.d0

&KELVIN_VOIGT etaH='GAUSSIAN' /
&DIST_GAUSSIAN length=1d6,282.094, ampli=0.1d0 /

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

&BC_DEF  tag = 1 , kind = 'SWFFLT' /
&BC_SWFFLT Dc=0.4d0, MuSHet='PWCONR', MuD=0.525d0,
           Tn=-120.d6,TtHet='PWCONR', IxOut=1,65,4 /
&DIST_PWCONR num=2/ # Static friction coefficient
     15d3
0.677d0  1d4
&DIST_PWCONR num=2/ # Initial shear stress
     1.5d3
81.6d6 70.d6

&BC_DEF  tag = 2 , kind = 'DT0TN0' /  
#&BC_DEF  tag = 2 , kind = 'ABSORB' /
#&BC_ABSORB side = 'R' /
 
&BC_DEF  tag = 3 , kind = 'ABSORB' /  
&BC_ABSORB side = 'U', stacey=T /

&BC_DEF  tag = 4 , kind = 'DT0TN0' /  

#---- Time scheme settings ----------------------
&TIME  TotalTime=16.d0, courant = 0.55d0, kind='leapfrog' /
#&TIME  TotalTime=16.d0, courant = 0.55d0, kind='newmark' /
#&TIME_NEWMARK rho=0.5d0 /
&TIME_NEWMARK alpha=1d0, beta=0d0, gamma=0.5d0 /

#----- Receivers ---------------------------------
#&REC_LINE number = 10 , first = 0d3,0.5d3, last = 15.d3,0.5d3, isamp=20 / 

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