<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal">Hi Spencer,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">>> I was just wondering if anyone has code to open and plot LLC90 input files (the ones that are in the LLC native grid like 'GRACE_jpl_rl05m_withland_noNaN_2008').<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The answer is definitely yes, and now there is an example of its use in the new tutorial “Loading LLC compact binary files”
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><a href="https://ecco-v4-python-tutorial.readthedocs.io/ECCO_v4_Loading_LLC_compact_binary_files.html">https://ecco-v4-python-tutorial.readthedocs.io/ECCO_v4_Loading_LLC_compact_binary_files.html</a><o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The key routine is ‘<span style="font-size:10.0pt;font-family:"Courier New"">read_llc_to_tiles</span>’. It reads the llc compact format binary files in the ‘and converts to a numpy ndarrays of dimension up to [nt, nk, 13 (tiles), llc x
llc].<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The tutorial also shows how to convert these ndarrays to xarray DataArrays.
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Another example of the use of ‘read_llc_to_tiles’ can be seen in the nice new tutorial that Jan Erik Tesdal put together on closing the heat budget:<o:p></o:p></p>
<p class="MsoNormal"><a href="https://ecco-v4-python-tutorial.readthedocs.io/ECCO_v4_Heat_budget_closure.html">https://ecco-v4-python-tutorial.readthedocs.io/ECCO_v4_Heat_budget_closure.html</a><o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Look for the lines related to the geothermal heat flux. The geothermal heat flux file is in the llc compact format.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Hope this helps!<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">-Ian<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black"><ecco-support-bounces@mit.edu> on behalf of Spencer Jones <spencerj@ldeo.columbia.edu><br>
<b>Reply-To: </b>"ECCO support list, wider membership" <ecco-support@mit.edu><br>
<b>Date: </b>Monday, August 3, 2020 at 1:52 PM<br>
<b>To: </b>"ecco-support@mit.edu" <ecco-support@mit.edu><br>
<b>Subject: </b>[EXTERNAL] [ecco-support] Reading input files for llc90<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Hi, <o:p></o:p></p>
<div>
<p class="MsoNormal">I was just wondering if anyone has code to open and plot LLC90 input files (the ones that are in the LLC native grid like 'GRACE_jpl_rl05m_withland_noNaN_2008'). I always have trouble plotting faces above face=7, and I seem to remember
that a couple of grid points are removed somewhere but I can't find that information anywhere now. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Python would be ideal, but matlab will be ok too if you can point me to the exact line that does this.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Thanks very much, <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Spencer<o:p></o:p></p>
</div>
</div>
</div>
</body>
</html>