/[MITgcm]/MITgcm/pkg/dic/dic_biotic_forcing.F
ViewVC logotype

Diff of /MITgcm/pkg/dic/dic_biotic_forcing.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.12 by stephd, Tue Nov 28 21:16:03 2006 UTC revision 1.13 by stephd, Tue Dec 12 22:37:28 2006 UTC
# Line 84  c  freefe                 :: iron not bo Line 84  c  freefe                 :: iron not bo
84        _RL  CAR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  CAR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
85        _RL  bioac(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  bioac(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
86        _RL  pflux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  pflux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
87          _RL  exportflux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
88        _RL  cflux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  cflux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
89  #ifdef ALLOW_FE  #ifdef ALLOW_FE
90        _RL  GFE(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  GFE(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
# Line 112  CEOP Line 113  CEOP
113             BIO_kar(i,j,k)=0.d0             BIO_kar(i,j,k)=0.d0
114             bioac(i,j,k)=0.d0             bioac(i,j,k)=0.d0
115             pflux(i,j,k)=0.d0             pflux(i,j,k)=0.d0
116               exportflux(i,j,k)=0.d0
117             cflux(i,j,k)=0.d0             cflux(i,j,k)=0.d0
118  #ifdef ALLOW_FE  #ifdef ALLOW_FE
119             GFE(i,j,k)=0.d0             GFE(i,j,k)=0.d0
# Line 153  c biological activity Line 155  c biological activity
155       I           myIter,myTime,myThid)       I           myIter,myTime,myThid)
156    
157  c flux of po4 from layers with biological activity  c flux of po4 from layers with biological activity
158         CALL PHOS_FLUX( bioac, pflux,         CALL PHOS_FLUX( bioac, pflux, exportflux,
159       &                    bi,bj,imin,imax,jmin,jmax,       &                    bi,bj,imin,imax,jmin,jmax,
160       &                    myIter,myTime,myThid)       &                    myIter,myTime,myThid)
161    
# Line 261  c save averages Line 263  c save averages
263       &                           OmegaC(i,j,k,bi,bj)*deltaTclock       &                           OmegaC(i,j,k,bi,bj)*deltaTclock
264              pfluxave(i,j,k,bi,bj)= pfluxave(i,j,k,bi,bj) +              pfluxave(i,j,k,bi,bj)= pfluxave(i,j,k,bi,bj) +
265       &                           pflux(i,j,k)*deltaTclock       &                           pflux(i,j,k)*deltaTclock
266                epfluxave(i,j,k,bi,bj)= epfluxave(i,j,k,bi,bj) +
267         &                           exportflux(i,j,k)*deltaTclock
268              cfluxave(i,j,k,bi,bj)= cfluxave(i,j,k,bi,bj) +              cfluxave(i,j,k,bi,bj)= cfluxave(i,j,k,bi,bj) +
269       &                           cflux(i,j,k)*deltaTclock       &                           cflux(i,j,k)*deltaTclock
270              if (k.eq.1) then              if (k.eq.1) then
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22