<div dir="ltr">I think this is a great suggestion.<div><br></div><div>I wonder if there is a way to directly calculate the tidal potentials within MITgcm itself. In my experience, the tidal forcing fields can be huge (hundreds of GB!). But they are essentially deterministic functions of a few orbital parameters. Calculating them "online" would save a lot of disk space and eliminate potential errors related to timing, coordinate systems, etc.</div><div><br></div><div>Joern Callies has a python code to calculate the potentials:</div><div><a href="https://naif.jpl.nasa.gov/naif/documentation.html">https://naif.jpl.nasa.gov/naif/documentation.html</a></div><div>This is based on the NASA SPICE toolkit, which has a Fortran API.</div><div><a href="https://naif.jpl.nasa.gov/naif/index.html">https://naif.jpl.nasa.gov/naif/index.html</a><br></div><div><br></div><div>An MITgcm tides package would depend on the SPICE library, in the same way mnc depends on the netCDF library.</div><div><br></div><div>It is clearly a significant amount of work. But it seems like the best long-term solution. I just thought I would raise the idea. </div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Oct 3, 2017 at 1:06 PM, Menemenlis, Dimitris (329C) <span dir="ltr"><<a href="mailto:Dimitris.Menemenlis@jpl.nasa.gov" target="_blank">Dimitris.Menemenlis@jpl.nasa.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Jean-Michel, this sounds very useful!<br>
Is Rui on MITgcm Devel list?<br>
It would be good to run your decisions listed below past him.<br>
They are definitely beyond my field of expertise.<br>
<br>
(Rui I owe you a response about tides soon.  Sorry right now I am a little overwhelmed.)<br>
<span class="HOEnZb"><font color="#888888"><br>
D.<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
> On Oct 3, 2017, at 11:47 AM, Jean-Michel Campin <<a href="mailto:jmc@mit.edu">jmc@mit.edu</a>> wrote:<br>
><br>
> Hi Dimitris and other,<br>
><br>
> Until now, when you run some global ocean with tides, the tidal forcing<br>
> was specified as an equivalent atmospheric pressure, presumably through pkg/exf.<br>
><br>
> There are few limitations with this approach:<br>
> 1) the exf atmospheric pressure field cannot be used for something else,<br>
> such as a more precise estimation of air density (#define EXF_CALC_ATMRHO).<br>
> 2) the tide geopotential is added to the ocean pressure field (totPhiHyd)<br>
> which prevent to use it inside the EOS (selectP_inEOS_Zc=2, default<br>
> with EOS: JMD95P, UNESCO, MDJWF and TEOS10).<br>
> 3) if one want to specify both an atmospheric loading and a tidal forcing,<br>
> the forcing input field become a mix, not practical nor very natural.<br>
> And this is also the case when coupled to atmospheric GCM.<br>
><br>
> Therefore, I propose to add an other 2-D forcing field to account for<br>
> tidal forcing. The most natural way would be a tidal geopotential (or<br>
> any time-dependent geopotential anomaly), units: m^2/s^2.<br>
><br>
> An other advantage is that it would be available in the same form for<br>
> Ocean in p-coords set-up as well as for atmospheric set-up.<br>
><br>
> The gradient of it would be added to horizontal momentum tendency as<br>
> part of the forcing (in S/R EXTERNAL_FORCING_U & EXTERNAL_FORCING_V)<br>
> which is different form current atmospheric loading, for instance regarding<br>
> the Adams-Bashforth (pressure gradient in or out AB depending on staggerTimeStep,<br>
> whereas forcing has it's own parameter: momForcingOutAB).<br>
><br>
> Any suggestion or comments ?<br>
><br>
> Cheers,<br>
> Jean-Michel<br>
><br>
> On Fri, Sep 29, 2017 at 04:41:26PM +0000, Menemenlis, Dimitris (329C) wrote:<br>
>> Jean-Michel and Atanas, two small coupling details that worry me:<br>
>><br>
>> 1. How do we include tides?  Right now for ocean only simulations<br>
>> we add hourly tidal potential forcing to atmospheric pressure.<br>
>> What is best way to do this in coupled model?<br>
>><br>
> ______________________________<wbr>_________________<br>
> MITgcm-devel mailing list<br>
> <a href="mailto:MITgcm-devel@mitgcm.org">MITgcm-devel@mitgcm.org</a><br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-devel" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-devel</a><br>
<br>
______________________________<wbr>_________________<br>
MITgcm-devel mailing list<br>
<a href="mailto:MITgcm-devel@mitgcm.org">MITgcm-devel@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-devel" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/<wbr>mailman/listinfo/mitgcm-devel</a><br>
</div></div></blockquote></div><br></div>