;
tit
FLAC Modelling Input File
;
config extra 4
grid  100   100 
gen 0.00e0,0.00e0 0.00e0,1.00e3 1.00e3,1.00e3 1.00e3,0.00e0 i=1,101 j=1,101
mark i=1
mark j=1
mark i = 101 
mark j = 101 
;
*************
model elastic
*************
;
tab  1   0 , 1000   500. , 1000   500. , 0 
tab  1   0 , 0   0 , 1000 
;
;
gen tab  1 
;
;
; Properties for each sub-regions: 
;  Basalt
prop dens= 2720  bulk=4.09e10 shear=2.45e10 reg= 9 , 32 
;  gnt
prop dens= 2720  bulk=4.09e10 shear=2.45e9 reg= 56 , 76 
;
*****
def inistr_phg
loop i (1,igp-1)
loop j (1,jgp-1)
jno = jgp - 1
If j = jno Then
hh1 = (abs(y(i,j)-y(i,j+1))+abs(y(i+1,j)-y(i+1,j+1)))*0.25
sss1 = density(i, j) * 9.8 * hh1
sss2 = sss1
Else
jno1 = j
jno2 = j + 1
hh1 = (abs(y(i,j)-y(i,j+1))+abs(y(i+1,j)-y(i+1,j+1)))*0.25
sss1 = density(i, j) * 9.8 * hh1
sss2 = sss1
loop jj (jno2,jgp-1)
aaa = abs(y(i, jj) - y(i, jj + 1))
bbb = abs(y(i + 1, jj) - y(i + 1, jj + 1))
ccc = 0.5 * (aaa + bbb)
sss3 = density(i, jj) * 9.8 * ccc
sss2 = sss2 + sss3
end_loop
end_if
syy(i, j) = -sss2
sxx(i, j) = -sss2
szz(i, j) = -sss2
end_loop
end_loop
end
inistr_phg
*
;
;
set large
;
;
fix x i=1
fix y i=1
fix x i= 101 
fix y i= 101 
fix x j=1
fix y j=1
fix x j= 101 
fix y j= 101 
ini xvel -0.025  var = 0.0, 0.05  i=1, 101  j=1, 101 
;
step 500
save melt1.sav ;    shear strain= 0.02 
step 500
save melt2.sav ;    shear strain= 0.04 
step 500
save melt3.sav ;    shear strain= 0.06 
step 500
save melt4.sav ;    shear strain= 0.08 
step 500
save melt5.sav ;    shear strain= 0.1 
step 500

return
stop