anchor.diff:      end if
     end do
   
  -  ! Recompute from loaded arrays
  -  oidpsv(1:napxo) = one/(one + dpsv(1:napxo))
  -  rvv(1:napxo)    = (ejv(1:napxo)*e0f)/(e0*ejfv(1:napxo))
  +  ! Recompute from loaded arrays (keep in sync with mod_particles)
  +  oidpsv(1:napxo)   = one/(one + dpsv(1:napxo))
  +  moidpsv(1:napxo)  = mtc(1:napxo)/(one + dpsv(1:napxo))
  +  omoidpsv(1:napxo) = ((one-mtc(1:napxo))*oidpsv(1:napxo))*c1e3
  +  rvv(1:napxo)      = (ejv(1:napxo)*e0f)/(e0*ejfv(1:napxo))
   
     ! Recompute the thick arrays
     if(ithick == 1) call synuthck
  
positive.code:       (dpsv1(j),   j=1,napxo),       &
        (ejv(j),     j=1,napxo),       &
        (ejfv(j),    j=1,napxo),       &
        (moidpsv(j), j=1,napxo),       &
        (omoidpsv(j),j=1,napxo),       &
        (nucm(j),    j=1,napxo),       &
        (mtc(j),     j=1,napxo),       &
        (naa(j),     j=1,napxo),       &
  
negative.code:   real(kind=fPrec), allocatable, private, save :: crdpsv1(:)    ! (npart)
    real(kind=fPrec), allocatable, private, save :: crejv(:)      ! (npart)
    real(kind=fPrec), allocatable, private, save :: crejfv(:)     ! (npart)
    real(kind=fPrec), allocatable, private, save :: craperv(:,:)  ! (npart,2)
    real(kind=fPrec), allocatable, private, save :: crnucm(:)     ! (npart)
    real(kind=fPrec), allocatable, private, save :: crmtc(:)      ! (npart)
  