--- MITgcm/pkg/dic/fe_chem.F 2008/04/18 19:32:12 1.14 +++ MITgcm/pkg/dic/fe_chem.F 2010/04/13 03:06:14 1.16 @@ -1,105 +1,107 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/fe_chem.F,v 1.14 2008/04/18 19:32:12 stephd Exp $ +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 $Name: $ #include "DIC_OPTIONS.h" -CStartOfInterFace +CBOP +C !ROUTINE: Fe_CHEM +C !INTERFACE: SUBROUTINE Fe_CHEM( - I bi,bj,iMin,iMax,jMin,jMax, + I bi,bj, iMin,iMax,jMin,jMax, I fe, freefe, I myIter, myThid ) -C /==========================================================\ -C | SUBROUTINE Fe_chem | -C | | -C | o Calculate L,FeL,Fe concentration | -C |==========================================================| + +C !DESCRIPTION: \bv +C *==========================================================* +C | SUBROUTINE Fe_CHEM +C | o Calculate L,FeL,Fe concentration +C *==========================================================* +C \ev + +C !USES: IMPLICIT NONE C == GLobal variables == #include "SIZE.h" -#include "DYNVARS.h" #include "EEPARAMS.h" #include "PARAMS.h" #include "GRID.h" #include "DIC_VARS.h" -C == Routine arguments == -C bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation -C results will be set. -C myThid - Instance number for this innvocation of CALC_GT - _RL freefe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) - _RL fe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) - INTEGER bi,bj,iMin,iMax,jMin,jMax - INTEGER myIter,myThid -#ifdef AD_SAFE - _RL thx, thy, theps -#endif -CEndOfInterface +C !INPUT/OUTPUT PARAMETERS: +C == Routine arguments == +C bi, bj :: current tile indices +C iMin,iMax,jMin,jMax :: Range of points for which calculation is performed. +C myThid :: my Thread Id number + INTEGER bi,bj + INTEGER iMin,iMax,jMin,jMax + _RL fe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr) + _RL freefe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr) + INTEGER myIter, myThid +CEOP #ifdef ALLOW_FE - - INTEGER I,J,K +C !LOCAL VARIABLES: + INTEGER i,j,k _RL lig, FeL _RL tmpfe +#ifdef AD_SAFE + _RL thx, thy, theps +#endif -CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc -CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc -CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc -CC -CC ADAPTED FROM PAYAL -CC -CC -CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc -CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc -CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc +CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC +C +C ADAPTED FROM PAYAL +C +CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C ligand balance in surface layer C in surface layer - DO j=jmin,jmax - DO i=imin,imax - DO k=1,Nr - IF (hFacC(i,j,k,bi,bj) .GT. 0. _d 0) THEN + DO k=1,Nr + DO j=jMin,jMax + DO i=iMin,iMax + IF ( maskC(i,j,k,bi,bj).GT.0. ) THEN #ifdef DIC_NO_NEG - tmpfe=max(0. _d 0 , fe (i,j,k,bi,bj)) + tmpfe =MAX( 0. _d 0 , fe(i,j,k) ) #else - tmpfe=fe (i,j,k,bi,bj) + tmpfe = fe(i,j,k) #endif - + C Ligand,FeL,Fe calculation lig=(-ligand_stab*tmpfe + & ligand_stab*ligand_tot-1. _d 0 - & +((ligand_stab*tmpfe + & +((ligand_stab*tmpfe & -ligand_stab*ligand_tot+1. _d 0)**2 & +4. _d 0*ligand_stab*ligand_tot)**0.5 _d 0 & )/(2. _d 0*ligand_stab) - + FeL = ligand_tot-lig - if (tmpfe.eq.0. _d 0) then - freefe(i,j,k,bi,bj) = tmpfe -FeL - else - freefe(i,j,k,bi,bj) = 0. _d 0 - endif + IF (tmpfe.NE.0. _d 0) THEN + freefe(i,j,k) = tmpfe -FeL + ELSE + freefe(i,j,k) = 0. _d 0 + ENDIF #ifdef MINFE #ifdef AD_SAFE - thx=freefe(i,j,k,bi,bj) + thx=freefe(i,j,k) thy=freefemax theps=1. _d -8 - freefe(i,j,k,bi,bj) = + freefe(i,j,k) = & ( 1. _d 0 - tanh((thx-thy)/theps) ) * thx/2.+ & ( 1. _d 0 + tanh((thx-thy)/theps) ) * thy/2. #else - freefe(i,j,k,bi,bj) = min(freefe(i,j,k,bi,bj),freefemax) + freefe(i,j,k) = MIN(freefe(i,j,k),freefemax) #endif - fe(i,j,k,bi,bj) = FeL+freefe(i,j,k,bi,bj) + fe(i,j,k) = FeL+freefe(i,j,k) #endif - END IF + ENDIF ENDDO ENDDO - ENDDO -c -#endif - RETURN - END + ENDDO + +#endif /* ALLOW_FE */ + RETURN + END