pkg/dic/fe_chem.F

C $Header: /u/gcmpack/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"

CStartOfInterFace
      SUBROUTINE Fe_CHEM(
     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     |==========================================================|
      IMPLICIT NONE

C     == GLobal variables ==
#include "SIZE.h"
#include "DYNVARS.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "DIC_BIOTIC.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

#ifdef ALLOW_FE

      INTEGER I,J,K
      _RL  lig, FeL

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
CC
CC ADAPTED FROM PAYAL
CC
CC
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
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

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)
 
              FeL = ligand_tot-lig
              freefe(i,j,k,bi,bj) = fe (i,j,k,bi,bj)-FeL
#ifdef MINFE
#ifdef AD_SAFE
              thx=freefe(i,j,k,bi,bj)
              thy=freefemax
              theps=1.d-8
              freefe(i,j,k,bi,bj) =
     &                 ( 1.d0 - tanh((thx-thy)/theps) ) * thx/2 +
     &                 ( 1.d0 + tanh((thx-thy)/theps) ) * thy/2

#else
              freefe(i,j,k,bi,bj) = min(freefe(i,j,k,bi,bj),freefemax)
#endif
              fe(i,j,k,bi,bj) = FeL+freefe(i,j,k,bi,bj)
#endif
          END IF  
         ENDDO
        ENDDO
       ENDDO    
c
#endif 
        RETURN
        END                      
1	stephd	1.9	C $Header: /u/gcmpack/MITgcm/pkg/dic/fe_chem.F,v 1.8 2005/12/16 21:07:53 stephd Exp $
2	jmc	1.4	C $Name: $
3
4	edhill	1.3	#include "DIC_OPTIONS.h"
5	stephd	1.1	#include "GCHEM_OPTIONS.h"
6
7			CStartOfInterFace
8			SUBROUTINE Fe_CHEM(
9			I bi,bj,iMin,iMax,jMin,jMax,
10			I fe, freefe,
11			I myIter, myThid )
12			C /==========================================================\
13			C \| SUBROUTINE Fe_chem \|
14			C \| \|
15			C \| o Calculate L,FeL,Fe concentration \|
16			C \|==========================================================\|
17			IMPLICIT NONE
18
19			C == GLobal variables ==
20			#include "SIZE.h"
21			#include "DYNVARS.h"
22			#include "EEPARAMS.h"
23			#include "PARAMS.h"
24			#include "GRID.h"
25			#include "DIC_BIOTIC.h"
26
27			C == Routine arguments ==
28			C bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
29			C results will be set.
30			C myThid - Instance number for this innvocation of CALC_GT
31			_RL freefe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
32			_RL fe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
33			INTEGER bi,bj,iMin,iMax,jMin,jMax
34			INTEGER myIter,myThid
35	stephd	1.9	#ifdef AD_SAFE
36			_RL thx, thy, theps
37			#endif
38	stephd	1.1	CEndOfInterface
39
40	jmc	1.5	#ifdef ALLOW_FE
41	stephd	1.1
42			INTEGER I,J,K
43			_RL lig, FeL
44
45			CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
46			CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
47			CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
48			CC
49			CC ADAPTED FROM PAYAL
50			CC
51			CC
52			CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
53			CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
54			CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
55
56			C ligand balance in surface layer
57			C in surface layer
58
59	stephd	1.8	DO j=jmin,jmax
60			DO i=imin,imax
61	stephd	1.1	DO k=1,nR
62			IF(hFacC(i,j,k,bi,bj) .gt. 0.0)THEN
63
64			C Ligand,FeL,Fe calculation
65
66			lig=(-ligand_stab*fe (i,j,k,bi,bj)+
67			& ligand_stab*ligand_tot-1
68			& +((ligand_stab*fe (i,j,k,bi,bj)
69			& -ligand_stabligand_tot+1)*2+4
70			& ligand_stabligand_tot)*0.5)/(2ligand_stab)
71
72			FeL = ligand_tot-lig
73			freefe(i,j,k,bi,bj) = fe (i,j,k,bi,bj)-FeL
74	stephd	1.6	#ifdef MINFE
75	stephd	1.9	#ifdef AD_SAFE
76			thx=freefe(i,j,k,bi,bj)
77			thy=freefemax
78			theps=1.d-8
79			freefe(i,j,k,bi,bj) =
80			& ( 1.d0 - tanh((thx-thy)/theps) ) * thx/2 +
81			& ( 1.d0 + tanh((thx-thy)/theps) ) * thy/2
82
83			#else
84	stephd	1.6	freefe(i,j,k,bi,bj) = min(freefe(i,j,k,bi,bj),freefemax)
85	stephd	1.9	#endif
86	stephd	1.6	fe(i,j,k,bi,bj) = FeL+freefe(i,j,k,bi,bj)
87			#endif
88	stephd	1.1	END IF
89			ENDDO
90			ENDDO
91			ENDDO
92			c
93	jmc	1.5	#endif
94	stephd	1.1	RETURN
95			END