Simplified model driver data and namelist files to test a General Lake Model particle tracking model, using Falling Creek Reservoir data.
1. 01_unstratified: a one-month run in January 2016, with a homogeneous temperature profile as the initial condition, constant air temperature, solar radiation, and humidity, and no precipitation, wind, or inflows/outflows.
Purpose: This scenario is designed to test vertical diffusion under still, unstratified conditions. Neutrally buoyant particles should slowly diffuse throughout the water column, similar to how a drop of dye spreads in a tub of still water.
2. 02_stratified: a one-month run in July-August 2015, with a stratified temperature profile as the initial condition, constant air temperature, solar radiation, and humidity, and no precipitation, wind, or inflows/outflows.
Purpose: This scenario is designed to test vertical diffusion under still, stratified conditions. Neutrally buoyant particles should accumulate in regions of low diffusivity (Kv), which typically occur at and below the thermocline.
3. 03_unstratified_wind: a one-month run in January 2016, with a homogeneous temperature profile, constant air temperature, solar radiation, and humidity, and no precipitation or inflows/outflows. This scenario includes artificial wind (six days no wind, then 1 day of wind, where each windy day is of increasing strength).
Purpose: This scenario is designed to test the movement of particles due to wind-driven mixing. In unstratified conditions, wind-driven mixing should cause particles to vertically diffuse more quickly over time in the surface layers, due to increases in turbulent diffusivity. Moreover, wind-driven mixing may lead to merging (and potentially post-wind-event splitting) of layers even in unstratified conditions. In this case, particles in layers that are merged should be redistributed throughout the newly created merged layer, and particles in a splitting layer should be reassigned to one of the two newly created split layers based on their elevation.
4. 04_stratified_wind: a one-month run in July-August 2015, with a stratified temperature profile, constant air temperature, solar radiation, and humidity, and no precipitation or inflows/outflows. This scenario includes artificial wind (six days no wind, then 1 day of wind, where each windy day is of increasing strength).
Purpose: This scenario is designed to test the movement of particles due to wind-driven mixing. In stratified conditions, wind should homogenize the temperature across surface layers, causing layers to merge. In this case, particles within the two layers that are merging should be randomly redistributed across the new merged layer. After a wind event, due to warm summer temperatures, layers should re-stratify and split. In this case, particles should be reassigned to the appropriate new layer ID based on their elevation.
5. 05_unstratified_inflow: a one-month run in January 2016, with a homogeneous temperature profile, constant air temperature, solar radiation, and humidity, and no precipitation or wind. This scenario includes one inflow/outflow of constant flow (cms) and temperature (the same as the water column). The lake starts half full and inflow is 5x outflow to fill the lake. Inflow is set to submerged and is colder than the water column so we know it will sink.
Purpose: This scenario is designed to test the movement of particles due to inflow discharge under weakly conditions, which includes (1) introduction of new particles through the inflow; and (2) changes in particle elevation due to corresponding water level elevation changes. Inflow to a lake in GLM may cause a layer to shift in elevation (raise or lower). In this case, the elevation of particles within that layer should be reassigned according to the newly assigned elevation of the layer. Under unstratified winter conditions, when the temperature of the inflow is colder than the lake and the inflow is submerged, the submerged inflow will remain submerged. Thus, particles above the inflow should have their height shifted upward as the layers above the inflow are shifted upwards.
6. 06_unstratified_outflow: a one-month run in January 2016, with a homogeneous temperature profile, constant air temperature, solar radiation, and humidity, and no precipitation or wind. This scenario includes one inflow/outflow of constant flow (cms) and temperature (the same as the water column). Outflow is 5x inflow to drain the lake. Outflow is set to submerged so we know it is withdrawing water from the bottom of the lake.
Purpose: This scenario is designed to test the movement of particles due to outflow under weakly stratified conditions, with a focus on changes in particle elevation due to corresponding water level elevation changes. Outflow from a lake in GLM may cause a layer to shift in elevation (raise or lower). In this case, the elevation of particles within that layer should be reassigned according to the newly assigned elevation of the layer. We expect particles above the submerged outflow to have their height shifted downward as the lake drains and the layers above the outflow are shifted downwards.
7. 07_stratified_inflow: a one-month run in July-August 2015, with a stratified temperature profile, constant air temperature, solar radiation, and humidity, and no precipitation or wind. This scenario includes one inflow/outflow of constant flow (cms) and temperature (about 15ºC cooler than the surface layer). The lake starts half full and inflow is 5x outflow to fill the lake. Inflow is set to submerged and is colder than the water column so we know it will sink.
Purpose: This scenario is designed to test the movement of particles due to inflow discharge under stratified conditions, which includes (1) introduction of new particles through the inflow; and (2) changes in particle elevation due to corresponding water level elevation changes. Inflow to a lake in GLM may cause a layer to shift in elevation (raise or lower). In this case, the elevation of particles within that layer should be reassigned according to the newly assigned elevation of the layer. Under stratified winter conditions, when the temperature of the inflow is cooler than the surface lake temperature, the surface inflow will plunge, leading to raising of layers above it and lowering of layers below it. Thus, particles introduced via the inflow should be assigned an elevation corresponding to the plunge elevation of the inflow. The expectation is that there will be more particles at the end of this simulation than at the beginning due to influx of particles from the inflow.
8. 08_stratified_outflow: a one-month run in July-August 2015, with a stratified temperature profile, constant air temperature, solar radiation, and humidity, and no precipitation or wind. This scenario includes one inflow/outflow of constant flow (cms) and temperature (about 15ºC cooler than the surface layer). Outflow is 5x inflow to drain the lake. Outflow is set to submerged so we know it is withdrawing water from the bottom of the lake.
Purpose: This scenario is designed to test the movement of particles due to outflow under stratified conditions, with a focus on changes in particle elevation due to corresponding water level elevation changes. Outflow from a lake in GLM may cause a layer to shift in elevation (raise or lower). In this case, the elevation of particles within that layer should be reassigned according to the newly assigned elevation of the layer. We expect particles above the submerged outflow to have their height shifted downward as the lake drains and the layers above the outflow are shifted downwards.
9. 09_unstratified_observed_wind_inflow: a one-month run in January 2016, with a homogeneous temperature profile and observed meteorology and inflow data (from the primary inflow only); outflow data are set to match inflow data to maintain the reservoir at full pond.
Purpose: this scenario is designed to test the movement of particles under realistic (variable) meteorology and inflow conditions applied to an unstratified lake.
10. 10_stratified_observed_wind_inflow: a one-month run in July-August 2015, with a stratified temperature profile and observed meteorology and inflow data (from the primary inflow only); outflow data are set to match inflow data to maintain the reservoir at full pond.
Purpose: this scenario is designed to test the movement of particles under realistic (variable) meteorology and inflow conditions applied to a stratified lake.
-
Download VSCode here.
-
Clone latest version of GLM-AED; see instructions here.
-
Install homebrew if you haven't; see instructions here.
-
Make sure you have the following libraries: gcc (fortran compiler); gd (graphics library); netcdf (to handle output); these can be installed with the command:
brew install LIBRARY_NAME
- Compile GLM. To do this you will need to navigate to the glm-source folder and build GLM:
cd glm-aed-ptm/glm-source
./build_glm.sh
If you run into an error and subsequently wish to compile GLM in debug mode, run the following:
clean_build.sh
./build_glm.sh --debug
- Then run GLM. You will need to point to both the folder where you have the .nml and driver files for your example as well as the folder containing the GLM executable you wish to run.
Navigate to directory that houses the model scenario you wish to run (be sure to check the filepath is correct for your directory):
cd RProjects/ptm-examples/1_unstratified
Run the GLM executable (be sure to check the filepath is correct for your executable):
/Users/MaryLofton/glm-aed-ptm/glm-source/GLM/glm
- If you edit the source code or download a new version of GLM, you need will to re-compile the model (Step 5) again before running it (Step 6)
Plots of the initial conditions and driver data for each scenario can be found in the plots folder, under the sub-folder for each scenario. An example IC_driver.png plot for Scenario 1 (1_unstratified) is shown below:
