/[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.17 by dfer, Fri Oct 26 21:08:13 2007 UTC revision 1.18 by dfer, Mon Oct 29 16:49:49 2007 UTC
# Line 71  C  SURO                   :: tendency of Line 71  C  SURO                   :: tendency of
71  C  GPO4                   :: tendency of PO4 due to biological productivity,  C  GPO4                   :: tendency of PO4 due to biological productivity,
72  C                            exchange with DOP pool and reminerization  C                            exchange with DOP pool and reminerization
73  C  CAR                    :: carbonate changes due to biological  C  CAR                    :: carbonate changes due to biological
74  C                             productivity and reminerization  C                             productivity and remineralization
75  C  BIOac                  :: biological productivity  C  BIOac                  :: biological productivity
76  C  pflux                  :: changes to PO4 due to flux and reminerlization  C  RDOP                   :: DOP sink due to remineralization
77  c  cflux                  :: carbonate changes due to flux and reminerlization  C  pflux                  :: changes to PO4 due to flux and remineralization
78  c  freefe                 :: iron not bound to ligand  C  CAR_S                  :: carbonate sink
79    C  cflux                  :: carbonate changes due to flux and remineralization
80    C  freefe                 :: iron not bound to ligand
81        _RL  GDIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  GDIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
82        _RL  GALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  GALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
83        _RL  GPO4(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  GPO4(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
# Line 85  c  freefe                 :: iron not bo Line 87  c  freefe                 :: iron not bo
87        _RL  SURO(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL  SURO(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
88        _RL  CAR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  CAR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
89        _RL  BIOac(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  BIOac(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
90          _RL  RDOP(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
91        _RL  pflux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  pflux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
92        _RL  exportflux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  exportflux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
93          _RL  CAR_S(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
94        _RL  cflux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  cflux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
95  #ifdef ALLOW_O2  #ifdef ALLOW_O2
96        _RL  GO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL  GO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
# Line 106  CEOP Line 110  CEOP
110         DO k=1,Nr         DO k=1,Nr
111           DO j=1-OLy,sNy+OLy           DO j=1-OLy,sNy+OLy
112            DO i=1-OLx,sNx+OLx            DO i=1-OLx,sNx+OLx
113               RDOP(i,j,k) =0. _d 0
114             GDIC(i,j,k) =0. _d 0             GDIC(i,j,k) =0. _d 0
115             GALK(i,j,k) =0. _d 0             GALK(i,j,k) =0. _d 0
116             GPO4(i,j,k) =0. _d 0             GPO4(i,j,k) =0. _d 0
# Line 115  CEOP Line 120  CEOP
120             pflux(i,j,k)   =0. _d 0             pflux(i,j,k)   =0. _d 0
121             exportflux(i,j,k)=0. _d 0             exportflux(i,j,k)=0. _d 0
122             cflux(i,j,k)   =0. _d 0             cflux(i,j,k)   =0. _d 0
123               CAR_S(i,j,k)   =0. _d 0
124  #ifdef ALLOW_O2  #ifdef ALLOW_O2
125             GO2(i,j,k)     =0. _d 0             GO2(i,j,k)     =0. _d 0
126  #endif  #endif
# Line 171  c flux of po4 from layers with biologica Line 177  c flux of po4 from layers with biologica
177       &                    bi,bj,imin,imax,jmin,jmax,       &                    bi,bj,imin,imax,jmin,jmax,
178       &                    myIter,myTime,myThid)       &                    myIter,myTime,myThid)
179    
180    C- Carbonate sink
181           DO k=1,Nr
182             DO j=jmin,jmax
183              DO i=imin,imax
184                 CAR_S(i,j,k)=BIOac(i,j,k)*R_CP*rain_ratio(i,j,bi,bj)*
185         &                    (1. _d 0-DOPfraction)
186              ENDDO
187             ENDDO
188           ENDDO
189    
190  c carbonate  c carbonate
191  #ifdef CAR_DISS  #ifdef CAR_DISS
192  c dissolution only below saturation horizon  c dissolution only below saturation horizon
# Line 188  c Line 204  c
204       &                    myIter,myTime,myThid)       &                    myIter,myTime,myThid)
205  #else  #else
206  c old OCMIP way  c old OCMIP way
207          CALL CAR_FLUX( BIOac, cflux,          CALL CAR_FLUX( CAR_S, cflux,
208       &                    bi,bj,imin,imax,jmin,jmax,       &                    bi,bj,imin,imax,jmin,jmax,
209       &                    myIter,myTime,myThid)       &                    myIter,myTime,myThid)
210  #endif  #endif
# Line 197  c add all tendencies for PO4, DOP, ALK, Line 213  c add all tendencies for PO4, DOP, ALK,
213         DO k=1,Nr         DO k=1,Nr
214           DO j=jmin,jmax           DO j=jmin,jmax
215            DO i=imin,imax            DO i=imin,imax
216               RDOP(i,j,k)= maskC(i,j,k,bi,bj)*KDOPRemin*PTR_DOP(i,j,k)
217    
218             GPO4(i,j,k)=-BIOac(i,j,k)+pflux(i,j,k)             GPO4(i,j,k)=-BIOac(i,j,k)+pflux(i,j,k) + RDOP(i,j,k)
      &         + maskC(i,j,k,bi,bj)*KDOPRemin*PTR_DOP(i,j,k)  
219    
220             car(i,j,k) =-BIOac(i,j,k)* R_CP*rain_ratio(i,j,bi,bj)*             car(i,j,k) = cflux(i,j,k) - CAR_S(i,j,k)
      &                (1. _d 0-DOPfraction)+cflux(i,j,k)  
221    
222             GDOP(i,j,k)=+BIOac(i,j,k)*DOPfraction             GDOP(i,j,k)=+BIOac(i,j,k)*DOPfraction - RDOP(i,j,k)
      &         - maskC(i,j,k,bi,bj)*KDOPRemin*PTR_DOP(i,j,k)  
223    
224             GALK(i,j,k)=+2. _d 0 *car(i,j,k)-R_NP*GPO4(i,j,k)             GALK(i,j,k)=+2. _d 0 *car(i,j,k)-R_NP*GPO4(i,j,k)
225    
# Line 219  c add all tendencies for PO4, DOP, ALK, Line 233  c add all tendencies for PO4, DOP, ALK,
233             endif             endif
234  #endif  #endif
235  #ifdef ALLOW_FE  #ifdef ALLOW_FE
236             GFE(i,j,k) =R_FeP*GPO4(i,j,k)             GFE(i,j,k) = R_FeP*GPO4(i,j,k)
237       &             -Kscav*freefe(i,j,k)       &                 -Kscav*freefe(i,j,k)
238  #endif  #endif
239            ENDDO            ENDDO
240           ENDDO           ENDDO
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22