anchor.diff:  SUBROUTINE Case_13(Total_char_freq, Total_harm_osc_freq, HR_params )
    
     USE heat_capacity_functions
  -   
   
     IMPLICIT NONE
  +  INTEGER :: Ulog
  +  COMMON /Fio/ Ulog
     
   ! Global Variables
     REAL(8), DIMENSION(7) :: CV_temps
  
positive.code:                               typical_frequencies, predicted_frequencies)
  
    IMPLICIT NONE
    INTEGER, INTENT(IN) :: Number_of_bonds
    REAL(8), DIMENSION(:), INTENT(IN) :: typical_frequencies
    REAL(8),DIMENSION(5*Number_of_bonds), INTENT(OUT) :: predicted_frequencies
    INTEGER                                      :: i,j, k
  
    IF (Number_of_bonds == 0 ) THEN
  !     WRITE (*,*) 'N = 0'
       ! predicted_frequencies has DIMENSION(0) so we can't assign anything to it
    
   ELSEIF (Number_of_bonds == 1) THEN
  
negative.code:    
  
    IMPLICIT NONE
    INTEGER :: Ulog
    COMMON /Fio/ Ulog
    
  ! Global Variables
    REAL(8), DIMENSION(7) :: CV_temps
  