--- MITgcm/pkg/dic/fe_chem.F	2007/12/06 22:24:34	1.11
+++ MITgcm/pkg/dic/fe_chem.F	2008/06/23 22:26:32	1.15
@@ -1,8 +1,7 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/fe_chem.F,v 1.11 2007/12/06 22:24:34 stephd Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/fe_chem.F,v 1.15 2008/06/23 22:26:32 stephd Exp $
 C $Name:  $
 
 #include "DIC_OPTIONS.h"
-#include "GCHEM_OPTIONS.h"
 
 CStartOfInterFace
       SUBROUTINE Fe_CHEM(
@@ -22,7 +21,7 @@
 #include "EEPARAMS.h"
 #include "PARAMS.h"
 #include "GRID.h"
-#include "DIC_BIOTIC.h"
+#include "DIC_VARS.h"
 
 C     == Routine arguments ==                
 C     bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
@@ -63,7 +62,7 @@
           IF (hFacC(i,j,k,bi,bj) .GT. 0. _d 0) THEN
 
 #ifdef DIC_NO_NEG
-              tmpfe=max(0. _d0 , fe (i,j,k,bi,bj))
+              tmpfe=max(0. _d 0 , fe (i,j,k,bi,bj))
 #else
               tmpfe=fe (i,j,k,bi,bj)
 #endif
@@ -77,7 +76,7 @@
      &            )/(2. _d 0*ligand_stab)
  
               FeL = ligand_tot-lig
-              if (tmpfe.eq.0. _d 0) then
+              if (tmpfe.ne.0. _d 0) then
                 freefe(i,j,k,bi,bj) = tmpfe -FeL
               else
                 freefe(i,j,k,bi,bj) = 0. _d 0