/[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.7 by stephd, Thu Jun 28 20:36:11 2012 UTC revision 1.8 by stephd, Fri Jun 29 20:41:59 2012 UTC
# Line 48  c ANNA pass extra variables if WAVEBANDS Line 48  c ANNA pass extra variables if WAVEBANDS
48       I                       phytoup, popuplocal, ponuplocal,       I                       phytoup, popuplocal, ponuplocal,
49       I                       pofeuplocal, psiuplocal,       I                       pofeuplocal, psiuplocal,
50       I                       PARlocal,Tlocal, Slocal,       I                       PARlocal,Tlocal, Slocal,
51         I                       pCO2local,
52       I                       freefelocal, inputFelocal,       I                       freefelocal, inputFelocal,
53       I                       bottom, dzlocal,       I                       bottom, dzlocal,
54       O                       Rstarlocal, RNstarlocal,       O                       Rstarlocal, RNstarlocal,
# Line 148  c zoo   = zooplankton Line 149  c zoo   = zooplankton
149           _RL PARlocal           _RL PARlocal
150           _RL Tlocal           _RL Tlocal
151           _RL Slocal           _RL Slocal
152             _RL pCO2local
153           _RL freefelocal           _RL freefelocal
154           _RL inputFelocal           _RL inputFelocal
155           _RL bottom           _RL bottom
# Line 247  c phytoplankton specific nutrient limita Line 249  c phytoplankton specific nutrient limita
249           _RL limit(npmax)           _RL limit(npmax)
250  c phytoplankton light limitation term  c phytoplankton light limitation term
251           _RL ilimit(npmax)           _RL ilimit(npmax)
252             _RL pCO2limit(npmax)
253           _RL ngrow(npmax)           _RL ngrow(npmax)
254           _RL grow(npmax)           _RL grow(npmax)
255           _RL PspecificPO4(npmax)           _RL PspecificPO4(npmax)
# Line 636  c using Platt-like equations with inhibi Line 639  c using Platt-like equations with inhibi
639  #endif  #endif
640  c ANNA endif  c ANNA endif
641    
642    c pCO2 limit - default to non
643             do np=1,npmax
644               pCO2limit(np)=1. _d 0
645    c          if (np.eq.6) then
646    c             pCO2limit(np)=1. _d 0 + (pCO2local-400. _d -6)/600 _d -6
647    c             pCO2limit(np)=max(pCO2limit(np),1. _d 0)
648    c             pCO2limit(np)=min(pCO2limit(np),2. _d 0)
649    c          endif
650               if (debug.eq.15) print*,'pco2limit',pCO2limit(np),pCO2local
651             enddo
652    
653    
654  ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc  ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
655  c Determine phytoplankton nutrient limitation as mimimum of  c Determine phytoplankton nutrient limitation as mimimum of
656  c P,N,Si,Fe. However N can be utilized in several forms, so  c P,N,Si,Fe. However N can be utilized in several forms, so
# Line 853  c ANNA endif Line 868  c ANNA endif
868  #endif  #endif
869    
870            do np = 1, npmax            do np = 1, npmax
871             pcm(np)=pcmax(np)*limit(np)*phytoTempFunction(np)             pcm(np)=pcmax(np)*limit(np)*phytoTempFunction(np)*
872         &                                         pCO2limit(np)
873  #ifdef DYNAMIC_CHL  #ifdef DYNAMIC_CHL
874             if (phyto(np).gt. 0. _d 0) then             if (phyto(np).gt. 0. _d 0) then
875               chl2c(np)=phychl(np)/(phyto(np)*R_PC(np))               chl2c(np)=phychl(np)/(phyto(np)*R_PC(np))
# Line 930  c         acclim(np)=max(acclim(np),chl2 Line 946  c         acclim(np)=max(acclim(np),chl2
946  c diagnostic version of the above that does not feed back to growth  c diagnostic version of the above that does not feed back to growth
947            ChlGeiderlocal = 0. _d 0            ChlGeiderlocal = 0. _d 0
948            do np = 1, npmax            do np = 1, npmax
949             tmppcm = mu(np)*limit(np)*phytoTempFunction(np)             tmppcm = mu(np)*limit(np)*phytoTempFunction(np)*
950         &                                       pCO2limit(np)
951             if (tmppcm.gt.0.d0) then             if (tmppcm.gt.0.d0) then
952               tmpchl2c = Geider_chl2cmax(np)/               tmpchl2c = Geider_chl2cmax(np)/
953       &         (1+(Geider_chl2cmax(np)*Geider_alphachl(np)*PARdaylocal)/       &         (1+(Geider_chl2cmax(np)*Geider_alphachl(np)*PARdaylocal)/
# Line 1290  c consumption - sum of phytoplankton con Line 1307  c consumption - sum of phytoplankton con
1307  c phospate uptake by each phytoplankton  c phospate uptake by each phytoplankton
1308  #ifndef GEIDER  #ifndef GEIDER
1309                 grow(np)=ngrow(np)*mu(np)*limit(np)*ilimit(np)*                 grow(np)=ngrow(np)*mu(np)*limit(np)*ilimit(np)*
1310       &                            phytoTempFunction(np)       &                            phytoTempFunction(np)*pCO2limit(np)
1311  #endif  #endif
1312  #ifdef GEIDER  #ifdef GEIDER
1313                 grow(np)=ngrow(np)*pcarbon(np)                 grow(np)=ngrow(np)*pcarbon(np)
# Line 1506  c%%% Line 1523  c%%%
1523  #ifdef DAR_DIAG_RSTAR  #ifdef DAR_DIAG_RSTAR
1524  #ifndef GEIDER  #ifndef GEIDER
1525                tmpgrow=ngrow(np)*mu(np)*ilimit(np)*                tmpgrow=ngrow(np)*mu(np)*ilimit(np)*
1526       &              phytoTempFunction(np)       &              phytoTempFunction(np)*pCO2limit(np)
1527  #endif  #endif
1528  #ifdef GEIDER  #ifdef GEIDER
1529                tmpgrow=grow(np)/limit(np)                tmpgrow=grow(np)/limit(np)
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22