anchor.diff:      goto 30
     end if
     if(nSplit == 0) goto 10 ! Skip empty lines
  -  
  +
     ! Look up the target collimator or family
     iFam  = -1
     iColl = -1
  
positive.code:     if(isFam) then
        do i=1,cdb_nColl
          if(cdb_cFamily(iFam) == iFam) then
            cdb_cSides(i) = iTemp
          end if
        end do
        write(lout,"(a,i0,a)") "COLLDB> Collimator family '"//trim(lnSplit(2))//"' is set as onesided (",iTemp,")"
      else
        cdb_cSides(iColl) = iTemp
        write(lout,"(a,i0,a)") "COLLDB> Collimator '"//trim(lnSplit(2))//"' is set as onesided (",iTemp,")"
      end if
  
    case default
      write(lerr,"(a)") "COLLDB> ERROR Unknown keyword '"//trim(lnSplit(1))//"' in SETTINGS section"
  
negative.code: 
  end subroutine cdb_writeDB_HDF5
  #endif
  
  ! ================================================================================================ !
  !  V.K. Berglyd Olsen, BE-ABP-HSS
  !  Created: 2019-09-02
  !  Updated: 2019-09-02
  !  Rotate a collimater by rAngle radians
  !  Note: does not check that zero <= rAngle <= twopi
  ! ================================================================================================ !
  subroutine cdb_rotateCollimator(collID, rAngle)
  
    use numerical_constants
  
    integer,          intent(in) :: collID
    real(kind=fPrec), intent(in) :: rAngle
  
    cdb_cRotation(collID) = cdb_cRotation(collID) + rAngle
    if(cdb_cRotation(collID) > twopi) then
      cdb_cRotation(collID) = cdb_cRotation(collID) - twopi
    else if(cdb_cRotation(collID) < zero) then
      cdb_cRotation(collID) = cdb_cRotation(collID) + twopi
    end if
  
  end subroutine cdb_rotateCollimator
  
  ! ================================================================================================ !
  !  V.K. Berglyd Olsen, BE-ABP-HSS
  