# Boundary Conditions

Boundary Conditions

## Horizontal Boundary Conditions

ROMS comes with a variety of boundary conditions, including open, closed, and periodic. See Marchesiello et al. (2001) for a more thorough exploration of the options. Some options require a value for the boundary points from either an included analytic expression (Functionals) or from an external NetCDF file. Here, represents the exterior value of a quantity .

This boundary condition is extremely simple and consists of setting the gradient of a field to zero at the edge. The outside value is set equal to the closest interior value.

### Wall boundary condition

ROMS now assumes a wall condition if no other boundary condition is chosen. This is a zero gradient condition for tracers and the surface elevation and zero flow for the normal velocity. For tangential velocities, the wall is treated as either no-slip or free-slip, depending on the value of gamma2 chosen by the user.

### Clamped boundary condition

Almost as simple is setting the boundary value to a known exterior value.

 $phi equals phi Superscript normal e normal x normal t$ (1)

### Flather boundary condition

For the normal component of the barotropic velocity, one option is to radiate out deviations from exterior values at the speed of the external gravity waves ( Flather (1976)):

 $u overbar equals u overbar Superscript normal e normal x normal t Baseline minus StartRoot StartFraction g Over upper D EndFraction EndRoot left-parenthesis zeta minus zeta Superscript normal e normal x normal t Baseline right-parenthesis$ (2)

The exterior values are often used to provide tidal boundary contitions to the barotropic mode. However, there are times when only the tidal elevation is known. A reduced physics option is available for estimating in that case.

### Chapman boundary condition

The corresponding condition for surface elevation was investigated by Chapman (1985), assuming all outgoing signals leave at the shallow-water wave speed of . This can be useful when using the Flather condition on the 2-D momentum equations.

 $StartFraction normal partial-differential zeta Over normal partial-differential t EndFraction equals plus-or-minus StartRoot g upper D EndRoot StartFraction normal partial-differential zeta Over normal partial-differential xi EndFraction$ (3)

The time derivative here can be handled either explicitly or implicitly. The model uses an implicit timestep, with the term being evaluated at the new timestep.

In realistic domains, open boundary conditions can be extremely difficult to get right. There are situations in which incoming flow and outgoing flow happen along the same boundary or even at different depths at the same horizontal location. Orlanski (1976) proposed a radiation scheme in which a local normal phase velocity is computed and used to radiate things out (if it is indeed going out). This works well for a wave propagating normal to the boundary, but has problems when waves approach the boundary at an angle. Raymond and Kuo (1984) have modified the scheme to account for propagation in all three directions. In ROMS, only the two horizontal directions are accounted for (with the recommended RADIATION_2D option):

 $StartFraction normal partial-differential phi Over normal partial-differential t EndFraction equals minus left-parenthesis phi Subscript xi Baseline StartFraction normal partial-differential phi Over normal partial-differential xi EndFraction plus phi Subscript eta Baseline StartFraction normal partial-differential phi Over normal partial-differential eta EndFraction right-parenthesis right-parenthesis$ (4)

where

 $StartLayout 1st Row 1st Column phi Subscript xi 2nd Column equals StartStartFraction upper F StartFraction normal partial-differential phi Over normal partial-differential xi EndFraction OverOver left-parenthesis StartFraction normal partial-differential phi Over normal partial-differential xi EndFraction right-parenthesis squared plus left-parenthesis StartFraction normal partial-differential phi Over normal partial-differential eta EndFraction right-parenthesis squared EndEndFraction 2nd Row 1st Column phi Subscript eta 2nd Column equals StartStartFraction upper F StartFraction normal partial-differential phi Over normal partial-differential eta EndFraction OverOver left-parenthesis StartFraction normal partial-differential phi Over normal partial-differential xi EndFraction right-parenthesis squared plus left-parenthesis StartFraction normal partial-differential phi Over normal partial-differential eta EndFraction right-parenthesis squared EndEndFraction 3rd Row 1st Column upper F 2nd Column equals minus StartFraction normal partial-differential phi Over normal partial-differential t EndFraction EndLayout$ (5)

These terms are evaluated at the closest interior point in a manner consistent with the time stepping scheme used. The phase velocities are limited so that the local CFL condition is satisfied. They are then applied to the boundary point using equation (4), again using a consistent time stepping scheme. Raymond and Kuo give the form used for centered differencing and a leapfrog time step while ROMS uses one-sided differences.

The radiation approach is appropriate for waves leaving the domain. A check is made to see which way the phase velocity is headed. If it is entering the domain, a zero gradient condition is applied unless the next option is also specified.