SPEED
Functions/Subroutines
MAKE_MECH_PROP_CASE_019.f90 File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine make_mech_prop_case_019 (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_019()

subroutine make_mech_prop_case_019 ( 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_019.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 if (vs30 .le. 300.0d0) then
42 vp30 = 1800
43 else
44 vp30 = 2000
45 endif
46
47
48 if ((depth.ge.0.0d0).and.(zs_all.ge.0.0d0)) then
49 ! + MATERIAL INSIDE THE BASIN
50
51 vs = vs30 + (2000.0d0-vs30)*(depth/1000.0d0)**(0.70d0)
52 vp = vp30 + (4500.0d0-vp30)*(depth/1000.0d0)**(0.70d0)
53 rho = 2000.d0 + 0.40d0*depth
54
55
56 if (depth .le. 50.0d0) then
57 qs = 20
58 elseif (depth .le. 200.0d0) then
59 qs = 50
60 elseif (depth .le. 500.0d0) then
61 qs = 100
62 else
63 qs = 150
64 endif
65
66 lambda = rho * (vp**2 - 2*vs**2)
67 mu = rho * vs**2
68 gamma = (3.1415*(2/3))/qs ! max freq of 2 Hz
69
70 else
71 ! + MATERIAL INSIDE THE BEDROCK
72
73 vs = 2000.d0
74 vp = 4500.d0
75 rho = 2400.d0
76 qs = 200.d0
77 lambda = rho * (vp**2 - 2*vs**2)
78 mu = rho * vs**2
79 gamma = (3.1415*(2/3))/qs
80 endif
81
82

Referenced by make_eltensor_for_cases().

Here is the caller graph for this function: