SPEED
Functions/Subroutines
MAKE_MECH_PROP_CASE_029.f90 File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine make_mech_prop_case_029 (rho, lambda, mu, gamma, qs, qp, xs, ys, zs, depth, zs_all, vs30, thickness, sub_tag_all)
 Makes not-honoring technique.
 

Function/Subroutine Documentation

◆ make_mech_prop_case_029()

subroutine make_mech_prop_case_029 ( real*8, intent(out)  rho,
real*8, intent(out)  lambda,
real*8, intent(out)  mu,
real*8, intent(out)  gamma,
real*8, intent(out)  qs,
real*8, intent(out)  qp,
real*8, intent(in)  xs,
real*8, intent(in)  ys,
real*8, intent(in)  zs,
real*8, intent(in)  depth,
real*8, intent(in)  zs_all,
real*8, intent(in)  vs30,
real*8, intent(in)  thickness,
integer*4  sub_tag_all 
)

Makes not-honoring technique.

Mechanical properties given node by node.

Definition at line 23 of file MAKE_MECH_PROP_CASE_029.f90.

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 real*8 :: vsini, vsfin, zini, zfin, rini, rfin
33
34
35 rho = 0.d0;
36 lambda = 0.d0;
37 mu = 0.d0;
38 gamma = 0.d0;
39 qs = 0.d0;
40 qp = 0.d0
41
42
43 if (vs30 .le. 300.0d0) then
44 vp30 = 1800
45 else
46 vp30 = 2000
47 endif
48
49
50 if ((depth.ge.0.0d0).and.(zs_all.ge.0.0d0)) then
51 ! + MATERIAL INSIDE THE BASIN
52
53 vs = vs30 + (2000.0d0-vs30)*(depth/1000.0d0)**(0.70d0)
54 vp = vp30 + (4500.0d0-vp30)*(depth/1000.0d0)**(0.70d0)
55 rho = 2000.d0 + 0.40d0*depth
56
57
58 if (depth .le. 50.0d0) then
59 qs = 20
60 elseif (depth .le. 200.0d0) then
61 qs = 50
62 elseif (depth .le. 500.0d0) then
63 qs = 100
64 else
65 qs = 150
66 endif
67
68 lambda = rho * (vp**2 - 2*vs**2)
69 mu = rho * vs**2
70 gamma = (3.1415*(2/3))/qs ! max freq of 2 Hz
71
72 else
73 ! + MATERIAL INSIDE THE BEDROCK
74
75 ! generic bedrock outcrop with Vs30 = 1500 m/s - from Cotton et al. 2006
76 if (depth .le. 1) then
77 vs = 1144
78 elseif(depth .le. 30) then
79 vsini = 1144
80 vsfin = 1696
81 zini = 1
82 zfin = 30
83 vs = vsini + (vsfin-vsini)*((depth-zini)/(zfin-zini))**0.50 !VS: S velocity in m/s
84 elseif(depth .le. 190) then
85 vsini = 1696
86 vsfin = 2381
87 zini = 30
88 zfin = 190
89 vs = vsini + (vsfin-vsini)*((depth-zini)/(zfin-zini))**0.50 !VS: S velocity in m/s
90 elseif(depth .le. 4000) then
91 vsini = 2381
92 vsfin = 3454
93 zini = 190
94 zfin = 4000
95 vs = vsini + (vsfin-vsini)*((depth-zini)/(zfin-zini))**0.50 !VS: S velocity in m/s
96 else
97 vs = 3440
98 endif
99
100 vp = 2.25*vs !VP: P velocity in m/s - (Poisson = 0.3 approx)
101
102
103 if (depth .le. 100) then
104 rini = 2200;
105 rfin = 2400;
106 zini = 0;
107 zfin = 100;
108 rho = rini + (rfin-rini)*((depth-zini)/(zfin-zini))**0.50 !RHO: MASS DENSITY in kg/m^3
109 elseif(depth .le. 1000) then
110 rini = 2400;
111 rfin = 2700;
112 zini = 100;
113 zfin = 1000;
114 rho = rini + (rfin-rini)*((depth-zini)/(zfin-zini))**0.50 !RHO: MASS DENSITY in kg/m^3
115 else
116 rho = 2700.0d0
117 endif
118
119
120
121 qs = 0.1d0*vs
122
123
124! VS = 2000.d0
125! VP = 4500.d0
126! rho = 2400.d0
127! qs = 200.d0
128
129 lambda = rho * (vp**2 - 2*vs**2)
130 mu = rho * vs**2
131 gamma = (3.1415*(2/3))/qs
132
133
134 endif
135
136

Referenced by make_eltensor_for_cases().

Here is the caller graph for this function: