need test case for rotating core-collapse progenitor?

The attached figure shows a (not-so-great) resolution test for a Z=0.001, 40Msun star initially rotating at ~0.6 critical which experiences chemically homogeneous evolution. The figure only shows the profile of specific angular momentum, density, and velocity of the inner 5Msun and the numbers in the legend are the number of zones at this particular snapshot.

This is a rather convoluted setup (damping spurious velocities in the layers not in sonic contact in the core and trying to add resolution for a bunch of gradients), but very similar to the setup [here](https://github.com/mathren/CHE_net_comparison) developed with input from @aurimontem . The point that worries me is the "spiky" infall velocity profile (more visible in the blue, but it's coming up in the orange and hidden a bit by the dashed linestyle).

I have been noticing this in many rotating models (even with less extreme rotation rates), but I just double checked that the test cases that go to core-collapse in r24.03.1 instead exhibit smooth velocity profiles, however they are also non-rotating.

I suspect this has to do with the rotation playing with the lagrangian hydrodynamics.
I would like to know if any of the developers has suggestions on how to fix this (@Debraheem maybe?), and if/once fixed, if there is interest in molding this into a test case for rotating core-collapse progenitors:


No MESA file or test entry point is named. Start by reviewing the CHE_net_comparison setup and reproducing the Z=0.001, 40Msun, ~0.6-critical rotating model, then compare its infall-velocity profile with the non-rotating r24.03.1 cases. Done means identifying the cause of the spikes and, if appropriate, adding a reproducible rotating core-collapse test case.

Written by the indexing model from the issue text.