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:   call alloc(part_indiv_cr,     npart, -c1m6,   "part_indiv_cr")
    call alloc(part_impact_cr,    npart, zero,    "part_impact_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")
  
negative.code:   call dealloc(part_indiv_cr,     "part_indiv_cr")
    call dealloc(part_impact_cr,    "part_impact_cr")
  
    call dealloc(cry_proc_cr,       "cry_proc_cr")
    call dealloc(cry_proc_prev_cr,  "cry_proc_prev_cr")
    call dealloc(cry_proc_tmp_cr,   "cry_proc_tmp_cr")
  
    call dealloc(cn_impact_cr,      "cn_impact")
    call dealloc(cn_absorbed_cr,    "cn_absorbed")
    call dealloc(caverage_cr,       "caverage")
  