SPEED
Functions/Subroutines
MAKE_MECH_PROP_CASE_028.f90 File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine make_mech_prop_case_028 (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_028()

subroutine make_mech_prop_case_028 ( 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,
real*8, intent(in)  sub_tag_all 
)

Makes not-honoring technique.

Mechanical properties given node by node.

Definition at line 24 of file MAKE_MECH_PROP_CASE_028.f90.

27
28 real*8, intent(out) :: rho, lambda, mu, gamma, qs, qp
29 real*8, intent(in) :: xs, ys, zs, depth, zs_all,&
30 vs30, thickness, sub_tag_all
31 real*8 :: ni, vs, vp, depth_real
32 real*8 :: vsini, vsfin, zini, zfin, rini, rfin
33
34 rho = 0.d0;
35 lambda = 0.d0;
36 mu = 0.d0;
37 gamma = 0.d0;
38 qs = 0.d0;
39 qp = 0.d0
40
41 if ((depth .ge. 0.0d0).and.(zs_all .ge. 0.0d0)) then
42
43 ! + material inside the basin
44 !! NEWER VS-RULE
45 if (depth .le. 150) then
46 vs = 281.64 + (2.0000*(548.33-281.64))/(1.0000+(15.0000/(depth+0.1000))**1.2900)
47 else
48 vs = 975
49 endif
50
51 vp = 1.855 * vs
52 rho = 1900 + vs / 1700 * 600
53 lambda = rho * (vp**2 - 2*vs**2)
54 mu = rho * vs**2
55 qs = 50
56 gamma = (3.1415*(2/3))/qs
57
58
59 else
60 ! + material inside the bedrock - first layer of crustal model
61 ! depth_real = zs
62
63 ! generic bedrock outcrop with Vs_min = 1000 m/s and Vs30 = 1200 m/s - from Cotton et al. 2006
64 if (depth .le. 1) then
65 vs = 1000
66 elseif(depth .le. 30) then
67 vsini = 1000
68 vsfin = 1421
69 zini = 1
70 zfin = 30
71 vs = vsini + (vsfin-vsini)*((depth-zini)/(zfin-zini))**0.50 !VS: S velocity in m/s
72 elseif(depth .le. 190) then
73 vsini = 1421
74 vsfin = 2217
75 zini = 30
76 zfin = 190
77 vs = vsini + (vsfin-vsini)*((depth-zini)/(zfin-zini))**0.50 !VS: S velocity in m/s
78 elseif(depth .le. 1000) then
79 vsini = 2217
80 vsfin = 2600
81 zini = 190
82 zfin = 1000
83 vs = vsini + (vsfin-vsini)*((depth-zini)/(zfin-zini))**0.50 !VS: S velocity in m/s
84 else
85 vs = 2600
86 endif
87
88 vp = sqrt(3.46)*vs !VP: P velocity in m/s - (Poisson = 0.3 approx)
89
90 if (depth .le. 100) then
91 rini = 2200;
92 rfin = 2500;
93 zini = 0;
94 zfin = 100;
95 rho = rini + (rfin-rini)*((depth-zini)/(zfin-zini))**0.50 !RHO: MASS DENSITY in kg/m^3
96 elseif(depth .le. 1000) then
97 rini = 2500;
98 rfin = 2840;
99 zini = 100;
100 zfin = 1000;
101 rho = rini + (rfin-rini)*((depth-zini)/(zfin-zini))**0.50 !RHO: MASS DENSITY in kg/m^3
102 else
103 rho = 2840.0d0
104 endif
105
106 qs = 0.1d0*vs
107 lambda = rho * (vp**2 - 2*vs**2)
108 mu = rho * vs**2
109 gamma = (3.1415*(2/3))/qs ! (Qs=Vs/10 f=0.67 Hz)
110
111
112 ! generic bedrock outcrop with Vs_min = 800 m/s and Vs30 = 900 m/s - SUPERSEDED
113 ! VS = 800.0d0 + (1700.0d0-800.0d0)*(0.001d0*Depth)**0.50 !VS: S velocity in m/s
114 ! VP = sqrt(3.46)*VS !VP: P velocity in m/s - (Poisson = 0.3 approx)
115 ! rho = 2200.0d0 + (2500.0d0-2200.0d0)*(0.001d0*Depth)**0.50 !RHO: MASS DENSITY in kg/m^3
116 ! qs = 0.1d0*VS
117 ! lambda = rho * (VP**2 - 2*VS**2)
118 ! mu = rho * VS**2
119 ! gamma = (3.1415*(2/3))/qs ! (Qs=Vs/10 f=0.67 Hz)
120
121 ! vs = 1700
122 ! vp = 3160
123 ! rho = 2500
124 ! lambda = rho * (vp**2 - 2*vs**2)
125 ! mu = rho * vs**2
126 ! qs = 200
127 ! gamma = (3.1415*(2/3))/qs
128 endif
129
130
131
132

Referenced by make_eltensor_for_cases().

Here is the caller graph for this function: