/[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.1 by stephd, Wed Jun 25 21:00:36 2003 UTC revision 1.9 by stephd, Tue Nov 28 21:16:03 2006 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    
7  CStartOfInterFace  CStartOfInterFace
# Line 20  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"
 #include "PTRACERS.h"  
26    
27  C     == Routine arguments ==                  C     == Routine arguments ==                
28  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 30  C     myThid - Instance number for this Line 32  C     myThid - Instance number for this
32        _RL  fe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL  fe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
33        INTEGER bi,bj,iMin,iMax,jMin,jMax        INTEGER bi,bj,iMin,iMax,jMin,jMax
34        INTEGER myIter,myThid        INTEGER myIter,myThid
35    #ifdef AD_SAFE
36          _RL thx, thy, theps
37    #endif
38  CEndOfInterface  CEndOfInterface
39    
40    #ifdef ALLOW_FE
41    
42        INTEGER I,J,K        INTEGER I,J,K
43        _RL  lig, FeL        _RL  lig, FeL
# Line 50  CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC Line 56  CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
56  C ligand balance in surface layer  C ligand balance in surface layer
57  C in surface layer  C in surface layer
58    
59         DO j=jMin,jMax         DO j=jmin,jmax
60          DO i=iMin,iMax          DO i=imin,imax
61           DO k=1,nR           DO k=1,nR
62            IF(hFacC(i,j,k,bi,bj) .gt. 0.0)THEN            IF(hFacC(i,j,k,bi,bj) .gt. 0.0)THEN
63    
# Line 65  C   Ligand,FeL,Fe calculation Line 71  C   Ligand,FeL,Fe calculation
71    
72                FeL = ligand_tot-lig                FeL = ligand_tot-lig
73                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
74    #ifdef MINFE
75    #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                  freefe(i,j,k,bi,bj) = min(freefe(i,j,k,bi,bj),freefemax)
85    #endif
86                  fe(i,j,k,bi,bj) = FeL+freefe(i,j,k,bi,bj)
87    #endif
88            END IF              END IF  
89           ENDDO           ENDDO
90          ENDDO          ENDDO
91         ENDDO             ENDDO    
92  c  c
93    #endif
94          RETURN          RETURN
95          END                                END                      
   
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22