SPEED
Functions/Subroutines
MAKE_MECH_PROP_CASE_027.f90 File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine make_mech_prop_case_027 (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_027()

subroutine make_mech_prop_case_027 ( 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 23 of file MAKE_MECH_PROP_CASE_027.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, sub_tag_all
30 real*8 :: ni, vs, vp, depth_real
31
32 rho = 0.d0;
33 lambda = 0.d0;
34 mu = 0.d0;
35 gamma = 0.d0;
36 qs = 0.d0;
37 qp = 0.d0;
38 vs = 0.d0;
39 vp = 0.d0;
40
41
42 if ((depth .ge. 0.0d0) .and. (zs_all .ge. 0.0d0)) then
43 ! Option #1
44 vs = 300.0d0 + 36.0d0*(depth)**0.43 !VS: S velocity in m/s
45 vp = sqrt(4.57)*vs !VP: P velocity in m/s (nu = 0.36)
46 rho = 1900.0d0 !RHO: MASS DENSITY in kg/m^3
47 lambda = rho * (vp**2 - 2*vs**2)
48 mu = rho * vs**2
49 gamma = 6.9813e-02 ! (Qs=Vs/10=30 f=0.67 Hz)
50 else
51 ! + MATERIAL INSIDE THE BEDROCK (Vs=1700m/s)
52 ! generic bedrock outcrop with Vs_min = 800 m/s and Vs30 = 900 m/s
53 vs = 800.0d0 + (1700.0d0-800.0d0)*(0.001d0*depth)**0.50 !VS: S velocity in m/s
54 vp = sqrt(3.46)*vs !VP: P velocity in m/s - (Poisson = 0.3 approx)
55 rho = 2200.0d0 + (2500.0d0-2200.0d0)*(0.001d0*depth)**0.50 !RHO: MASS DENSITY in kg/m^3
56 qs = 0.1d0*vs
57 lambda = rho * (vp**2 - 2*vs**2)
58 mu = rho * vs**2
59 gamma = (3.1415*(2/3))/qs ! (Qs=Vs/10 f=0.67 Hz)
60
61 ! Old bedrock outcrop
62 !lambda = 1.0514E+10
63 !mu = 7.2250E+09
64 !rho = 2500.0d0
65 !gamma = 2.0944E-02
66
67 endif
68

Referenced by make_eltensor_for_cases().

Here is the caller graph for this function: