Supplement 1. Fortran and Matlab files for conducting model calculations described in the main text.

File ListSupplement, original

calcB7_5curve.for (md5: 4d88e7313e582d7eb41e272f536467cd)
B7.inc (md5: 6b4d6683acd1765f93fce9fefd2b7ff3)
calcB7_5curve.exe (md5: c46092f613abb3cb6c0b69c735d132ef)
B_and_NPP_of_Trantime.m (md5: b02746500c7f393645567eb995a8997c)

Supplement, Revision 1 (posted 23 January 2013)

calcB7_5curve.for – corrected, (md5: b7233a525e671cf016281c86913aee91)
calcB7_5curve.exe – corrected, (md5: 42f4844387a0c2f5d9cfcdaa26c1d8a9)
B_and_NPP_of_Trantime.m – corrected, (md5: b31358b099d3038809795cacd1a9ff7a)

Description
The Fortran code “calcB7_5curve.for” and accompanying include file “B7.inc” contain the model used to test Hypothesis #1 in Lucas and Thompson (2012), i.e., to explore the relationship between phytoplankton biomass and productivity and water depth, and how that relationship changes as a function of benthic grazing rate. The include file, B7.inc, contains and describes all changeable model parameters. We include an executable version of the code “calcB7_5curve.exe”. However, if any changes are made to the include file, the code needs to be recompiled before running in order for changes to take effect. To produce the results presented in Lucas and Thompson (2012), the model was compiled with Compaq Visual Fortran Standard Edition 6.6.a. See Appendix and Table 1 in main manuscript and comments in source code for more detail.

The Matlab script “B_and_NPP_of_Trantime.m” contains the model used to test Hypothesis #2 in Lucas and Thompson (2012), i.e., to explore the relationship between phytoplankton biomass and productivity and transport time, and how that relationship changes as a function of benthic grazing rate. In this model, net phytoplankton growth rate, inflowing algal biomass concentration, zooplankton grazing rate, water depth, and cellular C:Chl ratio are specified by the user and assumed uniform and constant. In addition, a range of benthic grazing rates and transport times are specified. This model calculates steady-state habitat averaged phytoplankton biomass and productivity two ways: (1) via analytical expressions provided in the Appendix (see Eqs. 5 and 6) of Lucas and Thompson (2012); (2) via numerical integration along the streamwise habitat axis. This script also plots and compares results in various ways. See comments in script and Appendix and Table 1 in main manuscript for more detail.