anchor.diff:    output_file = 'rho_input'
     Uin = 12
     Uout = 22
  +  Ulog = 1
   
   ! 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 input file.
     OPEN (UNIT = Uin, FILE = input_file, STATUS = 'OLD', & 
  
positive.code:   CALL dqed ( Case_24_hd, mequa, nvars, mcon, ind, bl, bu, x, fj, ldfj, &
      fnorm, igo, iopt, ropt, iwork, work )
  
    WRITE ( *, '(a)' ) ' '
    WRITE ( *, '(a)' ) ' CASE 24'
    WRITE ( *, '(a)' ) ' '
    WRITE ( *, '(a,i6)' ) '  Output flag from DQED, IGO = ', igo
    WRITE ( *, '(a)' ) ' '
    WRITE ( *, '(a)' ) '  Computed X:'
    WRITE ( *, '(a)' ) ' '
    WRITE ( *, '(g14.8)' ) x(1:nvars)
    WRITE ( *, '(a)' ) ' '
    WRITE ( *, '(a,g14.6)' ) '  L2 norm of the residual, FNORM = ', fnorm
    WRITE ( *, '(a)' ) ' '
    WRITE ( *, '(a)' ) 'Expected X:'
   
    WRITE ( *, '(a)' ) ' '
    WRITE ( *, '(a)' ) 'DQED_PRB2'
    WRITE ( *, '(a)' ) '  Normal end of execution.'
  
    WRITE ( *, '(a)' ) ' '
    CALL timestamp ( )
  
    results = x
  
negative.code: SUBROUTINE assign_cases(N_rot, N_vib, Total_char_freq, Total_harm_osc_freq, HR_params )
  
  IMPLICIT NONE
    INTEGER :: Ulog
    COMMON /Fio/ Ulog
   INTEGER, INTENT(IN):: N_vib
   INTEGER, INTENT(IN) :: N_rot
  
  