advection scheme for MOMENTUM equations

Message

Barbara · #1 Unread post by **Barbara** » Sat Sep 22, 2007 2:56 pm

Hi all,

In roms2.2 the default advection for momentum equation is third-order upstream bias, but now the default horizontal and vertical advection is 4th-order centered.

The third-order upstream bias is no longer used in ROMS? Why?

Can we define UV_U3HADVECTION in /ROMS/Include/*.h?

Thanks!

jcwarner · #2 Unread post by **jcwarner** » Sat Sep 22, 2007 5:36 pm

the default horiz advection is still the 3rd order scheme. You need to look carefully in rhs3d.F.
try:
cd ROMS/Nonlinear
grep -nir '\#' rhs3d.F
and this will give you a list of the cpp options.

Barbara · #3 Unread post by **Barbara** » Sun Sep 23, 2007 4:31 pm

Thanks a lot!

I find it.

Then I think the following parts in the cppdefs.h in ~/ROMS/Include/ should be modified, so that the users won't get confused.

Code: Select all

** The following is short description of all available CPP options.          **
**                                                                           **
** OPTIONS associated with momentum equations:                               **
**                                                                           **
**   The default horizontal and vertical advection is 4th-order centered.    **
**   Use the splines vertical advection option is shallow, high vertical     **
**   resolution applications.                                                **
**                                                                           **
** UV_ADV              use to turn ON or OFF advection terms                 **
** UV_COR              use to turn ON or OFF Coriolis term                   **
** UV_C2ADVECTION      use to turn ON or OFF 2nd-order centered advection    **
** UV_C4ADVECTION      use to turn ON or OFF 4th-order centered advection    **
** UV_SADVECTION       use to turn ON or OFF splines vertical advection      **