/[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.13 by dfer, Mon Apr 7 20:31:16 2008 UTC revision 1.15 by stephd, Mon Jun 23 22:26:32 2008 UTC
# Line 62  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 76  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.13  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22