MPAS_Ocean_Julia

For sure another case for the "you doing it wrong" team at Nvidia. But the decision to add Julia/CUDA.jl here was not the best idea on a comparison with really interesting insights about Fortran vs. Julia. The CUDA part of the code is nearly a one-to-one copy of the MPI code. Especially one sentence[1] (there are some similar, comparable blocks) explains why the code not performs very well. [1] https://github.com/robertstrauss/MPAS_Ocean_Julia/blob/main/...

The authors seem to using the MPI package from Julia [1] so not something they did on their own. I didn't see any discussion on why MPI.jl is scales differently than Fortran-MPI. Something for the MPI.jl developers to investigate. [1] https://github.com/robertstrauss/MPAS_Ocean_Julia/blob/main/...