Changeset 655 for trunk/NEMO/TOP_SRC

Timestamp:

2007-04-18T11:34:10+02:00 (17 years ago)

Author:

opalod

Message:

nemo_v2_update_016:CE:Tuning of the iron concentration to a minimum level that is set to the detection limit

File:

• trunk/NEMO/TOP_SRC/SMS/p4zlim.F (modified) (2 diffs)

trunk/NEMO/TOP_SRC/SMS/p4zlim.F

@@ -45 +45 @@
 CC ==================
       INTEGER ji, jj, jk
-      REAL xlim1,xlim2,xlim3,xlim4
+      REAL xlim1,xlim2,xlim3,xlim4,zno3,zferlim
       REAL xconctemp,xconctemp2,xconctempn,xconctempn2
 C
@@ -52 +52 @@
 C  -------------------------------------
 C
-
-        trn(:,:,:,jpfer)=max(trn(:,:,:,jpfer),1.E-11)
+        do jk=1,jpk
+          do jj=1,jpj
+            do ji=1,jpi
+        zno3=trn(ji,jj,jk,jpno3)*1E6
+        zferlim=max(1.5E-11*(zno3/40)**2,3E-12)
+        zferlim=min(zferlim,1.5E-11)
+        trn(ji,jj,jk,jpfer)=max(trn(ji,jj,jk,jpfer),zferlim)
+            end do
+          end do
+        end do
 C
 C  Computation of a variable Ks for iron on diatoms