20240215 : Henrique F. Duarte

**********************************************
ORIGINAL WASSI CODE: WASSICB_WEB_2_1_ORIGINAL/
**********************************************

WaSSI Ecosystem Services Model - Version 2.1 (September 2013)

Credits: 
USDA Forest Service Eastern Forest Environmental Threat Assessment Center
Web application (WaSSIweb) developed with programming support from Praecipio Consulting and the USDA Forest Service Geospatial Technology and Application Center (GTAC) through a partnership between the Eastern Forest Environmental Threat Assessment Center
https://wassiweb.fs.usda.gov/

Contacts:
Model development and applications
* Ge Sun (Ge.Sun@usda.gov)
* Steve McNulty (Steven.McNulty@usda.gov)
Model functionality and use
* Peter Caldwell (Peter.V.Caldwell@usda.gov)
Technical issues
* Bill Christie (William.M.Christie@usda.gov)
* Danika Mosher (Danika.Mosher@usda.gov)

****************************************
MODIFIED WASSI CODE: WASSICB_WEB_2_1HFD/
****************************************

Citation:
Duarte, H.F., Kim, J.B., Sun, G., McNulty, S., & Xiao, J. (2024) Future impacts of climate and vegetation change on water yield in the conterminous United States. Under review for Journal of Hydrology.

WASSICB_WEB_2_1HFD is based on WASSICB_WEB_2_1_ORIGINAL, with modifications by Henrique F. Duarte:
* Source code was adapted for compilation with gfortran in a Linux environment
* Numerical instability fixes
* Added support for dynamic land cover input (annual area fraction and monthly LAI for each land cover type)
* The original WaSSI considers a distinct equation for PAET (potential actual evapotranspiration) for HUC8s in the water resources regions 1, 2 (only HUC8s above 40deg of latitude), 4, and 5 in the Northeastern U.S. that has more than 20% of forest cover (PET.for). To avoid inconsistencies over time (in the case that the total forest cover oscillates above and below 20%), the forest cover conditional was removed. The distinct PAET equation was applied to all HUC8s in the water resources regions 1, 2 (only HUC8s above 40deg of latitude), 4, and 5, regardless of forest cover.

### The modified code currently supports only the 'US' region and climate inputs of the type 'US_XX' (customized input files) ###

The original code is provided in WASSICB_WEB_2_1_ORIGINAL/ for comparison.