/[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.11 by stephd, Thu Dec 6 22:24:34 2007 UTC revision 1.15 by stephd, Mon Jun 23 22:26:32 2008 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "DIC_OPTIONS.h"  #include "DIC_OPTIONS.h"
 #include "GCHEM_OPTIONS.h"  
5    
6  CStartOfInterFace  CStartOfInterFace
7        SUBROUTINE Fe_CHEM(        SUBROUTINE Fe_CHEM(
# Line 22  C     == GLobal variables == Line 21  C     == GLobal variables ==
21  #include "EEPARAMS.h"  #include "EEPARAMS.h"
22  #include "PARAMS.h"  #include "PARAMS.h"
23  #include "GRID.h"  #include "GRID.h"
24  #include "DIC_BIOTIC.h"  #include "DIC_VARS.h"
25    
26  C     == Routine arguments ==                  C     == Routine arguments ==                
27  C     bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation  C     bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
# Line 63  C in surface layer Line 62  C in surface layer
62            IF (hFacC(i,j,k,bi,bj) .GT. 0. _d 0) THEN            IF (hFacC(i,j,k,bi,bj) .GT. 0. _d 0) THEN
63    
64  #ifdef DIC_NO_NEG  #ifdef DIC_NO_NEG
65                tmpfe=max(0. _d0 , fe (i,j,k,bi,bj))                tmpfe=max(0. _d 0 , fe (i,j,k,bi,bj))
66  #else  #else
67                tmpfe=fe (i,j,k,bi,bj)                tmpfe=fe (i,j,k,bi,bj)
68  #endif  #endif
# Line 77  C   Ligand,FeL,Fe calculation Line 76  C   Ligand,FeL,Fe calculation
76       &            )/(2. _d 0*ligand_stab)       &            )/(2. _d 0*ligand_stab)
77    
78                FeL = ligand_tot-lig                FeL = ligand_tot-lig
79                if (tmpfe.eq.0. _d 0) then                if (tmpfe.ne.0. _d 0) then
80                  freefe(i,j,k,bi,bj) = tmpfe -FeL                  freefe(i,j,k,bi,bj) = tmpfe -FeL
81                else                else
82                  freefe(i,j,k,bi,bj) = 0. _d 0                  freefe(i,j,k,bi,bj) = 0. _d 0
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22