90 |
lit=sfac(j) |
lit=sfac(j) |
91 |
#endif |
#endif |
92 |
#endif |
#endif |
93 |
|
|
94 |
|
#if (defined (ALLOW_SEAICE) && defined (USE_QSW)) |
95 |
|
c if using Qsw and seaice, then ice fraction is already |
96 |
|
c taken into account |
97 |
|
lit=lit |
98 |
|
#else |
99 |
|
#if (defined (ALLOW_THSICE) && defined (USE_QSW)) |
100 |
|
c if using Qsw and seaice, then ice fraction is already |
101 |
|
c taken into account |
102 |
|
lit=lit |
103 |
|
#else |
104 |
lit=lit*(1. _d 0 - FIce(i,j,bi,bj)) |
lit=lit*(1. _d 0 - FIce(i,j,bi,bj)) |
105 |
|
#endif |
106 |
|
#endif |
107 |
DO k=1,nlev |
DO k=1,nlev |
108 |
C$TAF STORE lit = bio_export |
C$TAF STORE lit = bio_export |
109 |
atten=(k0*drF(k)*hFacC(i,j,k,bi,bj)*.5 _d 0) |
atten=(k0*drF(k)*hFacC(i,j,k,bi,bj)*.5 _d 0) |