--- MITgcm/pkg/dic/fe_chem.F	2010/04/13 03:06:14	1.16
+++ MITgcm/pkg/dic/fe_chem.F	2012/08/22 00:40:56	1.17
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/fe_chem.F,v 1.16 2010/04/13 03:06:14 jmc Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/fe_chem.F,v 1.17 2012/08/22 00:40:56 jmc Exp $
 C $Name:  $
 
 #include "DIC_OPTIONS.h"
@@ -8,7 +8,8 @@
 C     !INTERFACE:
       SUBROUTINE Fe_CHEM(
      I           bi,bj, iMin,iMax,jMin,jMax,
-     I           fe, freefe,
+     I           fe,
+     O           freefe,
      I           myIter, myThid )
 
 C     !DESCRIPTION: \bv
@@ -55,6 +56,14 @@
 C
 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
 
+      DO k=1,Nr
+       DO j=1-OLy,sNy+OLy
+        DO i=1-OLx,sNx+OLx
+          freefe(i,j,k) = 0. _d 0
+        ENDDO
+       ENDDO
+      ENDDO
+
 C ligand balance in surface layer
 C in surface layer