/[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.23 by jmc, Tue Aug 4 18:33:07 2009 UTC revision 1.24 by dfer, Tue Dec 29 01:13:47 2009 UTC
# Line 30  C !USES: =============================== Line 30  C !USES: ===============================
30  #include "PARAMS.h"  #include "PARAMS.h"
31  #include "GRID.h"  #include "GRID.h"
32  #include "DIC_VARS.h"  #include "DIC_VARS.h"
33    #include "PTRACERS_SIZE.h"
34    #include "PTRACERS_PARAMS.h"
35    
36  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
37  C  myThid               :: thread number  C  myThid               :: thread number
# Line 259  C update Line 261  C update
261           DO j=jmin,jmax           DO j=jmin,jmax
262            DO i=imin,imax            DO i=imin,imax
263             PTR_DIC(i,j,k)=             PTR_DIC(i,j,k)=
264       &      PTR_DIC(i,j,k)+GDIC(i,j,k)*dTtracerLev(k)       &      PTR_DIC(i,j,k)+GDIC(i,j,k)*PTRACERS_dTLev(k)
265             PTR_ALK(i,j,k)=             PTR_ALK(i,j,k)=
266       &      PTR_ALK(i,j,k)+GALK(i,j,k)*dTtracerLev(k)       &      PTR_ALK(i,j,k)+GALK(i,j,k)*PTRACERS_dTLev(k)
267             PTR_PO4(i,j,k)=             PTR_PO4(i,j,k)=
268       &      PTR_PO4(i,j,k)+GPO4(i,j,k)*dTtracerLev(k)       &      PTR_PO4(i,j,k)+GPO4(i,j,k)*PTRACERS_dTLev(k)
269             PTR_DOP(i,j,k)=             PTR_DOP(i,j,k)=
270       &      PTR_DOP(i,j,k)+GDOP(i,j,k)*dTtracerLev(k)       &      PTR_DOP(i,j,k)+GDOP(i,j,k)*PTRACERS_dTLev(k)
271  #ifdef ALLOW_O2  #ifdef ALLOW_O2
272             PTR_O2(i,j,k)=             PTR_O2(i,j,k)=
273       &      PTR_O2(i,j,k)+GO2(i,j,k)*dTtracerLev(k)       &      PTR_O2(i,j,k)+GO2(i,j,k)*PTRACERS_dTLev(k)
274  #endif  #endif
275  #ifdef ALLOW_FE  #ifdef ALLOW_FE
276             PTR_FE(i,j,k)=             PTR_FE(i,j,k)=
277       &      PTR_FE(i,j,k)+GFE(i,j,k)*dTtracerLev(k)       &      PTR_FE(i,j,k)+GFE(i,j,k)*PTRACERS_dTLev(k)
278  #endif  #endif
279            ENDDO            ENDDO
280           ENDDO           ENDDO
Legend:
Removed from v.1.23  
changed lines
  Added in v.1.24
  ViewVC Help
Powered by ViewVC 1.1.22