Minor bugfix in p4zrem

src/TOP/PISCES/P4Z/p4zrem.F90

@@ -73,7 +73,7 @@ CONTAINS
       INTEGER  ::   ji, jj, jk
       REAL(wp) ::   zremik, zremikc, zremikn, zremikp, zsiremin, zfact 
       REAL(wp) ::   zsatur, zsatur2, znusil, znusil2, zdep, zdepmin, zfactdep
-      REAL(wp) ::   zbactfer, zonitr, zrfact2
+      REAL(wp) ::   zbactfer, zonitr
       REAL(wp) ::   zammonic, zoxyremc, zosil, ztem, zdenitnh4, zolimic
       CHARACTER (len=25) :: charout
       REAL(wp), DIMENSION(A2D(0),jpk) :: zdepbac, zolimi, zfacsi, zfacsib, zdepeff, zfebact
@@ -240,7 +240,6 @@ CONTAINS
       ENDIF
 
       IF( knt == nrdttrc ) THEN
-          zrfact2 = 1.e+3 * rfact2r  !  conversion from mol/l/kt to  mol/m3/s
           !
           ALLOCATE( zw3d(A2D(0),jpk) )  ;  zw3d(A2D(0),jpk) = 0._wp
           !
@@ -260,8 +259,8 @@ CONTAINS
           ENDIF
           !
           IF( iom_use( "DENIT" ) )  THEN ! Denitrification
-             denitr(:,:,jpk) = 0._wp  
-             CALL iom_put( "DENIT", denitr(:,:,:) * rdenit * rno3 * tmask(:,:,:) * zrfact2 )
+             zw3d(A2D(0),1:jpkm1) = denitr(A2D(0),1:jpkm1) * 1E+3 * rfact2r * rno3 * tmask(A2D(0),1:jpkm1)
+             CALL iom_put( "DENIT", zw3d )
           ENDIF
           !
           DEALLOCATE( zw3d )