Changeset 433 for trunk/NEMO/TOP_SRC/SMS/trcbio.F

Timestamp:

2006-04-10T17:40:29+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_044 : CT : update the passive tracers TOP component and the standard GYRE configuration

File:

• trunk/NEMO/TOP_SRC/SMS/trcbio.F (modified) (5 diffs)

trunk/NEMO/TOP_SRC/SMS/trcbio.F

@@ -212 +212 @@
 C 4. zooplankton liquide excretion
 C
-            zzoonh4 = tauzn * zzoo * fdoml 
-            zzoodom = tauzn * zzoo * (1-fdoml) 
-
+            zzoonh4 = tauzn * fzoolab * zzoo 
+            zzoodom = tauzn * (1 - fzoolab) * zzoo
+C
 C 5. mortality
 C
@@ -226 +226 @@
 C
             zzoobod = tmminz * zzoo * zzoo
-            fbod(ji,jj) = fbod(ji,jj) + zzoobod * fse3t(ji,jj,jk)
+            fbod(ji,jj) = fbod(ji,jj) 
+     $                 + (1-fdbod) * zzoobod * fse3t(ji,jj,jk)
+            zboddet = fdbod * zzoobod
 C
 C
@@ -232 +234 @@
 C
 C
-            zdetnh4 = taudn * fdoml * zdet
-            zdetdom = taudn * (1 - fdoml) * zdet
+            zdetnh4 = taudn * fdetlab * zdet
+            zdetdom = taudn * (1 - fdetlab) * zdet 
+
             zdomnh4 = taudomn * zdom
 C
+C flux added to express how the excess of nitrogen from
+C PHY, ZOO and DET to DOM goes directly to NH4 (flux of ajustment)
+            zdomaju = (1 - redf/reddom) * (zphydom + zzoodom + zdetdom)
 C
 C 7. Nitrification
@@ -254 +260 @@
             zno3a = - zno3phy + znh4no3
             znh4a = - znh4phy - znh4no3 + zphynh4 + zzoonh4 + zdomnh4
-     $          + zdetnh4
-            zdeta = zphydet + zzoodet  - zdetzoo - zdetnh4 - zdetdom
-            zdoma = zphydom + zzoodom + zdetdom - zdomnh4
+     $          + zdetnh4 + zdomaju
+            zdeta = zphydet + zzoodet  - zdetzoo - zdetnh4 - zdetdom +
+     $          zboddet
+            zdoma = zphydom + zzoodom + zdetdom - zdomnh4 - zdomaju
 C
 #if defined key_trc_diabio
@@ -296 +303 @@
      $          -zzoobod-zzoonh4-zzoodom) *ze3t(jk)
             trc2d(ji,jj,18)=trc2d(ji,jj,18)+zdetdom*ze3t(jk)
-
+c trend number 19 is in trcexp.F
             trc3d(ji,jj,jk,1)= zno3phy *86400     
             trc3d(ji,jj,jk,2)= znh4phy *86400