/[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.1 by jahn, Wed Apr 13 18:56:25 2011 UTC revision 1.4 by jahn, Mon Nov 7 21:55:11 2011 UTC
# Line 248  c tendencies Line 248  c tendencies
248         _RL do2l         _RL do2l
249         _RL dZooCl(nzmax)         _RL dZooCl(nzmax)
250  c air-sea fluxes  c air-sea fluxes
251         _RL flxCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)         _RL flxCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
252         _RL flxALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)         _RL flxALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
253         _RL flxO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)         _RL flxO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
254  #endif  #endif
255    
256         _RL tot_Nfix         _RL tot_Nfix
# Line 283  c Line 283  c
283  COJ for diagnostics  COJ for diagnostics
284             PParr(i,j,k) = 0. _d 0             PParr(i,j,k) = 0. _d 0
285             Nfixarr(i,j,k) = 0. _d 0             Nfixarr(i,j,k) = 0. _d 0
286    #ifdef DAR_DIAG_CHL
287               GeiderChlarr(i,j,k) = 0. _d 0
288               GeiderChl2Carr(i,j,k) = 0. _d 0
289               DoneyChlarr(i,j,k) = 0. _d 0
290               DoneyChl2Carr(i,j,k) = 0. _d 0
291               CloernChlarr(i,j,k) = 0. _d 0
292               CloernChl2Carr(i,j,k) = 0. _d 0
293    #endif
294  c ANNA_TAVE  c ANNA_TAVE
295  #ifdef WAVES_DIAG_PCHL  #ifdef WAVES_DIAG_PCHL
296             DO np=1,npmax             DO np=1,npmax
# Line 931  c            endif Line 939  c            endif
939  c  c
940  #ifdef IRON_SED_SOURCE  #ifdef IRON_SED_SOURCE
941  c only above minimum depth (continental shelf)  c only above minimum depth (continental shelf)
942               if (rF(k).lt.depthfesed) then               if (rF(k).gt.-depthfesed) then
943  c only if bottom layer  c only if bottom layer
944                 if (bottom.eq.1.0 _d 0) then                 if (bottom.eq.1.0 _d 0) then
945  #ifdef IRON_SED_SOURCE_VARIABLE  #ifdef IRON_SED_SOURCE_VARIABLE
# Line 966  c Line 974  c
974               picupl = PICl               picupl = PICl
975  c include surface forcing  c include surface forcing
976               if (k.eq.1) then               if (k.eq.1) then
977                ddicl  =  ddicl  + flxCO2(i,j,bi,bj)                ddicl  =  ddicl  + flxCO2(i,j)
978                dalkl  =  dalkl  + flxALK(i,j,bi,bj)                dalkl  =  dalkl  + flxALK(i,j)
979                do2l   =  do2l   + flxO2(i,j,bi,bj)                do2l   =  do2l   + flxO2(i,j)
980               endif               endif
981  #endif  #endif
982  c  c
# Line 1331  Coj            no Eu at surface (yet) Line 1339  Coj            no Eu at surface (yet)
1339  #ifdef ALLOW_CARBON  #ifdef ALLOW_CARBON
1340               if (k.eq.1) then               if (k.eq.1) then
1341                SURave(i,j,bi,bj)    =SURave(i,j,bi,bj)+                SURave(i,j,bi,bj)    =SURave(i,j,bi,bj)+
1342       &                              flxCO2(i,j,bi,bj)*dtplankton       &                              flxCO2(i,j)*dtplankton
1343                SURCave(i,j,bi,bj)    =SURCave(i,j,bi,bj)+                SURCave(i,j,bi,bj)    =SURCave(i,j,bi,bj)+
1344       &                              FluxCO2(i,j,bi,bj)*dtplankton       &                              FluxCO2(i,j,bi,bj)*dtplankton
1345                SUROave(i,j,bi,bj)   =SUROave(i,j,bi,bj)+                SUROave(i,j,bi,bj)   =SUROave(i,j,bi,bj)+
1346       &                              flxO2(i,j,bi,bj)*dtplankton       &                              flxO2(i,j)*dtplankton
1347                pCO2ave(i,j,bi,bj)   =pCO2ave(i,j,bi,bj)+                pCO2ave(i,j,bi,bj)   =pCO2ave(i,j,bi,bj)+
1348       &                              pCO2(i,j,bi,bj)*dtplankton       &                              pCO2(i,j,bi,bj)*dtplankton
1349                pHave(i,j,bi,bj)     =pHave(i,j,bi,bj)+                pHave(i,j,bi,bj)     =pHave(i,j,bi,bj)+
# Line 1463  c ANNA end TAVE Line 1471  c ANNA end TAVE
1471       &                         0,Nr,2,bi,bj,myThid )       &                         0,Nr,2,bi,bj,myThid )
1472  #endif  #endif
1473  #ifdef ALLOW_CARBON  #ifdef ALLOW_CARBON
1474          CALL DIAGNOSTICS_FILL( flxCO2(1-Olx,1-Oly,bi,bj), 'DICTFLX ',          CALL DIAGNOSTICS_FILL( flxCO2(1-Olx,1-Oly), 'DICTFLX ',
1475       &                         0,1,2,bi,bj,myThid )       &                         0,1,2,bi,bj,myThid )
1476          CALL DIAGNOSTICS_FILL( FluxCO2(1-Olx,1-Oly,bi,bj), 'DICCFLX ',          CALL DIAGNOSTICS_FILL( FluxCO2(1-Olx,1-Oly,bi,bj), 'DICCFLX ',
1477       &                         0,1,2,bi,bj,myThid )       &                         0,1,2,bi,bj,myThid )
1478          CALL DIAGNOSTICS_FILL( flxO2(1-Olx,1-Oly,bi,bj), 'DICOFLX ',          CALL DIAGNOSTICS_FILL( flxO2(1-Olx,1-Oly), 'DICOFLX ',
1479       &                         0,1,2,bi,bj,myThid )       &                         0,1,2,bi,bj,myThid )
1480          CALL DIAGNOSTICS_FILL( pCO2(1-Olx,1-Oly,bi,bj), 'DICPCO2 ',          CALL DIAGNOSTICS_FILL( pCO2(1-Olx,1-Oly,bi,bj), 'DICPCO2 ',
1481       &                         0,1,2,bi,bj,myThid )       &                         0,1,2,bi,bj,myThid )
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22