<div dir="ltr"><div>Hi Wenrui, </div><div><br></div><div>I hope you had a wonderful winter break and happy new year (belated!)</div><div><br></div><div>I tried your method, and I was able to get the 91-by-91-shaped tiles for the Arctic Ocean. Just like this - </div><div><br></div><img src="cid:ii_m65j6ifd2" alt="Screenshot 2025-01-20 at 4.00.49 PM.png" width="421" height="186" style="margin-right: 0px;"><br><div><br></div><div>However, one problem that I am facing is that the velocity map is itself a 90x90 matrix. So this code is showing an error in plotting. </div><div><br></div><div><img src="cid:ii_m65j9p8j4" alt="Screenshot 2025-01-20 at 4.03.18 PM.png" width="562" height="258"><br></div><img src="cid:ii_m65j7q8f3" alt="Screenshot 2025-01-20 at 4.01.45 PM.png" width="562" height="241"><div><br></div><div>Can you please suggest how can I pcolormesh them into the 91x91 grid system?</div><div><br><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">Your friend from the ECCO hackathon,<div>_________________</div><div><b>Ashfaq Ahmed</b></div><div>Ph.D. Candidate</div><div>Center for Fluid Mechanics</div><div>Brown University</div><div>Providence, Rhode Island, USA 02912</div></div></div></div><br></div></div><br><div class="gmail_quote gmail_quote_container"><div dir="ltr" class="gmail_attr">On Sat, Dec 7, 2024 at 9:04 AM Wenrui Jiang <<a href="mailto:wjiang33@jhu.edu">wjiang33@jhu.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg-2916281213894783649">
<div dir="ltr">
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Hi Ash,</div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
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.</div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
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 (<a href="https://matplotlib.org/3.1.1/api/_as_gen/matplotlib.pyplot.pcolormesh.html" id="m_-2916281213894783649LPlnk315041" target="_blank">https://matplotlib.org/3.1.1/api/_as_gen/matplotlib.pyplot.pcolormesh.html</a>).
I also recommend using shading = 'flat' in this configuration. If you are using ECCO,<br>
<br>
</div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
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 (<a href="https://macekuailv.github.io/seaduck/api_reference/apiref_topology.html#seaduck.topology.Topology.ind_tend" id="m_-2916281213894783649LPlnk696844" target="_blank">https://macekuailv.github.io/seaduck/api_reference/apiref_topology.html#seaduck.topology.Topology.ind_tend</a>).</div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Let me know if that works for you.</div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Your friend from the hackweek,</div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Wenrui<br>
<br>
</div>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
<br>
</div>
<div id="m_-2916281213894783649appendonsend"></div>
<hr style="display:inline-block;width:98%">
<div id="m_-2916281213894783649divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> ecco-support <<a href="mailto:ecco-support-bounces@mit.edu" target="_blank">ecco-support-bounces@mit.edu</a>> on behalf of <a href="mailto:ecco-support-bounces@mit.edu" target="_blank">ecco-support-bounces@mit.edu</a> <<a href="mailto:ecco-support-bounces@mit.edu" target="_blank">ecco-support-bounces@mit.edu</a>><br>
<b>Sent:</b> Wednesday, December 4, 2024 11:47 AM<br>
<b>To:</b> ECCO support list, wider membership <<a href="mailto:ecco-support@mit.edu" target="_blank">ecco-support@mit.edu</a>><br>
<b>Subject:</b> [ecco-support] Merging Arctic tiles</font>
<div> </div>
</div>
<div>
<p><strong><br>
</strong><strong><span style="background:red;color:white;font-family:"Times New Roman",serif;font-size:12pt"> External Email - Use Caution </span></strong></p>
<p><strong><span style="background:red;color:white;font-family:"Times New Roman",serif;font-size:12pt"><br>
</span></strong></p>
<p><strong><span style="background:red;color:white;font-family:"Times New Roman",serif;font-size:12pt"><br>
</span></strong></p>
<div>
<div dir="ltr">
<div>HI ECCO team, </div>
<div><br>
</div>
<div>I hope you guys are having a wonderful time.</div>
<div><br>
</div>
<div>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. </div>
<div><br>
</div>
<div>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:<br>
<br>
<img alt="Screenshot 2024-12-04 at 11.46.01 AM.png" width="515" height="239" style="margin-right: 0px;" src="cid:ii_19485828d88d8824bed1"><br>
</div>
<div><br>
</div>
<div>But it plotted this figure:<br>
<br>
</div>
<div><img alt="image.png" width="459" height="542" src="cid:ii_19485828d88cb971f162"><br>
</div>
<div><br>
</div>
<div>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. </div>
<div><br>
</div>
<div>
<div dir="ltr">
<div dir="ltr">Sincerely,
<div>_________________</div>
<div><b>Ashfaq Ahmed</b></div>
<div>Ph.D. Candidate</div>
<div>Center for Fluid Mechanics</div>
<div>Brown University</div>
<div>Providence, Rhode Island, USA 02912</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
_______________________________________________<br>
ecco-support mailing list<br>
<a href="mailto:ecco-support@mit.edu" target="_blank">ecco-support@mit.edu</a><br>
<a href="https://mailman.mit.edu/mailman/listinfo/ecco-support" rel="noreferrer" target="_blank">https://mailman.mit.edu/mailman/listinfo/ecco-support</a><br>
</div></blockquote></div>