Changeset 733 for tags

Timestamp:

2007-11-19T16:20:47+01:00 (17 years ago)

Author:

cetlod

Message:

update of iron aggregation computation, see ticket:14

File:

• tags/nemo_v2_3/NEMO/TOP_SRC/SMS/p4zrem.F (modified) (4 diffs)

tags/nemo_v2_3/NEMO/TOP_SRC/SMS/p4zrem.F

@@ -1 +1 @@
 
-CCC $Header$ 
+CCC $Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/SMS/p4zrem.F,v 1.8 2007/10/12 09:28:41 opalod Exp $ 
 CCC  TOP 1.0 , LOCEAN-IPSL (2005) 
 C This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
@@ -50 +50 @@
       REAL zsatur,zsatur2,znusil,zdepbac(jpi,jpj,jpk)
       REAL zlamfac,zstep,fesatur(jpi,jpj,jpk)
+      REAL ztempbac(jpi,jpj)
 C
 C      Time step duration for the biology
@@ -59 +60 @@
 C      --------------------------------------------------
 C
-        DO jk=1,12
-         zdepbac(:,:,jk)=min(0.7*(trn(:,:,jk,jpzoo)+2*trn(:,:,jk,jpmes))
-     &     ,4E-6)
+        DO jk=1,jpkm1
+          DO jj = 1, jpj
+            DO ji = 1, jpi
+         IF (fsdept(ji,jj,jk).lt.120.) THEN
+         zdepbac(ji,jj,jk)=min(0.7*(trn(ji,jj,jk,jpzoo)
+     &     +2*trn(ji,jj,jk,jpmes)),4E-6)
+         ztempbac(ji,jj)=zdepbac(ji,jj,jk)
+         ELSE
+         zdepbac(ji,jj,jk)=min(1.,120./fsdept(ji,jj,jk))
+     &      *ztempbac(ji,jj)
+         ENDIF
+            END DO
+          END DO
         END DO
-C
-C      Vertical decay of the bacterial activity
-C      ----------------------------------------
-C
-         do jk=13,jpk
-           do jj=1,jpj
-             do ji=1,jpi
-         zdepbac(ji,jj,jk)=min(1.,fsdept(ji,jj,12)/fsdept(ji,jj,jk))
-     &      *zdepbac(ji,jj,12)
-             end do
-           end do
-         end do
 
          DO jk = 1,jpkm1
@@ -253 +252 @@
          zlamfac=min(1.,zlamfac)
 #if ! defined key_trc_kriest
-         xlam1b=(80.*(trn(ji,jj,jk,jpdoc)+40E-6)+698.
+         xlam1b=(80.*(trn(ji,jj,jk,jpdoc)+35E-6)+698.
      &    *trn(ji,jj,jk,jppoc)+1.05E4*trn(ji,jj,jk,jpgoc))
-     &    *zdiss(ji,jj,jk)+1E-5*(1.-zlamfac)+xlam1*max(0.,
+     &    *zdiss(ji,jj,jk)+1E-4*(1.-zlamfac)+xlam1*max(0.,
      &    (trn(ji,jj,jk,jpfer)*1E9-1.))
 #else
-         xlam1b=(80.*(trn(ji,jj,jk,jpdoc)+40E-6)+698.
+         xlam1b=(80.*(trn(ji,jj,jk,jpdoc)+35E-6)+698.
      &    *trn(ji,jj,jk,jppoc))
-     &    *zdiss(ji,jj,jk)+1E-5*(1.-zlamfac)+xlam1*max(0.,
+     &    *zdiss(ji,jj,jk)+1E-4*(1.-zlamfac)+xlam1*max(0.,
      &    (trn(ji,jj,jk,jpfer)*1E9-1.))
 #endif
 
 
-         xaggdfe(ji,jj,jk)=xlam1b*zstep*0.76*(trn(ji,jj,jk,jpfer)
+         xaggdfe(ji,jj,jk)=xlam1b*zstep*0.5*(trn(ji,jj,jk,jpfer)
      &     -xfeequi)
 #    if defined key_off_degrad