Changeset 1289

Timestamp:

2009-02-03T15:18:40+01:00 (15 years ago)

Author:

cetlod

Message:

change the computation of xnegtr in PISCEs model, see ticket:318

File:

• trunk/NEMO/TOP_SRC/PISCES/p4zbio.F90 (modified) (4 diffs)

trunk/NEMO/TOP_SRC/PISCES/p4zbio.F90

@@ -59 +59 @@
       INTEGER, INTENT(in) :: kt, jnt
       INTEGER  ::  ji, jj, jk, jn
+      REAL(wp) ::  ztra
 #if defined key_kriest
       REAL(wp) ::  zcoef1, zcoef2
@@ -69 +70 @@
       !     OF PHYTOPLANKTON AND DETRITUS
 
-      xdiss(:,:,:) = 0.01
-
+      xdiss(:,:,:) = 1.
 !!gm the use of nmld should be better here?
-      DO jk = 1, jpkm1
+      DO jk = 2, jpkm1
          DO jj = 1, jpj
             DO ji = 1, jpi
-               IF( fsdepw(ji,jj,jk+1) .le. hmld(ji,jj) )   xdiss(ji,jj,jk) = 1.e0
+               IF( fsdepw(ji,jj,jk+1) > hmld(ji,jj) )   xdiss(ji,jj,jk) = 0.01
             END DO 
          END DO
@@ -92 +92 @@
       CALL p4z_meso ( kt, jnt )           ! mesozooplankton
 
-
-      !                                                      ! test if tracers concentrations fall below 0.
+      !                             ! test if tracers concentrations fall below 0.
       xnegtr(:,:,:) = 1.e0
       DO jn = jp_pcs0, jp_pcs1
@@ -99 +98 @@
             DO jj = 1, jpj
                DO ji = 1, jpi
-                  IF( ( trn(ji,jj,jk,jn) + tra(ji,jj,jk,jn) ) < 0.e0 )    xnegtr(ji,jj,jk) = 0.e0
-               END DO
+                  IF( ( trn(ji,jj,jk,jn) + tra(ji,jj,jk,jn) ) < 0.e0 ) THEN 
+                     ztra             = ABS(  ( trn(ji,jj,jk,jn) - rtrn ) &
+                                            / ( tra(ji,jj,jk,jn) + rtrn ) )
+                     xnegtr(ji,jj,jk) = MIN( xnegtr(ji,jj,jk),  ztra )
+                  ENDIF
+              END DO
             END DO
          END DO
       END DO
       !                                ! where at least 1 tracer concentration becomes negative
-      !                                ! all tracer tendancy are set to zero (i.e. trn = trb)
+      !                                ! 
       DO jn = jp_pcs0, jp_pcs1
          trn(:,:,:,jn) = trn(:,:,:,jn) + xnegtr(:,:,:) * tra(:,:,:,jn)