[MITgcm-support] Open Boundary tidal amplitude

Sat May 2 14:46:58 EDT 2020

Hi Stanislav and Jody,

Thank you both for your reply and helpful suggestions. I’ll look into it for a good solution.

All the best,

Weiyang

Weiyang Bao
Ph.D. Student, Physical Ocean Science and Engineering
School of Marine Science and Policy
University of Delaware
wbao at udel.edu
On May 1, 2020, 6:58 PM -0400, mitgcm-support-request at mitgcm.org, wrote:
> Send MITgcm-support mailing list submissions to
> mitgcm-support at mitgcm.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> or, via email, send a message with subject or body 'help' to
> mitgcm-support-request at mitgcm.org
>
> You can reach the person managing the list at
> mitgcm-support-owner at mitgcm.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of MITgcm-support digest..."
>
>
> Today's Topics:
>
> 1. Open boundary tidal amplitude (Weiyang Bao)
> 2. Re: Open boundary tidal amplitude (Stanislav Martyanov)
> 3. Re: Open boundary tidal amplitude (Jody Klymak)
> 4. SSH problems with nonlinFreeSurf and Star (Stanislav Martyanov)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Fri, 1 May 2020 12:42:09 -0400
> From: Weiyang Bao <wbao at udel.edu>
> To: mitgcm-support at mitgcm.org
> Subject: [MITgcm-support] Open boundary tidal amplitude
> Message-ID: <38112c6c-8ac6-4e24-8d4c-60baf519701f at Spark>
> Content-Type: text/plain; charset="utf-8"
>
> Hi MITgcm users,
>
> I have a question when trying to apply tidal forcing at the open boundary of my model domain, which is an idealized 3-D cuboid estuary. The prescribed tidal amplitude input for the OBCS package is velocity (unit of m/s), but I would like to specify the tidal signal using sea surface elevation (or sea level height) instead. Although this might be an old problem, I could find little discussion through previous posts. Could anyone please give me a hand? Thank you in advance!
>
> Best,
>
> Weiyang
>
>
> Weiyang Bao
> Ph.D. Student, Physical Ocean Science and Engineering
> School of Marine Science and Policy
> University of Delaware
> wbao at udel.edu
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200501/621fce01/attachment-0001.html>
>
> ------------------------------
>
> Message: 2
> Date: Fri, 1 May 2020 22:36:21 +0300
> From: Stanislav Martyanov <martyanov.sd at gmail.com>
> To: mitgcm-support at mitgcm.org
> Subject: Re: [MITgcm-support] Open boundary tidal amplitude
> Message-ID:
> <CANJ4tTN7BLAqTHG5c_27AkQQK3S0yWzmkFK0k9TTr1WNDrzwNA at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Hi Weiyang!
>
> According to our experience in solving the same problem with MITgcm, the
> only method implemented in the model is prescribing the tidal velocity at
> the open boundaries. We also tried setting the tidal level oscillations,
> but with no result.
>
> Also, you may find interesting the following thread:
> http://mailman.mitgcm.org/pipermail/mitgcm-support/2016-February/010259.html
>
> Push the 'Next message' every time to read the discussions. Also:
> http://mailman.mitgcm.org/pipermail/mitgcm-support/2016-February/010269.html
>
>
> Regards,
> Stanislav Martyanov
>
> ??, 1 ??? 2020 ?. ? 20:21, Weiyang Bao <wbao at udel.edu>:
>
> > Hi MITgcm users,
> >
> > I have a question when trying to apply tidal forcing at the open boundary
> > of my model domain, which is an idealized 3-D cuboid estuary. The
> > prescribed tidal amplitude input for the OBCS package is velocity (unit of
> > m/s), but I would like to specify the tidal signal using sea surface
> > elevation (or sea level height) instead. Although this might be an old
> > problem, I could find little discussion through previous posts. Could
> > anyone please give me a hand? Thank you in advance!
> >
> > Best,
> >
> > Weiyang
> >
> >
> > *Weiyang Bao*
> > Ph.D. Student, Physical Ocean Science and Engineering
> > School of Marine Science and Policy
> > University of Delaware
> > wbao at udel.edu
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> >
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200501/e0df3aec/attachment-0001.html>
>
> ------------------------------
>
> Message: 3
> Date: Fri, 1 May 2020 13:47:13 -0700
> From: Jody Klymak <jklymak at uvic.ca>
> To: mitgcm-support at mitgcm.org
> Subject: Re: [MITgcm-support] Open boundary tidal amplitude
> Message-ID: <AECFA7AC-D26B-4FBA-9E61-7F1DDFE9E80D at uvic.ca>
> Content-Type: text/plain; charset="utf-8"
>
> Hi Weiyang,
>
> Ken Hughes (OSU) implemented surface nudging, which is essentially what you want. I don?t know if he ever made it into a package, or if it is just a patch.
>
> Cheers, Jody
>
> > On May 1, 2020, at 9:42 AM, Weiyang Bao <wbao at udel.edu> wrote:
> >
> > Hi MITgcm users,
> >
> > I have a question when trying to apply tidal forcing at the open boundary of my model domain, which is an idealized 3-D cuboid estuary. The prescribed tidal amplitude input for the OBCS package is velocity (unit of m/s), but I would like to specify the tidal signal using sea surface elevation (or sea level height) instead. Although this might be an old problem, I could find little discussion through previous posts. Could anyone please give me a hand? Thank you in advance!
> >
> > Best,
> >
> > Weiyang
> >
> >
> > Weiyang Bao
> > Ph.D. Student, Physical Ocean Science and Engineering
> > School of Marine Science and Policy
> > University of Delaware
> > wbao at udel.edu
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200501/2bc20b55/attachment-0001.html>
>
> ------------------------------
>
> Message: 4
> Date: Sat, 2 May 2020 01:57:44 +0300
> From: Stanislav Martyanov <martyanov.sd at gmail.com>
> To: mitgcm-support at mitgcm.org
> Subject: [MITgcm-support] SSH problems with nonlinFreeSurf and Star
> Message-ID:
> <CANJ4tTN1w+gBGsFQe+u+_3Vt5S9DTn+r1avF9bLoBJt7utPwfw at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Hello everyone!
>
> Having tested my several MITgcm configurations (details here:
> http://mailman.mitgcm.org/pipermail/mitgcm-support/2020-April/012462.html)
> I encountered an issue with the model results. For unknown reasons the sea
> level has high oscillations and drops (up to minus 3-6 m) which can not be
> explained by the prescribed external forcing and OBCS (they are OK, SSH =
> plus/minus 1 m at the boundary). Also, delR=2 m in the upper layers, and
> minimum depth in the domain is 5 m. No tides, no sponge.
>
> At each time step the model used SURF_ADJUSTMENT...
>
> Below are some parameters from the data file used for the reference run (?
> 0):
> &PARM01
> hFacMin = 0.2,
> hFacMinDr = 0.5,
> hFacInf = 0.2,
> hFacSup = 2.0,
> vectorInvariantMomentum = .TRUE.,
> highOrderVorticity = .TRUE.,
> rigidLid = .FALSE.,
> implicitFreeSurface = .TRUE.,
> nonlinFreeSurf = 4,
> select_rStar = 0,
> exactConserv = .TRUE.,
> useRealFreshWaterFlux = .TRUE.,
> staggerTimeStep = .TRUE.,
> multiDimAdvection = .TRUE.,
> implicitViscosity = .TRUE.,
> implicitDiffusion = .TRUE.,
> viscAhGrid = 0.0,
> viscAhGridMax = 0.5,
> viscC2leith = 1.0,
> viscC2leithD = 1.0,
> useFullLeith = .TRUE.,
> viscAr = 1.0e-4,
> no_slip_sides = .FALSE.,
> no_slip_bottom = .TRUE.,
> bottomDragQuadratic = 0.0025,
> tempAdvScheme = 33,
> saltAdvScheme = 33,
> diffKhT = 0.0,
> diffKhS= 0.0,
> diffKrT = 1.0e-4,
> diffKrS = 1.0e-4,
> useSingleCpuIO = .TRUE.,
> &
> forcing_In_AB = .FALSE.,
> momDissip_In_AB = .FALSE.,
> cAdjFreq = -1.0,
> usingCurvilinearGrid = .TRUE.,
>
> Having spotted such erroneous results in STDOUT statistics, I started
> tuning some model's functionality to figure out what happened. Some
> conclusions for the subsequent runs:
> 1) I turn off the river runoff everywhere in OBCS. The problem remained.
> 2) + I removed eta boundary conditions from data.obcs. The problem
> remained.
> 3) + I set useOBCS = .FALSE. in data.pkg, thus closed all the open
> boundaries, suspecting them. But the problem remained.
> 4) I returned to the initial settings (run ? 0, see above), but set:
> nonlinFreeSurf
> = 0 and linFSConserveTr = .TRUE., and the problem vanished, no
> SURF_ADJUSTMENT in the STDOUT any more!
> 5) I returned to the initial settings (run ? 0, see above), but set:
> nonlinFreeSurf
> = 4 and select_rStar = 2, , and get the followng message from one of tiles:
> fail at i,j= 44 39 ; rStarFacC,H,eta = 0.199701 5.000000E+00
> -4.001496E+00
> WARNING: r*FacC < hFacInf at 1 pts : bi,bj,Thid,Iter= 1 1 1
> 350
> STOP in CALC_R_STAR : too SMALL rStarFac[C,W,S] !
>
> Could you please help me to figure out:
> a) If I want to use non-linear free suface but do not want to use rStar,
> may I set select_rStar=0 and nonlinFreeSurf=4 simultaniously? In the
> verification 'internal_wave' I found nonlinFreeSurf=3 and select_rStar=0.
> So, for select_rStar=0 the parameter nonlinFreeSurf=0 is not a must?
> b) When selecting nonlinFreeSurf=4 and select_rStar=2 the model crushes
> complaining about too SMALL rStarFac[C,W,S]. My vertical resolution is too
> small (2 m) or I should change the hFacMin or hFacDr somehow?
> c) Where does the key DISABLE_RSTAR_CODE come from? It is used in many
> checks throughout the code, but I could not find where it is defined.
> Should I change it if I do not want to use rStar? Or using just
> select_rStar=0 is enough?
> d) Does the parameter hFacInf works only in the case of nonlinFreeSurf !=
> 0, keeping the upper layer alive? What if I use nonlinFreeSurf= 0 and at
> some moment -eta > dr(1)? The upper layer thickness may become 0 when
> nonlinFreeSurf=0?
>
> Thanks!
>
> Regards,
> Stanislav Martyanov
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200502/a7f087a6/attachment.html>
>
> ------------------------------
>
> Subject: Digest Footer
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>
>
> ------------------------------
>
> End of MITgcm-support Digest, Vol 203, Issue 1
> **********************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200502/d81b88ba/attachment.html>