87 |
lit=sfac(j) |
lit=sfac(j) |
88 |
#endif |
#endif |
89 |
#endif |
#endif |
90 |
|
lit=lit*(1. _d 0 - FIce(i,j,bi,bj)) |
91 |
DO k=1,nlev |
DO k=1,nlev |
92 |
C$TAF STORE lit = bio_export |
C$TAF STORE lit = bio_export |
93 |
atten=(k0*drF(k)*hFacC(i,j,k,bi,bj)*.5) |
atten=(k0*drF(k)*hFacC(i,j,k,bi,bj)*.5) |
94 |
if (k.gt.1) atten=atten+(k0*drF(k-1) |
if (k.gt.1) atten=atten+(k0*drF(k-1) |
95 |
& *hFacC(i,j,k-1,bi,bj)*.5) |
& *hFacC(i,j,k-1,bi,bj)*.5) |
96 |
lit=lit*exp(-atten)*(1.d0-Fice(i,j,bi,bj)) |
lit=lit*exp(-atten) |
97 |
if (lit.lt.0.d0.or.lit.gt.150) then |
if (lit.lt.0.d0.or.lit.gt.150) then |
98 |
print*,'QQ lit', lit |
print*,'QQ lit', lit |
99 |
endif |
endif |