Changeset 1830 for branches/CMIP5_IPSL/NEMO/TOP_SRC/PISCES/p4zmeso.F90

Timestamp:

2010-04-12T15:03:51+02:00 (14 years ago)

Author:

cetlod

Message:

Computation of additional diagnostics for PISCES model ( under CPP key key_diaar5 )

• needed for AR5 outputs (vertical inventories, passive tracers at surface,... )
• new output file with suffix dbio_T

File:

• branches/CMIP5_IPSL/NEMO/TOP_SRC/PISCES/p4zmeso.F90 (modified) (3 diffs)

branches/CMIP5_IPSL/NEMO/TOP_SRC/PISCES/p4zmeso.F90

@@ -76 +76 @@
 #if defined key_trc_diaadd && defined key_trc_dia3d && defined key_iomput
       REAL(wp) :: zrfact2
-      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zw3d
+#if defined key_diaar5
+      REAL(wp), DIMENSION(jpi,jpj) ::   zpcalint
+#endif
 #endif
 
@@ -203 +205 @@
       END DO
       
+#if defined key_trc_dia3d
+      ! Total grazing ( grazing by microzoo is already computed in p4zmicro ) 
+      grazing(:,:,:) = grazing(:,:,:) + (  zgrazd  (:,:,:) + zgrazz  (:,:,:) + zgrazn(:,:,:) &
+                     &                   + zgrazpoc(:,:,:) + zgrazffe(:,:,:)  )
+#endif
+
 
       DO jk = 1,jpkm1
@@ -311 +319 @@
 #if defined key_trc_diaadd && defined key_trc_dia3d && defined key_iomput
       zrfact2 = 1.e3 * rfact2r
-      zw3d(:,:,:) = (     zgrazd(:,:,:) +   zgrazz(:,:,:) + zgrazn(:,:,:) &
-                    & + zgrazpoc(:,:,:) + zgrazffe(:,:,:)                 ) * zrfact2 * tmask(:,:,:)
-      IF( jnt == nrdttrc ) CALL iom_put( "Graz2" , zw3d )
-
-      zw3d(:,:,:) = prodcal(:,:,:) * zrfact2 * tmask(:,:,:)
-      IF( jnt == nrdttrc ) CALL iom_put( "Pcal"  , zw3d )
+      ! Total grazing of phyto by zoo
+      grazing(:,:,:) = grazing(:,:,:) * zrfact2 * tmask(:,:,:)
+      ! Calcite production
+      prodcal(:,:,:) = prodcal(:,:,:) * zrfact2 * tmask(:,:,:)
+      IF( jnt == nrdttrc ) then 
+         CALL iom_put( "GRAZ" , grazing  )  ! Total grazing of phyto by zooplankton
+         CALL iom_put( "PCAL" , prodcal  )  ! Calcite production
+      ENDIF
+#if defined key_diaar5
+      ! Vertically integrated calcite production
+      zpcalint(:,:) = 0.
+      DO jk = 1, jpkm1
+         zpcalint(:,:) = zpcalint(:,:) + prodcal(:,:,jk) * fse3t(:,:,jk)
+      ENDDO
+      IF( jnt == nrdttrc ) CALL iom_put( "INTPCAL", zpcalint )  ! Vertically integrated calcite production
+#endif
 #endif