anchor.diff:  70 FORMAT('Total computation time is: ',I2.2,' hrs, ',I2.2, ' min, ',I2.2, ' s, ',I4.4 ,' ms')  
   
   ! Close the debug log
  -  CLOSE(1)
  +  CLOSE(Ulog)
  +  
  +  WRITE(Uout,*) 'FRANKIE_COMPLETED_ONE_ITERATION'
  +  FLUSH(Uout)
  +  GOTO 1234
   
     STOP
   
  
positive.code:   Uout = 6
    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)
    
  
negative.code:   nu_high = 300.0
  
    ! Read everything from the input
    READ (Uin,*,END=9998) cp_data, data
   ! WRITE (0,*) "cp_data is ",cp_data
   ! WRITE (0,'(A,28I3)') "data is ",data
    
    ! call the subroutine which separates the data info and calculates
    ! how many RRHO frequencies it can determine from the structure
  