anchor.diff:    (part_indiv_cr(j)     ,j=1,npart), &
     (part_impact_cr(j)    ,j=1,npart), &
   
  +  (cry_proc_cr(j)       ,j=1,npart), &
  +  (cry_proc_prev_cr(j)  ,j=1,npart), &
  +  (cry_proc_tmp_cr(j)   ,j=1,npart), &
  +
     (cn_impact_cr(j)      ,j=1,cdb_nColl_cr), &
     (cn_absorbed_cr(j)    ,j=1,cdb_nColl_cr), &
     (caverage_cr(j)       ,j=1,cdb_nColl_cr), &
  
positive.code:   real(kind=fPrec), allocatable, private, save :: yineff_cr(:)  ! (npart)
    real(kind=fPrec), allocatable, private, save :: xpineff_cr(:) ! (npart)
    real(kind=fPrec), allocatable, private, save :: ypineff_cr(:) ! (npart)
  
    integer, private, save :: n_tot_absorbed_cr = 0
    integer, private, save :: n_absorbed_cr     = 0
  
negative.code:   call alloc(part_indiv_cr,     npart, -c1m6,   "part_indiv_cr")
    call alloc(part_impact_cr,    npart, zero,    "part_impact_cr")
  
    call alloc(cry_proc_cr,       npart, -1,      "cry_proc_cr")
    call alloc(cry_proc_prev_cr,  npart, -1,      "cry_proc_prev_cr")
    call alloc(cry_proc_tmp_cr,   npart, -1,      "cry_proc_tmp_cr")
  
    call alloc(cn_impact_cr,     cdb_nColl_cr, 0,    "cn_impact")
    call alloc(cn_absorbed_cr,   cdb_nColl_cr, 0,    "cn_absorbed")
    call alloc(caverage_cr,      cdb_nColl_cr, zero, "caverage")
  