Changeset 1082 for trunk/NEMO/OPA_SRC/ZDF/zdfkpp.F90

Timestamp:

2008-06-05T16:04:42+02:00 (16 years ago)

Author:

rblod

Message:

Bugs correction in KPP, see tickets #69 and #101

File:

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

trunk/NEMO/OPA_SRC/ZDF/zdfkpp.F90

@@ -350 +350 @@
                   ze3sqr = 1. / ( fse3w(ji,jj,jk) * fse3w(ji,jj,jk) )
                   ! Square of vertical shear  at interfaces
-                  zsh2   = 0.5 * 0.5 * ( zdku2 + zdkv2 ) * ze3sqr
+                  zsh2   = 0.5 * ( zdku2 + zdkv2 ) * ze3sqr
                   zri    = MAX( rn2(ji,jj,jk), 0. ) / ( zsh2 + epsln ) 
 #if defined key_c1d
@@ -461 +461 @@
             zBosol(ji,jj) = grav * zthermal * qsr(ji,jj)
             ! Non radiative surface buoyancy force
-            zBo   (ji,jj) = grav * zthermal * qns(ji,jj) -  grav * zhalin * emp(ji,jj)
+            zBo   (ji,jj) = grav * zthermal * qns(ji,jj) -  grav * zhalin * emps(ji,jj)
             ! Surface Temperature flux for non-local term
             wt0(ji,jj) = - ( qsr(ji,jj) + qns(ji,jj) )* ro0cpr * tmask(ji,jj,1)
             ! Surface salinity flux for non-local term
-            ws0(ji,jj) = - ( emp(ji,jj) * sn(ji,jj,1) * rcs ) * tmask(ji,jj,1)
+            ws0(ji,jj) = - ( emps(ji,jj) * sn(ji,jj,1) * rcs ) * tmask(ji,jj,1)
          ENDDO
       ENDDO