anchor.diff:    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)
     
  
positive.code: 70 FORMAT('Total computation time is: ',I2.2,' hrs, ',I2.2, ' min, ',I2.2, ' s, ',I4.4 ,' ms')  
  
  ! Close the debug log
    CLOSE(1)
  
    STOP
  
  
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
  