get_r Function

private function get_r(me, et) result(r)

Type Bound

body_eph_interface

Arguments

Type IntentOptional Attributes Name
class(body_eph_interface), intent(inout) :: me
real(kind=wp), intent(in) :: et

ephemeris time (sec)

Return Value real(kind=wp), dimension(3)

position vector


Calls

proc~~get_r~~CallsGraph proc~get_r body_eph_interface%get_r db1val db1val proc~get_r->db1val

Called by

proc~~get_r~~CalledByGraph proc~get_r body_eph_interface%get_r proc~get_r_splined jpl_ephemeris_splined%get_r_splined proc~get_r_splined->proc~get_r

Source Code

    function get_r(me,et) result(r)

        class(body_eph_interface),intent(inout) :: me
        real(wp),intent(in) :: et !! ephemeris time (sec)
        real(wp),dimension(3) :: r !! position vector

        integer :: iflag, i

        do i = 1, 3
            call db1val(et, 0, me%eph%tx(:,i), me%eph%nx, kx, me%eph%bcoef(:,i), r(i), iflag, &
                        me%inbvx, me%w0, extrap=.false.)
            if (iflag /= 0) then
                write(*,*) 'iflag = ', iflag
                error stop 'error calling get_r'
            end if
        end do

    end function get_r