A benchmark case for tide

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply

A benchmark case for tide

#1 Post by feroda »

In SVN repository, various cases are collected in the folder of '/test'.

I wanted to setup a 2-dimension run solely forced by tide. By looking through the 'test' folder, I did not find any case for reference.

I am wondering if it is physically sounding to have ROMS solely forced by tide. If it is:
1. Do zeta, ubar, vbar represent tidal amplitude and currents, respectively?
2. without any other forcing except tide, ROMS could keep running as far into the future as one
wants. Is this sort of tide prediction?

Anyone has done a 2-D case forced solely by tide, and would like to share the header file and input parameter file (.h and .in)?


User avatar
Posts: 141
Joined: Mon May 05, 2003 2:41 pm
Location: The University of Western Australia, Perth, Australia

Re: A benchmark case for tide

#2 Post by jivica »

I believe it is too simple to have it ;)
You can setup it easily without any grid or forcing, just using ana_*
For example use ana_grid.h to make rectangular basin, in that case you can increase resolution within ocean.in file where you specify number of cells, no need for recompile (specify some analytic function for bathymetry) then use for boundary again ana_fsobc.h to force elevation (and/or momentum) and that's it. Actually it is nice example of Taylor case when having Kelvin wave propagate and reflect in the channel giving amphidromic point and so on..
Just checked and all is setup for you;
have a look at the $ROMS/Functionals/ana_fsobc.h

you'll see:

Code: Select all

!  Free-surface open boundary conditions.
#if defined INLET_TEST
        DO i=IstrR,IendR
        END DO
      END IF
# endif
#elif defined KELVIN
      fac=1.0_r8                                ! zeta0
      omega=2.0_r8*pi/(12.42_r8*3600.0_r8)      ! M2 Tide period
# ifdef WEST_FSOBC
        DO j=JstrR,JendR
          val=fac*EXP(-GRID(ng)%f(Istr-1,j)*GRID(ng)%yp(Istr-1,j)/      &
     &                SQRT(g*GRID(ng)%h(Istr-1,j)))
        END DO
      END IF
# endif
# ifdef EAST_FSOBC

Post Reply