Tomás Chor

Thanks for checking into this. I haven't tested it myself, but it seems that indeed that `@warn` line is the cause of the error here, even though I never reach...

Also I see there's been some discussion about possibly moving this package to SpecialFunctions.jl in https://github.com/JuliaMath/SpecialFunctions.jl/pull/371. In my opinion moving `lambertw` there would result in more eyes on it for...

Unfortunately this example ```julia using Oceananigans grid = RectilinearGrid(size = (1, 1, 8), extent = (1,1,1)); model = NonhydrostaticModel(; grid, closure = ScalarDiffusivity(ν=1e-2)) set!(model, u=(x, y, z,) -> z) simulation...

> > Still pretty easy, but more work/time that I have right now. > > Reading the code this is just one more line. Can you be more specific about...

> Let me know what else is missing or if we just need a third review to merge. Tests still aren't passing.

> Is this still an issue? Apparently not anymore. I just tested @navidcy 's MWE are everything seems to work fine (on a Quadro GPU and CUDA 3.12 on `main`):...

> still an issue? Yes. I just tested this on `main`

It's common for this error to pop up when there's a file handle already open (maybe if there was an error the previous time you ran the code). @loganpknudsen if...

I fixed the tests and also started testing that `BackgroundField`s get written correctly. @glwagner If you're okay with it, I'm gonna merge this.

Ah, I see. Sounds like it wouldn't be trivial to add that support. I guess a workaround to avoid partially-averaged results when picking up would be to set the `Checkpointer`...