Naar aanleiding van de moeilijk te verwachten neerslagsoort voor morgenochtend, kwam Alwin met de opmerking dat je eigenlijk een algoritme zou moeten hebben om de neerslagsoort te bepalen uit de GFS-modeluitvoer. Toevallig liep ik al wat langer met dat idee rond en daarom heb ik het maar even een beetje uitgewerkt. Het is zomaar een voorstel en ik hoop op reacties. Het lijkt me mooi als we samen iets in elkaar kunnen zetten en dit eventueel later omzetten naar een GrADS script om het te testen!
Omdat ik op dit moment nog onvoldoende kennis heb van de GrADS programmeertaal (Fortran-achtig dacht ik), heb ik het algoritme uitgeschreven in C++/Java stijl en voorzien van wat commentaar.
Het is een eenvoudig model en puur gebaseerd op de accumulatie van de natteboltemperatuur op de niveaus waardoor de ontstane neerslag zal vallen. Het model is erg eenvoudig, maar kan wel al meerdere soorten neerslag onderscheiden, ook wordt rekening gehouden met ijzel, door te kijken naar de mogelijkheden van het opnieuw bevriezen van een sneeuwvlok die (volledig) gesmolten is.
Er moet wel nog veel gefinetuned worden en veel van de constanten die nu voorkomen zijn zeer grof geschat, maar het gaat nu puur om het idee.
Ik benieuwd naar jullie meningen! Ook als iemand een geheel andere methode in gedachten heeft hoor ik het graag.
Groeten,
Bas
ps. Ik heb zelf nog geen tijd gehad om wat literatuur over bepaling van neerslagsoort door te nemen. Zie het
bericht van Saskia voor bruikbare artikelen.
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
Tw is de natteboltemperatuur, deze moet uit de modeluitvoer worden gehaald. Het algoritme doorloopt vanaf het niveau waarop de neerslag ontstaat naar beneden alle niveaus tot de grond.
Quote selectie