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

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

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

revision 1.8 by stephd, Fri Jun 29 20:41:59 2012 UTC revision 1.9 by stephd, Thu Jul 26 18:01:22 2012 UTC
# Line 1456  c degradation of  CDOM - high when bleac Line 1456  c degradation of  CDOM - high when bleac
1456  #endif  #endif
1457  #endif  #endif
1458    
1459    #ifdef ALLOW_DENIT
1460               if (O2local.lt.O2crit) then
1461                  if (NO3local.lt.no3crit) then
1462    c no remineralization for N, P, Fe (not Si?)
1463                   DOPremin = 0. _d 0
1464                   DONremin = 0. _d 0
1465                   DOFeremin = 0. _d 0
1466                   preminP  = 0. _d 0
1467                   preminN  = 0. _d 0
1468                   preminFe =  0. _d 0
1469    #ifdef ALLOW_CDOM
1470                   preminP_cdom =  0. _d 0
1471                   preminN_cdom =  0. _d 0
1472                   preminFe_cdom =  0. _d 0
1473    #endif
1474    #ifdef ALLOW_CARBON
1475                  DOCremin = 0. _d 0
1476                  preminC  = 0. _d 0
1477    #ifdef ALLOW_CDOM
1478                  preminC_cdom =  0. _d 0
1479    #endif
1480    #endif
1481    
1482    #ifdef ALLOW_CDOM
1483    c degradation of  CDOM - high when bleached by light
1484               cdomp_degrd = reminTempFunction * cdomlocal
1485         &            *(cdombleach*min(PARlocal/PARcdom,1. _d 0) )
1486               cdomn_degrd = rnp_cdom * cdomp_degrd
1487               cdomfe_degrd= rfep_cdom * cdomp_degrd
1488    #ifdef ALLOW_CARBON
1489               cdomc_degrd = rcp_cdom  * cdomp_degrd
1490    #endif
1491    #endif
1492                  endif
1493               endif
1494    #endif
1495    c end denit caveats
1496    
1497  c chemistry  c chemistry
1498  c NH4 -> NO2 -> NO3 by bacterial action  c NH4 -> NO2 -> NO3 by bacterial action
1499              NO2prod = knita*( 1. _d 0-min(PARlocal/PAR0,1. _d 0) )              NO2prod = knita*( 1. _d 0-min(PARlocal/PAR0,1. _d 0) )
# Line 1635  c-ONLYNO3   dNO3dt =  - consumpNO3 +  pr Line 1673  c-ONLYNO3   dNO3dt =  - consumpNO3 +  pr
1673  #else  #else
1674       &                *(preminP  + DOPremin)       &                *(preminP  + DOPremin)
1675  #endif  #endif
1676                dNO3dt = dNO3dt - denit                dNO3dt = dNO3dt - (denit_no3/denit_np)*denit
1677                dNH4dt = dNH4dt -                  dNH4dt = dNH4dt -  
1678  #ifdef ALLOW_CDOM  #ifdef ALLOW_CDOM
1679       &                 (DONremin)       &                 (DONremin)
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22