--- MITgcm/pkg/dic/fe_chem.F 2003/10/09 04:19:19 1.3 +++ MITgcm/pkg/dic/fe_chem.F 2010/04/13 03:06:14 1.16 @@ -1,76 +1,107 @@ +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" -#include "GCHEM_OPTIONS.h" -#ifdef ALLOW_FE -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_BIOTIC.h" -#include "PTRACERS.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 -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 - - INTEGER I,J,K +#ifdef ALLOW_FE +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.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) ) +#else + tmpfe = fe(i,j,k) +#endif C Ligand,FeL,Fe calculation - - lig=(-ligand_stab*fe (i,j,k,bi,bj)+ - & ligand_stab*ligand_tot-1 - & +((ligand_stab*fe (i,j,k,bi,bj) - & -ligand_stab*ligand_tot+1)**2+4 - & *ligand_stab*ligand_tot)**0.5)/(2*ligand_stab) - + lig=(-ligand_stab*tmpfe + + & ligand_stab*ligand_tot-1. _d 0 + & +((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 - freefe(i,j,k,bi,bj) = fe (i,j,k,bi,bj)-FeL - END IF + 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) + thy=freefemax + theps=1. _d -8 + 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) = MIN(freefe(i,j,k),freefemax) +#endif + fe(i,j,k) = FeL+freefe(i,j,k) +#endif + ENDIF ENDDO ENDDO - ENDDO -c - RETURN - END -#endif + ENDDO + +#endif /* ALLOW_FE */ + RETURN + END