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