/[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.10 by stephd, Thu Oct 13 16:25:12 2005 UTC revision 1.11 by stephd, Fri Dec 16 21:07:53 2005 UTC
# Line 92  c  freefe                 :: iron not bo Line 92  c  freefe                 :: iron not bo
92         INTEGER I,J,k         INTEGER I,J,k
93         INTEGER nCALCITEstep         INTEGER nCALCITEstep
94  CEOP  CEOP
95           jmin=1
96           jmax=sNy
97           imin=1
98           imax=sNx
99    
100         DO k=1,Nr         DO k=1,Nr
101           DO j=1-OLy,sNy+OLy           DO j=1-OLy,sNy+OLy
# Line 177  c old OCMIP way Line 181  c old OCMIP way
181    
182  c add all tendencies for PO4, DOP, ALK, DIC  c add all tendencies for PO4, DOP, ALK, DIC
183         DO k=1,Nr         DO k=1,Nr
184           DO j=1-OLy,sNy+OLy           DO j=jmin,jmax
185            DO i=1-OLx,sNx+OLx            DO i=imin,imax
186             bio(i,j,k)=-bioac(i,j,k)+pflux(i,j,k)             bio(i,j,k)=-bioac(i,j,k)+pflux(i,j,k)
187       &         + maskC(i,j,k,bi,bj)*Kdopremin*PTR_DOP(i,j,k)       &         + maskC(i,j,k,bi,bj)*Kdopremin*PTR_DOP(i,j,k)
188             car(i,j,k)=-bioac(i,j,k)* R_cp*rain_ratio(i,j,bi,bj)*             car(i,j,k)=-bioac(i,j,k)* R_cp*rain_ratio(i,j,bi,bj)*
# Line 215  c add all tendencies for PO4, DOP, ALK, Line 219  c add all tendencies for PO4, DOP, ALK,
219    
220  C update  C update
221         DO k=1,Nr         DO k=1,Nr
222           DO j=1-OLy,sNy+OLy           DO j=jmin,jmax
223            DO i=1-OLx,sNx+OLx            DO i=imin,imax
224             PTR_DIC(i,j,k)=             PTR_DIC(i,j,k)=
225       &      PTR_DIC(i,j,k)+GDIC(i,j,k)*dTtracerLev(k)       &      PTR_DIC(i,j,k)+GDIC(i,j,k)*dTtracerLev(k)
226             PTR_ALK(i,j,k)=             PTR_ALK(i,j,k)=
# Line 247  c find free iron and get rid of insolubl Line 251  c find free iron and get rid of insolubl
251  #ifdef ALLOW_TIMEAVE  #ifdef ALLOW_TIMEAVE
252  c save averages  c save averages
253        DO k=1,Nr        DO k=1,Nr
254           DO j=1-OLy,sNy+OLy           DO j=jmin,jmax
255            DO i=1-OLx,sNx+OLx            DO i=imin,imax
256              BIOave(i,j,k,bi,bj)=BIOave(i,j,k,bi,bj)+              BIOave(i,j,k,bi,bj)=BIOave(i,j,k,bi,bj)+
257       &                          BIOac(i,j,k)*deltaTclock       &                          BIOac(i,j,k)*deltaTclock
258              CARave(i,j,k,bi,bj)=CARave(i,j,k,bi,bj)+              CARave(i,j,k,bi,bj)=CARave(i,j,k,bi,bj)+
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22