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) ) |