anchor.diff:    napx   = crnapx
     e0     = cre0
     e0f    = sqrt(e0**2-nucm0**2)
  -  beta0 = crbeta0
  +  beta0  = crbeta0
     brho   = crbrho
     nucmda = crnucmda
   
  
positive.code: 
      write(crlog,"(a)") "CR_CHECK>  * Tracking variables"
      flush(crlog)
      read(cr_pntUnit(nPoint),iostat=ioStat) crnumlcr,crnuml,crsixrecs,crbinrec,cril,cr_time,crnapxo, 
        crnapx,cre0,crbeta0,crbrho,crnucmda
      if(ioStat /= 0) cycle
  
      write(crlog,"(a)") "CR_CHECK>  * Particle arrays"
      flush(crlog)
      read(cr_pntUnit(nPoint),iostat=ioStat) 
        (crbinrecs(j), j=1,(crnapxo+1)/2), 
        (crnumxv(j),   j=1,crnapxo),       
        (crnnumxv(j),  j=1,crnapxo),       
        (crpartID(j),  j=1,crnapxo),       
        (crparentID(j),j=1,crnapxo),       
        (crpstop(j),   j=1,crnapxo),       
        (crxv1(j),     j=1,crnapxo),       
        (cryv1(j),     j=1,crnapxo),       
        (crxv2(j),     j=1,crnapxo),       
        (cryv2(j),     j=1,crnapxo),       
        (crsigmv(j),   j=1,crnapxo),       
        (crdpsv(j),    j=1,crnapxo),       
        (crdpsv1(j),   j=1,crnapxo),       
        (crejv(j),     j=1,crnapxo),       
        (crejfv(j),    j=1,crnapxo),       
        (crnucm(j),    j=1,crnapxo),       
        (crmtc(j),     j=1,crnapxo),       
        (crnaa(j),     j=1,crnapxo),       
        (crnzz(j),     j=1,crnapxo),       
        (crnqq(j),     j=1,crnapxo),       
        (crpdgid(j),   j=1,crnapxo),       
        (craperv(j,1), j=1,crnapxo),       
        (craperv(j,2), j=1,crnapxo),       
        (crllostp(j),  j=1,crnapxo)
      if(ioStat /= 0) cycle
  
  
negative.code:       write(crlog,"(a)") "CR_POINT>  * Particle arrays"
        flush(crlog)
      end if
      write(cr_pntUnit(nPoint),err=100) &
        (binrecs(j), j=1,(napxo+1)/2),  &
        (numxv(j),   j=1,napxo),        &
        (nnumxv(j),  j=1,napxo),        &
        (partID(j),  j=1,napxo),        &
        (parentID(j),j=1,napxo),        &
        (pstop(j),   j=1,napxo),        &
        (xv1(j),     j=1,napxo),        &
        (yv1(j),     j=1,napxo),        &
        (xv2(j),     j=1,napxo),        &
        (yv2(j),     j=1,napxo),        &
        (sigmv(j),   j=1,napxo),        &
        (dpsv(j),    j=1,napxo),        &
        (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),        &
        (nzz(j),     j=1,napxo),        &
        (nqq(j),     j=1,napxo),        &
        (pdgid(j),   j=1,napxo),        &
        (aperv(j,1), j=1,napxo),        &
        (aperv(j,2), j=1,napxo),        &
        (llostp(j),  j=1,napxo)
      flush(cr_pntUnit(nPoint))
  
  