eq2ec Subroutine

    pure subroutine eq2ec (pq, pe)

    implicit none

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

    pe (1) = pq (1)                              ! xec same as xeq
    pe (2) = s_cobl * pq (2) + s_sobl * pq (3)   ! yec
    pe (3) = - s_sobl * pq (2) + s_cobl * pq (3) ! zec
    pe (4) = pq (4)                              ! vxec same as vxeq
    pe (5) = s_cobl * pq (5) + s_sobl * pq (6)   ! vyec
    pe (6) = - s_sobl * pq (5) + s_cobl * pq (6) ! vzec

    end subroutine eq2ec