GET_NEAREST_NODE_PGM.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	get_nearest_node_pgm (nloc, xs, ys, zs, xt, yt, zt, nt, dist_min
	Computes the nearest node with respet to (xt,yt,zt) starting from a given elevation.

Function/Subroutine Documentation

get_nearest_node_pgm()

subroutine get_nearest_node_pgm	(	integer*4	nloc,
		real*8, dimension(nloc)	xs,
		real*8, dimension(nloc)	ys,
		real*8, dimension(nloc)	zs,
		real*8	xt,
		real*8	yt,
		real*8	zt,
		integer*4	nt,
		real*8	dist_min
	)

Computes the nearest node with respet to (xt,yt,zt) starting from a given elevation.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	nloc	number of local nodes
[in]	xs	vertex x-coordinate of local nodes
[in]	ys	vertex y-coordinate of local nodes
[in]	zs	vertex z-coordinate of local nodes
[in]	xt	x-coor current vertex
[in]	yt	y-coor current vertex
[in]	zt	z-coor current vertex
[in]	elevation	z-value from wich we start the search
[out]	nt	index of the nearest node from (xt,yt,zt)
[out]	dist_min	minimal distance from (xt,yt,zt)

Definition at line 35 of file GET_NEAREST_NODE_PGM.f90.

      
      implicit none
      
      integer*4 :: i, nloc, nt
 
      real*8 :: xt,yt,zt
      real*8 :: elevation
      real*8 :: dist,dist_min
 
      real*8, dimension(nloc) :: xs, ys, zs
      
      dist_min = 1.0d30
      nt = 1
 
      do i = 1,nloc
 
         if (zs(i) .gt. elevation) then
            
            dist = dsqrt((xs(i) - xt)*(xs(i) - xt) + (ys(i) - yt)*(ys(i) - yt) + (zs(i) - zt)*(zs(i) - zt)) 
 
            if (dist .lt. dist_min) then
               dist_min = dist; nt = i
            endif
 
        endif
 
      enddo
      
      return

Referenced by find_monitor_position(), and read_system_position().

Here is the caller graph for this function: