Castro icon indicating copy to clipboard operation
Castro copied to clipboard

zero out I_q when castro.drive_initial_convection resets state

Open zingale opened this issue 1 year ago • 0 comments

When we do SDC with castro.drive_initial_convection=1, every time we hit castro.drive_initial_convection_reinit_period we reset the fluid state, but keep the velocities. This means that in the next timestep, the energy generation will change discontinuously, so using the lagged I_q from the previous step across this reset is not a good idea. We should reset I_q in this case.

zingale avatar Dec 01 '22 19:12 zingale