/[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.16 by stephd, Mon Oct 15 16:02:58 2007 UTC revision 1.17 by dfer, Fri Oct 26 21:08:12 2007 UTC
# Line 93  C$TAF INIT bio_export = static, 10 Line 93  C$TAF INIT bio_export = static, 10
93             lit=lit*(1. _d 0 - FIce(i,j,bi,bj))             lit=lit*(1. _d 0 - FIce(i,j,bi,bj))
94             DO k=1,nlev             DO k=1,nlev
95  C$TAF STORE lit = bio_export  C$TAF STORE lit = bio_export
96               atten=(k0*drF(k)*hFacC(i,j,k,bi,bj)*.5)               atten=(k0*drF(k)*hFacC(i,j,k,bi,bj)*.5 _d 0)
97               if (k.gt.1) atten=atten+(k0*drF(k-1)               if (k.gt.1) atten=atten+(k0*drF(k-1)
98       &                         *hFacC(i,j,k-1,bi,bj)*.5)       &                         *hFacC(i,j,k-1,bi,bj)*.5 _d 0)
99               lit=lit*exp(-atten)               lit=lit*exp(-atten)
100               if (lit.lt.0.d0.or.lit.gt.350) then               IF (lit.LT.0. _d 0.OR.lit.GT.350. _d 0) THEN
101                   print*,'QQ lit', lit                   print*,'QQ lit', lit
102               endif               ENDIF
103    
104  #ifdef ALLOW_FE  #ifdef ALLOW_FE
105  #ifdef AD_SAFE  #ifdef AD_SAFE
106               thx = PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4)               thx = PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4)
107               thy = PTR_FE(i,j,k)/(PTR_FE(i,j,k)+KFE)               thy = PTR_FE(i,j,k)/(PTR_FE(i,j,k)+KFE)
108               theps = 1.d-6               theps = 1. _d -6
109               nutlimit= ( 1.d0 - tanh((thx-thy)/theps) ) * thx/2 +               nutlimit= ( 1. _d 0 - tanh((thx-thy)/theps) ) * thx/2. _d 0
110       &                 ( 1.d0 + tanh((thx-thy)/theps) ) * thy/2       &                +( 1. _d 0 + tanh((thx-thy)/theps) ) * thy/2. _d 0
111  #else  #else
112               nutlimit=   min(PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4),               nutlimit=   min(PTR_PO4(i,j,k)/(PTR_PO4(i,j,k)+KPO4),
113       &                        PTR_FE(i,j,k)/(PTR_FE(i,j,k)+KFE) )       &                        PTR_FE(i,j,k)/(PTR_FE(i,j,k)+KFE) )
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22