;
tit
FLAC Modelling Input File
;
config extra 4
grid  150   150 
gen 0.00e0,0.00e0 0.00e0,1.00e3 1.00e3,1.00e3 1.00e3,0.00e0 i=1,151 j=1,151
mark i=1
mark j=1
mark i = 151 
mark j = 151 
;
**********
model mohr
**********
;
tab  1   0 , 1000   500. , 1000   500.0 , 0 
tab  1   0 , 0   0 , 1000  
;
;
gen tab  1 
;
;
; Properties for each sub-regions: 
;  Basalt
pro dens= 2720  bulk=4.09e10 shear=2.45e9 reg= 15 , 60 
pro coh=3.00e7 tens=1.50e7 reg= 15 , 60 
pro fric=30 dil= 0  reg= 15 , 60 
;  granite
pro dens= 2720  bulk=4.09e10 shear=2.45e8 reg= 140 , 121 
pro coh=3.00e7 tens=1.50e7 reg= 140 , 121 
pro fric=30 dil= 0  reg= 140 , 121 
;
; this sets stress everywhere to 2kb and gravity is zero
ini sxx=-2e8
ini syy=-2e8
ini szz=-2e8
ini sxy=0e6
;
;
set gravity=0.0
;
set large
;
;
fix x i=1
fix y i=1
fix x i= 151 
fix y i= 151 
fix x j=1
fix y j=1
fix x j= 151 
fix y j= 151 

; Schaubs: make sure to put these commands after the fish function which computes your initial stress
hist 1 unbal nstep 1
set small
set mech on fl off
step 200
;ini   xd 0    yd 0
;ini xvel 0  yvel 0 
;set mech on fl on
;step 200  ; this can be increased if you think the unbalanced forces haven't bottomed out
;ini   xd 0    yd 0 
;ini xvel 0  yvel 0
;pl hold hist 1
set large
;initial velocity command should be here ie ini xvel 0.004 etc

; commands to calculate meanstress
define meanstress
while_stepping
loop i (1,izones)
 loop j (1,jzones)
   ex_1(i,j)=(sxx(i,j)+syy(i,j)+szz(i,j))/3.0
  end_loop
 end_loop
end
meanstress


ini xvel -0.0025  var = 0.0, 0.005  i=1, 151  j=1, 151 
;
step 100
save Inc2.sav
step 2000
save Inc2a.sav ;    shear strain= 0.1
step 2000
save Inc2b.sav ;    shear strain= 0.2 
step 2000
save Inc2c.sav ;    shear strain= 0.3 
step 2000
save Inc2d.sav ;    shear strain= 0.4 
step 2000
save Inc2e.sav ;    shear strain= 0.5 
step 2000
save Inc2f.sav ;    shear strain= 0.6 
step 2000
save Inc2g.sav ;    shear strain= 0.7 
step 2000
save Inc2h.sav ;    shear strain= 0.8 
step 2000
save Inc2i.sav ;    shear strain= 0.9 
step 2000
save Inc2j.sav ;    shear strain= 1 

return
stop