<div dir="ltr"><div dir="ltr"><div>Good morning Dani</div><div><br></div><div>May I ask why you are regridding the data for use by the EXF package?</div><div>My recollection is that it does on the fly grid interpolation, which is probably the best feature of EXF.<br></div><div>Here is what the documentation says:</div><div><a href="https://mitgcm.readthedocs.io/en/latest/phys_pkgs/exf.html">https://mitgcm.readthedocs.io/en/latest/phys_pkgs/exf.html</a></div><div>"An interpolation
routine provides on-the-fly interpolation of forcing fields an arbitrary
grid onto the model grid."</div><div>Check the EXF package CPP options that enable this interpolation also (link above).</div><div>You probably need the CAL package to get the timing right.<br></div><div><br></div><div>You may still need to convert to binary format (say, if not using netCDF/MNC package,).</div><div>If using the MDS package for IO, you still need swap bytes (it stuns me that MITgcm continues to insist in big-endian), <br></div><div>if needed, ensure that lat,lon,lev,time go in the correct order, create consistent data.exf  and data.cal namelists, etc.</div><div>I believe there are some Matlab routines in the MITgcm tarball that may help to convert to binary.<br></div><div><br></div><div>My two cents,</div><div>Gus Correa</div><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Feb 9, 2023 at 9:17 AM Dani Jones <<a href="mailto:dcjones.work@gmail.com">dcjones.work@gmail.com</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 dir="ltr">Hello MITgcm-support listers,<br><div><br></div><div>I am having trouble with the EXF package while trying to force a regional model configuration using external fields derived from ERA5. The diagnostics output from EXF shows horizontal stripes or distorted results.  </div><div><br></div><div>To set up my configuration, I followed these steps:  </div><div><ul><li style="margin-left:15px">Downloaded a subset of ERA5 using the cdsapi package </li><li style="margin-left:15px">Regridded the data onto my model grid using the xESMF package </li><li style="margin-left:15px">Ran MITgcm for a few timesteps and checked the EXF diagnostics output </li></ul></div><div>I have conducted various checks and tried different approaches to resolve the issue, including:  </div><div><ul><li style="margin-left:15px">Examined the contents of the ERA5 data </li><li style="margin-left:15px">Examined the binary files for any problems </li><li style="margin-left:15px">Tried "on the fly" interpolation with EXF and repeatedly checked the data.exf parameters </li><li style="margin-left:15px">Checked the precision and endian-ness throughout the entire process </li><li style="margin-left:15px">Created binary files using Matlab and Python </li><li style="margin-left:15px">Updated to the latest version of MITgcm and recompiled with the -D_BYTESWAPIO flag </li><li style="margin-left:15px">Tried every permutation of the dimensions for the binary files </li></ul></div><div>Despite these efforts, the binary files appear to be fine, and it seems that MITgcm/EXF is not reading them as expected. For further information and details, I have documented this issue on GitHub with plots from the diagnostics package and links to relevant files: </div><div><br></div><div><a href="https://github.com/so-wise/so-wise-gyre/issues/34" target="_blank">https://github.com/so-wise/so-wise-gyre/issues/34</a>  </div><div><br></div><div>Thanks in advance for any assistance, guidance, or suggestions that you can provide.  </div><div><br></div><div>Best regards, </div><div>Dan(i)<br><br></div><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div style="color:rgb(0,0,0);font-family:Calibri,Arial,Helvetica,sans-serif;font-size:small"><div style="margin:0px;padding:0px;border:0px none;font-stretch:inherit;line-height:inherit;vertical-align:baseline">--<br><b>Dr Dan(i) Jones</b> • Physical Oceanographer </div><div style="margin:0px;padding:0px;border:0px none;font-stretch:inherit;line-height:inherit;vertical-align:baseline"><a href="https://www.bas.ac.uk/team/science-teams/oceans/" title="https://www.bas.ac.uk/team/science-teams/oceans/" style="margin:0px;padding:0px;border:0px none;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;vertical-align:baseline" target="_blank">Polar Oceans Team</a><span style="margin:0px;padding:0px;border:0px none;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;vertical-align:baseline;color:inherit"><span style="margin:0px;padding:0px;border:0px none;font:inherit;vertical-align:baseline;color:inherit"> </span>• </span><a href="https://www.bas.ac.uk/project/ai" title="https://www.bas.ac.uk/project/ai" style="margin:0px;padding:0px;border:0px none;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;vertical-align:baseline" target="_blank">BAS AI Lab</a><span style="margin:0px;padding:0px;border:0px none;font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;line-height:inherit;vertical-align:baseline;color:inherit"> • </span><a href="https://www.bas.ac.uk/" title="https://www.bas.ac.uk/" style="margin:0px;padding:0px;border:0px none;font:inherit;vertical-align:baseline" target="_blank">British Antarctic Survey </a></div><div style="margin:0px;padding:0px;border:0px none;font-stretch:inherit;line-height:inherit;vertical-align:baseline"><a href="http://www.danjonesocean.com/" title="http://www.danjonesocean.com/" style="margin:0px;padding:0px;border:0px none;font:inherit;vertical-align:baseline" target="_blank">Website</a> • <a href="https://twitter.com/DanJonesOcean" title="https://twitter.com/DanJonesOcean" style="margin:0px;padding:0px;border:0px none;font:inherit;vertical-align:baseline" target="_blank">Twitter</a> • they/them</div></div></div></div></div></div></div></div></div></div></div></div></div>
_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</blockquote></div></div>