pump1.dat ; properties ; limestone (E = 7e10 v=0.25) pro bul 4.6667e10 shear 2.8e10 coh 2.0e7 ten 2e6 prop dens 2500 prop porosity 0.3 perm 1e-10 prop dil 2 fric 30 ; impermeable layer def depth depth1 = -17.2e3 depth2 = -17.8e3 loop i (1,100) loop j (1,100) if y(i,j) < depth1 then if y(i,j) > depth2 then k11(i,j) = 1.0e-15 k22(i,j) = 1.0e-15 friction(i,j) = 35 end_if end_if end_loop end_loop end depth * ; fault as a zone with ubiquitous joints pro jfric 5 jcoh 5e6 ja 50.63 jtens 2e6 i=35,66 j=50,51 prop perm 1e-9 i=35,66 j=50,51 water bulk 2e9 dens 1000