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

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

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

revision 1.13 by stephd, Thu Aug 9 19:51:18 2007 UTC revision 1.23 by dfer, Mon Apr 7 20:31:16 2008 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "DIC_OPTIONS.h"  #include "DIC_OPTIONS.h"
 #include "GCHEM_OPTIONS.h"  
5    
6  CBOP  CBOP
7  C !ROUTINE: BIO_EXPORT  C !ROUTINE: BIO_EXPORT
8    
9  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
10        SUBROUTINE BIO_EXPORT( PTR_PO4 ,        SUBROUTINE BIO_EXPORT( PTR_PO4 ,
11  #ifdef ALLOW_FE  #ifdef ALLOW_FE
12       I           PTR_FE,       I           PTR_FE,
13  #endif    #endif
14       I           bioac,       I           bioac,
15       I           bi,bj,imin,imax,jmin,jmax,       I           bi,bj,imin,imax,jmin,jmax,
16       I           myIter,myTime,myThid)       I           myIter,myTime,myThid)
17    
18  c !DESCRIPTION:  c !DESCRIPTION:
19  C  Calculate biological activity and export                  C  Calculate biological activity and export
20    
21  C !USES: ===============================================================  C !USES: ===============================================================
22        IMPLICIT NONE        IMPLICIT NONE
# Line 23  C !USES: =============================== Line 25  C !USES: ===============================
25  #include "EEPARAMS.h"  #include "EEPARAMS.h"
26  #include "PARAMS.h"  #include "PARAMS.h"
27  #include "GRID.h"  #include "GRID.h"
28  #include "DIC_ABIOTIC.h"  #include "DIC_VARS.h"
 #include "DIC_BIOTIC.h"  
29  #ifdef USE_QSW  #ifdef USE_QSW
30  #include "FFIELDS.h"  #include "FFIELDS.h"
31  #endif  #endif
# Line 59  C  i,j,k                  :: loop indice Line 60  C  i,j,k                  :: loop indice
60        _RL sfac(1-OLy:sNy+OLy)        _RL sfac(1-OLy:sNy+OLy)
61        _RL lit, atten        _RL lit, atten
62        _RL nutlimit        _RL nutlimit
63          _RL tmpfe, tmppo4
64  #ifdef AD_SAFE  #ifdef AD_SAFE
65        _RL thx, thy, theps        _RL thx, thy, theps
66  #endif  #endif
# Line 80  C FORTRAN-77 with know max of nlev Line 82  C FORTRAN-77 with know max of nlev
82  C$TAF INIT bio_export = static, 10  C$TAF INIT bio_export = static, 10
83  #ifdef READ_PAR  #ifdef READ_PAR
84             lit=PAR(i,j,bi,bj)             lit=PAR(i,j,bi,bj)
85  #else  #elif (defined USE_QSW)
 #ifdef USE_QSW  
86             lit=-parfrac*Qsw(i,j,bi,bj)*maskC(i,j,1,bi,bj)             lit=-parfrac*Qsw(i,j,bi,bj)*maskC(i,j,1,bi,bj)
87  #else  #else
88             lit=sfac(j)             lit=sfac(j)
89  #endif  #endif
90  #endif  
91               IF ( .NOT. QSW_underice ) THEN
92    c if using Qsw but not seaice/thsice or coupled, then
93    c ice fraction needs to be taken into account
94                  lit=lit*(1. _d 0 - FIce(i,j,bi,bj))
95               ENDIF
96    c
97             DO k=1,nlev             DO k=1,nlev
98  C$TAF STORE lit = bio_export  C$TAF STORE lit = bio_export
99               atten=(k0*drF(k)*hFacC(i,j,k,bi,bj)*.5)               atten=(k0*drF(k)*hFacC(i,j,k,bi,bj)*.5 _d 0)
100               if (k.gt.1) atten=atten+(k0*drF(k-1)               if (k.gt.1) atten=atten+(k0*drF(k-1)
101       &                         *hFacC(i,j,k-1,bi,bj)*.5)       &                         *hFacC(i,j,k-1,bi,bj)*.5 _d 0)
102               lit=lit*exp(-atten)*(1.d0-Fice(i,j,bi,bj))               lit=lit*exp(-atten)
103               if (lit.lt.0.d0.or.lit.gt.150) then               IF (lit.LT.0. _d 0.OR.lit.GT.350. _d 0) THEN
104                   print*,'QQ lit', lit                   print*,'QQ lit', lit
105               endif               ENDIF
106    
107    #ifdef DIC_NO_NEG
108                 tmppo4=max(0. _d 0, PTR_PO4(i,j,k))
109                 lit=max(0. _d 0,lit)
110    #else
111                 tmppo4=PTR_PO4(i,j,k)
112    #endif
113    
114  #ifdef ALLOW_FE  #ifdef ALLOW_FE
115    #ifdef DIC_NO_NEG
116                 tmpfe=max(0. _d 0,PTR_FE(i,j,k))
117    #else
118                 tmpfe=PTR_FE(i,j,k)
119    #endif
120  #ifdef AD_SAFE  #ifdef AD_SAFE
121               thx = PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4)               thx = tmppo4/(tmppo4+KPO4)
122               thy = PTR_FE(i,j,k)/(PTR_FE(i,j,k)+KFE)               thy = tmpfe/(tmpfe+KFE)
123               theps = 1.d-6  c            thx = PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4)
124               nutlimit= ( 1.d0 - tanh((thx-thy)/theps) ) * thx/2 +  c            thy = PTR_FE(i,j,k)/(PTR_FE(i,j,k)+KFE)
125       &                 ( 1.d0 + tanh((thx-thy)/theps) ) * thy/2               theps = 1. _d -6
126                 nutlimit= ( 1. _d 0 - tanh((thx-thy)/theps) ) * thx/2. _d 0
127         &                +( 1. _d 0 + tanh((thx-thy)/theps) ) * thy/2. _d 0
128  #else  #else
129               nutlimit=   min(PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4),               nutlimit=   min(tmppo4/(tmppo4+KPO4),
130       &                        PTR_FE(i,j,k)/(PTR_FE(i,j,k)+KFE) )       &                        tmpfe/(tmpfe+KFE) )
131  #endif  #endif
132  #else  #else
133               nutlimit=    PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4)               nutlimit=    tmppo4/(tmppo4+KPO4)
134  #endif  #endif
135    
136               bioac(i,j,k)=alpha(i,j,bi,bj)*               bioac(i,j,k)=alpha(i,j,bi,bj)*
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.23
  ViewVC Help
Powered by ViewVC 1.1.22