do while (time_step_rem .gt. 1.e-15)
& time_step_rem )
time_step_rem = time_step_rem - time_step_loc
_RL time_step_full, time_step_rem
time_step_rem = time_step_loc