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:   DEALLOCATE(Total_predicted_freq )
    DEALLOCATE( HR_params )
  
    RETURN
  9998 CONTINUE
    WRITE(0,*) 'End of input stream.'
    ! Close the input/output file units
    CLOSE(Uin)
    CLOSE(Uout)
    CLOSE(Ulog)
    STOP
  
  END SUBROUTINE Calculate_RRHO_HR_params
  
  
  