[MITgcm-support] Resources for generating curvilinear grid for MITgcm (displaced-pole Arctic)

[Sam Phillips]

Hi all,

I'm setting up a barotropic, wind-forced Arctic Ocean simulation on a
displaced-pole curvilinear grid (pole at -40E, 75N over Greenland), NX=1530
NY=520 at ~10.7km resolution, single vertical level just as a baseline.

I'm running into what appear to be grid generation problems - my mitgrid
metrics show implausible values (dxF and dxC with minimums around 22m on a
nominally 10.7km grid, dyC constant everywhere, rAz reaching zero) that I
believe are causing immediate cg2d divergence.

Before going further with debugging, I'd like to go back to basics and
generate a proper mitgrid file. I'm currently using a Python-based workflow
but I'm open to anything. A few specific questions:
1. Is there a recommended tool or workflow for generating displaced-pole
curvilinear grids for MITgcm?
2. Are there example mitgrid files or grid generation scripts for Arctic or
regional curvilinear configurations that I could use as a reference?
Googling hasn't proven very helpful.
3. Is there anything specific to watch out for when the displaced pole
falls over land (Greenland in my case) in terms of how the grid metrics
behave near the seam?

Any pointers to scripts, papers, or configurations that have worked for
similar domains would be very helpful.

Thanks,
Sam
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20260222/13664b92/attachment.html>