Hydrobricks

Fixed

- Fixing the computation of the potential radiation for the temperature index melt model.

Added

- Adding a function to compute the units connectivity.
- Adding the cast shadow computation for the potential radiation.

Added

- Adding a new functionality to create a BehaviourLandCoverChange object from shapefiles of glacier extents.
- Adding a function to compute potential clear-sky direct solar radiation as defined by Hock (1999).
- Adding a temperature index melt model based on Hock (1999).
- Adding the possibility to discretize the hydrological units based on the mean annual potential radiation used in Hock (1999).
- Adding an aspect-related degree day factor.
- Adding a reference prediction based on a bootstrapping approach of the observed discharge.
- Adding the option to define a parameter shared by different processes (e.g., snow and ice melt).
- Allow calibrating a single parameter set for multiple catchments.
- Adding a Results class to simplify outputs analysis.
- Adding maps plotting and spatio-temporal animations of the model outputs.
- Handling yearly recurrent forcing with a 'day_of_year' approach.
- Moving the model structure definition from C++ to Python.

Changed

- Transitioning to a dictionary-based structure definition.
- Moving the parameters generation out of the Socont class.
- Moving the process parameters and forcing definitions to respective classes.
- Changing workflow in load_from_csv.
- Restructuring examples.
- An error is raised if some parameter values are missing.
- The parameter estimation and evaluation functions are restructured to a new trainer.py file.
- Adding the date/time to the log file name.

Fixed

- Socont quick discharge has been fixed... For real, this time.
- Fixing the SpotpySetup initialization.
- Fixing issue with hydro unit initialization in catchment discretization.
- Fixing a bug that occurred when the study area outline comprised multiple polygons.
- Adapting the FixLandCoverFractionsTotal() function to also work with cases where the hydro unit is fully glaciated at the beginning of the Behaviour time series.

Breaking changes

- The csv files defining the hydro units (e.g., elevation bands) properties have now to include the units of the values (e.g., m for elevation) below the column title.

Added

- Adding the computation of the slope and aspect from the DEM.
- Transferring slope and aspect values to the hydrological model (e.g., for Socont quick discharge using the unit slope).
- Any hydro unit property can now be extracted from data and provided to the hydrological model (flexible generic approach).
- Adding the computation of the mean latitude and longitude for each hydro unit and allow using it for the PET computation.
- Adding an extraction of the forcing using a faster weighted approach.
- Addition of an example script for the catchment data preparation.
- Allow loading unit ids from raster.
- Adding a hydro unit discretization by elevation and aspect (separately or combined; with an example script).

Changed

- Refactoring catchment attributes parsing.
- The csv files defining the hydro units (e.g., elevation bands) properties have now to include the units of the values (e.g., m for elevation) below the column title.

Fixed

- Socont quick discharge has been fixed.

Fixed

- Fixing an issue with the data container shape when using the spatialization from gridded data.

Breaking changes

- Many changes in the Forcing class:
- load_from_csv() was renamed to load_station_data_from_csv().
- define_spatialization() was renamed to spatialize_from_station_data() and is only meant for spatialization from station data.
- correct_station_data() was added and is to be used for applying a correction factor, for example.
- spatialize_from_gridded_data() was added to load data from gridded netCDF files.
- compute_pet() was added and uses the pyet package.
- The operations are not performed immediately, but only applied when apply_operations() is called, which is done automatically before the model run or before saving the forcing data to a netcdf file.
- The Catchment class is now part of the main module.

Added

- Adding the spatialization from gridded data with spatialize_from_gridded_data().
- Adding PET computation using the pyet library. New function: compute_pet().
- The forcing data can now be loaded from the dumped netCDF file with load_from().

Changed

- Refactoring of the Forcing and TimeSeries classes to separate 1D and 2D data.
- Enums are used to specify the Forcing variables.
- The Catchment class is now part of the main module.
- The Catchment class can generate a geotiff of the hydro unit ids with save_unit_ids_raster()

Fixed

- Fixing an issue with the elevation range condition in the elevation bands creation.