# NOT RUN {
library(tealeaves)
cs <- make_constants()
ep <- make_enviropar()
lp <- make_leafpar()
ep$T_sky <- ep$T_sky(ep)
T_leaf <- set_units(298.15, K)
energy_balance(T_leaf, lp, ep, cs, FALSE, TRUE, TRUE)
# }
Run the code above in your browser using DataLab