Changeset 1163

Timestamp:

2008-07-08T17:11:10+02:00 (16 years ago)

Author:

ctlod

Message:

trunk: remove extra coefficient zavdt when computing zavds scalar in zdfddm.F90, see ticket: #226

File:

• trunk/NEMO/OPA_SRC/ZDF/zdfddm.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/ZDF/zdfddm.F90

@@ -163 +163 @@
                zrr = zrr * zrr
                zavfs = avts / ( 1 + zrr*zrr*zrr ) * zmsks(ji,jj) *zmskf(ji,jj)
-               zavft = 0.7 * zavfs / rrau(ji,jj,jk)
+               zavft = 0.7 * zavfs * zinr
                ! diffusive layering
                zavdt = 1.3635e-6 * EXP(4.6*EXP(-0.54*(zinr-1.) ) )   &
@@ -169 +169 @@
                zavds = zavdt * zmsks(ji,jj)   &
                      * ( (1.85 * rrau(ji,jj,jk) - 0.85 ) * zmskd3(ji,jj)   &
-                         + zavdt * 0.15 * rrau(ji,jj,jk) * zmskd2(ji,jj)  )
+                       +  0.15 * rrau(ji,jj,jk)          * zmskd2(ji,jj)  )
                ! add to the eddy viscosity coef. previously computed
                avs (ji,jj,jk) = avt(ji,jj,jk) + zavfs + zavds