199 |
DO i=imin,imax |
DO i=imin,imax |
200 |
|
|
201 |
GPO4(i,j,k)=-BIOac(i,j,k)+pflux(i,j,k) |
GPO4(i,j,k)=-BIOac(i,j,k)+pflux(i,j,k) |
202 |
& + maskC(i,j,k,bi,bj)*Kdopremin*PTR_DOP(i,j,k) |
& + maskC(i,j,k,bi,bj)*KDOPRemin*PTR_DOP(i,j,k) |
203 |
|
|
204 |
car(i,j,k) =-BIOac(i,j,k)* R_CP*rain_ratio(i,j,bi,bj)* |
car(i,j,k) =-BIOac(i,j,k)* R_CP*rain_ratio(i,j,bi,bj)* |
205 |
& (1. _d 0-DOPfraction)+cflux(i,j,k) |
& (1. _d 0-DOPfraction)+cflux(i,j,k) |
206 |
|
|
207 |
GDOP(i,j,k)=+BIOac(i,j,k)*DOPfraction |
GDOP(i,j,k)=+BIOac(i,j,k)*DOPfraction |
208 |
& - maskC(i,j,k,bi,bj)*Kdopremin*PTR_DOP(i,j,k) |
& - maskC(i,j,k,bi,bj)*KDOPRemin*PTR_DOP(i,j,k) |
209 |
|
|
210 |
GALK(i,j,k)=+2. _d 0 *car(i,j,k)-R_NP*GPO4(i,j,k) |
GALK(i,j,k)=+2. _d 0 *car(i,j,k)-R_NP*GPO4(i,j,k) |
211 |
|
|