Skip to content

New Mesh: EC15to60kmL64E3SMv2r01: Draft of new EC mesh with refinement around Greenland#522

Closed
xylar wants to merge 3 commits intoMPAS-Dev:ocean/developfrom
xylar:refined_NA_arctic_test
Closed

New Mesh: EC15to60kmL64E3SMv2r01: Draft of new EC mesh with refinement around Greenland#522
xylar wants to merge 3 commits intoMPAS-Dev:ocean/developfrom
xylar:refined_NA_arctic_test

Conversation

@xylar
Copy link
Collaborator

@xylar xylar commented Apr 17, 2020

Mesh long name: EC15to60kmL64E3SMv2r01

Adds a new EC15to60 test case with refinement around Greenland

This is just do document a mesh @maltrud and I are playing around with that might eventually inform a future EC, low res mesh for v2.

@xylar
Copy link
Collaborator Author

xylar commented Apr 17, 2020

Plots of resolution

cellWidthGlobal

arctic

south_atlantic

Comment on lines 12 to 15
Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Probably need to be adjusted

Comment on lines +3 to +10
Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Probably need to be adjusted

Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@vanroekel, this would also be where to add other namelist changes. Hit the "Insert a suggestion" icon and edit the code as you would like.

@xylar
Copy link
Collaborator Author

xylar commented Apr 17, 2020

The new EC grid can be found here:

/lustre/scratch4/turquoise/xylar/test_EC15to60/ocean/global_ocean/EC15to60/

I'm still spinning it up (and probably doing so very inefficiently) but this should be a good starting point for you to play around with the new mesh.

For now (for comparison with the EC60to30 mesh), we're running with 60, not 64 vertical levels:

nVertLevels = 60 ;
nCells = 220321 ;
nEdges = 671300 ;
nVertices = 450694 ;

@xylar xylar changed the title Proposed EC15to60kmL64E2r01 mesh Proposed EC15to60kmL60E2r01 mesh Apr 17, 2020
@vanroekel
Copy link
Contributor

@xylar or @mark-petersen before this develops too much further, can we socialize this mesh with water cycle as well, especially if there are discussions within cryosphere of this potentially being the default v2 low res mesh? I don't think anyone outside cryosphere has been a part of this discussion for developing this mesh. Looking at this, I like the cell count, but am worried the equatorial region is too coarse. I'm concerned that the transition zone will harm tropical instability wave fidelity and potentially degrade ENSO fidelity as well.

@maltrud
Copy link
Contributor

maltrud commented Apr 17, 2020

@vanroekel i've been fully involved with this, even safe to say i've been driving it. i want to emphasize that it is not intended to be the default v2 mesh, but if it looks good it could be a prototype. my main purpose is to see if a ec60to30-like grid with ~15km resolution in the subpolar gyre gives us a more realistic AMOC. i didn't notice if the latest has the narrower equatorial band--if it does, that could be interesting to test, as well.

@xylar xylar requested review from lconlon and vanroekel April 18, 2020 06:34
@xylar
Copy link
Collaborator Author

xylar commented Apr 18, 2020

@vanroekel, I can't find Chris Golaz's GitHub username. Feel free to mention him or add him as a reviewer.

@xylar xylar requested a review from milenaveneziani April 18, 2020 06:40
@xylar xylar changed the title Proposed EC15to60kmL60E2r01 mesh Early potential draft of a new EC15to60kmL60E2r01 mesh Apr 18, 2020
@xylar xylar added the For discussion PRs and Issues that are open for discussion and feedback label Apr 18, 2020
@xylar
Copy link
Collaborator Author

xylar commented Apr 18, 2020

@mark-petersen's argument for tighter resolution around the equator was based on:

  1. a discussion with FESOM developers at AWI that suggested this and
  2. The following figure suggesting that this should be sufficient to capture the major equatorial current:
    image

However, I'm sure discussion is welcome, and would also help us with revisions to our SO refined mesh (#518)

@xylar xylar force-pushed the refined_NA_arctic_test branch 2 times, most recently from f6ae664 to d0627a2 Compare April 18, 2020 12:12
@vanroekel
Copy link
Contributor

@xylar perhaps I'm not clear from the figures about the equatorial region. Is the 30 km region 7.5degrees or 10 or something else? The figure you post suggests we should use 10, which I think is perfect, but one of your figures, I think on the SO mesh one seem to show a transition at 7.5. Could you let me know what your cutoff is? Thanks

@vanroekel
Copy link
Contributor

one other question. Would you like me to suggest any parameter settings for this mesh? I saw a few spots saying values should be changed. I'm happy to help. One thing I want to make sure of is that we use the config_eddy_resolution_taper = 'ramp' and not constant.

Comment on lines 32 to 35
Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@vanroekel, the edge of the 30-km region around the equator is at +/- 7.5 degrees (latPosEq). This could be increased to 10.0. I would probably also increase latWidthEq to 4.0.

Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@mark-petersen had suggested that the counter current might not need high resolution, which is why we had halved the high-res region from 30 degrees total to 15 degrees total. It sounds like @vanroekel thinks it needs high res.

Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I updated this as stated above.

@xylar
Copy link
Collaborator Author

xylar commented Apr 19, 2020

