Changeset 2528 for trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/trcrst_pisces.F90

Timestamp:

2010-12-27T18:33:53+01:00 (13 years ago)

Author:

rblod

Message:

Update NEMOGCM from branch nemo_v3_3_beta

File:

• trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/trcrst_pisces.F90 (modified) (3 diffs, 1 prop)

trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/trcrst_pisces.F90

@@ -21 +21 @@
    USE iom
    USE trcdta
+   USE lib_mpp
+   USE lib_fortran
 
    IMPLICIT NONE
@@ -118 +120 @@
       IF(lwp)  WRITE(numout,*)
 
-      IF( cp_cfg == "orca" .AND. .NOT. lk_trc_c1d ) THEN      ! ORCA condiguration (not 1D) !
+      IF( cp_cfg == "orca" .AND. .NOT. lk_c1d ) THEN      ! ORCA condiguration (not 1D) !
          !                                                    ! --------------------------- !
          ! set total alkalinity, phosphate, nitrate & silicate
 
-         zalksum = 0.e0
-         zpo4sum = 0.e0
-         zno3sum = 0.e0
-         zsilsum = 0.e0
-         DO jk = 1, jpk
-            DO jj = 1, jpj
-               DO ji = 1, jpi
-                  zvol = cvol(ji,jj,jk)
-#  if defined key_off_degrad
-                  zvol = zvol * facvol(ji,jj,jk)
-#  endif
-                  zalksum = zalksum + trn(ji,jj,jk,jptal) * zvol
-                  zpo4sum = zpo4sum + trn(ji,jj,jk,jppo4) * zvol
-                  zno3sum = zno3sum + trn(ji,jj,jk,jpno3) * zvol
-                  zsilsum = zsilsum + trn(ji,jj,jk,jpsil) * zvol
-               END DO
-            END DO
-         END DO
-         IF( lk_mpp )   CALL mpp_sum( zalksum )     ! sum over the global domain
-         IF( lk_mpp )   CALL mpp_sum( zpo4sum )     ! sum over the global domain
-         IF( lk_mpp )   CALL mpp_sum( zno3sum )     ! sum over the global domain
-         IF( lk_mpp )   CALL mpp_sum( zsilsum )     ! sum over the global domain
          zarea   = 1. / areatot * 1.e6
-         zalksum = zalksum * zarea
-         zpo4sum = zpo4sum * zarea / 122.
-         zno3sum = zno3sum * zarea / 7.6
-         zsilsum = zsilsum * zarea
+# if defined key_degrad
+         zalksum = glob_sum( trn(:,:,:,jptal) * cvol(:,:,:) * facvol(:,:,:) ) * zarea
+         zpo4sum = glob_sum( trn(:,:,:,jppo4) * cvol(:,:,:) * facvol(:,:,:) ) * zarea / 122.
+         zno3sum = glob_sum( trn(:,:,:,jpno3) * cvol(:,:,:) * facvol(:,:,:) ) * zarea / 7.6
+         zsilsum = glob_sum( trn(:,:,:,jpsil) * cvol(:,:,:) * facvol(:,:,:) ) * zarea
+# else
+         zalksum = glob_sum( trn(:,:,:,jptal) * cvol(:,:,:)  ) * zarea
+         zpo4sum = glob_sum( trn(:,:,:,jppo4) * cvol(:,:,:)  ) * zarea / 122.
+         zno3sum = glob_sum( trn(:,:,:,jpno3) * cvol(:,:,:)  ) * zarea / 7.6
+         zsilsum = glob_sum( trn(:,:,:,jpsil) * cvol(:,:,:)  ) * zarea
+# endif
 
          IF(lwp) WRITE(numout,*) '       TALK mean : ', zalksum
@@ -263 +250 @@
 #if defined key_dtatrc
       ! Restore close seas values to initial data
-      CALL trc_dta( nittrc000 ) 
+      CALL trc_dta( nit000 ) 
       DO jn = 1, jptra
          IF( lutini(jn) ) THEN