92 |
INTEGER I,J,k |
INTEGER I,J,k |
93 |
INTEGER nCALCITEstep |
INTEGER nCALCITEstep |
94 |
CEOP |
CEOP |
95 |
|
jmin=1 |
96 |
|
jmax=sNy |
97 |
|
imin=1 |
98 |
|
imax=sNx |
99 |
|
|
100 |
DO k=1,Nr |
DO k=1,Nr |
101 |
DO j=1-OLy,sNy+OLy |
DO j=1-OLy,sNy+OLy |
181 |
|
|
182 |
c add all tendencies for PO4, DOP, ALK, DIC |
c add all tendencies for PO4, DOP, ALK, DIC |
183 |
DO k=1,Nr |
DO k=1,Nr |
184 |
DO j=1-OLy,sNy+OLy |
DO j=jmin,jmax |
185 |
DO i=1-OLx,sNx+OLx |
DO i=imin,imax |
186 |
bio(i,j,k)=-bioac(i,j,k)+pflux(i,j,k) |
bio(i,j,k)=-bioac(i,j,k)+pflux(i,j,k) |
187 |
& + maskC(i,j,k,bi,bj)*Kdopremin*PTR_DOP(i,j,k) |
& + maskC(i,j,k,bi,bj)*Kdopremin*PTR_DOP(i,j,k) |
188 |
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)* |
219 |
|
|
220 |
C update |
C update |
221 |
DO k=1,Nr |
DO k=1,Nr |
222 |
DO j=1-OLy,sNy+OLy |
DO j=jmin,jmax |
223 |
DO i=1-OLx,sNx+OLx |
DO i=imin,imax |
224 |
PTR_DIC(i,j,k)= |
PTR_DIC(i,j,k)= |
225 |
& PTR_DIC(i,j,k)+GDIC(i,j,k)*dTtracerLev(k) |
& PTR_DIC(i,j,k)+GDIC(i,j,k)*dTtracerLev(k) |
226 |
PTR_ALK(i,j,k)= |
PTR_ALK(i,j,k)= |
251 |
#ifdef ALLOW_TIMEAVE |
#ifdef ALLOW_TIMEAVE |
252 |
c save averages |
c save averages |
253 |
DO k=1,Nr |
DO k=1,Nr |
254 |
DO j=1-OLy,sNy+OLy |
DO j=jmin,jmax |
255 |
DO i=1-OLx,sNx+OLx |
DO i=imin,imax |
256 |
BIOave(i,j,k,bi,bj)=BIOave(i,j,k,bi,bj)+ |
BIOave(i,j,k,bi,bj)=BIOave(i,j,k,bi,bj)+ |
257 |
& BIOac(i,j,k)*deltaTclock |
& BIOac(i,j,k)*deltaTclock |
258 |
CARave(i,j,k,bi,bj)=CARave(i,j,k,bi,bj)+ |
CARave(i,j,k,bi,bj)=CARave(i,j,k,bi,bj)+ |