Binnen het NWM kan het LSM model gedraaid worden voor verschillende scenario’s. Hiervoor zullen bovenstroomse en laterale instroom, neerslag en verdamping en benedenstroomse waterstanden moeten worden aangepast. De afvoerreeksen voor Lobith en Maas zijn afkomstig vanuit het HBV project (Sperna Weiland, 2015), de afvoerreeksen voor de verschillende scenario’s voor de kleinere modelranden Vecht, Nier, Swalm, Dommel, Bovenmark, Roer en Tongelreep zijn berekend m.b.v. een script (voor details zie https://repos.deltares.nl/repos/Rekenfaciliteit_Deltamodel/trunk/Deltamodel_SA/scripts/ZW_make_timeseries/script/README.txt) dat de percentuele verschillen voor de Maas (tussen elk scenario en de referentiereeks) toepast op de genoemde rivieren (Sperna Weiland, 2016). Vervolgens kan het LSM model zelf gedraaid worden.
Voor de referentiesituatie worden gemeten debieten voor de bovenstroomse invoer (Lobith, Eijsden, Roer, Niers, Swalm, Bovenmark, Dommel, Tongelreep en Overijsselse Vecht) en gemeten waterstanden voor benedenstroomse zeeranden gebruikt (Haringvliet, Maasmond, Bathsluis, Stavenisse, Noordwijk, Lauwersoog, Nieuwe Statenzijl, Den Oever, Kornwerderzand, Harlingen, Vlissingen, Roompot Buiten, Brouwershavensche Gat, Scheveningen, IJmuiden Buitenhaven, Den Helder, en Delfzijl).. Deze zijn als csv bestanden beschikbaar en worden met de randvoorwaarden ingelezen. Laterale debieten voor het LSM zijn afkomstig uit het LHM, welke daar berekend worden vanuit de DM districten. Aan het eind van de NHI workflow worden alle relevante fluxen in Centrale Delft-FEWS database geïmporteerd.
De hieronder staande stappen worden uitgevoerd binnen de LSM workflow. Bij de model run workflow kan gekozen worden tussen standaard LSM en LSM Light (LT). Meer informatie over deze workflows is te vinden in sectie "Workflows en ModuleInstances Zoetwater".
- Import KNMI grid data via OPeNDAP;
- LSM converteren meteo;
- LSM importeren randvoorwaarden;
- LSM (LT) model run.
Wanneer om wat voor reden dan ook een nieuwe LSM Sobek Module aangemaakt wordt dan dient gecheckt te worden of, in het bestand SETTINGS.DAT, de variable PeriodFromEvent op -1 is gezet.
In de Sobek model-schematisatie worden voor de TRIGGERS.def en de CONTROLL.def template bestanden gebruikt; hierin staan periodieke tabellen. In deze template tabellen zijn de jaren met YYYY gemarkeerd. Deze worden per berekening via een het tcl-script (\bin\replaceYYYYinControlTriggersTables.tcl) aangepast door de ‘tag’ YYYY te vervangen met het juiste jaar overeenkomend met het jaar dat doorgerekend wordt. Bij de inbouw van een nieuwe LSM versie dienen daarom template versies van de TRIGGERS.def en de CONTROLL.def (waar in de periodieke tabellen YYYY wordt gebruikt) aangeleverd te worden, zodat FEWS het juiste jaar hier in weg kan schrijven bij aanvang van de modelberekening.
Model invoer voor het LTM Light
De ZW_LSMLTSobek taak roept het LSM Light model aan. De resultaten hiervan, die nodig zijn voor het draaien van de waterkwaliteits module LTM Light, worden aan het eind van de modelrun in de NWM FEWS database geïmporteerd. Het betreft een viertal *.his bestanden die in FEWS als mapstack worden geïmporteerd:
Reachvol.his # volumes van de reaches
Reachflw.his # waterbeweging van de reaches
Nodesvol.his # volumes van de nodes
Qlat.his # laterale debieten
Er wordt een conversie-tool aangeroepen om de *.his bestanden om te zetten naar mapstacks. Dit zijn een reeks *.asc bestanden met een bijbehorende mapstack.xml en *.table file.
De conversie tool kan zowel *.his bestanden omzetten naar ascii-GRID bestanden (TOGRID) als omgekeerd (TOHIS).
Er moeten 5 argumenten meegegeven worden:
- TOGRID of TOHIS
- Locatie en naam van HIS-bestand
- Locatie folder waar uitvoer (*.asc + mapstacks.xml) weggeschreven moet worden
- Naam van grid (deze is in grids.xml gedefinieerd en moet overeenkomen met de dimensie van de *.his file, dat wil zeggen 1 bij x, waarbij x het aantal locaties in de desbetreffende his file is)
- Timezone in gmt (bijv. 1.0; alleen indien TOGRID is gekozen bij eerste argument)
In een convert_mapstacks.inp bestand staat een vertaal tabel van externe parameters (zoals aanwezig in de HIS bestanden) en interne parameters (welke wordt gebruikt voor de naamgeving van de ascii bestanden). Parameters die in meer dan één HIS bestand voorkomen moeten naar unieke parameters gemapped worden zodat ze in FEWS van elkaar onderscheiden kunnen worden.
Figuur 1 De parameters die per *.his bestand worden geïmporteerd.
Allereerst wordt een batchbestand aangeroepen om de convert_mapstacks.exe en convert_mapstacks.inp vanuit de juiste folders in de workdirectory te plaatsen. Vervolgens wordt per his-bestand de executable aangeroepen met de 5 argumenten. In de uitvoer folder worden de ascii-grids en de bijbehorende *mapstack.xml weggeschreven. In deze *mapstack.xml file staan de start en eindtijd van de ascii´s (daarbij komt *0000.asc overeen met de start tijd en *0001.asc de eerstvolgende tijdstap etc.), zodat deze door FEWS met behulp van de importMapstackActivity kunnen worden geimporteerd in de database.
Bij de conversie van *.his naar *.asc (TOGRID) worden naast de *.asc en het *mapstack.xml bestanden ook een *.table files aangemaakt, hierin staan zowel het aantal locaties als de parameternamen uit de his file. Deze file wordt ook gebruikt bij het terugzetten van de *.asc files naar *.his files (TOHIS) en dienen bij die transformatie in de folder te staan van waaruit de transformatie wordt aangeroepen (in veel gevallen de ´work´ dir).
Naast de dynamische *.his bestanden zijn ook nog de bijbehorende *.hia bestanden en een reachdimSTATIC.his bestand als invoer nodig voor LTM Light model.
