pkg/dic/fe_chem.F

C $Header: /u/gcmpack/MITgcm/pkg/dic/fe_chem.F,v 1.5 2004/10/18 16:01:13 jmc 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"
#include "PTRACERS_SIZE.h"
#include "PTRACERS.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

#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
              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                      
1	C $Header: /u/gcmpack/MITgcm/pkg/dic/fe_chem.F,v 1.5 2004/10/18 16:01:13 jmc Exp $
2	C $Name: $
3
4	#include "DIC_OPTIONS.h"
5	#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	#include "PTRACERS_SIZE.h"
27	#include "PTRACERS.h"
28
29	C == Routine arguments ==
30	C bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
31	C results will be set.
32	C myThid - Instance number for this innvocation of CALC_GT
33	_RL freefe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
34	_RL fe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
35	INTEGER bi,bj,iMin,iMax,jMin,jMax
36	INTEGER myIter,myThid
37	CEndOfInterface
38
39	#ifdef ALLOW_FE
40
41	INTEGER I,J,K
42	_RL lig, FeL
43
44	CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
45	CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
46	CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
47	CC
48	CC ADAPTED FROM PAYAL
49	CC
50	CC
51	CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
52	CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
53	CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
54
55	C ligand balance in surface layer
56	C in surface layer
57
58	DO j=jMin,jMax
59	DO i=iMin,iMax
60	DO k=1,nR
61	IF(hFacC(i,j,k,bi,bj) .gt. 0.0)THEN
62
63	C Ligand,FeL,Fe calculation
64
65	lig=(-ligand_stab*fe (i,j,k,bi,bj)+
66	& ligand_stab*ligand_tot-1
67	& +((ligand_stab*fe (i,j,k,bi,bj)
68	& -ligand_stabligand_tot+1)*2+4
69	& ligand_stabligand_tot)*0.5)/(2ligand_stab)
70
71	FeL = ligand_tot-lig
72	freefe(i,j,k,bi,bj) = fe (i,j,k,bi,bj)-FeL
73	#ifdef MINFE
74	freefe(i,j,k,bi,bj) = min(freefe(i,j,k,bi,bj),freefemax)
75	fe(i,j,k,bi,bj) = FeL+freefe(i,j,k,bi,bj)
76	#endif
77	END IF
78	ENDDO
79	ENDDO
80	ENDDO
81	c
82	#endif
83	RETURN
84	END