getvec Subroutine

public subroutine getvec(unitv)

THIS SUBROUTINE ALLOWS THE USER TO RETRIEVE THE LAST COMPUTED POSITION ON THE SKY AS A UNIT VECTOR.

 UNITV  = UNIT VECTOR TOWARD LAST COMPUTED POSITION ON THE
          SKY, IN THE COORDINATE SYSTEM USED FOR THAT
          POSITION (OUT)

Arguments

Type IntentOptional Attributes Name
real(kind=wp), dimension(3) :: unitv

Calls

proc~~getvec~~CallsGraph proc~getvec novas_module::getvec setvec setvec proc~getvec->setvec

Source Code

subroutine getvec ( unitv )

real(wp) :: unitv, p, pos, r, dsqrt
dimension unitv(3), p(3), pos(3)
save p

r = dsqrt ( p(1)**2 + p(2)**2 + p(3)**2 )

do j = 1, 3
    unitv(j) = p(j) / r
end do

return

entry setvec ( pos )
!
!     THIS ENTRY STORES THE LAST COMPUTED POSITION ON THE SKY.
!
!          POS    = VECTOR TOWARD LAST COMPUTED POSITION ON THE
!                   SKY, IN THE COORDINATE SYSTEM USED FOR THAT
!                   POSITION (IN)

do j = 1, 3
    p(j) = pos(j)
end do

return

end subroutine getvec