24 xs, ys, zs, Depth, zs_all,&
25 vs30, thickness, sub_tag_all)
27 real*8,
intent(out) :: rho, lambda, mu, gamma, qs, qp
28 real*8,
intent(in) :: xs, ys, zs, depth, zs_all,&
30 integer*4 :: sub_tag_all
31 real*8 :: ni, vs, vp, depth_real
42 if ((depth .ge. 0.0d0) .and. (zs_all.ge.0.0d0))
then
49 vs = 300.0d0 + 30.0d0*(depth)**0.67
50 vp = 600.0d0 + 30.0d0*(depth)**0.77
51 rho = 1900.0d0 + 1.25*depth
52 lambda = rho * (vp**2 - 2*vs**2)
subroutine make_mech_prop_case_006(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.