one other question. Would you like me to suggest any parameter settings for this mesh? I saw a few spots saying values should be changed. I'm happy to help. One thing I want to make sure of is that we use the config_eddy_resolution_taper = 'ramp' and not constant.

It doesn't hurt to change these parameters here for the spin-up but it's probably more important for the E3SM pull request when that happens. I'll point out where these changes would go and you can use the GitHub suggestion feature to essentially edit the code here if you want.

@xylar xylar force-pushed the refined_NA_arctic_test branch 2 times, most recently from 6686056 to 8ff2fa2 Compare April 19, 2020 11:03
@xylar xylar force-pushed the refined_NA_arctic_test branch from edf02aa to d624acc Compare April 19, 2020 16:38
@xylar xylar changed the title Early potential draft of a new EC15to60kmL60E2r01 mesh EC15to60kmL60E2r01: Early potential draft of a new EC mesh with refinement around Greenland Apr 20, 2020
@lconlon
Copy link

lconlon commented Apr 20, 2020

I think this version looks pretty good, subpolar gyre-wise. What is the max GM that you're using?

@xylar
Copy link
Collaborator Author

xylar commented Apr 20, 2020

What is the max GM that you're using?

I don't think we've got far enough to decide that yet, and would certainly welcome input. @maltrud, @milenaveneziani, etc. what is the default?

@vanroekel
Copy link
Contributor

I'd imagine we would use a value of 600, which is the default for other EC6030 Gcases, but it is something we could and maybe should test.

@xylar xylar marked this pull request as ready for review April 20, 2020 16:46
@xylar xylar changed the title EC15to60kmL60E2r01: Early potential draft of a new EC mesh with refinement around Greenland EC15to60kmE2r01: Early potential draft of a new EC mesh with refinement around Greenland Apr 20, 2020
@lconlon
Copy link

lconlon commented Apr 20, 2020

I'd imagine we would use a value of 600, which is the default for other EC6030 Gcases, but it is something we could and maybe should test.

Sounds good. Happy to help with any testing if needed.

@milenaveneziani
Copy link
Contributor

We used a max GM of 600 for the ARRM meshes, but in those simulations we don't have 3d GM.

@pwolfram pwolfram added the mesh label Apr 20, 2020
@xylar xylar force-pushed the refined_NA_arctic_test branch from d624acc to 13e3190 Compare April 22, 2020 15:49
@xylar
Copy link
Collaborator Author

xylar commented May 21, 2020

@maltrud, is testing progressing on this mesh at all? Or are you still stuck? No rush, since you were the one that requested this one but figured we could put in some results here if we wanted.

At some point, we probably want to merge this just to document this version of the mesh, even if it isn't our final E3SM v2 EC mesh.

@vanroekel
Copy link
Contributor

@maltrud could you summarize your testing on this mesh? I seem to remember you saying we won't want this for v2, but don't remember the reason. As @xylar mentions in a comment above we should merge this to document, but it would be great to have a summary of what we might want to change for an r02 version of this mesh.

@xylar
Copy link
Collaborator Author

xylar commented Aug 26, 2020

It sounds like this hadn't made it far enough in testing that we can use it for E3SM v2. Nevertheless, I would like to have this rebased and merged for documentation purposes.

@xylar xylar force-pushed the refined_NA_arctic_test branch from 13e3190 to 7777730 Compare August 26, 2020 16:42
@xylar
Copy link
Collaborator Author

xylar commented Aug 26, 2020

Actually, I don't want to merge until #563 and #577

@xylar xylar changed the title EC15to60kmE2r01: Early potential draft of a new EC mesh with refinement around Greenland New Mesh: EC15to60kmL64E3SMv2r01: Draft of new EC mesh with refinement around Greenland Aug 26, 2020
@xylar xylar mentioned this pull request Nov 3, 2020
mark-petersen added a commit that referenced this pull request Nov 5, 2020
New Mesh: EC15to60E2r4 #740

Long Name: EC15to60kmL60E3SMv2r04

MPAS Eddy Closure mesh for E3SM version2 with enhanced resolution around
the equator (30 km), South pole (35 km), Greenland (15 km ), 60-km
resolution at mid latitudes, and 60 vertical levels. Bathymetry is from
GEBCO 2019, combined with BedMachineAntarctica around Antarctica. The
default initial condition is from EN4 1900, but PHC is also an option.

The motivation for this new mesh is enhanced resolution around
Greenland, similar to the EC15to60E2r1 mesh (#522), that should
hopefully increase the strength of the AMOC, which was weak in E3SM v1
and is even weaker in E3SM v2.
@xylar
Copy link
Collaborator Author

xylar commented Nov 6, 2020

Unfortunately, this branch was never merged and COMPASS has moved on to the point where it is not feasible to merge now.

I have created a branch on the new COMPASS repo to preserve this branch in case we need to recreate the mesh: https://github.com/MPAS-Dev/compass/tree/EC15to60E2r1

@xylar xylar closed this Nov 6, 2020
@xylar xylar deleted the refined_NA_arctic_test branch November 6, 2020 12:30
@matthewhoffman matthewhoffman added COMPASS For discussion PRs and Issues that are open for discussion and feedback in progress Ocean labels Mar 17, 2021
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

COMPASS For discussion PRs and Issues that are open for discussion and feedback in progress Ocean

Projects

None yet

Development

Successfully merging this pull request may close these issues.

8 participants