Fixed bug where the original heliocentric position was not properly s…

…aved when oblateness is turned on
MintonGroup · Jun 15, 2021 · b7d25dc · b7d25dc
1 parent 63ad803
commit b7d25dc
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/rmvs/rmvs_step.f90 b/src/rmvs/rmvs_step.f90
@@ -222,7 +222,7 @@ subroutine rmvs_interp_in(pl, cb, dt, outer_index, config)
          xtmp(:, :) = pl%inner(0)%x(:, :)
          vtmp(:, :) = pl%inner(0)%v(:, :)
          if (config%loblatecb) then
-            allocate(xh_original,mold=pl%xh)
+            allocate(xh_original,source=pl%xh)
             pl%xh(:, :) = xtmp(:, :) ! Temporarily replace heliocentric position with inner substep values to calculate the oblateness terms
             call pl%obl_acc(cb)
             pl%inner(0)%aobl(:, :) = pl%aobl(:, :) ! Save the oblateness acceleration on the planet for this substep