--- MITgcm/pkg/dic/fe_chem.F	2004/10/18 16:01:13	1.5
+++ MITgcm/pkg/dic/fe_chem.F	2005/12/16 21:07:53	1.8
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/fe_chem.F,v 1.5 2004/10/18 16:01:13 jmc Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/fe_chem.F,v 1.8 2005/12/16 21:07:53 stephd Exp $
 C $Name:  $
 
 #include "DIC_OPTIONS.h"
@@ -55,8 +55,8 @@
 C ligand balance in surface layer
 C in surface layer
 
-       DO j=jMin,jMax
-        DO i=iMin,iMax
+       DO j=jmin,jmax
+        DO i=imin,imax
          DO k=1,nR
           IF(hFacC(i,j,k,bi,bj) .gt. 0.0)THEN
 
@@ -70,6 +70,10 @@
  
               FeL = ligand_tot-lig
               freefe(i,j,k,bi,bj) = fe (i,j,k,bi,bj)-FeL
+#ifdef MINFE
+              freefe(i,j,k,bi,bj) = min(freefe(i,j,k,bi,bj),freefemax)
+              fe(i,j,k,bi,bj) = FeL+freefe(i,j,k,bi,bj)
+#endif
           END IF  
          ENDDO
         ENDDO