question-mark
Stuck on an issue?

Lightrun Answers was designed to reduce the constant googling that comes with debugging 3rd party libraries. It collects links to all the places you might be looking at while hunting down a tough bug.

And, if you’re still stuck at the end, we’re happy to hop on a call to see how we can help out.

I’m trying to regrid the cmip5 data from the ocean model NEMO, which used a tripolar grid (still curvilinear). It seems none of the bilinear, patch and conservative method works. I thought bilinear and patch do not require grid corners but do they require some other specific grid configurations? The conservative method doesn’t work because I can’t figure out the corners. The corners in the data are specified as 4 vertices but the orientations seem not constant. I mean in the tripolar grid the four faces of the grid are not necessarily oriented east-west or north-south, what is the orientation that is acceptable? The notebook is uploaded here. https://github.com/JianghuiDu/cmip5

Issue Analytics

  • State:open
  • Created 5 years ago
  • Comments:13 (6 by maintainers)

github_iconTop GitHub Comments

1reaction
JianghuiDucommented, Sep 7, 2018

Thanks! Finally figured out the orientation convention.

0reactions
bekozicommented, Sep 7, 2018

CF convention, as far as I know, follows from SCRIP which ESMF adheres to closely: http://cfconventions.org/Data/cf-conventions/cf-conventions-1.7/cf-conventions.html#cell-boundaries. And, yes, cell boundaries should defined in a CCW manner in ESMF.

I put a short notebook together to demonstrate how the corner conversion functions I linked to are used. This was for my own benefit as well. ocgis corners may be out of compliance with CF…never confirmed it was compliant necessarily (hence the lack of convention in the function names). You may find you need to adjust indexing. In case it’s useful, this is where the (m, n, 4) conversion to ESMF corners is used in ocgis-esmpy.

I guess that get_esmf_corners_from_ocgis_corners() assumes Fortran-ordering (lon, lat)?

It does produce Fortran ordering, but it only works off a single coordinate array. Need to call the function twice for lat and lon corners.

Read more comments on GitHub >

github_iconTop Results From Across the Web

ORCA family: global ocean with tripolar grid
The ORCA grid is a tripolar is based on the semi-analytical method of Madec and Imbard [1996]. It allows to construct a global...
Read more >
CM2.X model ocean and sea ice grids ... - GFDL Data Portal
Why the tripolar grid? To avoid a singularity at the North Pole when running the model, the CM2.X ocean and sea ice model...
Read more >
The tripolar grid, often used in ocean modeling. Polar ...
The tripolar grid, often used in ocean modeling. Polar singularities are placed over land and excluded from the simulation. Source publication.
Read more >
ORCA family: global ocean with tripolar grid - Forge
ORCA family: global ocean with tripolar grid ... Figure: ORCA mesh conception. The departure from an isotropic Mercator grid start poleward of 20....
Read more >
Spatial selection with tripolar ACCESS-OM2 grid
The ACCESS-OM2 model collection utilises a tripolar grid north of 65N to avoid a point of convergence over the ocean at the north...
Read more >

github_iconTop Related Medium Post

No results found

github_iconTop Related StackOverflow Question

No results found

github_iconTroubleshoot Live Code

Lightrun enables developers to add logs, metrics and snapshots to live code - no restarts or redeploys required.
Start Free

github_iconTop Related Reddit Thread

No results found

github_iconTop Related Hackernoon Post

No results found

github_iconTop Related Tweet

No results found

github_iconTop Related Dev.to Post

No results found

github_iconTop Related Hashnode Post

No results found