anchor.diff:    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(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)
   
  
positive.code:   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
  
negative.code:   ALLOCATE(  Alcohol_pred_freq(  Alcohol  * 2 ) )
    ALLOCATE(  Ether_pred_freq(    Ether    * 1 ) )
    ALLOCATE(  ROOH_pred_freq(     ROOH     * 4 ) )
    ALLOCATE(  ROOR_pred_freq(     ROOR     * 2 ) )
    ALLOCATE(  peroxy_pred_freq(   peroxy   * 2 ) )
  
    ALLOCATE(  Rings_pred_freq(    Rings    * 2 ) )
  