Ik vermoed dat je inderdaad waarden die buiten het interval 0-1 vallen moet afronden naar 0 of 1 naargelang waar ze erbuiten vallen.
Ik heb net toevallig een ingescande PDF gevonden met daarin het programma dat Bergman heeft gemaakt op bassi van de initiële Slob methode: http://publicaties.minienm.nl/documenten/het-programma-voor-berekening-van-zonneschijnduur-uit-globale-st
Als ik daar door de code kijk, zie ik nu dat daar waardes <0 en >1 worden vervangen door resp. 0 en 1. Dus ik ga dit aanpassen in mijn script.
Er bestaan meerdere algoritmes om de zonneschijnduur te bepalen. Zelf heb ik ook eenvoudiger versie in gebruik omdat mijn oude VP enkel de gemiddelde straling per 10 minuten doorgeeft. Op mijn site wordt het berekend met een harde grens per 10 minuten: is de straling meer dan 40% van de theoretisch mogelijke waarde (bovenaan de atmosfeer) en meer dan 20 W/m², dan schijnt de zon, anders niet. Op maandbasis komt deze zeer goed overeen met de omliggende KNMI stations, en op dagbasis komt het meestal ook goed overeen. Dagen met dikkere (middel)hoge bewolking durven wel eens redelijk af te wijken.
Ik heb ook ooit een vrij eenvoudig algoritme gevonden in een studie welke goed overeen (ook op dagbasis) komt met de moderne meettechniek van directe zonnestraling met drempel op 120 W/m². Die heb ik wel niet meteen bij de hand, want die staat thuis op de PC.
Interessant. Ik ben erg benieuwd naar een eenvoudig algoritme. Met name om met terugwerkende kracht de zonneschijnduur te berekenen aan de hand van alleen de gemiddelde globale straling, datum/tijd en GPS locatie. Want voordat ik de Meteobridge had, had ik ook geen Gmin en Gmax waardes. Maar ik heb wel een database met Ggem die teruggaat tot 2003
De waardes die het KNMI opgeeft (en bij uitbreiding de techniek die jij implementeert), proberen de oude campbell stokes metingen te benaderen. De moderne meettechniek op basis van directe straling geeft lagere waardes.
Zover ik uit het document heb begrepen, heeft Schipper het Slob algoritme aangepast en directe straling als referentie gebruikt uit een grotere dataset:
"Uitgangspunt is nog steeds het algoritme van Slob, echter is bij het algoritme Schipper gebruik gemaakt van een grotere dataset (Jan-95 tot dec-01) en is de directe straling als referentie gekozen. "
Maar dit is natuurlijk ook al 20 jaar oud. Wellicht er weer nieuwe inzichten zijn??