[ecco-support] Merging Arctic tiles

[Wenrui Jiang]

Hi Ash,

For pcolormesh requires X,Y, and data. Say the data is 90*90. The x,y you are using are also 90*90. This is basically telling matplotlib the center location of the box to fill in color, but it does not specify where the edges and corners of the box are.

Another way to do this is to give it X,Y with shape (91,91) that specified the corner (XG,YG) of each box. For more explanation, see here (https://matplotlib.org/3.1.1/api/_as_gen/matplotlib.pyplot.pcolormesh.html). I also recommend using shading = 'flat' in this configuration. If you are using ECCO,

I have attached the X and Y I use for plotting with dimension (13,91,91). Basically, if you replace your x,y with those files (x_g = np.load('plotting_xg,npy'); x = x_g[tile]), the bondary will disappear. If you are not using LLC90, you can generate something like this using seaduck.Topology .ind_tend (https://macekuailv.github.io/seaduck/api_reference/apiref_topology.html#seaduck.topology.Topology.ind_tend).

Let me know if that works for you.

Your friend from the hackweek,
Wenrui
[https://res.cdn.office.net/assets/mail/file-icon/png/generic_16x16.png]plotting_xg.npy<https://livejohnshopkins-my.sharepoint.com/:u:/g/personal/wjiang33_jh_edu/EQPGX1gSfKpDnruEdySmW_oBis9bHskr3YptNd0vtFEsoA>
[https://res.cdn.office.net/assets/mail/file-icon/png/generic_16x16.png]plotting_yg.npy<https://livejohnshopkins-my.sharepoint.com/:u:/g/personal/wjiang33_jh_edu/EUaVbzyaTFBKvVnwrWfu_hkB9HZeRO7L1aGyFNcgejmuow>

________________________________
From: ecco-support <ecco-support-bounces at mit.edu> on behalf of ecco-support-bounces at mit.edu <ecco-support-bounces at mit.edu>
Sent: Wednesday, December 4, 2024 11:47 AM
To: ECCO support list, wider membership <ecco-support at mit.edu>
Subject: [ecco-support] Merging Arctic tiles


      External Email - Use Caution



HI ECCO team,

I hope you guys are having a wonderful time.

I am writing this email to solve a tricky problem that I am facing as I am trying to merge the tiles to project a complete Arctic Ocean sea ice velocity map.

We know that tiles 2, 5, 6, 7, and 10 have contents of Arctic (e.g., that go beyond 65°N). I was writing this code:

[Screenshot 2024-12-04 at 11.46.01 AM.png]

But it plotted this figure:

[image.png]

While this figure is almost what I want but I was wondering if there is any way to fill the gaps between the tile. If anyone knows, please let me know.

Sincerely,
_________________
Ashfaq Ahmed
Ph.D. Candidate
Center for Fluid Mechanics
Brown University
Providence, Rhode Island, USA 02912
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mit.edu/pipermail/ecco-support/attachments/20241204/f0a38ed5/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screenshot 2024-12-04 at 11.46.01?AM.png
Type: image/png
Size: 131951 bytes
Desc: Screenshot 2024-12-04 at 11.46.01?AM.png
URL: <http://mailman.mit.edu/pipermail/ecco-support/attachments/20241204/f0a38ed5/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 228396 bytes
Desc: image.png
URL: <http://mailman.mit.edu/pipermail/ecco-support/attachments/20241204/f0a38ed5/attachment-0003.png>