Getting scipy's solve_ivp function to execute intermediate task

Question

I am solving a fluid-structure interaction problem and I need to modify solve_ivp in order to execute a intermediate task after each time step. The algorithm should be somet