From 2dc254c372220bfc4cb8a566ae199f3bd2ac141a Mon Sep 17 00:00:00 2001 From: David A Minton Date: Thu, 8 Jul 2021 08:54:22 -0400 Subject: [PATCH] Fixed bad index in planetocentric to heliocentric conversion --- src/rmvs/rmvs_step.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/rmvs/rmvs_step.f90 b/src/rmvs/rmvs_step.f90 index c4535d884..6fe920863 100644 --- a/src/rmvs/rmvs_step.f90 +++ b/src/rmvs/rmvs_step.f90 @@ -592,8 +592,8 @@ subroutine rmvs_end_planetocentric(system) tp%xh(j, tpind(1:pl%nenc(i))) = tpenci%xh(j,1:pl%nenc(i)) + pl%inner(NTPHENC)%x(j, i) tp%vh(j, tpind(1:pl%nenc(i))) = tpenci%vh(j,1:pl%nenc(i)) + pl%inner(NTPHENC)%v(j, i) end do - tp%lperi(tpind(1:pl%nenc(i))) = tpenci%lperi(tpind(1:pl%nenc(i))) - tp%plperP(tpind(1:pl%nenc(i))) = tpenci%plperP(tpind(1:pl%nenc(i))) + tp%lperi(tpind(1:pl%nenc(i))) = tpenci%lperi(1:pl%nenc(i)) + tp%plperP(tpind(1:pl%nenc(i))) = tpenci%plperP(1:pl%nenc(i)) deallocate(pl%planetocentric(i)%tp) deallocate(cbenci%inner) do inner_index = 0, NTPHENC