anchor.diff:    ROOR      = data(26)
     Peroxy    = data(27)
   
  +! cyclic species C-H stretch
  +  Rings     = data(28)
  +
   Bond_degeneracy =  9*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 + 6*RsCH2r + & 
                      3*RdCHr + 4*RsCHrsR + 2*CdCrsR + 2*OdCrsR + 3*RsCrsR2 + & 
  -                   4*Alcohol + 1*Ether + 5*ROOH + 3*ROOR + 3*Peroxy
  +                   4*Alcohol + 1*Ether + 5*ROOH + 3*ROOR + 3*Peroxy + &
  +                   2*Rings
   
   WRITE(*,*) 'degeneracy = ', Bond_degeneracy
     
  
positive.code:   INTEGER :: Ketene, CtCsR, RsCHsR2, CdCsR2, Ketone, RsCsR3
    INTEGER :: RsCH2r, RdCHr, RsCHrsR, CdCrsR, OdCrsR, RsCrsR2
    INTEGER :: Alcohol, Ether, ROOH, ROOR, Peroxy
  
  ! All the arrays corresponding to the characteristic frequencies for each bond
    REAL(8), ALLOCATABLE, DIMENSION(:) :: RsCH3_pred_freq, RdCH2_pred_freq 
  
negative.code:   CALL  Bond_avg_5x1(ROOH,  ROOH_typ_freq, ROOH_pred_freq)
    CALL  Bond_avg_3x1(ROOR,  ROOR_typ_freq, ROOR_pred_freq)
    CALL  Bond_avg_3x1(peroxy,  peroxy_typ_freq, peroxy_pred_freq)
    CALL  Bond_avg_2x1(Rings,  Rings_typ_freq, Rings_pred_freq)
  
  !------------------------------------------------------------------------------
  !  END calculating the bond degeneracy, sump cp bonds, and char. freqs.
  