/[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.11 by stephd, Wed Mar 27 20:51:06 2013 UTC revision 1.13 by stephd, Thu Mar 28 18:14:17 2013 UTC
# Line 1040  c sum all palatability*phyto and find ph Line 1040  c sum all palatability*phyto and find ph
1040  #endif  #endif
1041             do np=1,npmax             do np=1,npmax
1042              tmpz=max(0. _d 0,(allphyto(nz)-phygrazmin) )              tmpz=max(0. _d 0,(allphyto(nz)-phygrazmin) )
1043              grazphy(np,nz)=grazemax(nz)*              grazphy(np,nz)=grazemax(nz)*zooTempFunction(nz)*
1044  #ifdef SER_GRAZ  #ifdef SER_GRAZ
1045  c as in Vallina et al, 2011  c as in Vallina et al, 2011
1046       &          (((palat(np,nz)*phyto(np)/allphyto(nz))*phyto(np))/       &          (((palat(np,nz)*phyto(np)/allphyto(nz))*phyto(np))/
# Line 1337  c phospate uptake by each phytoplankton Line 1337  c phospate uptake by each phytoplankton
1337                 if (debug.eq.14) print*,'grow', grow(np), pcarbon(np)                 if (debug.eq.14) print*,'grow', grow(np), pcarbon(np)
1338  #ifdef DYNAMIC_CHL  #ifdef DYNAMIC_CHL
1339  c geider 97 for dChl/dt (source part) Eq. 3  c geider 97 for dChl/dt (source part) Eq. 3
1340                 if (acclim(np).gt. 0. _d 0) then                 if (acclim(np).gt. 0. _d 0.and.
1341         &                                 alpha_I(np).gt. 0. _d 0) then
1342                  rhochl(np)=chl2cmax(np) *                  rhochl(np)=chl2cmax(np) *
1343       &                      (grow(np)/(alpha_I(np)*acclim(np)) )       &                      (grow(np)/(alpha_I(np)*acclim(np)) )
1344                 else                 else
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22