--- MITgcm/pkg/dic/fe_chem.F	2003/10/02 18:30:07	1.1.2.1
+++ MITgcm/pkg/dic/fe_chem.F	2005/10/13 16:25:12	1.6
@@ -1,3 +1,6 @@
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/fe_chem.F,v 1.6 2005/10/13 16:25:12 stephd Exp $
+C $Name:  $
+
 #include "DIC_OPTIONS.h"
 #include "GCHEM_OPTIONS.h"
 
@@ -20,6 +23,7 @@
 #include "PARAMS.h"
 #include "GRID.h"
 #include "DIC_BIOTIC.h"
+#include "PTRACERS_SIZE.h"
 #include "PTRACERS.h"
 
 C     == Routine arguments ==                
@@ -32,6 +36,7 @@
       INTEGER myIter,myThid
 CEndOfInterface
 
+#ifdef ALLOW_FE
 
       INTEGER I,J,K
       _RL  lig, FeL
@@ -65,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                      
-