model mohr
prop dens 1400 bulk 1.5e9 shear 0.9e9 coh 0.7e6 ten 1.2e6 fric 24 ran group MEI
prop dens 1500 bulk 0.38e9 shear 2.35e8 coh 0.19e6 ten 0.01e6 fric 17 ran group DC
prop dens 2810 bulk 15.3e9 shear 8.6e9 coh 16.9e6 ten 8.7e6 fric 39 ran group XSY
prop dens 2740 bulk 13.1e9 shear 7.2e9 coh 8.9e6 ten 14.5e6 fric 33 ran group ZSY
prop dens 1560 bulk 7.6e9 shear 4.7e9 coh 3.8e6 ten 2.6e6 fric 35 ran group NY
prop dens 2620 bulk 9.7e9 shear 5.6e9 coh 4.5e6 ten 3.3e6 fric 28 ran group SZNY
prop dens 2690 bulk 11.8e9 shear 6.5e9 coh 5.7e6 ten 12.8e6 fric 30 ran group FSY
fix x y z ran z 0
;fix x ran x 0
;fix x ran x 200
;fix y ran y 0
;fix y ran y 100
ini szz -14.5e6 grad 0 0 2.5e4 ran z 0 80
ini sxx -17.4e6 grad 0 0 3e4
ini syy -17.4e6 grad 0 0 3e4
app szz -12.5e6 ran z 79.9 80.1
app sxx -17.4e6 grad 0 0 3e4 ran x 0
app sxx -17.4e6 grad 0 0 3e4 ran x 200
app syy -17.4e6 grad 0 0 3e4 ran y 0
app syy -17.4e6 grad 0 0 3e4 ran y 100
set grav 0 0 -10
solve