READ_SISM.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	read_sism (xipo, yipo, zipo, x1, y1, z1, x2, y2, z2, x3, y3, z3, nnod, xs, ys, zs, num_ns, sour_ns, i_sism, dist_sour_ns, nl_sism, max_num_ns, loc_n_num, nn_loc, pos_sour_nx, pos_sour_ny, pos_sour_nz)
	Generates seismic triangular faults.

Function/Subroutine Documentation

read_sism()

subroutine read_sism	(	real*8	xipo,
		real*8	yipo,
		real*8	zipo,
		real*8	x1,
		real*8	y1,
		real*8	z1,
		real*8	x2,
		real*8	y2,
		real*8	z2,
		real*8	x3,
		real*8	y3,
		real*8	z3,
		integer*4	nnod,
		real*8, dimension(nn_loc)	xs,
		real*8, dimension(nn_loc)	ys,
		real*8, dimension(nn_loc)	zs,
		integer*4	num_ns,
		integer*4, dimension(max_num_ns,nl_sism)	sour_ns,
		integer*4	i_sism,
		real*8, dimension(max_num_ns,nl_sism)	dist_sour_ns,
		integer*4	nl_sism,
		integer*4	max_num_ns,
		integer*4, dimension(nn_loc)	loc_n_num,
		integer*4	nn_loc,
		real*8, dimension(max_num_ns,nl_sism)	pos_sour_nx,
		real*8, dimension(max_num_ns,nl_sism)	pos_sour_ny,
			pos_sour_nz
	)

Generates seismic triangular faults.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	Xipo	Hipocenter x-coordinate
[in]	Yipo	Hipocenter y-coordinate
[in]	Zipo	Hipocenter z-coordinate
[in]	X1	x-coordinate of the triangular fault 1-node
[in]	Y1	y-coordinate of the triangular fault 1-node
[in]	Z1	z-coordinate of the triangular fault 1-node
[in]	X2	x-coordinate of the triangular fault 2-node
[in]	Y2	y-coordinate of the triangular fault 2-node
[in]	Z2	z-coordinate of the triangular fault 2-node
[in]	X3	x-coordinate of the triangular fault 3-node
[in]	Y3	y-coordinate of the triangular fault 3-node
[in]	Z3	z-coordinate of the triangular fault 3-node
[in]	xs	x-coord local spectral nodes
[in]	ys	y-coord local spectral nodes
[in]	zs	z-coord local spectral nodes
[in]	num_ns	number of seismic nodes
[in]	i_sism	column index for filling sour_ns
[in]	nl_sism	number of seismic loads
[in]	nn_loc	number of local nodes
[in]	nnod	number of local nodes
[in]	loc_n_num	local node numeration
[in]	max_num	max number of nodes for the specific triangular fault
[out]	sour_ns	info about seismic nodes
[out]	dist_sour_ns	distance of seismic nodes in sour_ns

Definition at line 48 of file READ_SISM.f90.

 
 
      implicit none
      
      integer*4 :: i_sism,isn,i,nn_loc
      integer*4 :: nnod,node_sism,num_ns,nl_sism,max_num_ns
      
      integer*4, dimension(nn_loc) :: loc_n_num(nn_loc)
 
      integer*4, dimension(max_num_ns,nl_sism) :: sour_ns
 
      real*8 :: xipo,yipo,zipo,x1,y1,z1,x2,y2,z2,x3,y3,z3,x4,y4,z4,ux,uy,uz,vx,vy,vz,tol
      real*8 :: p1x,p1y,p1z,p2x,p2y,p2z 
      real*8 :: epsilon, twopi, anglesum, costheta,m1,m2
 
      real*8, dimension(nn_loc) :: xs,ys,zs
      real*8, dimension(4) :: x,y,z
 
      real*8, dimension(max_num_ns,nl_sism) :: dist_sour_ns
      real*8, dimension(max_num_ns,nl_sism) :: pos_sour_nx,pos_sour_ny,pos_sour_nz
 
      node_sism = 0
 
      tol = 6.28
      epsilon = 1.d-8
      twopi = 6.283185307179586476925287
                    
      x(1) = x1
      y(1) = y1
      z(1) = z1
 
      x(2) = x2
      y(2) = y2
      z(2) = z2 
      
      x(3) = x3
      y(3) = y3
      z(3) = z3
 
      x(4) = x1
      y(4) = y1
      z(4) = z1
 
 
      do isn = 1,nnod
 
          anglesum = 0.0d0
          do i = 1,3
                   p1x = x(i) - xs(isn)
                  p1y = y(i) - ys(isn)
                  p1z = z(i) - zs(isn)
                  p2x = x(i+1) - xs(isn)
                  p2y = y(i+1) - ys(isn)
                  p2z = z(i+1) - zs(isn)
 
                  m1 = dsqrt(p1x*p1x + p1y*p1y + p1z*p1z)
                  m2 = dsqrt(p2x*p2x + p2y*p2y + p2z*p2z)
      
                  if ((m1*m2).le.epsilon) then
                        anglesum = twopi
                  else
                        costheta = (p1x*p2x + p1y*p2y + p1z*p2z)/(m1*m2)
                        anglesum = anglesum + dacos(costheta)
                  endif
          enddo
 
          if (anglesum.ge.tol) then
 
              node_sism = node_sism + 1
              sour_ns(node_sism,i_sism) = loc_n_num(isn) 
              dist_sour_ns(node_sism,i_sism) = dsqrt((xipo - xs(isn))**2 +(yipo - ys(isn))**2 +(zipo - zs(isn))**2)
              pos_sour_nx(node_sism,i_sism) = xs(isn)
              pos_sour_ny(node_sism,i_sism) = ys(isn)              
              pos_sour_nz(node_sism,i_sism) = zs(isn)
 
          endif
 
      enddo
 
     
      return