Changeset 7617 for branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90

Timestamp:

2017-01-27T18:17:01+01:00 (7 years ago)

Author:

aumont

Message:

update diagnostics + changes in quota code

File:

• branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90 (modified) (5 diffs)

branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90

@@ -89 +89 @@
       REAL(wp), POINTER, DIMENSION(:,:,:) :: zprorcan, zprorcad, zprofed, zprofen, zpronewn, zpronewd
       REAL(wp), POINTER, DIMENSION(:,:,:) :: zmxl_fac, zmxl_chl
+#if defined key_ligand
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: zpligprod, zpligprod2
+#endif
       !!---------------------------------------------------------------------
       !
@@ -98 +101 @@
       CALL wrk_alloc( jpi, jpj, jpk, zmxl_fac, zmxl_chl )
       CALL wrk_alloc( jpi, jpj, jpk, zprorcan, zprorcad, zprofed, zprofen, zpronewn, zpronewd )
+#if defined key_ligand
+      CALL wrk_alloc( jpi, jpj, jpk, zpligprod, zpligprod2 )
+#endif
       !
       zprorcan(:,:,:) = 0._wp ; zprorcad(:,:,:) = 0._wp ; zprofed (:,:,:) = 0._wp
@@ -345 +351 @@
 #if defined key_ligand
                  tra(ji,jj,jk,jplgw) = tra(ji,jj,jk,jplgw) + zdocprod * ldocp  - zfeup * plig(ji,jj,jk) * lthet
+                 zpligprod(ji,jj,jk) = zdocprod * ldocp
+                 zpligprod2(ji,jj,jk) = zfeup * plig(ji,jj,jk) * lthet
 #endif
                  tra(ji,jj,jk,jpoxy) = tra(ji,jj,jk,jpoxy) + o2ut * ( zproreg + zproreg2) &
@@ -424 +432 @@
               CALL iom_put( "TPBFE"  , zw3d )
           ENDIF
+#if defined key_ligand
+          IF( iom_use( "LPRODP" ) )  THEN
+              zw3d(:,:,:) = zpligprod(:,:,:) * 1e9 * zfact * tmask(:,:,:)
+              CALL iom_put( "LPRODP"  , zw3d )
+          ENDIF
+          IF( iom_use( "LDETP" ) )  THEN
+              zw3d(:,:,:) = zpligprod2(:,:,:) * 1e9 * zfact * tmask(:,:,:)
+              CALL iom_put( "LDETP"  , zw3d )
+          ENDIF
+#endif
           IF( iom_use( "INTPPPHY" ) .OR. iom_use( "INTPPPHY2" ) ) THEN  
              zw2d(:,:) = 0.
@@ -495 +513 @@
      CALL wrk_dealloc( jpi, jpj, jpk, zmxl_fac, zmxl_chl )
      CALL wrk_dealloc( jpi, jpj, jpk, zprorcan, zprorcad, zprofed, zprofen, zpronewn, zpronewd )
+#if defined key_ligand
+     CALL wrk_dealloc( jpi, jpj, jpk, zpligprod, zpligprod2 )
+#endif
      !
      IF( nn_timing == 1 )  CALL timing_stop('p4z_prod')