Benchmarking scripts for immersed boundaries and immersed boundary conditions

[glwagner]

We need to benchmark various flavors of ImmersedBoundaryGrid including

• GridFittedBoundary for immersed boundaries that are both function-based (computed on the fly) and array-based (precomputed)
• GridFittedBottom for the same (I think for this case, precomputation is always wise. That's not the case for GridFittedBoundary though. The benchmarks will hopefully show this.)
• Various immersed boundary conditions implemented in Support boundary conditions on immersed boundaries #2437

These should be implemented for both NonhydrostaticModel and HydrostaticFreeSurfaceModel.

[francispoulin]

If we are doing this, could we also try doing immersed coastal boundaries for ShallowWaterModel? It shouldn't be any harder and would be nice to make sure everything works there as well.

[glwagner]

I'm closing this issue because I'm judging that it's not of current, timely relevance to Oceananigans development. If you would like to make it a higher priority or if you think the issue was closed in error please feel free to re-open.