Fixed bad index in planetocentric to heliocentric conversion

MintonGroup · Jul 8, 2021 · 2dc254c · 2dc254c
1 parent 7d42193
commit 2dc254c
Showing 1 changed file with 2 additions and 2 deletions.
diff --git 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