SPEED
Functions/Subroutines
MAKE_MECH_PROP_CASE_012.f90 File Reference

Go to the source code of this file.

Functions/Subroutines

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

Function/Subroutine Documentation

◆ make_mech_prop_case_012()

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

Makes not-honoring technique.

Mechanical properties given node by node.

Definition at line 23 of file MAKE_MECH_PROP_CASE_012.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, f_distance
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
43 if ((depth.ge.0.0d0).and.(zs_all.ge.0.0d0)) then
44 ! + MATERIAL INSIDE THE BASIN
45 if (depth .le. 150.0d0) then
46 vs = 300.d0
47 vp = 1500.d0
48 rho = 1800.d0
49 lambda = rho * (vp**2 - 2*vs**2)
50 mu = rho * vs**2
51 qs = 0.1*vs;
52 gamma = (3.1415*(2.d0/3.d0))/qs !hy: fpeak = 2/3 Hz
53
54 elseif(depth .gt. 150.d0 .and. depth .le. f_distance) then
55 vs = 300.d0 + 10.d0*(depth-150.d0)**0.5
56 vp = 1500.d0 + 10.d0*(depth-150.d0)**0.5
57 rho = 1800.d0 + 6.d0*(depth-150.d0)**0.5
58 lambda = rho * (vp**2 - 2*vs**2)
59 mu = rho * vs**2
60 qs = 0.1*vs;
61 gamma = (3.1415*(2.d0/3.d0))/qs
62
63 else
64 vs = 800.d0 + 15.d0*(depth-f_distance)**0.5
65 vp = 2000.d0 + 15.d0*(depth-f_distance)**0.5
66 rho = 2100.d0 + 4.d0*(depth-f_distance)**0.5
67 lambda = rho * (vp**2 - 2*vs**2)
68 mu = rho * vs**2
69 qs = 0.1*vs
70 gamma = (3.1415*(2/3))/qs
71
72 endif
73
74 else
75 ! + MATERIAL INSIDE THE BEDROCK
76 depth_real = abs(zs)
77 if (depth_real .le. 1000.0d0) then
78 vs = 1200.d0
79 vp = 2300.d0
80 rho = 2100.d0
81 lambda = rho * (vp**2 - 2*vs**2)
82 mu = rho * vs**2
83 gamma = (3.1415*(2/3))/(150.d0)
84 !if(damping_type .eq. 2) then
85 ! qs=120.d0; qp=230.d0;
86 !endif
87 elseif (depth_real.le.3000.0d0) then
88 vs = 2100.d0
89 vp = 3500.d0
90 rho = 2200.d0
91 lambda = rho * (vp**2 - 2*vs**2)
92 mu = rho * vs**2
93 gamma = (3.1415*(2/3))/(200.d0)
94 !if(damping_type .eq. 2) then
95 ! qs=200.d0; qp=400.d0;
96 !endif
97
98 elseif (depth_real.le.6000.0d0) then
99 vs = 2750.d0
100 vp = 4750.d0
101 rho = 2400.d0
102 lambda = rho * (vp**2 - 2*vs**2)
103 mu = rho * vs**2
104 gamma = (3.1415*(2/3))/(250.d0)
105 !if(damping_type .eq. 2) then
106 ! qs=250.d0; qp=500.d0;
107 !endif
108 else
109 vs = 3670.d0
110 vp = 6340.d0
111 rho = 2800.d0
112 lambda = rho * (vp**2 - 2*vs**2)
113 mu = rho * vs**2
114 gamma = (3.1415*(2/3))/(350.d0)
115 !if(damping_type .eq. 2) then
116 ! qs=350.d0; qp=700.d0;
117 !endif
118
119 endif
120 endif
121
122
123
124

Referenced by make_eltensor_for_cases().

Here is the caller graph for this function: