/[MITgcm]/MITgcm_contrib/darwin2/pkg/monod/monod_forcing.F
ViewVC logotype

Diff of /MITgcm_contrib/darwin2/pkg/monod/monod_forcing.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.14 by jahn, Thu Oct 25 15:58:24 2012 UTC revision 1.16 by jahn, Fri Dec 27 17:29:00 2013 UTC
# Line 892  c taking from igsm dic run - fed through Line 892  c taking from igsm dic run - fed through
892  c               pCO2local=surfaceForcingT(i,j,bi,bj)  c               pCO2local=surfaceForcingT(i,j,bi,bj)
893  c or from darwin carbon module  c or from darwin carbon module
894  #ifdef ALLOW_CARBON  #ifdef ALLOW_CARBON
895    #ifdef pH_3D
896                   pCO2local=pCO2(i,j,k,bi,bj)
897    #else
898                 pCO2local=pCO2(i,j,bi,bj)                 pCO2local=pCO2(i,j,bi,bj)
899    #endif
900  #else  #else
901                 pCO2local=280. _d -6                 pCO2local=280. _d -6
902  #endif  #endif
# Line 1549  Coj            no Eu at surface (yet) Line 1553  Coj            no Eu at surface (yet)
1553       &                              FluxCO2(i,j,bi,bj)*dtplankton       &                              FluxCO2(i,j,bi,bj)*dtplankton
1554                SUROave(i,j,bi,bj)   =SUROave(i,j,bi,bj)+                SUROave(i,j,bi,bj)   =SUROave(i,j,bi,bj)+
1555       &                              flxO2(i,j)*dtplankton       &                              flxO2(i,j)*dtplankton
1556                 endif
1557    #ifdef pH_3D
1558                  pCO2ave(i,j,k,bi,bj)   =pCO2ave(i,j,k,bi,bj)+
1559         &                              pCO2(i,j,k,bi,bj)*dtplankton
1560                  pHave(i,j,k,bi,bj)     =pHave(i,j,k,bi,bj)+
1561         &                              pH(i,j,k,bi,bj)*dtplankton
1562    #else
1563                 if (k.eq.1) then
1564                pCO2ave(i,j,bi,bj)   =pCO2ave(i,j,bi,bj)+                pCO2ave(i,j,bi,bj)   =pCO2ave(i,j,bi,bj)+
1565       &                              pCO2(i,j,bi,bj)*dtplankton       &                              pCO2(i,j,bi,bj)*dtplankton
1566                pHave(i,j,bi,bj)     =pHave(i,j,bi,bj)+                pHave(i,j,bi,bj)     =pHave(i,j,bi,bj)+
1567       &                              pH(i,j,bi,bj)*dtplankton       &                              pH(i,j,bi,bj)*dtplankton
1568               endif               endif
1569  #endif  #endif
1570    #endif
1571            endif              endif  
1572  c end if hFac>0  c end if hFac>0
1573    
# Line 1699  c ANNA end TAVE Line 1712  c ANNA end TAVE
1712       &                         0,1,2,bi,bj,myThid )       &                         0,1,2,bi,bj,myThid )
1713          CALL DIAGNOSTICS_FILL( flxO2(1-Olx,1-Oly), 'DICOFLX ',          CALL DIAGNOSTICS_FILL( flxO2(1-Olx,1-Oly), 'DICOFLX ',
1714       &                         0,1,2,bi,bj,myThid )       &                         0,1,2,bi,bj,myThid )
1715    #ifdef pH_3D
1716            CALL DIAGNOSTICS_FILL( pCO2(1-Olx,1-Oly,1,bi,bj), 'DICPCO2 ',
1717         &                         0,Nr,2,bi,bj,myThid )
1718            CALL DIAGNOSTICS_FILL( pH(1-Olx,1-Oly,1,bi,bj), 'DICPHAV ',
1719         &                         0,Nr,2,bi,bj,myThid )
1720    #else
1721          CALL DIAGNOSTICS_FILL( pCO2(1-Olx,1-Oly,bi,bj), 'DICPCO2 ',          CALL DIAGNOSTICS_FILL( pCO2(1-Olx,1-Oly,bi,bj), 'DICPCO2 ',
1722       &                         0,1,2,bi,bj,myThid )       &                         0,1,2,bi,bj,myThid )
1723          CALL DIAGNOSTICS_FILL( pH(1-Olx,1-Oly,bi,bj), 'DICPHAV ',          CALL DIAGNOSTICS_FILL( pH(1-Olx,1-Oly,bi,bj), 'DICPHAV ',
1724       &                         0,1,2,bi,bj,myThid )       &                         0,1,2,bi,bj,myThid )
1725    #endif
1726  #endif /* ALLOW_CARBON */  #endif /* ALLOW_CARBON */
1727         ENDIF         ENDIF
1728  #endif /* ALLOW_DIAGNOSTICS */  #endif /* ALLOW_DIAGNOSTICS */
# Line 1715  c determine iron partitioning  - solve f Line 1735  c determine iron partitioning  - solve f
1735  c  c
1736  #ifdef ALLOW_TIMEAVE  #ifdef ALLOW_TIMEAVE
1737  c save averages  c save averages
1738         do k=1,nR           dar_timeave(bi,bj) = dar_timeave(bi,bj) + dtplankton
          dar_timeave(bi,bj,k)=dar_timeave(bi,bj,k)  
      &                         +dtplankton  
1739  #ifdef ALLOW_CARBON  #ifdef ALLOW_CARBON
1740           dic_timeave(bi,bj,k)=dic_timeave(bi,bj,k)           dic_timeave(bi,bj) = dic_timeave(bi,bj) + dtplankton
      &                         +dtplankton  
1741  #endif  #endif
        enddo  
1742  #endif  #endif
1743  c  c
1744  c -----------------------------------------------------  c -----------------------------------------------------
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22