24 xs, ys, zs, Depth, zs_all,&
25 vs30, thickness, sub_tag_all, f_distance)
27 real*8,
intent(out) :: rho, lambda, mu, gamma, qs, qp
28 real*8,
intent(in) :: xs, ys, zs, depth, zs_all,&
29 vs30, thickness, f_distance
30 integer*4 :: sub_tag_all
31 real*8 :: ni, vs, vp, depth_real
43 if ((depth.ge.0.0d0).and.(zs_all.ge.0.0d0))
then
45 if (depth .le. 150.0d0)
then
49 lambda = rho * (vp**2 - 2*vs**2)
52 gamma = (3.1415*(2.d0/3.d0))/qs
54 elseif(depth .gt. 150.d0 .and. depth .le. f_distance)
then
55 vs = 300.d0 + 10.d0*(depth-150.d0)**0.5
56 vp = 1500.d0 + 10.d0*(depth-150.d0)**0.5
57 rho = 1800.d0 + 6.d0*(depth-150.d0)**0.5
58 lambda = rho * (vp**2 - 2*vs**2)
61 gamma = (3.1415*(2.d0/3.d0))/qs
64 vs = 800.d0 + 15.d0*(depth-f_distance)**0.5
65 vp = 2000.d0 + 15.d0*(depth-f_distance)**0.5
66 rho = 2100.d0 + 4.d0*(depth-f_distance)**0.5
67 lambda = rho * (vp**2 - 2*vs**2)
70 gamma = (3.1415*(2/3))/qs
77 if (depth_real .le. 1000.0d0)
then
81 lambda = rho * (vp**2 - 2*vs**2)
83 gamma = (3.1415*(2/3))/(150.d0)
87 elseif (depth_real.le.3000.0d0)
then
91 lambda = rho * (vp**2 - 2*vs**2)
93 gamma = (3.1415*(2/3))/(200.d0)
98 elseif (depth_real.le.6000.0d0)
then
102 lambda = rho * (vp**2 - 2*vs**2)
104 gamma = (3.1415*(2/3))/(250.d0)
112 lambda = rho * (vp**2 - 2*vs**2)
114 gamma = (3.1415*(2/3))/(350.d0)
subroutine make_mech_prop_case_012(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all, f_distance)
Makes not-honoring technique.