anchor.diff:    call alloc(crdpsv1,    npart_new,    zero,    "crdpsv1")
     call alloc(crejv,      npart_new,    zero,    "crejv")
     call alloc(crejfv,     npart_new,    zero,    "crejfv")
  -  call alloc(crmoidpsv,  npart_new,    one,     "crmoidpsv")
  -  call alloc(cromoidpsv, npart_new,    zero,    "cromoidpsv")
     call alloc(crnucm,     npart_new,    nucm0,   "crnucm")
     call alloc(crmtc,      npart_new,    one,     "crmtc")
     call alloc(crnaa,      npart_new,    aa0,     "crnaa")
  
positive.code:       (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),       &
  
negative.code:       (dpsv1(j),   j=1,napxo),       &
        (ejv(j),     j=1,napxo),       &
        (ejfv(j),    j=1,napxo),       &
        (nucm(j),    j=1,napxo),       &
        (mtc(j),     j=1,napxo),       &
        (naa(j),     j=1,napxo),       &
  