/[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.2 by stephd, Wed Jul 9 19:59:18 2003 UTC revision 1.3 by stephd, Mon Oct 6 20:11:10 2003 UTC
# Line 147  c add all tendencies for PO4, DOP, ALK, Line 147  c add all tendencies for PO4, DOP, ALK,
147       &         - maskC(i,j,k,bi,bj)*Kdopremin*PTR_DOP(i,j,k)       &         - maskC(i,j,k,bi,bj)*Kdopremin*PTR_DOP(i,j,k)
148             GALK(i,j,k)=+2.d0*car(i,j,k)-R_NP*bio(i,j,k)             GALK(i,j,k)=+2.d0*car(i,j,k)-R_NP*bio(i,j,k)
149             GDIC(i,j,k)=car(i,j,k)+R_CP*bio(i,j,k)             GDIC(i,j,k)=car(i,j,k)+R_CP*bio(i,j,k)
 cQQchange start  
150             if (PTR_O2(i,j,k).gt.o2crit) then             if (PTR_O2(i,j,k).gt.o2crit) then
151               GO2(i,j,k)=R_OP*bio(i,j,k)               GO2(i,j,k)=R_OP*bio(i,j,k)
152             else             else
153               GO2(i,j,k)=0.d0               GO2(i,j,k)=0.d0
154             endif             endif
 cQQchange end  
155  #ifdef ALLOW_FE  #ifdef ALLOW_FE
156             GFE(i,j,k)=R_FeP*bio(i,j,k)             GFE(i,j,k)=R_FeP*bio(i,j,k)
157       &             -Kscav*freefe(i,j,k)       &             -Kscav*freefe(i,j,k)
# Line 190  C update Line 188  C update
188             PTR_FE(i,j,k)=             PTR_FE(i,j,k)=
189       &      PTR_FE(i,j,k)+GFE(i,j,k)*deltaTtracer       &      PTR_FE(i,j,k)+GFE(i,j,k)*deltaTtracer
190  #endif  #endif
 cQQchange start  
 cQQ        PTR_O2(i,j,k)=max(PTR_O2(i,j,k),O2crit)  
 cQQchange end  
191            ENDDO            ENDDO
192           ENDDO           ENDDO
193         ENDDO         ENDDO
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22