Changeset 3946 for branches/2012

Timestamp:

2013-06-28T11:22:56+02:00 (11 years ago)

Author:

smasson

Message:

NEMOGCM_v3_4_STABLE: major bugfix in UBS, see #118

File:

• branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv_ubs.F90 (modified) (3 diffs)

branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv_ubs.F90

@@ -28 +28 @@
    PRIVATE
 
-   REAL(wp), PARAMETER :: gamma1 = 1._wp/3._wp  ! =1/4 quick      ; =1/3  3rd order UBS
-   REAL(wp), PARAMETER :: gamma2 = 1._wp/8._wp  ! =0   2nd order  ; =1/8  4th order centred
+   REAL(wp), PARAMETER :: gamma1 = 1._wp/3._wp   ! =1/4 quick      ; =1/3  3rd order UBS
+   REAL(wp), PARAMETER :: gamma2 = 1._wp/32._wp  ! =0   2nd order  ; =1/32 4th order centred
 
    PUBLIC   dyn_adv_ubs   ! routine called by step.F90
@@ -57 +57 @@
       !!                       = 1/3  3rd order Upstream biased scheme
       !!                gamma2 = 0    2nd order finite differencing 
-      !!                       = 1/8  4th order finite differencing
+      !!                       = 1/32 4th order finite differencing
       !!      For stability reasons, the first term of the fluxes which cor-
       !!      responds to a second order centered scheme is evaluated using  
@@ -64 +64 @@
       !!      before velocity (forward in time). 
       !!      Default value (hard coded in the begining of the module) are 
-      !!      gamma1=1/3 and gamma2=1/8.
+      !!      gamma1=1/3 and gamma2=1/32.
       !!
       !! ** Action : - (ua,va) updated with the 3D advective momentum trends