Skip to content

A heat and mass balance model, written in fortran (90).

Notifications You must be signed in to change notification settings

rniemeyer07/Reservoir_Model

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

Two-Layer Reservoir Model

  • A simple two-layer heat and mass balance reservoir model, written in fortran (90)
  • Authors: Ryan Niemeyer, Yifan Cheng, John Yearsley, Yixin Mao, and Bart Nijssen
  • University of Washington - [Computational Hydrology Lab] (http://uw-hydro.github.io/) ([GitHub site] (https://github.com/UW-Hydro))

####Model Overview: This repository is public source code for a simple two-layer reservoir model. The model does a simple mass and energy balance for a reservoir with two distinct layers: the epilimnion (surface) and hypolimnion (below surface). Inflow is partitioned based on the calculated densities of the incoming water, epilimnion, and hypolimnion. The net surface energy is calculated with VIC-derived incoming energy and an energy subroutine from RBM. The core subroutine (reservoir_subroutine.f90) calculates the change in temperature of each layer based on net change in energy from advection, surface energy (epilimnoin only), and diffusion (see conceptual model below). A more complete documentation of model development and equations can be found here. The model allows for measured releases from a spillway (epilimnion) and penstock/sluiceway (hypolimnion) to be uploaded. The model code comes with no guarantees, expressed or implied, as to suitability, completeness, accuracy, or any claim you would like to make.

####Conceptual model of stratified reservoir: ####Conceptual model of two-layer mass and energy fluxes*:

*parameters (i.e. As, Kz, etc.) are defined in the "Documentation" document

####File descriptions:

  • Block_Energy.f90: parameters for energy subroutine (originated from RBM Block)
  • Block_Reservoir.f90: parameters for the model
  • Block_Flow.f90: parameters for the "flow_subroutine.f90" (continuity/mass balance equation)
  • reservoir.f90: main program that calls each subroutine
  • density_subroutine.f90: calculates density of incoming river flow, epilimnion, and hypolimnion
  • flow_subroutine.f90: continuity/mass balance equation - calculates the flow into hypolimnion and hypolimnion (based on the density) and calculates change in volume of each layer
  • Energy.f90: RBM's energy subroutine to calculate surface energy exchange
  • reservoir_subroutine.f90: calculate transfer of energy between layers through advection, diffusion, and surface energy
  • input_file: file path of line 1- comments, line 2- flow in to reservoir, line 3-flow out reservoir, line 4-energy flux (from VIC), line 5-measured stream temperatures, line 6-measured releases
  • reservoir_file: file with information about each reservoir to be modeled. line 1-column names, line 2 to n-rows of data. Start node is upstream point (where water enters the reservoir), end node is downstream point (where water exits the reservoir)
  • makefile: make file with code to compile all the fortran code and to remove the compiled fortran code

####Model code flow diagram: alt text

####To run the model: Note: This must be run from the terminal while in the folder with all the fortran code. You must have gfortran installed on your machine. This also assumes input_file and reservoir_file are ready and all files for the input_file are present.

  1. enter "make clean" - remove any old compiled fortran code
  2. enter "make" - compile the fortran code
  3. enter "./reservoir input_file reservoir_file"

About

A heat and mass balance model, written in fortran (90).

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published