/[MITgcm]/MITgcm/pkg/dic/fe_chem.F
ViewVC logotype

Diff of /MITgcm/pkg/dic/fe_chem.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by stephd, Mon Oct 6 20:11:10 2003 UTC revision 1.7 by stephd, Mon Dec 12 19:07:36 2005 UTC
# Line 1  Line 1 
1  #include "CPP_OPTIONS.h"  C $Header$
2    C $Name$
3    
4    #include "DIC_OPTIONS.h"
5  #include "GCHEM_OPTIONS.h"  #include "GCHEM_OPTIONS.h"
6    
 #ifdef ALLOW_FE  
7  CStartOfInterFace  CStartOfInterFace
8        SUBROUTINE Fe_CHEM(        SUBROUTINE Fe_CHEM(
9       I           bi,bj,iMin,iMax,jMin,jMax,       I           bi,bj,iMin,iMax,jMin,jMax,
# Line 21  C     == GLobal variables == Line 23  C     == GLobal variables ==
23  #include "PARAMS.h"  #include "PARAMS.h"
24  #include "GRID.h"  #include "GRID.h"
25  #include "DIC_BIOTIC.h"  #include "DIC_BIOTIC.h"
26    #include "PTRACERS_SIZE.h"
27  #include "PTRACERS.h"  #include "PTRACERS.h"
28    
29  C     == Routine arguments ==                  C     == Routine arguments ==                
# Line 33  C     myThid - Instance number for this Line 36  C     myThid - Instance number for this
36        INTEGER myIter,myThid        INTEGER myIter,myThid
37  CEndOfInterface  CEndOfInterface
38    
39    #ifdef ALLOW_FE
40    
41        INTEGER I,J,K        INTEGER I,J,K
42        _RL  lig, FeL        _RL  lig, FeL
# Line 51  CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC Line 55  CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
55  C ligand balance in surface layer  C ligand balance in surface layer
56  C in surface layer  C in surface layer
57    
58         DO j=jMin,jMax         DO j=1-OLy,sNy+OLy
59          DO i=iMin,iMax          DO i=1-OLx,sNx+OLx
60           DO k=1,nR           DO k=1,nR
61            IF(hFacC(i,j,k,bi,bj) .gt. 0.0)THEN            IF(hFacC(i,j,k,bi,bj) .gt. 0.0)THEN
62    
# Line 66  C   Ligand,FeL,Fe calculation Line 70  C   Ligand,FeL,Fe calculation
70    
71                FeL = ligand_tot-lig                FeL = ligand_tot-lig
72                freefe(i,j,k,bi,bj) = fe (i,j,k,bi,bj)-FeL                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              END IF  
78           ENDDO           ENDDO
79          ENDDO          ENDDO
80         ENDDO             ENDDO    
81  c  c
82    #endif
83          RETURN          RETURN
84          END                                END                      
 #endif  
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22