anchor.diff:    READ (Uin,*,END=9998) cp_data, data
    ! WRITE (0,*) "cp_data is ",cp_data
    ! WRITE (0,'(A,28I3)') "data is ",data
  -  
  +
  +! The program may pause at the previous line until it is given more data.
  +! Now that we have some data to start processing,
  +! we can start the clock and open the log file.
  +
  +  ! Feel free to cut this.  I use it to determine how long the code runs.  
  +  CALL  date_and_time(date, time, zone, start_value)
  +
  +  ! Open the debug log file
  +  ! use /dev/null to silence
  +  open (UNIT=Ulog, FILE='frankie.log')
  +    
     ! call the subroutine which separates the data info and calculates
     ! how many RRHO frequencies it can determine from the structure
     call read_bonds(data, bond_info,  degeneracy)
  
positive.code:   INTEGER :: Uout
    INTEGER :: Ulog
    COMMON /Fio/ Ulog
    
    
    ! Fluffy stuff  FEEL FREE TO REMOVE
    CHARACTER(10) :: time
    CHARACTER(8) :: date
    CHARACTER(5) :: zone
    INTEGER, DIMENSION(8) :: start_value, end_value
  
  ! Should be getting these unit numbers from iso_fortran_env module
  ! See http://stackoverflow.com/questions/8508590/stderr-in-fortran90
  
negative.code:   INTEGER :: Uin  ! input file unit number
    INTEGER :: Uout ! output file unit number
  
    ! These variables are just used to calculate the running time. FEEL FREE TO REMOVE
    CHARACTER(10) :: time
    CHARACTER(8) :: date
    CHARACTER(5) :: zone
    INTEGER, DIMENSION(8) :: start_value, end_value
    
  
    ! These variables are set to the common block.
    ! They are required by the fitting routine within DQED
    COMMON N_vib, N_rot, nu_low, nu_high, nu_mid, cp_difference, CV_temps
  
  