The function calculates chill units using the Utah model (Richardson
et al, 1974). This model is characterized by differential weighting of
temperature ranges, including negative weights for temperatures
above 15.9°C. This model recognizes that different temperatures
vary in effectiveness in accumulating chill as well as a negative
influence of high temperatures on previously accumulated chill.
Chill Units (Utah or Anderson model) perform better than chill hours
for a wider range of climates, and it could be considered as the
'reference' method nowadays, but it is ill-suited for warm or
Mediterranean conditions. To date, Chill portions is the best
existing model for most growing regions, so chill fulfilment
should preferably be calculated using that method, especially
when transferring varieties from one region to another, or in
studies on climate change.