-
Notifications
You must be signed in to change notification settings - Fork 7
use TripolarGrid for Oceananigans sim #1409
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
base: main
Are you sure you want to change the base?
Conversation
c04aad4 to
9194823
Compare
|
@juliasloan25 how do I run the coupled simulation ? Is there a way to run it outside the caltech cluster (eg on a GPU system we have at MIT)? |
05c6433 to
845571c
Compare
d8b8872 to
7b17877
Compare
Sorry I just saw this - I get a lot of Github notifications so they get missed often. We have a docs page on running the coupled simulation. It shows how to run interactively or from the command line. You should be able to run it locally or on a different cluster, and it'll pick up the device to run on automatically (i.e. if a GPU is available it will run on GPU). Let me know if there are problems. |
022d2e7 to
c452c22
Compare
d04aa25 to
11addc7
Compare
11addc7 to
f5c917d
Compare
Purpose
Use a TripolarGrid instead of a capped lat/lon grid when constructing an Oceananigans simulation.
closes #1382
Note:
Oceananigans.Fields.interpolatedoesn't work for the TripolarGrid. We use this to remap Oceananigans fields from the OC grid to the ClimaCore field, so we can't switch to TripolarGrid until Oceananigans has a method that can regrid the TripolarGrid. This should be added in CliMA/ClimaOcean.jl#463 - requires Oceananigans v0.99.2.Note 2: Oceananigans is now using XESMF.jl to regrid for their TripolarGrid. This PR uses XESMF bilinear regridding to go from the Oceananigans grid to the ClimaCore boundary space. We continue to use ClimaCore.Remapping to go from ClimaCore to Oceananigans because it uses more information about the spectral element space than XESMF is able to. See this link for information about the limitations of using XESMF with spectral element spaces.
To-do
TripolarGridinstead ofLatitudeLongitudeGridInterfacer.get_remapper_to_ccthat retrieves the remapping object to get from a component model's grid to ClimaCore (this isnothingfor components that use ClimaCore, but for other components it's necessary to have access to the remapper object when remap is called)Additional notes
TripolarGridconstructs aOrthogonalSphericalShellGrid, which also uses lat/lon, so we can still use the same ClimaCore remapper to go from ClimaCore to Oceananigansset_from_extrinsic_vectors!