SPEED
Functions/Subroutines
MAKE_MECH_PROP_CASE_014.f90 File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine make_mech_prop_case_014 (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_014()

subroutine make_mech_prop_case_014 ( 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_014.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 ((depth .ge. 0.0d0) .and. (zs_all .ge. 0.0d0)) then
41
42 ! + MATERIAL INSIDE THE ALLUVIAL BASIN
43 if (depth .le. 100.0d0) then
44 vs = 300.d0
45 vp = 520.d0
46 rho = 2200.d0
47 lambda = rho * (vp**2 - 2.d0*vs**2)
48 mu = rho * vs**2
49 gamma = (3.1415*(3.d0/3.d0))/(30.d0)
50
51 elseif (depth .le. 200.0d0) then
52 vs = 400.d0
53 vp = 700.d0
54 rho = 2300.d0
55 lambda = rho * (vp**2 - 2.d0*vs**2)
56 mu = rho * vs**2
57 gamma = (3.1415*(2.d0/3.d0))/(40.d0)
58
59 elseif (depth .le. 500.0d0) then
60 vs = 500.d0
61 vp = 850.d0
62 rho = 2400.d0
63 lambda = rho * (vp**2 - 2.d0*vs**2)
64 mu = rho * vs**2
65 gamma = (3.1415*(2.d0/3.d0))/(50.d0)
66
67 else
68 vs = 1000.d0
69 vp = 1700.d0
70 rho = 2400.d0
71 lambda = rho * (vp**2 - 2.d0*vs**2)
72 mu = rho * vs**2
73 gamma = (3.1415*(2.d0/3.d0))/(100.d0)
74 endif
75
76 else
77 ! + MATERIAL INSIDE THE BEDROCK (Vs=1500m/s)
78 vs = 1500.d0
79 vp = 2800.d0
80 rho = 2400.d0
81 lambda = rho * (vp**2 - 2.d0*vs**2)
82 mu = rho * vs**2
83 gamma = (3.1415*(3.d0/3.d0))/(150.d0)
84 endif
85
86

Referenced by make_eltensor_for_cases().

Here is the caller graph for this function: