Changeset 14786 for NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zrem.F90

Timestamp:

2021-05-05T09:09:10+02:00 (3 years ago)

Author:

aumont

Message:

Various bug fixes in PISCES

File:

• NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zrem.F90 (modified) (3 diffs)

NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zrem.F90

@@ -97 +97 @@
          zdep = MAX( hmld(ji,jj), heup_01(ji,jj) )
          zdepbac(ji,jj,jk) = 0.6 * ( MAX(0.0, tr(ji,jj,jk,jpzoo,Kbb) + tr(ji,jj,jk,jpmes,Kbb) ) * 1.0E6 )**0.6 * 1.E-6
-         IF( gdept(ji,jj,jk,Kmm) < zdep ) THEN
+         IF( gdept(ji,jj,jk,Kmm) >= zdep ) THEN
             zdepmin = MIN( 1., zdep / gdept(ji,jj,jk,Kmm) )
             zdepeff(ji,jj,jk) = zdepeff(ji,jj,jk) * zdepmin**0.3
@@ -114 +114 @@
          ! -----------------------------------------------------
          zolimic = zremikc * ( 1.- nitrfac(ji,jj,jk) ) * tr(ji,jj,jk,jpdoc,Kbb) 
-         zolimic = MIN( ( tr(ji,jj,jk,jpoxy,Kbb) - rtrn ) / o2ut, zolimic ) 
+         zolimic = MAX(0., MIN( ( tr(ji,jj,jk,jpoxy,Kbb) - rtrn ) / o2ut, zolimic ) ) 
          zolimi(ji,jj,jk) = zolimic
 
@@ -155 +155 @@
          &         / ( 1.+ emoy(ji,jj,jk) ) * ( 1. + fr_i(ji,jj) * emoy(ji,jj,jk) ) 
          zdenitnh4 = nitrif * xstep * tr(ji,jj,jk,jpnh4,Kbb) * nitrfac(ji,jj,jk)
-         zdenitnh4 = MIN(  ( tr(ji,jj,jk,jpno3,Kbb) - rtrn ) / rdenita, zdenitnh4 ) 
+         zdenitnh4 = MAX(0., MIN(  ( tr(ji,jj,jk,jpno3,Kbb) - rtrn ) / rdenita, zdenitnh4 ) )
          ! Update of the tracers trends
          ! ----------------------------