From b7d25dc8303fb2819586f38415b0e38910e51314 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Tue, 15 Jun 2021 16:53:48 -0400 Subject: [PATCH] Fixed bug where the original heliocentric position was not properly saved when oblateness is turned on --- src/rmvs/rmvs_step.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/rmvs/rmvs_step.f90 b/src/rmvs/rmvs_step.f90 index 906b511c1..10de44f33 100644 --- 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