ec2eq Subroutine

private pure subroutine ec2eq(pe, pq)

converts cartesian heliocentric j2000 ecliptic to equatorial

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in), dimension(6) :: pe

ecliptic
real(kind=wp), intent(out), dimension(6) :: pq

equatorial

Called by

proc~~ec2eq~~CalledByGraph proc~ec2eq standish_module::ec2eq proc~standish_rv_func standish_module::standish_ephemeris%standish_rv_func proc~standish_rv_func->proc~ec2eq proc~standish_module_test standish_module::standish_module_test proc~standish_module_test->proc~standish_rv_func

Source Code

    pure subroutine ec2eq (pe, pq)

    implicit none

    real(wp), dimension(6), intent (in) :: pe  !! ecliptic
    real(wp), dimension(6), intent (out) :: pq !! equatorial

    pq (1) = pe (1)                             ! xeq same as xec
    pq (2) = s_cobl * pe (2) - s_sobl * pe (3)  ! yeq
    pq (3) = s_sobl * pe (2) + s_cobl * pe (3)  ! zeq
    pq (4) = pe (4)                             ! vxeq same as vxec
    pq (5) = s_cobl * pe (5) - s_sobl * pe (6)  ! vyeq
    pq (6) = s_sobl * pe (5) + s_cobl * pe (6)  ! vzeq

    end subroutine ec2eq