anchor.diff:    nu_high = 300.0
   
     ! Read everything from the input
  -  READ (Uin,*) cp_data, data
  +  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
  
positive.code: ! This is the command that calls the main program.
    CALL Calculate_RRHO_HR_params(Uin, Uout)
    
  ! Close the input/output file units
    CLOSE(Uin)
    CLOSE(Uout)
  
  ! Ditto for the next five lines.  Cut them.
    CALL  date_and_time(date, time, zone, end_value)
    WRITE(Ulog,*) ''
  
negative.code:   Uout = 6
    Ulog = 1
  
  1234 CONTINUE
  ! Feel free to cut this.  I use it to determine how long the code runs.  
    CALL  date_and_time(date, time, zone, start_value)
    
  