SPEED
Functions/Subroutines
MAKE_MECH_PROP_CASE_005.f90 File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine make_mech_prop_case_005 (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_005()

subroutine make_mech_prop_case_005 ( 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_005.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
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 ! + MATERIAL INSIDE THE ALLUVIAL BASIN - 2nd Layer
47 elseif (sub_tag_all.eq.2) then
48 rho = 2100
49 lambda = 6.289500e+09
50 mu = 2.572500e+08
51 gamma = 1.196797e-01
52
53 ! + MATERIAL INSIDE THE ALLUVIAL BASIN - 3rd Layer
54 elseif (sub_tag_all.eq.3) then
55 rho = 2200
56 lambda = 1.189100e+10
57 mu = 9.295000e+08
58 gamma = 6.444293e-02
59
60 ! + MATERIAL INTO THE BEDROCK
61 elseif (sub_tag_all.eq.4) then
62 if (zs.ge.-3000.0) then
63 vs = 0.4100*(-zs) + 2190.0000 !VS: S velocity in m/s
64 vp = 0.8100*(-zs) + 3690.0000 !VP: P velocity in m/s
65 rho = 0.0680*(-zs) + 2532.0000 !RHO: MASS DENSITY in kg/m^3
66 lambda = rho * (vp**2 - 2*vs**2)
67 mu = rho * vs**2
68 gamma = 1.6111e-02
69 else
70 vs = 0.0050*(-zs) + 3405.0000 !VS: S velocity in m/s
71 vp = 0.0050*(-zs) + 6105.0000 !VP: P velocity in m/s
72 rho = 0.0040*(-zs) + 2724.0000 !RHO: MASS DENSITY in kg/m^3
73 lambda = rho * (vp**2 - 2*vs**2)
74 mu = rho * vs**2
75 gamma = 1.6111e-02
76 endif
77
78 ! + MATERIAL INSIDE THE ALLUVIAL BASIN - 5th Layer (nu modified)
79 elseif (sub_tag_all.eq.5) then
80 rho = 2100
81 lambda = 1.260000e+08
82 mu = 8.400000e+07
83 gamma = 2.094395e-01
84
85 ! + MATERIAL INSIDE THE ALLUVIAL BASIN - 6th Layer (nu modified)
86 elseif (sub_tag_all.eq.6) then
87 rho = 2100
88 lambda = 3.858750e+08
89 mu = 2.572500e+08
90 gamma = 1.196797e-01
91
92 ! + MATERIAL INSIDE THE ALLUVIAL BASIN - 7th Layer (nu modified)
93 else
94 rho = 2200
95 lambda = 1.394250e+09
96 mu = 9.295000e+08
97 gamma = 6.444293e-02
98 endif
99
100

Referenced by make_eltensor_for_cases().

Here is the caller graph for this function: