Disv1dGeom.f90

Go to the documentation of this file.

module disv1dgeom
 
  use kindmodule, only: dp, i4b
  implicit none
  private
  public :: calcdist, line_unit_vector
 
contains
 
  !> @brief Calculate distance between two vertices
  !<
  function calcdist(vertices, ivert1, ivert2) result(dist)
    ! -- dummy
    real(dp), dimension(:, :), intent(in) :: vertices
    integer(I4B), intent(in) :: ivert1
    integer(I4B), intent(in) :: ivert2
    real(dp) :: dist, xdist, ydist
    ! -- local
    !
    ! -- calc distance
    xdist = abs(vertices(1, ivert1) - vertices(1, ivert2))
    ydist = abs(vertices(2, ivert1) - vertices(2, ivert2))
    dist = sqrt(xdist * xdist + ydist * ydist)
  end function calcdist
 
  !> @brief Calculate distance between two vertices
  !!
  !! Calculate the vector components (xcomp, ycomp, and zcomp)
  !! for a line defined by two points, (x0, y0, z0), (x1, y1, z1). Also
  !! return the magnitude of the original vector, vmag.
  !!
  !<
  subroutine line_unit_vector(x0, y0, x1, y1, &
                              xcomp, ycomp, vmag)
    real(dp), intent(in) :: x0
    real(dp), intent(in) :: y0
    real(dp), intent(in) :: x1
    real(dp), intent(in) :: y1
    real(dp), intent(out) :: xcomp
    real(dp), intent(out) :: ycomp
    real(dp) :: dx, dy, vmag
    !
    dx = x1 - x0
    dy = y1 - y0
    vmag = sqrt(dx**2 + dy**2)
    xcomp = dx / vmag
    ycomp = dy / vmag
  end subroutine line_unit_vector
 
end module disv1dgeom