Changeset 3441 for trunk/NEMOGCM

Timestamp:

2012-08-07T13:12:15+02:00 (12 years ago)

Author:

beppe

Message:

Minor bug in closea when key_mpp_rep is not active. Avoid mpp_sum of the
global area (surf(jpncs+1)) cause it is already computed using glob_area.

File:

• trunk/NEMOGCM/NEMO/OPA_SRC/DOM/closea.F90 (modified) (4 diffs)

trunk/NEMOGCM/NEMO/OPA_SRC/DOM/closea.F90

@@ -193 +193 @@
       COMPLEX(wp)         ::   ctmp 
       REAL(wp), DIMENSION(jpncs) ::   zfwf   ! 1D workspace
+      REAL(wp) :: zdum
       !!----------------------------------------------------------------------
       !
@@ -228 +229 @@
                END DO 
             END DO 
-            IF( lk_mpp )   CALL mpp_sum ( surf, jpncs+1 )       ! mpp: sum over all the global domain
+            IF( lk_mpp )   CALL mpp_sum ( surf, jpncs )       ! mpp: sum over all the global domain
          ENDIF
 
@@ -243 +244 @@
       ENDIF
       !                                                   !--------------------!
+      zdum = glob_sum(emp(:,:)*e1e2t(:,:))
+      WRITE(numout,*) 'kt=', kt, ' EMP BEFORE=', zdum
       !                                                   !  update emp, emps  !
       zfwf = 0.e0_wp                                      !--------------------!
@@ -343 +346 @@
       CALL lbc_lnk( emps, 'T', 1._wp )
       !
+      zdum = glob_sum(emp(:,:)*e1e2t(:,:))
+      WRITE(numout,*) 'kt=', kt, ' EMP BEFORE=', zdum
       IF( nn_timing == 1 )  CALL timing_stop('sbc_clo')
       !