Transform a position (and optionally velocity) vector from LVLH to IJK.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=wp), | intent(in), | dimension(3) | :: | rt_ijk |
Target IJK absolute position vector [km] |
|
real(kind=wp), | intent(in), | dimension(3) | :: | vt_ijk |
Target IJK absolute position vector [km] |
|
real(kind=wp), | intent(in), | dimension(3) | :: | dr_lvlh |
Chaser LVLH position vector relative to target [km] |
|
real(kind=wp), | intent(in), | dimension(3) | :: | dv_lvlh |
Chaser LVLH position vector relative to target [km] |
|
real(kind=wp), | intent(out), | dimension(3) | :: | r_ijk |
Chaser IJK absolute position vector [km] |
|
real(kind=wp), | intent(out), | optional, | dimension(3) | :: | v_ijk |
Chaser IJK absolute position vector [km] |
subroutine from_lvlh_to_ijk_rv(rt_ijk,vt_ijk,dr_lvlh,dv_lvlh,r_ijk,v_ijk) implicit none real(wp),dimension(3),intent(in) :: rt_ijk !! Target IJK absolute position vector [km] real(wp),dimension(3),intent(in) :: vt_ijk !! Target IJK absolute position vector [km] real(wp),dimension(3),intent(in) :: dr_lvlh !! Chaser LVLH position vector relative to target [km] real(wp),dimension(3),intent(in) :: dv_lvlh !! Chaser LVLH position vector relative to target [km] real(wp),dimension(3),intent(out) :: r_ijk !! Chaser IJK absolute position vector [km] real(wp),dimension(3),intent(out),optional :: v_ijk !! Chaser IJK absolute position vector [km] real(wp),dimension(3,3) :: c real(wp),dimension(3,3) :: cdot if (present(v_ijk)) then call from_lvlh_to_ijk(rt_ijk,vt_ijk,c=c,cdot=cdot) !chaser = target + delta: r_ijk = rt_ijk + matmul( c, dr_lvlh ) v_ijk = vt_ijk + matmul( cdot, dr_lvlh ) + matmul( c, dv_lvlh ) else call from_lvlh_to_ijk(rt_ijk,vt_ijk,c=c) r_ijk = rt_ijk + matmul( c, dr_lvlh ) end if end subroutine from_lvlh_to_ijk_rv