lanl / phoebus

Phifty One Ergs Blows Up A Star
BSD 3-Clause "New" or "Revised" License
32 stars 0 forks source link

Pre-normalized torus B field can trigger MHD floors #149

Open brryan opened 1 year ago

brryan commented 1 year ago

I just ran into this weird bug -- the q we calculate is used to calculate a fiducial B field for the torus which is then normalized. But that normalization occurs AFTER an ApplyFloors call. So the nonsense B field values can trigger MHD floors that mess up the fluid variables before the B field is normalized.