Het script kijkt niet naar een specifieke temperatuur.
Als uitgangspunt neemt het script deze kaarten.
Daar gaat een raster over heen om de kleur van de pixel te vinden waar De Bilt ligt.
De hele kaart met raster ziet er als volgt uit:
De Bilt is te vinden op het pixel dat ligt op 578 / 394
Helemaal ingezoomd ziet dat er zo uit, waarbij de rode X markering de locatie van De Bilt aan geeft.
Op de gebruikte kaarten wordt de T2m getoond in blokken van 3°C
Om die reden zie je bij elk tijdstip altijd een twee waarden staan, zoals bijvoorbeeld op +96 uit Oper GFS0212: +96 | 06-08-25 | 12z | 18 | 21
Vervolgens kijkt het script naar de dag als geheel en pakt daar de laagste categorie uit voor de Tn en de hoogste categorie voor de Tx:
+84 | 06-08-25 | 00z | 12 | 15
+87 | 06-08-25 | 03z | 12 | 15
+90 | 06-08-25 | 06z | 12 | 15
+93 | 06-08-25 | 09z | 18 | 21
+96 | 06-08-25 | 12z | 18 | 21
+99 | 06-08-25 | 15z | 18 | 21
+102 | 06-08-25 | 18z | 15 | 18
+105 | 06-08-25 | 21z | 12 | 15 || 06-08-25 - Tn: 12 | 15 (12,5 °C | 8,5 °C) , Tx: 18 | 21 (20,5 °C)
Vroeger ging het script er vanuit de de Tx rond 14:00 uur zou vallen en voor de Tn was ook een vast tijdstip.
Hoewel dat op veel dagen op gaat, kan het zijn dat een Tx al vroeg in de ochtend valt als bijvoorbeeld de warmte verdreven wordt en dus is dat niet meer en geldt dat ook voor het zoeken naar de Tn.
Vervolgens zie je op de laatste regel van die dag een samenvatting met achtereenvolgens: de datum - laagste Tn categorie (gecorrigeerde Tn | gecorrigeerde Tn indien stralingscondities) , hoogste Tx categorie (gecorrigeerde Tx)
Omdat de categories altijd blokken van 3°C zijn, ligt de Tn en Tx 1,5°C hoger dan de laagste waarde van hun categorie.
Of anders gezegd, de beide waarden van de categorie worden bij elkaar op geteld en gedeeld door 2.
Zodoende komt er voor 6 augustus uit de DMO: 18+21=39 / 2 = 19,5°C
Doordat de standaard afwijking van de hele dataset voor de Tx 1°C is (7,5 jaar aan data), is de verwachting dat de Tx uit komt op 20,5°C
Let wel, dat is voor die ene Oper, wat uiteindelijk door mij in Excel verwerkt gaat worden en dan deel uit maakt van 15 Opers die een verwachting (temperatuur categorie) hebben gedaan voor betreffende dag.
Die 15 waarnemingen worden bij elkaar opgeteld en gedeeld door 15 en daarmee heb je dan de gemiddelde DMO voor die dag.
Voor de 6e staan er nog geen 15 waarnemingen in, maar momenteel staat de gemiddelde DMO op 21,8°C
Daar gaat dan de correctie voor Tx bij en dan kom je op dit moment uit op een meest waarschijnlijke Tx van 22,8°C
Die waarde is terug te vinden bij de gele lijn in onderstaande grafiek (daar staat 22,7 bij maar dat heeft te maken met afrondingsverschillen).
Omdat de ene maand de andere niet is, heb ik in dit hele systeem sinds afgelopen juni ook nog een correctie toegevoegd waarbij de gemiddelde DMO gecorrigeerd wordt met de verzamelde correcties voor betreffende maand.
Die geeft op dit moment voor 6 augustus een meest waarschijnlijke Tx van 22,4°C
