MAKE_MECH_PROP_CASE_001.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	make_mech_prop_case_001 (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_001()

subroutine make_mech_prop_case_001	(	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_001.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
 
     vs  = 300.0 + 19.0 * sqrt(depth)        !VS: S velocity in m/s
     vp  = 1450.0 + 1.2 * depth              !VP: P velocity in m/s 
     rho = 2140.0 + 0.125 * depth            !RHO: MASS DENSITY in kg/m^3
     lambda = rho * (vp**2 - 2*vs**2)
     mu = rho * vs**2
     gamma = 6.2832e-02
   
                                                                                     

Referenced by make_eltensor_for_cases().

Here is the caller graph for this function: