--- MITgcm/pkg/dic/fe_chem.F 2003/10/09 04:19:19 1.3 +++ MITgcm/pkg/dic/fe_chem.F 2005/12/16 21:07:53 1.8 @@ -1,7 +1,9 @@ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/fe_chem.F,v 1.8 2005/12/16 21:07:53 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, @@ -21,6 +23,7 @@ #include "PARAMS.h" #include "GRID.h" #include "DIC_BIOTIC.h" +#include "PTRACERS_SIZE.h" #include "PTRACERS.h" C == Routine arguments == @@ -33,6 +36,7 @@ INTEGER myIter,myThid CEndOfInterface +#ifdef ALLOW_FE INTEGER I,J,K _RL lig, FeL @@ -51,8 +55,8 @@ C ligand balance in surface layer C in surface layer - DO j=jMin,jMax - DO i=iMin,iMax + DO j=jmin,jmax + DO i=imin,imax DO k=1,nR IF(hFacC(i,j,k,bi,bj) .gt. 0.0)THEN @@ -66,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