anchor.diff:    c_length   = cdb_cLength(icoll)
     c_offset   = cdb_cOffset(icoll)
     c_tilt(:)  = cdb_cTilt(:,icoll)
  +  cRot       = cos_mb(c_rotation)
  +  sRot       = sin_mb(c_rotation)
     if(firstrun .and. dowrite_amplitude) then
       call coll_computeStats
     end if
  
positive.code:   real(kind=fPrec) nsig,c_length,jawLength,jawAperture,jawOffset,jawTilt(2),x_Dump,xpDump,y_Dump,   &
      ypDump,s_Dump,xmax,ymax,calc_aperture,zpj,xmax_pencil,ymax_pencil,xmax_nom,ymax_nom,            &
      nom_aperture,scale_bx,scale_by,c_tilt(2),c_offset,c_aperture,c_rotation,cry_bendangle,cry_tilt, &
      cry_tilt0
  
    call time_startClock(time_clockCOLL)
  
  
negative.code: 
    ! Assign nominal OR design beta functions for later
    if(do_nominal) then
      bx_dist = (cdb_cBx(icoll)*scale_bx)/scale_bx0
      by_dist = (cdb_cBy(icoll)*scale_by)/scale_by0
    else
      bx_dist = (tbetax(ie)*scale_bx)/scale_bx0
      by_dist = (tbetay(ie)*scale_by)/scale_by0
    end if
  
    ! Write beam ellipse at selected collimator
  