<div dir="ltr"><div dir="ltr">Hi Oliver,<div><br></div><div>You've found the error! Well spotted. By changing `field_lat_inc = 0.25` to `field_lat_inc = 249*0.25` throughout, where Ny=249 is the number of latitude points in my ERA5 subset, I was able to get sensible-looking EXF diagnostics output. <br><br>I am still a bit surprised that using a regridded dataset didn't fix my issue, though. I was thinking that I could regrid the ERA5 data onto the model grid, comment out the interpolation-related grid parameters in `data.exf`, and use the default values. Oddly, this didn't work. I might do some more tests on this later, but for now I'm happy to have something that runs. <br><br>Thanks again!<br><br>Best wishes,<br>Dan(i)</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Feb 9, 2023 at 3:34 PM Oliver Jahn <<a href="mailto:jahn@mit.edu">jahn@mit.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">Hi Dan,<br>
<br>
I noticed one thing in your data.exf: there is only 1 value for vwind_lat_inc which is an array. You have to provide vwind_nlat-1 values. (This is needed for things like Gaussian grids but definitely could use more coverage in the docs!)<br>
<br>
Oliver<br>
<br>
<br>
On 2023-02-09 09:16, Dani Jones wrote:<br>
> Hello MITgcm-support listers,<br>
> <br>
> I am having trouble with the EXF package while trying to force a <br>
> regional model configuration using external fields derived from ERA5. <br>
> The diagnostics output from EXF shows horizontal stripes or distorted <br>
> results.<br>
> <br>
> To set up my configuration, I followed these steps:<br>
> <br>
> * Downloaded a subset of ERA5 using the cdsapi package<br>
> * Regridded the data onto my model grid using the xESMF package<br>
> * Ran MITgcm for a few timesteps and checked the EXF diagnostics output <br>
> <br>
> I have conducted various checks and tried different approaches to <br>
> resolve the issue, including:<br>
> <br>
> * Examined the contents of the ERA5 data<br>
> * Examined the binary files for any problems<br>
> * Tried "on the fly" interpolation with EXF and repeatedly checked the<br>
> data.exf parameters<br>
> * Checked the precision and endian-ness throughout the entire process<br>
> * Created binary files using Matlab and Python<br>
> * Updated to the latest version of MITgcm and recompiled with the<br>
> -D_BYTESWAPIO flag<br>
> * Tried every permutation of the dimensions for the binary files <br>
> <br>
> Despite these efforts, the binary files appear to be fine, and it seems <br>
> that MITgcm/EXF is not reading them as expected. For further information <br>
> and details, I have documented this issue on GitHub with plots from the <br>
> diagnostics package and links to relevant files:<br>
> <br>
> <a href="https://github.com/so-wise/so-wise-gyre/issues/34" rel="noreferrer" target="_blank">https://github.com/so-wise/so-wise-gyre/issues/34</a> <br>
> <<a href="https://github.com/so-wise/so-wise-gyre/issues/34" rel="noreferrer" target="_blank">https://github.com/so-wise/so-wise-gyre/issues/34</a>><br>
> <br>
> Thanks in advance for any assistance, guidance, or suggestions that you <br>
> can provide.<br>
> <br>
> Best regards,<br>
> Dan(i)<br>
> <br>
> --<br>
> *Dr Dan(i) Jones* • Physical Oceanographer<br>
> Polar Oceans Team <<a href="https://www.bas.ac.uk/team/science-teams/oceans/" rel="noreferrer" target="_blank">https://www.bas.ac.uk/team/science-teams/oceans/</a>>• <br>
> BAS AI Lab <<a href="https://www.bas.ac.uk/project/ai" rel="noreferrer" target="_blank">https://www.bas.ac.uk/project/ai</a>> • British Antarctic Survey <br>
> <<a href="https://www.bas.ac.uk/" rel="noreferrer" target="_blank">https://www.bas.ac.uk/</a>><br>
> Website <<a href="http://www.danjonesocean.com/" rel="noreferrer" target="_blank">http://www.danjonesocean.com/</a>> • Twitter <br>
> <<a href="https://twitter.com/DanJonesOcean" rel="noreferrer" target="_blank">https://twitter.com/DanJonesOcean</a>> • they/them<br>
> <br>
> _______________________________________________<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>
_______________________________________________<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>