/[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.5 by jmc, Tue Oct 28 22:19:53 2003 UTC revision 1.6 by stephd, Thu Feb 12 16:11:46 2004 UTC
# Line 1  Line 1 
 C modified for external_forcing_DIC.F  August 1999  
 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC  
 c  
 c modified swd Oct 01 and Feb 02, for use as package for c40_patch1  
 c modified to use with c44 and ptracers: swd May 2002  
 c modified to have carbonate and biological influences: swd June 2002  
 C  
   
1  #include "DIC_OPTIONS.h"  #include "DIC_OPTIONS.h"
2  #include "GCHEM_OPTIONS.h"  #include "GCHEM_OPTIONS.h"
3    
4  CStartOfInterFace  CBOP
5    C !ROUTINE: DIC_BIOTIC_FORCING
6    
7    C !INTERFACE: ==========================================================
8        SUBROUTINE DIC_BIOTIC_FORCING( PTR_DIC, PTR_ALK, PTR_PO4,        SUBROUTINE DIC_BIOTIC_FORCING( PTR_DIC, PTR_ALK, PTR_PO4,
9       &                            PTR_DOP, PTR_O2,       &                            PTR_DOP, PTR_O2,
10  #ifdef ALLOW_FE  #ifdef ALLOW_FE
# Line 18  CStartOfInterFace Line 13  CStartOfInterFace
13       &                            bi,bj,imin,imax,jmin,jmax,       &                            bi,bj,imin,imax,jmin,jmax,
14       &                             myIter,myTime,myThid)       &                             myIter,myTime,myThid)
15    
16  C     /==========================================================\  C !DESCRIPTION:
17  C     | SUBROUTINE DIC_BIOTIC_FORCING                                   |  C updates all the tracers for the effects of air-sea exchange, biological
18  C     | o Calculate the changes to DIC through air-sea  fluxes   |    c activity and remineralization
 C     |   carbonates, and biology (parameterizes by PO4)         |  
 C     |==========================================================|  
       IMPLICIT NONE  
19    
20  C     == GLobal variables ==  C !USES: ===============================================================
21          IMPLICIT NONE
22  #include "SIZE.h"  #include "SIZE.h"
23  #include "DYNVARS.h"  #include "DYNVARS.h"
24  #include "EEPARAMS.h"  #include "EEPARAMS.h"
# Line 34  C     == GLobal variables == Line 27  C     == GLobal variables ==
27  #include "DIC_BIOTIC.h"  #include "DIC_BIOTIC.h"
28  #include "DIC_ABIOTIC.h"  #include "DIC_ABIOTIC.h"
29    
30  C     == Routine arguments ==  C !INPUT PARAMETERS: ===================================================
31    C  myThid               :: thread number
32    C  myIter               :: current timestep
33    C  myTime               :: current time
34    C  PTR_DIC              :: dissolced inorganic carbon
35    C  PTR_ALK              :: alkalinity
36    C  PTR_PO4              :: phosphate
37    c  PTR_DOP              :: dissolve organic phosphurous
38    c  PTR_O2               :: oxygen
39    C  PTR_FE               :: iron
40        INTEGER myIter        INTEGER myIter
41        _RL myTime        _RL myTime
42        INTEGER myThid        INTEGER myThid
# Line 50  C     == Routine arguments == Line 52  C     == Routine arguments ==
52    
53  #ifdef ALLOW_PTRACERS  #ifdef ALLOW_PTRACERS
54  #ifdef DIC_BIOTIC  #ifdef DIC_BIOTIC
55  C     == Local variables ==  
56    C !LOCAL VARIABLES: ====================================================
57    C  i,j,k                  :: loop indices
58    C  G*                     :: tendency term for the tracers
59    C  SURA                   :: tendency of alkalinity due to freshwater
60    C  SURC                   :: tendency of DIC due to air-sea exchange
61    C                            and virtual flux
62    C  SURO                   :: tendency of O2 due to air-sea exchange
63    C  BIO                    :: tendency of PO4 due to biological productivity,
64    C                            exchange with DOP pool and reminerization
65    C  CAR                    :: carbonate changes due to biological
66    C                             productivity and reminerization
67    C  bioac                  :: biological productivity
68    C  pflux                  :: changes to PO4 due to flux and reminerlization
69    c  cflux                  :: carbonate changes due to flux and reminerlization
70    c  freefe                 :: iron not bound to ligand
71        _RL  GDIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  GDIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
72        _RL  GALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  GALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
73        _RL  GPO4(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  GPO4(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
# Line 69  C     == Local variables == Line 86  C     == Local variables ==
86        _RL  freefe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  freefe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
87  #endif  #endif
88         INTEGER I,J,k         INTEGER I,J,k
89    CEOP
90    
91         DO k=1,Nr         DO k=1,Nr
92           DO j=1-OLy,sNy+OLy           DO j=1-OLy,sNy+OLy
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22