Joseph Kump

> I suggest adding a test to Oceananigans for `Simulation` if we want to use it Although it wasn't breaking here, I also want to add test coverage for momentum...

Updated error: [error_0121.txt](https://github.com/user-attachments/files/18497443/error_0121.txt)

A reduced script reproducing the error can be found here: https://github.com/DJ4Earth/Enzymanigans.jl/blob/jlk9/reduce-for-active-arg-0122/dynamical_core/autodiff_double_gyre.jl

Oceananigans has been eliminated from the above script. It now replicates the AssertionError with Enzyme as its only dependency.

[jl_Y3qLCvQdWt.mlir.txt](https://github.com/user-attachments/files/19539853/jl_Y3qLCvQdWt.mlir.txt)

I tried `closure = nothing`, using `0` or an empty field instead of `0.0` in `no_slip_bc`, and omitting the boundary conditions altogether. Each of those still produced this error. I...

Realize this is on the back burner until after 4/15, but I think I isolated the cause of this error. I created a similar Oceananigans + Reactant script here: [https://github.com/DJ4Earth/Enzymanigans.jl/blob/jlk9/with-simulation/dynamical_core/reactant_attempt.jl](https://github.com/DJ4Earth/Enzymanigans.jl/blob/jlk9/with-simulation/dynamical_core/reactant_attempt.jl)...

This specific error output is no longer reproduced with the current versions of Reactant / Oceananigans, but we now get a "DimensionMismatch" error. I've made a new issue for it...

Enzyme is definitely affecting `mesh.Edges.dcEdge`. The code block you posted includes the `@show grad` statement in `gradient_normSq`, which I already found makes the change to `dcEdge` disappear. ``` using MOKA:...

Reduced the example in [https://github.com/jlk9/MPAS-Ocean.jl/tree/bug-reduce](https://github.com/jlk9/MPAS-Ocean.jl/tree/bug-reduce) further. Now it doesn't actually use the MOKA.jl package, instead using a minimal mesh struct as a fill-in (found in `bugMesh.jl`), and since the datafile...