|
1 |
| -# Structure |
2 |
| -- analysis: Tools for analyzing divergences between rates at each depth, diversity of rates, |
3 |
| - and the profiles for rates and concs. |
| 1 | +# Mystic |
| 2 | + |
| 3 | +An inferred biomass biogeochemical model of Mystic Lake. |
| 4 | + |
| 5 | +- by: Scott W Olesen ( [email protected]) |
| 6 | +- development page: http://github.com/swo/mystic |
| 7 | +- main page: http://almlab.mit.edu/mystic.html |
| 8 | + |
| 9 | +## Dependencies |
| 10 | + |
| 11 | +- Matlab (developed with version 8) |
| 12 | +- Python (developed with version 2.7) |
| 13 | + |
| 14 | +## Installation |
| 15 | + $ git clone http://github.com/swo/mystic.git |
| 16 | + |
| 17 | +## File structure |
| 18 | +- analysis: Tools for analyzing divergences between rates at each depth, diversity of rates, and the profiles for rates and concs. |
4 | 19 | - bin: Main matlab scripts
|
5 |
| -- calibration: Tools for calibrating the model parameters using the measured rates and/or |
6 |
| - OTU data. |
7 |
| -- doc: Documentation of the model. |
8 |
| -- interactive: A script for running the model with easily-adjustable parameters and |
9 |
| - displaying the output. |
10 |
| -- sensitivity\_analysis: Tools for analyzing the sensitivity of the rate profiles to the |
11 |
| - values for each of the input parameters. |
12 |
| -- lake.cfg: Configuration for different tools, with lists of the names for processes and |
13 |
| - some default values from which to generate scripts for the various tools. |
14 |
| -- set\_matlab\_path.sh: Convenience tool for setting the matlab path when running scripts |
15 |
| - from the terminal. |
16 |
| - |
17 |
| -# Usage |
18 |
| - |
19 |
| -You should be able to run the sensitivity analysis in parallel. For running on 20 nodes on coyote, |
20 |
| - |
21 |
| -``` |
22 |
| -cat submit_0.sh | ssub -n 20 |
23 |
| -``` |
| 20 | +- calibration: Tools for calibrating the model parameters using the measured rates and/or OTU data. |
| 21 | +- interactive: A script for running the model with easily-adjustable parameters and displaying the output. |
| 22 | +- sensitivity\_analysis: Tools for analyzing the sensitivity of the rate profiles to the values for each of the input parameters. |
| 23 | +- lake.cfg: Configuration for different tools, with lists of the names for processes and some default values from which to generate scripts for the various tools. |
| 24 | +- set\_matlab\_path.sh: Convenience tool for setting the matlab path when running scripts from the terminal. |
| 25 | + |
| 26 | +## Usage |
| 27 | + |
| 28 | +### Interactive |
| 29 | +To run the model interactively, |
| 30 | +- Update `lake.cfg` with the desired parameters. |
| 31 | +- Under interactive, use `write_default_values_script.py` to create `run_interactive_defaults.m`, which supplies the values in `lake.cfg` to `interactive.m`. |
| 32 | +- In Matlab, run `run_interactive.m`. This requires the `bin/` folder to be on Matlab's path. If running in the terminal, you can use the `set_matlab_path.py` script to do this. |
| 33 | +- The output of the simulation is stored in `concs_history` and `rates_history` variables. You can write them to files using `bin/write_data_to_files.m`. |
| 34 | + |
| 35 | +### Analysis |
| 36 | +To make a plot of the Jensen-Shannon divergences between the composition of rates at each depth, run the `analysis/divergenes/pipeline.sh`. |
| 37 | + |
| 38 | +To compare the simulated Shannon diversity of composition of rates with some observed diversity of rates, run `analysis/diversity/get_diversity_depth.m`. |
| 39 | + |
| 40 | +To make timecourses and end-time plots for all the rates and chemical concentrations, run `analysis/profiles/pipeline.sh`. |
| 41 | + |
| 42 | + |
0 commit comments