<div dir="ltr"><div>Thank you Stanislav, Ryan, Yu-Kun, Jody and Jean-Michel for the help! I think these are all very helpful suggestions. As some of you mentioned EXF, I want to ask a follow up question about this package: what kind of input file it's expecting for tidal forcing? I couldn't find an answer to this from the MITGCM's user manual. <br></div><div><br></div><div>Thanks,</div><div>Qiyu<br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Dec 4, 2021 at 12:00 PM <<a href="mailto:mitgcm-support-request@mitgcm.org">mitgcm-support-request@mitgcm.org</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">Send MITgcm-support mailing list submissions to<br>
<a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=1dA1yYQQgG8k61OS2EZM79DOLvhNmYYbHS3E-FoYqfs&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=1dA1yYQQgG8k61OS2EZM79DOLvhNmYYbHS3E-FoYqfs&e=</a> <br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:mitgcm-support-request@mitgcm.org" target="_blank">mitgcm-support-request@mitgcm.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:mitgcm-support-owner@mitgcm.org" target="_blank">mitgcm-support-owner@mitgcm.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of MITgcm-support digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Re: Tides in channel model (Jody Klymak)<br>
2. Re: Tides in channel model (Jean-Michel Campin)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Fri, 3 Dec 2021 18:13:26 +0100<br>
From: Jody Klymak <<a href="mailto:jklymak@uvic.ca" target="_blank">jklymak@uvic.ca</a>><br>
To: <a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a><br>
Subject: Re: [MITgcm-support] Tides in channel model<br>
Message-ID: <<a href="mailto:2743C27F-C1CB-4AF0-82E2-D0E698F5FB5A@uvic.ca" target="_blank">2743C27F-C1CB-4AF0-82E2-D0E698F5FB5A@uvic.ca</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
I hesitate to share this because it is undocumented, but it works fine:<br>
<br>
<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_jklymak_MITgcm-5Fbtforcing&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=LJ77rs2ArlyEBb2sDDFurZt3ZkLqzPnP0k7HBucDYD4&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_jklymak_MITgcm-5Fbtforcing&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=LJ77rs2ArlyEBb2sDDFurZt3ZkLqzPnP0k7HBucDYD4&e=</a> <<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_jklymak_MITgcm-5Fbtforcing&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=LJ77rs2ArlyEBb2sDDFurZt3ZkLqzPnP0k7HBucDYD4&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_jklymak_MITgcm-5Fbtforcing&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=LJ77rs2ArlyEBb2sDDFurZt3ZkLqzPnP0k7HBucDYD4&e=</a> ><br>
<br>
You should definitely adapt to your needs, - this one forces in y by fU for a mean, and x by a tidal body force?<br>
<br>
<br>
Cheers, Jody<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_pipermail_mitgcm-2Dsupport_attachments_20211203_4c4372ae_attachment-2D0001.html&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=pwRweddP-SliGI-TQblrzu1JdjC4zqGlx7v3PcsFNlo&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_pipermail_mitgcm-2Dsupport_attachments_20211203_4c4372ae_attachment-2D0001.html&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=pwRweddP-SliGI-TQblrzu1JdjC4zqGlx7v3PcsFNlo&e=</a> ><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Fri, 3 Dec 2021 12:28:28 -0500<br>
From: Jean-Michel Campin <<a href="mailto:jmc@mit.edu" target="_blank">jmc@mit.edu</a>><br>
To: <a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a><br>
Subject: Re: [MITgcm-support] Tides in channel model<br>
Message-ID: <<a href="mailto:20211203172828.GB46701@ocean.mit.edu" target="_blank">20211203172828.GB46701@ocean.mit.edu</a>><br>
Content-Type: text/plain; charset=us-ascii<br>
<br>
Hi Ryan,<br>
<br>
Since the tidal forcing is just variation of gravity, and gravity can be expressed<br>
as gradient of a potential, the tidal forcing can be provided as a tidal potential.<br>
This is how it's done in MITgcm, except that it's only a 2-D forcing field (neglecting vertical<br>
gradient of tidal geopotential is generally OK and fits better with hydrostatic approximation<br>
and, to a lesser extent with thin-shell approximation).<br>
<br>
It has been implemented in pkg/exf (back in October 2017):<br>
#define EXF_ALLOW_TIDES , variable "tidePot", parameters "tidePot*" (e.g.,<br>
tidePotFile, tidePotStartdate1 ...).<br>
<br>
And I am copying below the content of a message that got sent to some MITgcm email list<br>
back at that time, regarding the advantage of separating tidal geopotential from <br>
atmospheric pressure loading.<br>
<br>
Cheers,<br>
Jean-Michel<br>
<br>
-------------------------<br>
> Tue Oct 3 14:47:22 EDT 2017<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>
<br>
On Fri, Dec 03, 2021 at 11:46:58PM +0800, ?????? wrote:<br>
> Hi Ryan,<br>
> <br>
> <br>
> I have recently tried to add a body forcing to mimic the M2 tide<br>
> in a channel model (in 2D x-z plane). It is quite easy to implement<br>
> this by adding some code snippets in 'apply_forcing.F' (see attached).<br>
> I found these codes somewhere on the web (but do not remember now).<br>
> Its effect can be seen in this animation:<br>
> <a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_miniufo_InternalWaveBreaking&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=2DjJqdD4A8I_dMn9gxusC4B_Vi9JphQk2P8imYs7KBY&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_miniufo_InternalWaveBreaking&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=2DjJqdD4A8I_dMn9gxusC4B_Vi9JphQk2P8imYs7KBY&e=</a> <br>
> <br>
> <br>
> Notice that it is applied to the first 5 tiles (bi.LT.5). Hope this helps.<br>
> <br>
> <br>
> <br>
> Best regards <br>
> <br>
> Yu-Kun Qian (??????) <br>
> Center for Monsoon and Environment Research <br>
> Department of Atmospheric Sciences<br>
> School of Environmental Science and Engineering <br>
> Sun Yat-sen University <br>
> No. 135 Xingang West Road, Haizhu District <br>
> Guangzhou, 510275, P.R. China <br>
> Tel; 020-84115227 <br>
> Email: <a href="mailto:qianyk@mail3.sysu.edu.cn" target="_blank">qianyk@mail3.sysu.edu.cn</a> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> ------------------ ???????? ------------------<br>
> ??????: "Ryan Abernathey"<<a href="mailto:ryan.abernathey@gmail.com" target="_blank">ryan.abernathey@gmail.com</a>>;<br>
> ????????: 2021??12??3??(??????) ????10:21<br>
> ??????: "MITgcm-support"<<a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a>>; <br>
> <br>
> ????: Re: [MITgcm-support] Tides in channel model<br>
> <br>
> <br>
> <br>
> Thanks for the reply Stanislav! <br>
> <br>
> In this case, Qiyu is working with a zonally reentrant channel (similar to <a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__mitgcm.readthedocs.io_en_latest_examples_reentrant-5Fchannel_reentrant-5Fchannel.html&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=MXz3_QUcZCQkgP_w05JV9FpDCPDmvpD7RCBnLOaitnU&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=https-3A__mitgcm.readthedocs.io_en_latest_examples_reentrant-5Fchannel_reentrant-5Fchannel.html&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=MXz3_QUcZCQkgP_w05JV9FpDCPDmvpD7RCBnLOaitnU&e=</a> ), so there is no open boundary at which to apply tidal forcing. So although it is a small domain, it is more similar to a global model in this respect.<br>
> <br>
> <br>
> My understanding was that the most common way to force tides in global models was via a specified surface pressure field in the EXF package. (E.g. that's how the LLC4320 simulations were forced.) But I also know there has been a lot of recent work in this area, so maybe there is a simpler option.<br>
> <br>
> <br>
> Has anyone tried specifying a body force in analytic form using a custom code file? Where would be the right place in the code to insert such a forcing?<br>
> <br>
> <br>
> Best,<br>
> Ryan<br>
> <br>
> <br>
> On Fri, Dec 3, 2021 at 4:18 AM Stanislav Martyanov <<a href="mailto:martyanov.sd@gmail.com" target="_blank">martyanov.sd@gmail.com</a>> wrote:<br>
> <br>
> Hello Qiyu,<br>
> the easiest way is to generate a simple grid with constant dx (and/or dy) with one or several open boundaries specified in OBCS_OPTIONS.h (#define ALLOW_OBCS_NORTH, etc). Then, you will need to configure the open boundary conditions (if you specify induced tide) and external forcing (if you want to set geopotantial): i.e. #define OBCS_ALLOW_TIDES in OBCS_OPTIONS.h, and #define EXT_ALLOW_TIDES in EXT_OPTIONS.h<br>
> <br>
> <br>
> And, of course, you need to set the corresponding filenames with the tidal forcing in data.obcs (and in data.exf, if you want geopotential). Note that only the normal component of tidal current at the open boundaries is currently implemented in the model.<br>
> <br>
> <br>
> You also should get familiar with the OBCS package where tides can be specified (as open boundary conditions), and with the EXF package where tides can be added through geopotental.<br>
> <br>
> <br>
> Best,<br>
> Stanislav<br>
> <br>
> ????, 1 ??????. 2021 ??. ?? 21:50, Qiyu Xiao <<a href="mailto:qx344@nyu.edu" target="_blank">qx344@nyu.edu</a>>:<br>
> <br>
> Dear MITgcm users/developers,<br>
> <br>
> <br>
> If I want to run a channel model with tides, what is the best way to do that? What functions/variables should I change or is there any package for this purpose? <br>
> <br>
> <br>
> <br>
> Thanks,<br>
> Qiyu<br>
> <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="https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=1dA1yYQQgG8k61OS2EZM79DOLvhNmYYbHS3E-FoYqfs&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=1dA1yYQQgG8k61OS2EZM79DOLvhNmYYbHS3E-FoYqfs&e=</a> <br>
> <br>
> <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="https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=1dA1yYQQgG8k61OS2EZM79DOLvhNmYYbHS3E-FoYqfs&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=1dA1yYQQgG8k61OS2EZM79DOLvhNmYYbHS3E-FoYqfs&e=</a> <br>
<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="https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=1dA1yYQQgG8k61OS2EZM79DOLvhNmYYbHS3E-FoYqfs&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=1dA1yYQQgG8k61OS2EZM79DOLvhNmYYbHS3E-FoYqfs&e=</a> <br>
<br>
<br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<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="https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=1dA1yYQQgG8k61OS2EZM79DOLvhNmYYbHS3E-FoYqfs&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=5BDRdcOUIuLA0nIj7UBsNw&m=y1zd1j4BNqJV66mCcYi5svvmGr5gN1TgbfQ3HsncVPj6iFSQzQVtU6RZ1_TgDGnk&s=1dA1yYQQgG8k61OS2EZM79DOLvhNmYYbHS3E-FoYqfs&e=</a> <br>
<br>
<br>
------------------------------<br>
<br>
End of MITgcm-support Digest, Vol 222, Issue 3<br>
**********************************************<br>
</blockquote></div></div>