Changeset 10788 for NEMO/trunk/src/TOP/PISCES/P4Z/p4zsed.F90

Timestamp:

2019-03-21T12:15:14+01:00 (5 years ago)

Author:

cetlod

Message:

trunk: Some bugfixes at ice/biogeochmeical interfaces, see ticket #2262

File:

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

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

@@ -53 +53 @@
       INTEGER  ::  ji, jj, jk, ikt
       REAL(wp) ::  zrivalk, zrivsil, zrivno3
-      REAL(wp) ::  zwflux, zfminus, zfplus
-      REAL(wp) ::  zlim, zfact, zfactcal
+      REAL(wp) ::  zwflux, zlim, zfact, zfactcal
       REAL(wp) ::  zo2, zno3, zflx, zpdenit, z1pdenit, zolimit
       REAL(wp) ::  zsiloss, zcaloss, zws3, zws4, zwsc, zdep
@@ -104 +103 @@
                zdep    = rfact2 / e3t_n(ji,jj,1)
                zwflux  = fmmflx(ji,jj) / 1000._wp
-               zfminus = MIN( 0._wp, -zwflux ) * trb(ji,jj,1,jpfer) * zdep
-               zfplus  = MAX( 0._wp, -zwflux ) * icefeinput * zdep
-               zironice(ji,jj) =  zfplus + zfminus
+               zironice(ji,jj) =  MAX( -0.99 * trb(ji,jj,1,jpfer), -zwflux * icefeinput * zdep )
             END DO
          END DO