Opened 12 years ago

Closed 10 years ago

#233 closed Bug (fixed)

bottom friction issue...

Reported by: jchanut Owned by: nemo
Priority: low Milestone:
Component: OCE Version: release-2
Severity: Keywords:
Cc:

Description

It seems that there is a conflict between bottom friction and the way the surface pressure gradient is computed in NEMO. This is the case whatever the method used for the free surface (explicit, filtered or time split) and as long an implicit scheme is used for vertical diffusion (ln_zdfexp=.FALSE.).

In this case indeed, the bottom stress is computed implicitely through the dyn_zdf routine but before the call to dyn_spg so that the contribution of the pressure gradient is missing in the evaluation of the "next" bottom velocities. As a consequence the tendency due to the bottom stress is wrong.

In the explicit free surface case, moving the call to dyn_zdf after dyn_spg would solve the problem.

In other cases, it is more difficult to fix it, except if one use an explicit form of the bottom friction (at least for the barotropic contribution). It is likely that this would be much less stable (by doing this, my model blows up, but I probably made a mistake). If one would like to keep an implicit form of bottom friction, a possible conflict with the implicit computation of the filtering force will also arise…

In the time splitting case, a solution could be to explicitely add the bottom stress in the barotropic equations.

Jerome

Commit History (2)

ChangesetAuthorTimeChangeLog
1671rblod2009-10-16T12:09:46+02:00

Add missing pieces related to previous commit on bottom friction …, see ticket #233

1662rblod2009-10-14T18:51:06+02:00

Correction of major bug on bottom friction, see ticket #233

Change History (2)

comment:1 Changed 11 years ago by rblod

At this time we add an explicit computation with control on the coefficient


Update of the documentation follows
Note: tke_old had not been updated

Rachid

comment:2 Changed 10 years ago by gm

  • Resolution set to fixed
  • Status changed from new to closed
Note: See TracTickets for help on using tickets.