--- MITgcm/pkg/dic/fe_chem.F 2008/04/07 20:31:16 1.13 +++ MITgcm/pkg/dic/fe_chem.F 2008/06/23 22:26:32 1.15 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/fe_chem.F,v 1.13 2008/04/07 20:31:16 dfer 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" @@ -62,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 @@ -76,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