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 |
1797 |
enddo |
enddo |
1798 |
#endif |
#endif |
1799 |
|
|
1800 |
|
#endif /* ALLOW_CARBON */ |
1801 |
|
|
1802 |
#ifdef ALLOW_CDOM |
#ifdef ALLOW_CDOM |
1803 |
dcdomdt = totphy_cdomp + totzoo_cdomp + preminP_cdom |
dcdomdt = totphy_cdomp + totzoo_cdomp + preminP_cdom |
1804 |
& -cdomp_degrd |
& -cdomp_degrd |
1805 |
#endif |
#endif |
1806 |
|
|
|
#endif |
|
|
|
|
1807 |
if (debug.eq.10) print*,'dDOPdt', dDOPdt |
if (debug.eq.10) print*,'dDOPdt', dDOPdt |
1808 |
if (debug.eq.10) print*,'dpopdt',dpopdt |
if (debug.eq.10) print*,'dpopdt',dpopdt |
1809 |
if (debug.eq.10) print*,'dDONdt',dDONdt |
if (debug.eq.10) print*,'dDONdt',dDONdt |