I am using a third order Runge-Kutta for time integration. Do I apply boundary conditions at the end of each sub-iteration or at the end? Shuo

I have used boundary conditions at each sub-iteration. You would rather use it for sure.