213 |
DO k=1,Nr |
DO k=1,Nr |
214 |
DO j=jmin,jmax |
DO j=jmin,jmax |
215 |
DO i=imin,imax |
DO i=imin,imax |
216 |
RDOP(i,j,k)= maskC(i,j,k,bi,bj)*KDOPRemin*PTR_DOP(i,j,k) |
#ifdef DIC_NO_NEG |
217 |
|
RDOP(i,j,k)= MAX(maskC(i,j,k,bi,bj)*KDOPRemin*PTR_DOP(i,j,k) |
218 |
|
& ,0. _d 0) |
219 |
|
#else |
220 |
|
RDOP(i,j,k)= maskC(i,j,k,bi,bj)*KDOPRemin*PTR_DOP(i,j,k) |
221 |
|
#endif |
222 |
GPO4(i,j,k)=-BIOac(i,j,k)+pflux(i,j,k) + RDOP(i,j,k) |
GPO4(i,j,k)=-BIOac(i,j,k)+pflux(i,j,k) + RDOP(i,j,k) |
223 |
|
|
224 |
car(i,j,k) = cflux(i,j,k) - CAR_S(i,j,k) |
car(i,j,k) = cflux(i,j,k) - CAR_S(i,j,k) |