SPEED
Functions/Subroutines
MAKE_MECH_PROP_CASE_015.f90 File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine make_mech_prop_case_015 (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_015()

subroutine make_mech_prop_case_015 ( 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_015.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
41
42 if ((depth .ge. 0.0d0).and.(zs_all .ge. 0.0d0)) then
43 ! + MATERIAL INSIDE THE BASIN
44 vs = 100.d0 + 10.d0 * depth**(0.6d0) !100.d0 + 10.d0 * Depth**(0.6d0)
45 vp = dsqrt(10.d0)*vs
46 rho = 1530.d0 + 0.1d0*depth**(0.54d0)
47 lambda = rho * (vp**2 - 2*vs**2)
48 mu = rho * vs**2
49 qs = 0.1*vs
50 gamma = (3.1415*(2/3))/qs
51
52 else
53 ! + MATERIAL INSIDE THE BEDROCK
54 depth_real = zs
55 if (depth_real .ge. -500.0d0) then
56 vs = 1000.d0
57 vp = 1800.d0
58 rho = 2300.d0
59 lambda = rho * (vp**2 - 2*vs**2)
60 mu = rho * vs**2
61 qs = 0.1*vs
62 gamma = (3.1415*(2.d0/3.d0))/qs
63
64 elseif (depth_real .le. -500.d0 .and. depth_real .ge. -1000.0d0) then
65 vs = 1700.d0
66 vp = 3160.d0
67 rho = 2500.d0
68 lambda = rho * (vp**2 - 2*vs**2)
69 mu = rho * vs**2
70 qs = 0.1*vs
71 gamma = (3.1415*(2.d0/3.d0))/qs
72
73 else
74 vs = 2600.d0;
75 vp = 4830.d0;
76 rho = 2840.d0;
77 lambda = rho * (vp**2 - 2*vs**2)
78 mu = rho * vs**2
79 qs = 0.1*vs
80 gamma = (3.1415*(2.d0/3.d0))/qs
81
82
83 endif
84 endif
85
86
87
88

Referenced by make_eltensor_for_cases().

Here is the caller graph for this function: