Changeset 1521

Timestamp:

2009-07-22T12:48:06+02:00 (15 years ago)

Author:

ctlod

Message:

wrong scale factors, see ticket: #365

File:

• trunk/NEMO/OPA_SRC/BDY/bdyvol.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/BDY/bdyvol.F90

@@ -73 +73 @@
       INTEGER  ::   ji, jj, jk, jb, jgrd
       INTEGER  ::   ii, ij
-      REAL(wp) ::   zubtpecor, z_cflxemp, ztranst
+      REAL(wp) ::   zubtpecor, z_cflxemp, ztranst, zraur
       !!-----------------------------------------------------------------------------
 
@@ -85 +85 @@
       ! -----------------------------------------------------------------------
       z_cflxemp = 0.e0
-      DO jj = 1, jpj
-         DO ji = 1, jpi
-            z_cflxemp = z_cflxemp + emp(ji,jj) * bdytmask(ji,jj) * tmask_i(ji,jj) / rauw * e1v(ji,jj) * e2u(ji,jj)
-         END DO
-      END DO
+      zraur = 1.e0 / rauw
+      z_cflxemp = SUM ( emp(:,:) * bdytmask(:,:) * e1t(:,:) * e2t(:,:) * zraur )
       IF( lk_mpp )   CALL mpp_sum( z_cflxemp )   ! sum over the global domain