TF2A Subroutine

  1. astro_module.f90
  2. astro_module
  3. TF2A

public subroutine TF2A(s, ihour, imin, sec, rad, j)

Convert hours, minutes, seconds to radians.

Status: support routine.

Notes

  1. If the s argument is a string, only the leftmost character is used and no warning status is provided.

  2. The result is computed even if any of the range checks fail.

  3. Negative IHOUR, IMIN and/or SEC produce a warning status, but the absolute value is used in the conversion.

  4. If there are multiple errors, the status value reflects only the first, the smallest taking precedence.

History

  • IAU SOFA revision: 2013 December 2

Arguments

TypeIntentOptionalAttributesName
character(len=1), intent(in) :: s

sign: '-' = negative, otherwise positive
integer, intent(in) :: ihour

hours
integer, intent(in) :: imin

minutes
real(kind=wp), intent(in) :: sec

seconds
real(kind=wp), intent(out) :: rad

angle in radians
integer, intent(out) :: j
  • 0 = OK
  • 1 = IHOUR outside range 0-23
  • 2 = IMIN outside range 0-59
  • 3 = SEC outside range 0-59.999...

Contents

Source Code

TF2A

Source Code

    subroutine TF2A ( s, ihour, imin, sec, rad, j )

    implicit none

    character(len=1),intent(in) :: s !! sign:  '-' = negative, otherwise positive
    integer,intent(in) :: ihour !! hours
    integer,intent(in) :: imin !! minutes
    real(wp),intent(in) :: sec !! seconds
    real(wp),intent(out) :: rad !! angle in radians
    integer,intent(out) :: j !! status:
                             !! * 0 = OK
                             !! * 1 = IHOUR outside range 0-23
                             !! * 2 = IMIN outside range 0-59
                             !! * 3 = SEC outside range 0-59.999...

    real(wp) :: w

    !  Preset the status.
    j = 0

    !  Validate seconds, minutes, hours.
    if ( sec<0.0_wp .or. sec>=60.0_wp ) j=3
    if ( imin<0 .or. imin>59 ) j=2
    if ( ihour<0 .or. ihour>23 ) j=1

    !  Compute the angle.
    w = ( 60.0_wp * ( 60.0_wp * real( abs(ihour), wp ) + &
                                real( abs(imin), wp ) ) + &
                                      abs(sec) ) * ds2r

    !  Apply the sign.
    if ( s == '-' ) w = -w

    !  Return the result.
    rad = w

    end subroutine TF2A