GET_DIME_SISM.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	get_dime_sism (xipo, yipo, zipo, x1, y1, z1, x2, y2, z2, x3, y3, z3, nnod, xs, ys, zs, node_sism, nn_loc, mpi_id, ind, loc_n_n
	Computes local number of nodes where seismic source is imposed.

Function/Subroutine Documentation

get_dime_sism()

subroutine get_dime_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	node_sism,
		integer*4	nn_loc,
		integer*4	mpi_id,
		integer*4	ind,
			loc_n_n
	)

Computes local number of nodes where seismic source is imposed.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	Xipo	x-coordinate of Hypocenter
[in]	Yipo	y-coordinate of Hypocenter
[in]	Zipo	z-coordinate of Hypocenter
[in]	X1	vertex coordinate (x) of the triangular fault
[in]	Y1	vertex coordinate (y) of the triangular fault
[in]	Z1	vertex coordinate (z) of the triangular fault
[in]	X2	vertex coordinate (x) of the triangular fault
[in]	Y2	vertex coordinate (y) of the triangular fault
[in]	Z2	vertex coordinate (z) of the triangular fault
[in]	X3	vertex coordinate (x) of the triangular fault
[in]	Y3	vertex coordinate (y) of the triangular fault
[in]	Z3	vertex coordinate (z) of the triangular fault
[in]	nnod	number of local nodes (according to nodedomain.mpi)
[in]	xs	x-coordinate of spectral nodes
[in]	ys	y-coordinate of spectral nodes
[in]	zs	z-coordinate of spectral nodes
[in]	nn_loc	Total number of local nodes
[in]	loc_n_num	numeration of local nodes
[in]	mpi_id	id of MPI process
[in]	ind	index dummy
[out]	node_sism	local number of seismic source

Definition at line 47 of file GET_DIME_SISM.f90.

         
      implicit none
      
      integer*4 :: isn,i,nn_loc,mpi_id,ind
      integer*4 :: nnod,node_sism
      
      integer*4, dimension(nn_loc) :: loc_n_num
 
      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(4) :: x,y,z
      real*8, dimension(nn_loc) :: xs,ys,zs
 
      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
        endif
 
      enddo
 
 
      return