83 |
C$TAF INIT bio_export = static, 10 |
C$TAF INIT bio_export = static, 10 |
84 |
#ifdef READ_PAR |
#ifdef READ_PAR |
85 |
lit=PAR(i,j,bi,bj) |
lit=PAR(i,j,bi,bj) |
86 |
#else |
#elif (defined USE_QSW) |
|
#ifdef USE_QSW |
|
87 |
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) |
88 |
#else |
#else |
89 |
lit=sfac(j) |
lit=sfac(j) |
90 |
#endif |
#endif |
|
#endif |
|
91 |
|
|
92 |
if (QSW_underice) then |
IF ( .NOT. QSW_underice ) THEN |
93 |
c if using Qsw and seaice/thsice or coupled, then ice |
c if using Qsw but not seaice/thsice or coupled, then |
94 |
c fraction is already taken into account |
c ice fraction needs to be taken into account |
|
lit=lit |
|
|
else |
|
95 |
lit=lit*(1. _d 0 - FIce(i,j,bi,bj)) |
lit=lit*(1. _d 0 - FIce(i,j,bi,bj)) |
96 |
endif |
ENDIF |
97 |
c |
c |
98 |
DO k=1,nlev |
DO k=1,nlev |
99 |
C$TAF STORE lit = bio_export |
C$TAF STORE lit = bio_export |