Changeset 3446 for branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmeso.F90

Timestamp:

2012-08-10T13:13:55+02:00 (12 years ago)

Author:

cetlod

Message:

branch:2012/dev_r3438_LOCEAN15_PISLOB : 2nd step new PISCES updates from Olivier, see ticket #972

File:

• branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmeso.F90 (modified) (4 diffs)

branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmeso.F90

@@ -81 +81 @@
       REAL(wp) :: zgrazn, zgrazpoc, zgraznf, zgrazf
       REAL(wp) :: zgrazfff, zgrazffe
+      REAL(wp) :: zrfact2
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: zgrazing
       CHARACTER (len=25) :: charout
-      REAL(wp) :: zrfact2
       !!---------------------------------------------------------------------
       !
       IF( nn_timing == 1 )  CALL timing_start('p4z_meso')
       !
-
+      IF( ln_diatrc .AND. lk_iomput )  CALL wrk_alloc( jpi, jpj, jpk, zgrazing )
+      !
       DO jk = 1, jpkm1
          DO jj = 1, jpj
@@ -145 +147 @@
 
               ! Total grazing ( grazing by microzoo is already computed in p4zmicro )
-              grazing(ji,jj,jk) = grazing(ji,jj,jk) + zgraztot
+              IF( ln_diatrc .AND. lk_iomput )  zgrazing(ji,jj,jk) = zgraztot
+
               !    Mesozooplankton efficiency
               !    --------------------------
@@ -213 +216 @@
       END DO
       !
-      IF( ln_diatrc .AND. lk_iomput ) THEN
+      IF( ln_diatrc .AND. lk_iomput .AND. jnt == nrdttrc ) THEN
          zrfact2 = 1.e3 * rfact2r
-         grazing(:,:,:) = grazing(:,:,:) * zrfact2 * tmask(:,:,:)   ! Total grazing of phyto by zoo
-         prodcal(:,:,:) = prodcal(:,:,:) * zrfact2 * tmask(:,:,:)   ! Calcite production
-         IF( jnt == nrdttrc ) THEN
-            CALL iom_put( "GRAZ" , grazing  )  ! Total grazing of phyto by zooplankton
-            CALL iom_put( "PCAL" , prodcal  )  ! Calcite production
-         ENDIF
+         CALL iom_put( "GRAZ2", zgrazing(:,:,:) * zrfact2 * tmask(:,:,:) )  ! Total grazing of phyto by zooplankton
+         CALL iom_put( "PCAL" , prodcal(:,:,:)  * zrfact2 * tmask(:,:,:) )  ! Calcite production
       ENDIF
       !
@@ -228 +227 @@
         CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm)
       ENDIF
+      !
+      IF( ln_diatrc .AND. lk_iomput )  CALL wrk_dealloc( jpi, jpj, jpk, zgrazing )
       !
       IF( nn_timing == 1 )  CALL timing_stop('p4z_meso')