1337 |
if (debug.eq.14) print*,'grow', grow(np), pcarbon(np) |
if (debug.eq.14) print*,'grow', grow(np), pcarbon(np) |
1338 |
#ifdef DYNAMIC_CHL |
#ifdef DYNAMIC_CHL |
1339 |
c geider 97 for dChl/dt (source part) Eq. 3 |
c geider 97 for dChl/dt (source part) Eq. 3 |
1340 |
if (acclim(np).gt. 0. _d 0) then |
if (acclim(np).gt. 0. _d 0.and. |
1341 |
|
& alpha_I(np).gt. 0. _d 0) then |
1342 |
rhochl(np)=chl2cmax(np) * |
rhochl(np)=chl2cmax(np) * |
1343 |
& (grow(np)/(alpha_I(np)*acclim(np)) ) |
& (grow(np)/(alpha_I(np)*acclim(np)) ) |
1344 |
else |
else |