Changeset 10222 for NEMO/trunk/src/TOP/PISCES/P4Z/p4zrem.F90

Timestamp:

2018-10-25T11:42:23+02:00 (6 years ago)

Author:

aumont

Message:

update of the sediment module + enhancement, bug correction in PISCES generating O2 negative values

File:

• NEMO/trunk/src/TOP/PISCES/P4Z/p4zrem.F90 (modified) (2 diffs)

NEMO/trunk/src/TOP/PISCES/P4Z/p4zrem.F90

@@ -116 +116 @@
                   zammonic = zremik * nitrfac(ji,jj,jk) * trb(ji,jj,jk,jpdoc)
                   denitr(ji,jj,jk)  = zammonic * ( 1. - nitrfac2(ji,jj,jk) )
-                  zoxyrem           = zammonic *        nitrfac2(ji,jj,jk)
+                  denitr(ji,jj,jk)  = MIN( ( trb(ji,jj,jk,jpno3) - rtrn ) / rdenit, denitr(ji,jj,jk) )
+                  zoxyrem           = zammonic - denitr(ji,jj,jk)
                   !
                   zolimi (ji,jj,jk) = MAX( 0.e0, zolimi (ji,jj,jk) )
@@ -189 +190 @@
                &         / ( 1.+ emoy(ji,jj,jk) ) * ( 1. + fr_i(ji,jj) * emoy(ji,jj,jk) ) 
                zdenitnh4 = nitrif * xstep * trb(ji,jj,jk,jpnh4) * nitrfac(ji,jj,jk)
+               zdenitnh4 = MIN(  ( trb(ji,jj,jk,jpno3) - rtrn ) / rdenita, zdenitnh4 ) 
                ! Update of the tracers trends
                ! ----------------------------