anchor.diff:    bond_info(23) = ROOR 
     bond_info(24) = Peroxy
     bond_info(25) = Rings
  -
  -  ! RMG will probably double-count the number of ROOR bonds.
  -  ! Assuming that is the case, divide by 2:
  -  IF (ROOR/2==0) THEN
  -     !it is an odd number, so do nothing
  -  ELSE
  -     bond_info(23) = bond_info(23) / 2
  -  ENDIF
   
   
   END SUBROUTINE read_bonds
  
positive.code:   CALL  Bond_avg_2x1(Alcohol,  Alcohol_typ_freq, Alcohol_pred_freq)
    CALL  Bond_avg_1x1(Ether,  Ether_typ_freq, Ether_pred_freq)
    CALL  Bond_avg_4x1(ROOH,  ROOH_typ_freq, ROOH_pred_freq)
    CALL  Bond_avg_3x1(ROOR,  ROOR_typ_freq, ROOR_pred_freq)
    CALL  Bond_avg_2x1(peroxy,  peroxy_typ_freq, peroxy_pred_freq)
    CALL  Bond_avg_2x1(Rings,  Rings_typ_freq, Rings_pred_freq)
  
  
negative.code: ! cyclic species C-H stretch
    Rings     = data(28)
  
    ! RMG will probably double-count the number of ROOR bonds.
    ! Assuming that is the case, divide by 2:
    IF (ROOR/2==0) THEN
       !it is an odd number, so do nothing
    ELSE
       ROOR = ROOR / 2
    ENDIF
  
  Bond_degeneracy =  8*RsCH3 + 5*RdCH2 + 3*CtCH + 7*RSCH2sR + 5*CdCHsR + &
                     5*Aldehyde + 3*Cumulene + 3*Ketene + 2*CtCsR + & 
                     6*RsCHsR2 + 4*CdCsR2 + 4*Ketone + 5*RsCsR3 + 5*RsCH2r + & 
                     3*RdCHr + 4*RsCHrsR + 2*CdCrsR + 2*OdCrsR + 3*RsCrsR2 + & 
                     2*Alcohol + 1*Ether + 4*ROOH + 2*ROOR + 2*Peroxy + &
                     2*Rings
  
  WRITE(*,*) 'degeneracy = ', Bond_degeneracy
  