Functions/Subroutines
subroutine	make_mech_prop_case_091 (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_091()

subroutine make_mech_prop_case_091	(	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(inout)	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_091.f90.

                      
real*8, intent(out)   :: rho, lambda, mu, gamma, qs, qp
real*8, intent(in)    :: xs, ys, zs, zs_all,&
                         vs30,thickness
integer*4             :: sub_tag_all                                        
real*8, intent(inout) :: depth                                  
real*8                :: ni, vs, vp, depth_real                                                 
                      
rho    = 0.d0;
lambda = 0.d0;
mu     = 0.d0;
gamma  = 0.d0;
qs     = 0.d0;
qp     = 0.d0
 
if (depth .lt. 0.0d0) depth = 0.0d0
 
if (depth .lt. 1000.d0) then
     vs = 700.d0 + 39*(depth)**(0.5)
     vp  = vs*1.8;
     rho = 1960.d0 + 11*(depth)**(0.5) 
else
     vs = 1930.d0
     vp  = vs*1.8;
     rho = 2300.d0 
endif
 
lambda = rho * (vp**2.d0 - 2.d0*vs**2.d0);
mu = rho * vs**2;
qs = 0.1*vs;
qp = 0.1*vp;
gamma = 4.d0*datan(1.d0)/qs;
 
 
                                                     

Referenced by make_eltensor_for_cases().

Here is the caller graph for this function:

Functions/Subroutines

Function/Subroutine Documentation

◆ make_mech_prop_case_091()