anchor.diff:        (crdpsv1(j),   j=1,crnapxo),       &
         (crejv(j),     j=1,crnapxo),       &
         (crejfv(j),    j=1,crnapxo),       &
  -      (crmoidpsv(j), j=1,crnapxo),       &
  -      (cromoidpsv(j),j=1,crnapxo),       &
         (crnucm(j),    j=1,crnapxo),       &
         (crmtc(j),     j=1,crnapxo),       &
         (crnaa(j),     j=1,crnapxo),       &
  
positive.code:     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 the thick arrays
    if(ithick == 1) call synuthck
  
negative.code:   call alloc(crdpsv1,    npart_new,    zero,    "crdpsv1")
    call alloc(crejv,      npart_new,    zero,    "crejv")
    call alloc(crejfv,     npart_new,    zero,    "crejfv")
    call alloc(crnucm,     npart_new,    nucm0,   "crnucm")
    call alloc(crmtc,      npart_new,    one,     "crmtc")
    call alloc(crnaa,      npart_new,    aa0,     "crnaa")
  