search

pschmidtwalter / LWFBrook90R

Run the LWF-Brook90 hydrological model within R.
https://pschmidtwalter.github.io/LWFBrook90R/
11 stars 7 forks source link

Issue with interception parameters in setparam_LWFB90 #25

Closed rhabel closed 4 years ago

rhabel commented 4 years ago

Hello, there's a bug in the default settings of setparam_LWFB90():

The interception parameter's defaults are described as

frintlai | Intercepted fraction of rain per unit LAI. Default: 0.06 frintsai | Intercepted fraction of rain per unit SAI. Default: 0.04 fsintlai | Intercepted fraction of snow per unit LAI. Default: 0.06 fsintsai | Intercepted fraction of snow per unit SAI. Default: 0.04

However, according to the Brook90 documentation (http://www.ecoshift.net/brook/b90doc.html - Parameters and variables): FRINTL & FRINTS are 0.06 and FSINTL and FSINTS are 0.04

"FRINTL and FRINTS (Fixed parameters) - intercepted fraction of rain per unit LAI and per unit SAI respectively, dimensionless. See also FSINTL. FRINTL and FRINTS are both fixed at 0.06. ... FSINTL and FSINTS (Fixed parameters) - intercepted fraction of snow per unit LAI and per unit SAI respectively, dimensionless. See also FRINTL. FSINTL and FSINTS are both fixed at 0.04.

Best regards, Raphael