anchor.diff:  CONTAINS
   
   SUBROUTINE read_bonds(data, bond_info, bond_degeneracy)
  +
  +  INTEGER :: Ulog
  +  COMMON /Fio/ Ulog
   
     INTEGER :: RsCH3, RdCH2, CtCH, RSCH2sR, CdCHsR, Aldehyde,Cumulene
     INTEGER :: Ketene, CtCsR, RsCHsR2, CdCsR2, Ketone, RsCsR3
  
positive.code: SUBROUTINE Case_25(Total_char_freq, Total_harm_osc_freq, HR_params )
  
    USE heat_capacity_functions
     
  
    IMPLICIT NONE
  
  ! Global Variables
    REAL(8), DIMENSION(7) :: CV_temps
  
negative.code: 
    CALL dqed ( Case_16_hd, mequa, nvars, mcon, ind, bl, bu, x, fj, ldfj, &
      fnorm, igo, iopt, ropt, iwork, work )
    WRITE(Ulog, '(a)' ) ' '
    WRITE(Ulog, '(a)' ) ' CASE 16 Baby Yea!'
    WRITE(Ulog, '(a)' ) ' '
    WRITE(Ulog, '(a,i6)' ) '  Output flag from DQED, IGO = ', igo
    WRITE(Ulog, '(a)' ) ' '
    WRITE(Ulog, '(a)' ) '  Computed X:'
    WRITE(Ulog, '(a)' ) ' '
    WRITE(Ulog, '(g14.8)' ) x(1:nvars)
    WRITE(Ulog, '(a)' ) ' '
    WRITE(Ulog, '(a,g14.6)' ) '  L2 norm of the residual, FNORM = ', fnorm
    WRITE(Ulog, '(a)' ) ' '
    WRITE(Ulog, '(a)' ) 'Expected X:'
   
    WRITE(Ulog, '(a)' ) ' '
    WRITE(Ulog, '(a)' ) 'DQED_PRB2'
    WRITE(Ulog, '(a)' ) '  Normal end of execution.'
  
    WRITE(Ulog, '(a)' ) ' '
    CALL timestamp ( )
  
    results = x
  