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
41 if( vs30 .lt. 325.d0)
then
42 if ( depth .le. 160.d0)
then
43 vs = 250.d0 + 43.d0*depth**(0.5d0);
44 vp = 700.d0 + 45.d0*depth**(0.5d0);
46 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
49 gamma = 4.d0*datan(1.d0)/qs;
50 elseif(depth .le. 2000.d0)
then
51 vs = 800.d0 + 37.13d0*(depth-160d0)**(0.5)
53 rho = 1800 + 12.92d0*(depth-160d0)**(0.5);
54 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
57 gamma = 4.d0*datan(1.d0)/qs;
59 vs = 1350.d0 + 23.33*(depth)**(0.5);
61 rho = 2100 + 5.69*(depth)**(0.5);
62 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
65 gamma = 4.d0*datan(1.d0)/qs;
68 elseif (vs30 .lt. 500.d0)
then
70 if ( depth .le. 80.d0)
then
71 vs = 325.d0 + 30.74*depth**(0.5);
72 vp = 800.d0 + 42 *depth**(0.5);;
74 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
77 gamma = 4.d0*datan(1.d0)/qs;
78 elseif(depth .le. 120.d0)
then
79 vs = 600.d0 + 31.62*(depth-80.d0)**(0.5);
80 vp = 1175.d0 + 26.72*(depth-80.d0)**(0.5);
82 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
85 gamma = 4.d0*datan(1.d0)/qs;
86 elseif(depth .le. 250.d0)
then
87 vs = 800.d0 + 40.75*(depth-120.d0)**(0.5);
89 rho = 1850.d0 + 9.64*(depth-120.d0)**(0.5);
90 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
93 gamma = 4.d0*datan(1.d0)/qs;
94 elseif(depth .le. 2000.d0)
then
95 vs = 700.d0 + 38.14*(depth-30.d0)**(0.5);
97 rho = 1960.d0 + 9.43*(depth-250.d0)**(0.5);
98 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
101 gamma = 4.d0*datan(1.d0)/qs;
104 vs = 1350.d0 + 23.33*(depth)**(0.5);
106 rho = 2100 + 5.69*(depth)**(0.5);
107 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
110 gamma = 4.d0*datan(1.d0)/qs;
113 elseif (vs30 .lt. 700.d0)
then
115 if ( depth .le. 50.d0)
then
116 vs = 500.d0 + 42.42*depth**(0.5);
117 vp = 900.d0 + 42.42*depth**(0.5);
119 lambda = rho * (vp**2 - 2*vs**2);
122 gamma = 4.d0*datan(1.d0)/qs;
123 elseif ( depth .le. 250.d0)
then
124 vs = 800.d0 + 33.1*(depth-50.d0)**(0.5);
126 rho = 1900.d0 + 4.89*(depth-50.d0)**(0.5);
127 lambda = rho * (vp**2 - 2*vs**2);
130 gamma = 4.d0*datan(1.d0)/qs;
132 elseif(depth .le. 2000.d0)
then
133 vs = 700.d0 + 38.14*(depth-30.d0)**(0.5);
135 rho = 1960.d0 + 9.43*(depth-250.d0)**(0.5);
136 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
139 gamma = 4.d0*datan(1.d0)/qs;
141 vs = 1350.d0 + 23.33*(depth)**(0.5);
143 rho = 2100 + 5.69*(depth)**(0.5);
144 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
147 gamma = 4.d0*datan(1.d0)/qs;
150 elseif (vs30 .lt. 900.d0)
then
152 if ( depth .le. 4000.d0)
then
153 vs = 700.d0 + 37.9*(depth)**(0.5)
155 rho = 1960.d0 + 8.885*(depth)**(0.5)
156 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
159 gamma = 4.d0*datan(1.d0)/qs;
161 vs = 1350.d0 + 23.33*(depth)**(0.5);
163 rho = 2100 + 5.69*(depth)**(0.5);
164 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
167 gamma = 4.d0*datan(1.d0)/qs;
170 elseif (vs30 .lt. 1350.d0)
then
172 if ( depth .le. 2000.d0)
then
173 vs = 900.d0 + 33.38 * (depth)**(0.5);
175 rho = 2050.d0 + 215.1*(depth*0.001)**(0.5);
176 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
179 gamma = 4.d0*atan(1.d0)/qs;
181 vs = 1350.d0 + 23.33*(depth)**(0.5);
183 rho = 2100 + 5.69*(depth)**(0.5);
184 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
187 gamma = 4.d0*datan(1.d0)/qs;
191 vs = 1350.d0 + 23.33*(depth)**(0.5);
193 rho = 2100 + 5.69*(depth)**(0.5);
194 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
197 gamma = 4.d0*datan(1.d0)/qs;
247 if(dabs(xs - 576059.d0) .le. 2000.d0)
then
251 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
254 gamma = 4.d0*atan(1.d0)/qs;
256 elseif(dabs(xs - 740948.d0) .le. 2000.d0)
then
260 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
263 gamma = 4.d0*atan(1.d0)/qs;
265 elseif(dabs(ys - 4602206.d0) .le. 2000.d0)
then
269 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
272 gamma = 4.d0*atan(1.d0)/qs;
274 elseif(dabs(ys - 4502679.d0) .le. 2000.d0)
then
278 lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
281 gamma = 4.d0*atan(1.d0)/qs;