SPEED
MAKE_MECH_PROP_CASE_070.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_070(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 if (sub_tag_all.eq.1) then
41 !rho = 2100
42 !lambda = 4.557000E+09
43 !mu = 8.400000E+07
44 !gamma = 2.094395E-01
45
46 vs = 300.0d0 + 36.0d0*(depth)**0.43 !VS: S velocity in m/s
47 vp = sqrt(4.57)*vs !VP: P velocity in m/s (nu = 0.36)
48 rho = 1900.0d0 !RHO: MASS DENSITY in kg/m^3
49 lambda = rho * (vp**2 - 2*vs**2)
50 mu = rho * vs**2
51 gamma = 6.9813e-02 ! (Qs=Vs/10=30)
52
53
54
55
56
57 !! + MATERIAL INSIDE THE ALLUVIAL BASIN - 2nd Layer
58 !elseif (sub_tag_all.eq.2) then
59 ! rho = 2100
60 ! lambda = 6.289500E+09
61 ! mu = 2.572500E+08
62 ! gamma = 1.196797E-01
63
64 !! + MATERIAL INSIDE THE ALLUVIAL BASIN - 3rd Layer
65 !elseif (sub_tag_all.eq.3) then
66 ! rho = 2200
67 ! lambda = 1.189100E+10
68 ! mu = 9.295000E+08
69 ! gamma = 6.444293E-02
70
71 ! + MATERIAL INTO THE BEDROCK
72 elseif (sub_tag_all.eq.2) then
73
74
75 vs = 800.0d0 + (1700.0d0-800.0d0)*(0.001d0*depth)**0.50
76 vp = sqrt(3.46)*vs
77 rho = 2200.0d0 ! + (2500.0d0-2200.0d0)*(0.001d0*Depth)**0.50
78 ! + MATERIAL INSIDE THE BEDROCK (Vs=1700m/s)
79 !lambda = 1.0514E+10
80 !mu = 7.2250E+09
81 !rho = 2500.0d0
82 lambda = rho * (vp**2 - 2*vs**2)
83 mu = rho * vs**2
84 gamma = 2.0944e-02
85 else
86 vs=1700 !m/s
87 vp=3160
88 rho=2500
89 lambda = 1.0514e+10
90 mu = 7.2250e+09
91 rho = 2500.0d0
92 gamma = 2.0944e-02
93
94
95
96
97
98
99
100
101 ! if (zs.ge.-3000.0) then
102 ! VS = 0.4100*(-zs) + 2190.0000 !VS: S velocity in m/s
103 ! VP = 0.8100*(-zs) + 3690.0000 !VP: P velocity in m/s
104 ! rho = 0.0680*(-zs) + 2532.0000 !RHO: MASS DENSITY in kg/m^3
105 ! lambda = rho * (VP**2 - 2*VS**2)
106 ! mu = rho * VS**2
107 ! gamma = 1.6111E-02
108 ! else
109 ! VS = 0.0050*(-zs) + 3405.0000 !VS: S velocity in m/s
110 ! VP = 0.0050*(-zs) + 6105.0000 !VP: P velocity in m/s
111 ! rho = 0.0040*(-zs) + 2724.0000 !RHO: MASS DENSITY in kg/m^3
112 ! lambda = rho * (VP**2 - 2*VS**2)
113 ! mu = rho * VS**2
114 ! gamma = 1.6111E-02
115
116 endif
117
118 !! + MATERIAL INSIDE THE ALLUVIAL BASIN - 5th Layer (nu modified)
119 !elseif (sub_tag_all.eq.5) then
120 ! rho = 2100
121 ! lambda = 1.260000E+08
122 ! mu = 8.400000E+07
123 ! gamma = 2.094395E-01
124
125 !! + MATERIAL INSIDE THE ALLUVIAL BASIN - 6th Layer (nu modified)
126 !elseif (sub_tag_all.eq.6) then
127 ! rho = 2100
128 ! lambda = 3.858750E+08
129 ! mu = 2.572500E+08
130 ! gamma = 1.196797E-01
131
132 !! + MATERIAL INSIDE THE ALLUVIAL BASIN - 7th Layer (nu modified)
133 !else
134 ! rho = 2200
135 ! lambda = 1.394250E+09
136 ! mu = 9.295000E+08
137 ! gamma = 6.444293E-02
138 !endif
139
140
141 end subroutine make_mech_prop_case_070
make_mech_prop_case_070
subroutine make_mech_prop_case_070(rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
Makes not-honoring technique.
Definition MAKE_MECH_PROP_CASE_070.f90:26