--- MITgcm/pkg/dic/fe_chem.F	2004/07/13 18:03:31	1.4
+++ MITgcm/pkg/dic/fe_chem.F	2005/12/12 19:07:36	1.7
@@ -1,10 +1,9 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/fe_chem.F,v 1.4 2004/07/13 18:03:31 jmc Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/fe_chem.F,v 1.7 2005/12/12 19:07:36 stephd Exp $
 C $Name:  $
 
 #include "DIC_OPTIONS.h"
 #include "GCHEM_OPTIONS.h"
 
-#ifdef ALLOW_FE
 CStartOfInterFace
       SUBROUTINE Fe_CHEM(
      I           bi,bj,iMin,iMax,jMin,jMax,
@@ -37,6 +36,7 @@
       INTEGER myIter,myThid
 CEndOfInterface
 
+#ifdef ALLOW_FE
 
       INTEGER I,J,K
       _RL  lig, FeL
@@ -55,8 +55,8 @@
 C ligand balance in surface layer
 C in surface layer
 
-       DO j=jMin,jMax
-        DO i=iMin,iMax
+       DO j=1-OLy,sNy+OLy
+        DO i=1-OLx,sNx+OLx
          DO k=1,nR
           IF(hFacC(i,j,k,bi,bj) .gt. 0.0)THEN
 
@@ -70,11 +70,15 @@
  
               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
        ENDDO    
 c
+#endif 
         RETURN
         END                      
-#endif