Changes in version 0.8.3 o Rd manual files improved with respect to readability. Changes in version 0.8.2 (2022-08-22) o Rd files re-generated using latest version of roxygen2 in response to CRAN request. Changes in version 0.8.1 (2022-05-01) o code examples in vignette updated o example in wind.profile() coded more efficiently Changes in version 0.8.0 o new function optimum.temperature() added by Alison Bennett. Calculates optimum temperature of GPP (or any other flux) using boundary line analysis. Changes in version 0.7.6 o new function roughness.length.heat() added. Calculates roughness length for heat (z0h). o new output z0h added to function aerodynamic.conductance(). New test function for aerodynamic conductance. o bug fix in stability.correction() related to NA handling. New test function for stability.correction() o new roxygen2 version Changes in version 0.7.5 o growing season filter within filter.data() adjusted so that it can handle time series that start or end with incomplete days (i.e. time series that start/end in the middle of the day). Changes in version 0.7.4 o Solver functions within dew.point() and wetbulb.temp() not called if any of the inputs is NA, NaN, or Inf. In this case the function returns NA (and not the upper bound as in previous versions). Changes in version 0.7.3 o Factor 0.93 in wetbulb.temp documented and added to list of constants as 'Le067', defined as the Lewis number to the power of 0.67. Changes in version 0.7.2 o z0m as input to Choudhury_1988 and Su_2001 functions for boundary layer conductance. Avoids the implicit assumptions that z0m is constant over the entire time period. If z0m is not provided as input, the assumption of a constant z0m is still made and z0m is estimated from the roughness.parameters (method='wind_profile') function. o wind.profile() function rewritten: heights (argument z) for which wind speed is calculated must have the same length as the other input data (or of length 1). If wind speed at multiple height needs to be calculated, this needs to be done outside the function (e.g. as a for loop as in the function example). Changes in version 0.7.1 (2019-06-08) o package rebuilt with knitr version 1.23 (fixed encoding issue) Changes in version 0.7.0 (2019-03-20) o new functions added: extraterrestrial.radiation() and potential.radiation() o vignette updated o function e.to.rH() added to unit conversions o function kg.to.mol() added to unit conversions (not yet in use) o solar_constant and H2Omol added to bigleaf.constants o CITATION file added Changes in version 0.6.5 (2018-07-04) o Vignette "Tutorial on bigleaf" added Changes in version 0.6.4 o bug fix to energy.closure(): G and S were previously always ignored in the calculation of the intercept and slope. Changes in version 0.6.3 o rownames corrected in the example datasets (starting at 1 now) o second part of output in filter.data() function is only shown if at least one variable is actually filtered. o default values for several arguments changed (e.g. "Ga" to "Ga_h") to ensure internal consistency Changes in version 0.6.2 o hard-coded conversion factors replaced by constants (e.g. g2kg, Pa2kPa, umol2mol, etc.) which were added to the bigleaf.constants() list. Changes in version 0.6.1 (2018-04-21) o default names for 'formulation' arguments changed for clarity in surface.conductance() and decoupling(). "PenmanMonteith" is changed to "Penman-Monteith", "FluxGradient" to "Flux-Gradient", and "JarvisMcNaughton_1986" to "Jarvis&McNaughton_1986". Changes in version 0.6.0 o function reference.ET() is now covered by potential.ET(). If 'approach=Penman-Monteith', potential ET is calculated according to the Penman-Monteith equation (as previously in reference.ET()). Argument 'Gs_pot' (formerly 'Gs_ref') is given in mol m-2 s-1 and defaults to 0.6 mol m-2 s-1. Changes in version 0.5.2 (2018-04-10) o Gb output in aerodynamic.conductance and Gb functions renamed to Gb_h to clarify that values correspond to heat transfer. Same for Rb and kB. Changes in version 0.5.1 o results of additional quantities (for which Sc is provided) are expressed as conductances instead of resistances in aerodynamic.conductance() o 'data' included as argument in radiometric.surface.temp() o calculation of 'nr_data' accounts for missing values in Gamma if model=Leuning (function stomatal.slope()) Changes in version 0.5.0 o wind.profile() function restructured. z0m can be provided or calculated from the logarithmic wind profile equation o reflectance of downwelling longwave radiation is now accounted for in the calculation of radiometric surface temperature (function radiometric.surface.temp()) o variable 'LW_down' added to dataset DE_Tha_Jun_2014 o argument 'calc.surface.CO2' removed from intercellular.CO2() Changes in version 0.4.6 o important bugfix concerning several functions: all required arguments are passed to nested functions o argument 'N' in aerodynamic.conductance() now defaults to 2 as in Gb.Su() Changes in version 0.4.5 o variables 'LW_up' (upward longwave radiation) and 'Reco' (ecosystem respiration) were added to all three example datasets o bug fix to 'surface.conditions' and 'pressure.from.elevation' regarding the 'Esat.formulation' argument. In 'Esat.formulation', this argument is not implemented and the default formulation is taken for all calculations. Changes in version 0.4.4 o 'reference.ET': argument Gs renamed to Gs_ref Changes in version 0.4.3 o Function 'bigleaf.constants' restructured in a way that allows changing the constants for each function call Changes in version 0.4.2 o constants (a,b,c) from Allen et al. 1998 added as option to 'Esat.slope' (pull request from Rémi Vezy) o Argument 'Esat.formula' added to every function that calculates esat or slope of esat to allow consistency among functions. o 'light.response': argument "..." added to the nls function o 'stomatal.slope': argument "..." within each call of nls in the function Changes in version 0.4.1 o aerodynamic.conductance: Ga_CO2 added to function output Changes in version 0.4.0 o 'aerodynamic.conductance': Rb can be calculated for other quantities if the respective Schmidt number is provided. I.e. new arguments "Sc" and "Sc_name" as the value of the Schmidt number, and the name of the quantity for which Sc is provided, respectively. o 'bigleaf.constants': Prandtl number (Pr) and Schmidt number for CO2 (Sc_CO2) added o 'light.response': bug fix: +Reco replaced by -Reco (sign was reversed) Changes in version 0.3.2 o rH.to.VPD: if statement vectorized (thanks to Rémi Vezy) o radiometric.surface.temp: Output renamed ("." replaced by "_") Changes in version 0.3.1 o Monin.Obukhov.length: default arguments added