Hello Everyone,

I am trying to set a initial condition of geostrophically balanced current in a channel.

I set velocity field first and calculate the temperature field through geostrophic balance and linear relationship between the density and temperature as salinity is set to be zero. However, strong internal wave generated and break the geostrophic field. Anyone have a suggestion on how to set a geostrophic flow for initial condition? Thank you!

## Set geostrphically balanced initial current

### Re: Set geostrphically balanced initial current

It is usually works the other way around: density field + free surface (via level of no motion or some other idea of this sort) ==> pressure field ==> pressure gradient ==> velocity fieldI set velocity field first and calculate the temperature field through geostrophic balance

### Re: Set geostrphically balanced initial current

Thank you very much! it works well now!

shchepet wrote:It is usually works the other way around: density field + free surface (via level of no motion or some other idea of this sort) ==> pressure field ==> pressure gradient ==> velocity fieldI set velocity field first and calculate the temperature field through geostrophic balance

### Re: Set geostrphically balanced initial current

There is one more thing to add to this topic specifically related to the fact that you may

want not just geostrophic balance, but

ROMS uses C-grid, which means that discrete

subject to 4-point horizontal averaging, so the discrete geostrophic balance exists in a bit

awkward sense: <4-point-averaged u> = -(1/f) * (1/rho0) * dp/dy

and <4-point-averaged v> = +(1/f) * (1/rho0) * dp/dx where p is defined at rho-points.

The problem is that imagine that p contains checker-board mode, so dp/dy and dp/dx would

contain it as well because the differencing is only over 1 dx and 1 dy in each direction.

However 4-pint averaged velocity components cannot admit it, so there is no way to balance.

To address this one have to ensure that checkerboard mode should be excluded entirely from

the pressure and tracer fields defined at RHO points. This leads to a rather simple recipe

how to generate analytical geostrophically-balanced initial conditions in ROMS (or any other

C-gridded model):

1. Define temporal T,S (or density field) and free surface at horizontal

place them into scratch arrays.

2. Compute U and V by vertically integrating density computed from these temporal field

(still defined at vorticity points) and differencing the resultant pressure field (naturally

differencing over 1 dx in each direction).

3. Horizontally average T,S and vorticity fields defined in (1) using 4-point averaging, so

the outcome in at RHO-points. Accept them as the final fields.

want not just geostrophic balance, but

*the geostrophic balance as it would be felt by the*

discrete model, i.e., ROMS.discrete model

ROMS uses C-grid, which means that discrete

*u*s and*v*s in Coriolis terms aresubject to 4-point horizontal averaging, so the discrete geostrophic balance exists in a bit

awkward sense: <4-point-averaged u> = -(1/f) * (1/rho0) * dp/dy

and <4-point-averaged v> = +(1/f) * (1/rho0) * dp/dx where p is defined at rho-points.

The problem is that imagine that p contains checker-board mode, so dp/dy and dp/dx would

contain it as well because the differencing is only over 1 dx and 1 dy in each direction.

However 4-pint averaged velocity components cannot admit it, so there is no way to balance.

To address this one have to ensure that checkerboard mode should be excluded entirely from

the pressure and tracer fields defined at RHO points. This leads to a rather simple recipe

how to generate analytical geostrophically-balanced initial conditions in ROMS (or any other

C-gridded model):

1. Define temporal T,S (or density field) and free surface at horizontal

*vorticity points*andplace them into scratch arrays.

2. Compute U and V by vertically integrating density computed from these temporal field

(still defined at vorticity points) and differencing the resultant pressure field (naturally

differencing over 1 dx in each direction).

3. Horizontally average T,S and vorticity fields defined in (1) using 4-point averaging, so

the outcome in at RHO-points. Accept them as the final fields.

### Re: Set geostrphically balanced initial current

shchepet wrote:There is one more thing to add to this topic specifically related to the fact that you may

want not just geostrophic balance, butthe geostrophic balance as it would be felt by the, i.e., ROMS.

discrete model

ROMS uses C-grid, which means that discreteus andvs in Coriolis terms are

subject to 4-point horizontal averaging, so the discrete geostrophic balance exists in a bit

awkward sense: <4-point-averaged u> = -(1/f) * (1/rho0) * dp/dy

and <4-point-averaged v> = +(1/f) * (1/rho0) * dp/dx where p is defined at rho-points.

The problem is that imagine that p contains checker-board mode, so dp/dy and dp/dx would

contain it as well because the differencing is only over 1 dx and 1 dy in each direction.

However 4-pint averaged velocity components cannot admit it, so there is no way to balance.

Thank you! really helpful

To address this one have to ensure that checkerboard mode should be excluded entirely from

the pressure and tracer fields defined at RHO points. This leads to a rather simple recipe

how to generate analytical geostrophically-balanced initial conditions in ROMS (or any other

C-gridded model):

1. Define temporal T,S (or density field) and free surface at horizontalvorticity pointsand

place them into scratch arrays.

2. Compute U and V by vertically integrating density computed from these temporal field

(still defined at vorticity points) and differencing the resultant pressure field (naturally

differencing over 1 dx in each direction).

3. Horizontally average T,S and vorticity fields defined in (1) using 4-point averaging, so

the outcome in at RHO-points. Accept them as the final fields.