SPEED
MAKE_MECH_PROP_CASE_006.f90
Go to the documentation of this file.
1! Copyright (C) 2012 The SPEED FOUNDATION
2! Author: Ilario Mazzieri
3!
4! This file is part of SPEED.
5!
6! SPEED is free software; you can redistribute it and/or modify it
7! under the terms of the GNU Affero General Public License as
8! published by the Free Software Foundation, either version 3 of the
9! License, or (at your option) any later version.
10!
11! SPEED is distributed in the hope that it will be useful, but
12! WITHOUT ANY WARRANTY; without even the implied warranty of
13! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14! Affero General Public License for more details.
15!
16! You should have received a copy of the GNU Affero General Public License
17! along with SPEED. If not, see <http://www.gnu.org/licenses/>.
18
19
21
22
23 subroutine make_mech_prop_case_006(rho, lambda, mu, gamma, qs, qp, & !outputs
24 xs, ys, zs, Depth, zs_all,&
25 vs30, thickness, sub_tag_all)
26
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
30 integer*4 :: sub_tag_all
31 real*8 :: ni, vs, vp, depth_real
32
33 rho = 0.d0;
34 lambda = 0.d0;
35 mu = 0.d0;
36 gamma = 0.d0;
37 qs = 0.d0;
38 qp = 0.d0
39
40
41
42 if ((depth .ge. 0.0d0) .and. (zs_all.ge.0.0d0)) then
43 ! Option #1
44 ! VS = 300.0d0 + 6.50d0*Depth !VS: S velocity in m/s
45 ! VP = 600.0d0 + 12.0d0*Depth !VP: P velocity in m/s
46 ! rho = 1900.0d0 + 1.25*Depth !RHO: MASS DENSITY in kg/m^3
47
48 ! Option #2
49 vs = 300.0d0 + 30.0d0*(depth)**0.67 !VS: S velocity in m/s
50 vp = 600.0d0 + 30.0d0*(depth)**0.77 !VP: P velocity in m/s
51 rho = 1900.0d0 + 1.25*depth !RHO: MASS DENSITY in kg/m^3
52 lambda = rho * (vp**2 - 2*vs**2)
53 mu = rho * vs**2
54 gamma = 3.1416e-02 ! (Qs=100)
55 else
56 ! + MATERIAL INSIDE THE BEDROCK (Vs=1500m/s)
57 lambda = 1.0350e+10
58 mu = 5.1750e+09
59 rho = 2300.0d0
60 gamma = 2.0944e-02
61 endif
62
63 end subroutine make_mech_prop_case_006
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.