Ik heb zelf manueel even wat berekeningen gedaan voor de situatie die je aanhaalde rond de middag:
G0 = 1143 = 1322*sin(59.84)
Turbiditeitsfactor (Tl) = 4.36
Ggrens = G0 * (0.22 + exp(-4.36/(0.9+9.4*sin(59.84))) = 957
Gmax = 861, dus Gmax > 0.4*G0 (=457)
Gmin = 851, dus Gmin < Ggrens.
Gmax < Ggrens, dus moet de fractionele berekening gebruikt worden.
In scenario d is de turbiditeitsfactor 13
fractie zon = ((Ggem/G0) - 1.27*Gmin/G0)/exp(-13/(0.9+9.4*sin(h)) < 0, dus te herleiden naar 0.
Hoewel ik een lagere grenswaarde vind, is die naar mijn mening nog steeds hoog. 957 W/m² bereiken we nooit bij helder weer. 900 is zowat het hoogst haalbare. Het resultaat "geen zon" kan niet kloppen.
Bij het algoritme Slob kom ik uit op Ggrens = 720, en dat lijkt me wel realistisch.
Ggrens = G0 * [0,3 + exp {- TL/(0,9 + 9,4sinuszonshoogte)}] met TL = 10
Quote selectie