Skip to content

Conversation

@pwolfram
Copy link
Contributor

This adds an additional option that allows direct forcing of the tidal boundary, instead
of use of a thickness tendency to apply the tidal boundary.

The test case has been distinguished into these two separate options:

  • thickness_source that uses restoring to apply the tidal boundary condition via a thickness tendency
  • direct that sets the layer thicknesses to obtain a desired sea surface height

@pwolfram
Copy link
Contributor Author

Testing

thickness_source

tidalcomparison

direct

tidalcomparison

Note that the direct approach has a clamped boundary that perfectly respects the applied tidal forcing, whereas the thickness_source approach has the lag as noted in #206

@pwolfram
Copy link
Contributor Author

@mark-petersen, note I only implemented this for RK4 because at present because of #217, although we could probably extend this to split explicit without too much trouble. Given that a primary need for this is to support testing of wetting and drying, which is only implemented for RK4 at present, I'm not sure we need to support split explicit to start.

@pwolfram
Copy link
Contributor Author

Hm, this is interesting @mark-petersen. When I rebased onto ocean/develop the thickness_source approach improved in quality...

tidalcomparison

Did you update the thickness tendency terms anytime recently? The phase lag is decreased, which is great!

@pwolfram pwolfram force-pushed the coastal/clamped_tidal_BC branch from ea4b726 to 29f40d0 Compare June 18, 2019 14:28
This ensures that exact tidal boundary can be obtained
@pwolfram pwolfram force-pushed the coastal/clamped_tidal_BC branch from 29f40d0 to 7fd8c26 Compare June 18, 2019 14:49
@pwolfram
Copy link
Contributor Author

@mark-petersen, this also is needed by the drying slope cases too and it works when incorporated. I'm ready for it to be merged pending your checks (using the drying slope prototype as a half-verification). Note the tidal amplitude is forced on the left side and results are as expected too using the direct approach, much better than the thickness_source:

dryingslopecomparison
tidalcomparison

Copy link
Contributor

@mark-petersen mark-petersen left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@pwolfram if you are happy with this, go ahead and merge. I will do full testing and bfb comparisons when I merge coastal into ocean/develop.

@pwolfram pwolfram merged commit 7fd8c26 into MPAS-Dev:ocean/coastal Jun 21, 2019
pwolfram added a commit that referenced this pull request Jun 21, 2019
Adds option to directly force with tidal boundary

This ensures that exact tidal boundary can be obtained via
the new 'config_tidal_forcing_type'='direct' option, which is
a direct application of the free surface.
@pwolfram pwolfram deleted the coastal/clamped_tidal_BC branch June 21, 2019 19:29
@pwolfram
Copy link
Contributor Author

Thanks @mark-petersen!

mark-petersen added a commit that referenced this pull request Oct 19, 2019
This PR is an accumulation of PRs into the ocean/coastal branch. We are
bringing them in at once for efficiency:
 #285, #289, #284, #295, #310, #311, #312, #335, #354, #356, #358, #359,
 #365, #371
mark-petersen added a commit that referenced this pull request Oct 19, 2019
This PR is an accumulation of PRs into the ocean/coastal branch. We are
bringing them in at once for efficiency.
 #285, #289, #284, #295, #310, #311, #312, #335, #354, #356, #358, #359,
 #365, #371
ashwathsv pushed a commit to ashwathsv/MPAS-Model that referenced this pull request Jul 21, 2020
…stal

Adds option to directly force with tidal boundary

This ensures that exact tidal boundary can be obtained via
the new 'config_tidal_forcing_type'='direct' option, which is
a direct application of the free surface.
ashwathsv pushed a commit to ashwathsv/MPAS-Model that referenced this pull request Jul 21, 2020
mark-petersen pushed a commit to mark-petersen/MPAS-Model that referenced this pull request Jan 11, 2021
…stal

Adds option to directly force with tidal boundary

This ensures that exact tidal boundary can be obtained via
the new 'config_tidal_forcing_type'='direct' option, which is
a direct application of the free surface.
mark-petersen added a commit to mark-petersen/MPAS-Model that referenced this pull request Jan 11, 2021
caozd999 pushed a commit to caozd999/MPAS-Model that referenced this pull request Jan 14, 2021
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants