SIMONA B&O. Release Notes Export

SIMONA B&O Release Notes Export 2005-02

Rapport Project SIMONA B&O Versie 1.0

Registratienummer P4246-R11

Status Definitief

Datum 25 november 2005

Samengesteld door MX.Systems

In opdracht van RIKZ

Gecontroleerd door M. Elorche (RIKZ)

Datum controle november 2005

Release Notes Export 2005-02

MX . Systems Patrijsweg 32 2289 EX Rijswijk Nederland Tel.: 31(0)70-30 73 737 Fax: 31(0)70-30 73 700 e-mail: [email protected]

Rapport Project SIMONA B&O


Inhoud 1

Inleiding ....................................................................................................3

2

Informatie voor gebruikers.....................................................................4

2.1 2.2

Intro2005-02, Ontwikkelproject Droogv01 fase 3 ...............................................4 Intro2005-03, Onderhoudsintroductie................................................................12 2.2.1 M04036 - Transport-histories in Triwaq..................................................12 2.2.2 M04040 - Domeinverdeling door boxinvoer ...........................................12 2.2.3 M05030 - Installatie-identificatie ............................................................12 2.2.4 M05057 - Waterstanden op enclosure .....................................................12 2.2.5 M05063 - Cross-secties met aflopende coördinaten................................12 2.2.6 P05014 - Foutje in klassewaarden-array..................................................12 2.2.7 P05016 - Keyword Series bij Time_and_Values.....................................12 2.2.8 W04003 - User’s Guide Couple in User’s Guide Waqua ........................13 2.2.9 W05006 - Indeling van R-code onder Roughcombination ......................13 Intro2005-04, GSC: Generalized Spherical Coordinates ...................................14 2.3.1 M05069 - Opname Generalized Spherical Coordinates ..........................14 Intro2005-05, Onderhoudsintroductie................................................................15 2.4.1 M05045 - Backtracking optie in Simpar..................................................15 2.4.2 M05074 - Langere enclosures mogelijk ..................................................15 2.4.3 M05083 - PVM vervangen door MPI; nieuw stopcriterium Waqua .........15 2.4.4 M05091 - Nieuwe condities dynamische barriersturing ..........................15 2.4.5 P05017 - Verdubbelde waarden in Rrsbar ...............................................16 2.4.6 W05001 - Templatebeschrijving en overbodige return Waqview ...........16 2.4.7 W05003 - Kilometreringen bij paskruizen Waqview ..............................16 2.4.8 W05007 - Nieuwe optie minimumwaarden.............................................16

2.3 2.4

3

Informatie voor programmeurs ...........................................................17

3.1

Intro2005-03.......................................................................................................17 3.1.1 P04007 - Getting Started en User’s Guide Dirsys ...................................17 3.1.2 P05013 - Foutje in Sireft..........................................................................17 Intro2005-05.......................................................................................................17 3.2.1 W05007 - Nieuwe optie minimumwaarden.............................................17

3.2

4

Simona onder Windows XP ..................................................................18

4.1 4.2

Het Simona pakket .............................................................................................18 De Ipw................................................................................................................18

5

Sidonia.....................................................................................................19

5.1

Onderhoud..........................................................................................................19 5.1.1 Sid014 - Realisatie van documentatie......................................................19 5.1.2 Sid020 - Variabele root-path.............................................................19 5.1.3 Sid026 - Koppen boven datakolommen min/max waarden .....................19 5.1.4 Sid027 - Stationsnaam bij min/max waarden ..........................................19 5.1.5 Sid030 - Padnamen Gui en Getdata.........................................................19 5.1.6 Sid031 - Inlezen van een Ascii uitvoerreeks in Excel .............................19 5.1.7 Sid032 - Parameter -d gis in commandoregel..........................................20 5.1.8 Sid034 - Runscript getdata.bat........................................................20 5.1.9 Sid037 - Verschillen tussen Arcview en Waqview..................................20

©

MX, datum: 25 november 2005, versie 1.0

pagina 1 van 26



5.2

Nieuwe ontwikkelingen .....................................................................................20 5.2.1 Ontsluiten informatie m.b.t. dynamische barriers....................................20 5.2.2 Het kunnen ophalen van de initiële tijdstap voor veldwaarden ...............21 5.2.3 Ontsluiten van Slib3D SDS-files .............................................................21 5.2.4 Ontsluiten van wind SDS-files ................................................................22 5.2.5 Implementatie van de verbeterde afvoerpotentiaal ..................................22 5.2.6 Het kunnen ophalen van verschillen tussen twee SDS-files ....................22

6

Configuratie............................................................................................24

6.1 6.2 6.3 6.4

Platforms ............................................................................................................24 Configuratie-eisen..............................................................................................24 Producten ...........................................................................................................24 Bestanden...........................................................................................................25

7

Versienummers ......................................................................................26

©


pagina 2 van 26


1


Inleiding

Deze Release Notes horen bij Simona Export2005-02 (november 2005). Dit document bevat de belangrijkste wijzigingen die sinds de vorige export (Export2005-01) in Simona zijn doorgevoerd. Deze wijzigingen betreffen de introducties 2005-02 tot en met 2005-05, inclusief opnames van nieuwe functionaliteiten: •

Intro2005-02

• • •

Intro2005-03 Intro2005-04 Intro2005-05

Ontwikkelproject Droogv01, fase 3: “Verbetering droogvallen en onderlopen in Waqua/Triwaq” Onderhoudsintroductie (B&O wijzigingen) Ontwikkelproject GSC: Generalized Spherical Coordinates Onderhoudsintroductie (B&O wijzigingen)

Verder wordt melding gemaakt van de onderhoudswijzigingen en nieuwe ontwikkelingen in Sidonia, en de stand van zaken met betrekking tot de Windows XP versie van Simona en de Ipw. Alle documenten waarnaar in dit document wordt verwezen, kunnen op aanvraag beschikbaar gesteld worden. Voor ondersteuning van Simona kan men t/m 31 december 2005 terecht bij de Simona B&O Groep bij MX.Systems: 070 - 3073 733 • Servicedesk, [email protected] Voor ondersteuning van Sidonia zijn de contactpersonen t/m 31 december 2005: • Mustapha Elorche RIKZ, 070-3114249 • Joost Beckers RIKZ, 070-3114251 Advies voor nieuwe ontwikkelingen: Hierbij willen wij erop wijzen dat uitbreidingen van bepaalde onderdelen van Simona zoals lokale arrays (intgda, reagda, intloc, enz.), lokale data structuren (LDS-en), foutmeldingenbestanden, enz., vooraf door ontwikkelaars moeten worden gereserveerd bij de Simona B&O Groep, om “overboekingen” te voorkomen.

©


pagina 3 van 26


2


Informatie voor gebruikers

In dit hoofdstuk worden de voor de gebruiker belangrijkste wijzigingen opgesomd.

2.1 Intro2005-02, Ontwikkelproject Droogv01 fase 3 In het droogvalproject van 2004 heeft VORtech Computing met assistentie van het WL | Delft Hydraulics voor Rijkswaterstaat / RIKZ een nieuwe versie van WAQUA / TRIWAQ gerealiseerd waarin verbeteringen zijn aangebracht met betrekking tot droogvallen en onderlopen. Bovendien zijn WAQUA en TRIWAQ in dit project voor een belangrijk deel geüniformeerd en afgestemd op het OMS-project. Het project is in drie fases verdeeld: • In fase 1 is een detailontwerp voor de beoogde wijzigingen geschreven (rapport TR04-02). • In fase 2 zijn informatica-technische aanpassingen voor WAQUA / TRIWAQ geïmplementeerd (opgenomen in Export2005-01). • In fase 3 zijn aanpassingen aan de numerieke algoritmes van WAQUA / TRIWAQ geïmplementeerd (opgenomen in onderhavige Export2005-02). Nadere informatie is te vinden in de volgende documenten: Technisch Rapport TR04-02, versie 1.1, augustus 2004 "Detailontwerp verbetering droogvallen en onderlopen in WAQUA/TRIWAQ" dr. ir. E.A.H. Vollebregt, dr. ir. B. van ’t Hof (VORtech Computing, Delft) ir. J.A.Th.M. van Kester (WL | Delft Hydraulics) Technisch Rapport TR04-04, versie 0.5, september 2004 "Testen van verbeterde algoritmen voor droogvallen en onderlopen in WAQUA/TRIWAQ" dr. ir. E.A.H. Vollebregt, dr. ir. B. van 't Hof (VORtech Computing, Delft) ir. J.A.Th.M. van Kester (WL | Delft Hydraulics) Acceptatietesten ‘Verbetering droogvallen en onderlopen in WAQUA/TRIWAQ’, RIKZ/KW/2005.125W, september 2005. “User’s Guide WAQUA, General Information”, § 3.2.1.2 en § 3.6.2. “User’s Guide WAQPRE” Nieuwe keywords onder MESH / BATHYMetry / GLOBAL: § 2.6.5 DPD_GIVEN, DPS_GIVEN, METH_DPS en METH_DPUv Nieuwe keywords onder FLOW / PROBlem / DRYING: § 2.8.1.4 UPWIND_zeta, THRES_UV_flooding, THRES_WL_flooding en CHECK_WL Memo EV/M04.100, versie 1.0, 10 december 2004 "Overzicht van nieuwe droogvalopties voor gebruikers van WAQUA / TRIWAQ" Erik de Goede (WL | Delft Hydraulics) Edwin Vollebregt en Bas van ’t Hof (VORtech Computing) De hiernavolgende tekst is voornamelijk ontleend aan het laatstgenoemde memo. Samenvatting

©


pagina 4 van 26



In het memo EV/M04.100 en het acceptatietestrapport RIKZ/KW/2005.125W wordt een overzicht gegeven van de nieuwe opties van WAQUA / TRIWAQ met betrekking tot de diepteschematisatie en met betrekking tot droogvallen en onderlopen. Deze nieuwe opties zijn gerealiseerd in het droogvalproject 2004 dat door VORtech Computing in samenwerking met WL | Delft Hydraulics is uitgevoerd. Het memo heeft tot doel de WAQUA / TRIWAQ gebruiker meer inzicht te geven in de droogval- en onderloopopties die er nu zijn. Er wordt aangegeven welke combinaties van opties voor zover nu bekend aan te raden zijn. Tenslotte wordt aangegeven hoe de nieuwe opties in bestaande modellen kunnen worden ingevoerd. Het acceptatietestrapport heeft tot doel om duidelijkheid te verschaffen over de mate waarin de nieuwe geïmplementeerde functionaliteiten omtrent diepteschematisatie en droogvalprocedure voldoen aan de wensen c.q. eisen van de opdrachtgever. Verder had Rikz ook behoefte om ervaring op te doen met de nieuwe functionaliteiten teneinde praktische adviezen aan gebruikers te kunnen geven. Mogelijkheden voor de diepteschematisatie in WAQUA / TRIWAQ Inleiding In het project "Verbeteren droogvallen in onderlopen in WAQUA / TRIWAQ'' (hierna te noemen project "DROOGV01'') zijn de mogelijkheden voor de diepteschematisatie flink uitgebreid. In het memo worden de verschillende (oude en nieuwe) mogelijkheden op een rij gezet. Verder wordt aangegeven welke combinaties de voorkeur verdienen en hoe bestaande modellen geconverteerd kunnen worden naar de nieuwe opties. Diepteschematisatie in WAQUA / TRIWAQ Het specificeren van een diepteschematisatie in WAQUA / TRIWAQ is gerelateerd aan het rekenrooster dat gebruikt wordt. In WAQUA / TRIWAQ wordt een zogeheten "gestaggerd rooster'' toegepast. Dit betekent dat waterstands-, snelheids- en dieptepunten op verschillende, ten opzichte van elkaar verschoven locaties, zijn gedefinieerd. Het gestaggerde rooster van WAQUA wordt geïllustreerd in Figuur 1.

©


pagina 5 van 26



Figuur 1: Overzicht van het gestaggerde rooster in WAQUA / TRIWAQ Bij het specificeren van een diepteschematisatie in WAQUA/TRIWAQ spelen de volgende onderdelen een rol: 1. Waar worden de dieptes gespecificeerd? 2. Welk algoritme wordt gebruikt om de diepte in waterstandspunten te berekenen indien deze niet gegeven zijn (optie DPD_GIVEN)? 3. Welk algoritme wordt gebruikt om de diepte in snelheidspunten te berekenen?

Figuur 2: Definitieschets voor opgegeven dieptes D1, D2, D3 en D4 in dieptepunten voor de berekening van de diepte in een waterstandspunt Oude opties voor de diepteschematisatie De beschikbare opties in WAQUA / TRIWAQ waren voordat het project DROOGV01 werd uitgevoerd als volgt:

©


pagina 6 van 26



1. Dieptes werden opgegeven in dieptepunten. 2. Voor de diepte in waterstandspunten (keywoord IDRYFLAG) waren er vier opties, namelijk 0, 1, 2 en 3. Deze komen min of meer overeen met de nieuwe opties METH_DPS = 'MAX_DPUV', 'MEAN_DPD' en 'MIN_DPUV'. Deze nieuwe opties worden in het memo uitgebreid toegelicht. 3. Voor de diepte in snelheidspunten werd altijd het gemiddelde genomen van de twee naburige dieptepunten. Dit komt overeen met de nieuwe optie METH_DPU = 'MEAN_DPD'. Als we aannemen dat de vier bodemhoogtes in de hoekpunten van een roostercel D1, D2, D3 en D4 zijn conform de definities in Figuur 2, en dat de bodemhoogtes positief zijn in benedenwaartse richting t.o.v. de referentiehoogte, dan zijn de bovengenoemde drie opties voor diepte in waterstandspunten te definiëren als: •

IDRYFLAG = 0 (nieuwe optie METH_DPS = 'MAX_DPUV'):

 D1 + D 2 D 2 + D 4 D 4 + D 3 D 3 + D1  DPS = D = D ς = max  , , ,  2 2 2 2   Dit is het zogenaamde "maximum criterium'' waarbij de diepte in het midden van een cel gelijk is aan de grootste diepte van de vier zijdes van de cel. •

IDRYFLAG = 1 (nieuwe optie METH_DPS = 'MEAN_DPD'):

DPS = D ς =

D1 + D 2 + D 3 + D 4 4

Dit is het zogenaamde "gemiddelde criterium''. •

IDRYFLAG = 2 (nieuwe optie METH_DPS = 'MIN_DPUV'):

 D1 + D 2 D 2 + D 4 D 4 + D 3 D 3 + D1  DPS = D ς = min  , , ,  2 2 2 2   Dit is het zogenaamde "minimum criterium''. Hierbij is de diepte in een cel gelijk aan de kleinste diepte van de vier zijdes van de cel. •

In WAQUA / TRIWAQ vóór het project DROOGV01 waren er al drie opties mogelijk voor de berekening van de diepte in de waterstandspunten. In project DROOGV01 is een vierde optie toegevoegd, namelijk: METH_DPS = 'MAX_DPD'

DPS = Dς = max ( D1, D 2, D3, D 4) Dit is een "maximum criterium" waarbij naar de diepte in de hoekpunten in plaats van de zijden van een cel wordt gekeken. Met deze optie wordt een bodem gegenereerd die dieper is dan die van het maximale criterium dat al bestond ('MAX_DPUV'). Deze optie is toegevoegd omdat dit het "maximum criterium" is dat in Delft3D-FLOW wordt toegepast.

©


pagina 7 van 26



Als de eerste drie opties van boven naar beneden worden gevolgd, dan wordt de bodem in de middens van een roostercel steeds ondieper. De gebruikte formule voor de berekening voor de diepte in snelheidspunten was voorafgaand aan het project DROOGV01 niet instelbaar voor gebruikers. De gebruikte formule is:

DPU = D U =

D2 + D4 2

Deze formule correspondeert met de nieuwe invoeroptie METH_DPUV = 'MEAN_DPD'. Voorafgaand aan het project Droogval2004 kon de gebruiker dus eigenlijk kiezen uit drie methodes voor de bodemschematisatie, IDRYFLAG = 0, 1 en 2. De waarde IDRYFLAG = 3 levert dezelfde bodemhoogtes op als IDRYFLAG = 1. Gebruik van deze optie zorgt ervoor dat er in de simulatie minder droogvalcontroles worden uitgevoerd. Uitbreiding van opties voor de diepte in snelheidspunten In het project DROOGV01 zijn extra opties voor de diepte in snelheidspunten toegevoegd. Met deze opties worden discrepanties tussen de dieptes in waterstandspunten en snelheidspunten voorkomen, waardoor het gedrag van onderlopen en droogvallen gelijkmatiger en derhalve beter wordt. De volgende drie opties zijn toegevoegd voor de berekening van de diepte in een snelheidspunt: •

METH_DPUV = 'MIN_DPS' met

(

)

(

DPU = D U = min D mς , n , D mς +1, n , DPV = D V = min D mς , n , Dmς , n +1

•

METH_DPUV = 'MEAN_DPS' met

DPU = D = U

•

)

Dmς ,n + Dmς +1,n 2

, DPV = D = V

Dmς , n + Dmς ,n +1 2

METH_DPUV = 'MAX_DPS' met

(

)

(

DPU = D U = max Dmς , n , Dmς +1, n , DPV = D V = max Dmς ,n , Dmς ,n +1

)

Bij deze drie opties is de diepte in een snelheidspunt gelijk aan respectievelijk de minst diepe, het gemiddelde en de meest diepe van de waardes in de tweede naburige waterstandspunten.

©


pagina 8 van 26



Uitbreiding van locaties voor invoerdieptes In het project DROOGV01 is de optie ingebouwd dat bij invoer (dus via de "siminp"file) de dieptewaarden ook in waterstandspunten opgegeven kunnen worden. Dus naast de optie van het specificeren van dieptes in dieptepunten (DPD_GIVEN) is er nu een tweede optie mogelijk met dieptes in waterstandspunten (DPS_GIVEN). Een belangrijk voordeel van de laatste methode is dat de komberging van een model veel directer kan worden aangestuurd c.q. beïnvloed. Voor de komberging is namelijk van belang bij welke waterstand een roostercel uit de berekening genomen wordt, en dit wordt bepaald door de diepte in het midden van de cel. Deze kan met de nieuwe optie direct gespecificeerd worden in de siminp-file, in plaats van dat ze door interpolatie wordt bepaald. Een ander voordeel van de nieuwe methode is smalle geultjes met een celbreedte weergegeven kunnen worden, zie Figuur 4. In Figuur 3 is aangegeven hoe in de oude situatie een geul gemodelleerd werd.

Figuur 3: Natte doorsnede op basis van een "gemiddelde aanpak" (diepte in snelheidspunt is het gemiddelde van de dieptes in naastgelegen dieptepunten). Deze oude optie heet DPD_GIVEN.

Figuur 4: Natte doorsnede op basis van een "tegelbenadering" (diepte in snelheidspunt is het minimum van de dieptes in naastgelegen waterstandspunten). Deze nieuwe optie heet DPS_GIVEN

©


pagina 9 van 26



In de siminp-file wordt een diepteschematisatie in dieptepunten aangegeven met het keywoord DPD_GIVEN en voor de nieuwe optie van dieptes in waterstandspunten wordt het keywoord DPS_GIVEN gebruikt. Indien de dieptewaarden in waterstandspunten gedefinieerd worden, dan is de berekening van dieptewaarden in waterstandspunten (zie keywoord METH_DPS) uiteraard niet meer van toepassing. Dieptewaarden in dieptepunten zijn ook niet meer nodig. Ze worden wel berekend in WAQUA voor bijvoorbeeld de postprocessing. De dieptewaarden in dieptepunten hebben echter geen invloed op de modelresultaten. In WAQUA / TRIWAQ is het voldoende als de dieptes in het waterstandspunt en in de en -snelheidspunten bekend zijn. Uitbreiding van de opties voor het droogvalalgoritme In het project DROOGV01 zijn vele aanpassingen aan de algoritmes voor het droogvallen en onderlopen in WAQUA / TRIWAQ geïmplementeerd. Alleen wanneer dit in de gebruikersinvoer is terug te zien zal dat worden besproken in dit rapport. Voor de overige werkzaamheden in het project DROOGV01 wordt verwezen naar Technisch Rapport TR04-04 (zie bovengenoemde verwijzing). Er zijn drie bestaande opties (namelijk DUPWND, DEPCRIT en IDRYFLAG) omgezet naar nieuwe keywoorden. De functionaliteit is iets generieker geworden. Zo is het keywoord DEPCRIT min of meer opgesplitst in twee keywoorden (namelijk in THRES_WL_FLOODING en THRES_UV_FLOODING). Hierdoor kan nu zowel in waterstands- als in snelheidspunten een droogvaldrempel opgegeven worden. Tabellen 1 en 2 bevatten een overzicht van de oude en de nieuwe keywoorden. Voor een meer gedetailleerde beschrijving van deze oude en nieuwe keywoorden zij verwezen naar paragraaf 2.1 van bovengenoemd memo. Tabel 1: Overzicht van oude en nieuwe keywoorden voor droogvallen

oude keywoorden DUPWND DEPCRIT IDRYFLAG Tabel 2

nieuwe keywoorden UPWIND_ZETA THRES_WL_FLOODING, THRES_UV_FLOODING CHECK_WL

Oude opties in relatie tot de nieuwe keywords

equivalente nieuwe optie

Oude optie

waterstand in snelheidspunten (METH_DPUV)

waterstand in waterstandspunten (METH_DPS)

DPD_GIVEN DPD_GIVEN DPD_GIVEN DPD_GIVEN

MEAN_DPD MEAN_DPD MEAN_DPD MEAN_DPD

MAX_DPUV MEAN_DPD MIN_DPUV MEAN_DPD

DPD_GIVEN

MEAN_DPD

MAX_DPD

waarde idryflag bodempunten 0 1 2 3 (Max-optie Delft3D)

©


pagina 10 van 26



Opgemerkt wordt dat bij nieuw op te zetten modelschematisaties in het algemeen aanbevolen wordt een diepteschematisatie in waterstandspunten (DPS_GIVEN) te gebruiken. Een dergelijke aanpak biedt voordelen boven de oorspronkelijke aanpak. De komberging kan directer worden aangestuurd en er treden minder artefacten door droogvallen en onderlopen op. Verder is een koppeling met morfologie eenvoudiger te realiseren. Ten slotte kunnen impulsbehoudende schema's die worden overwogen voor WAQUA gemakkelijker worden gerealiseerd wanneer een diepteschematisatie in waterstandspunten wordt gebruikt. Aan de gebruiker wordt het volgende advies gegeven voor gebruik van de nieuwe opties: 1. De belangrijkste verbetering van de nieuwe opties is te verwachten door de dieptes in de snelheidspunten (METH_DPUV) te laten uitrekenen op basis van de diepte in de waterstandspunten (dus MIN_DPS, MEAN_DPS of MAX_DPS maar niet MEAN_DPD). 2. De juiste keuze uit de 3 bovengenoemde goede opties (METH_DPUV= MIN_DPS, MEAN_DPS of MAX_DPS) hangt van de praktijksituatie af. Veelal zal de optie MIN_DPS goede resultaten geven maar niet altijd (zgn. ‘tegelaanpak’) . Vaak is ook de optie MEAN_DPS goed en soms ook de optie MAX_DPS. 3. Om verbeteringen te zien moet ook DEPCRIT worden verlaagd in geval deze waarden heeft van orde 5cm of groter. In geval van de nieuwe opties leidt deze verlaging niet tot instabiel modelgedrag wat bij een aantal (met name estuariummodellen) vroeger wel het geval was. 4. De modellen moeten veelal opnieuw worden afgeregeld bij gebruik van de nieuwe opties, in het bijzonder wordt de komberging van het model veranderd bij keuze van nieuwe opties.

©


pagina 11 van 26



2.2 Intro2005-03, Onderhoudsintroductie 2.2.1

M04036 - Transport-histories in Triwaq

In Triwaq werden bij de berekening van transport-histories stroomsnelheidswaarden gebruikt van een halve rekenstap oud. Dit in tegenstelling tot Waqua, waarin de meest recente waarden gebruikt worden. Deze onvolkomenheid in Triwaq is hersteld. De gebruiker moet rekening houden met (soms flinke) verschillen in het Advectief Constituent Transport (ATR) op constituent stations.

2.2.2

M04040 - Domeinverdeling door boxinvoer

Naar aanleiding van een foutsituatie bij de domeinverdeling voor een parallelle rekensom door middel van de box-invoermethode is geconstateerd dat de Algemene Preprocessor, die de siminp-invoerfiles verwerkt, een fout bevat. Om te voorkomen dat deze fout vaker optreedt, is in de User’s Guides van alle betrokken applicaties de volgende waarschuwing met een uitleg van het probleem opgenomen: “If the last keyword block in the input file contains a sequential keyword, the Simona application independent preprocessor is not able to check the correctness of the block. This can result in incorrect processing of the input file!”

2.2.3

M05030 - Installatie-identificatie

Bij foutsituaties die aangemeld worden bij de beheer en onderhoud groep is het handig om te weten met welke installatie van de betrokken applicatie is gewerkt. Het versienummer alleen biedt daartoe niet voldoende informatie. Daarom wordt nu ook het pad “$SIMONADIR”, dat aangeeft in welke directory de gebruikte executable staat, naar de messagefile geschreven. Dit pad maakt duidelijk met welke programmainstallatie gewerkt is. NB Voorlopig is dit alleen geïmplementeerd in Waqpre en Waqpro.

2.2.4

M05057 - Waterstanden op enclosure

Door Waqview werden ten onrechte waterstanden afgebeeld in cellen aan de zuidrand van het model die toebehoren aan de grid enclosure. Deze fout is hersteld, zodat nu alleen nog waterstanden en andere grootheden in het feitelijke rekengebied getoond worden.

2.2.5

M05063 - Cross-secties met aflopende coördinaten

Waqview beeldde ten onrechte alleen cross-secties af als die gedefinieerd waren met oplopende m- en/of n-coördinaten. Deze fout is nu hersteld, zodat nu alle cross-secties naar behoren getoond worden.

2.2.6

P05014 - Foutje in klassewaarden-array

Er bleek een foutje te zitten in de verwerking van de klassewaarden (keyword Class_limits in de Waqpre invoer). De waarde van de positie 1 (variabele Nclas) in het array Class_changes_sep werd door Waqpro verkeerd berekend. Deze fout is nu hersteld.

2.2.7

P05016 - Keyword Series bij Time_and_Values

Tijdens het onderzoek van een probleem met een model met uniforme wind is geconstateerd dat de Algemene Preprecessor, die de Siminp-invoerfiles verwerkt, een bepaalde situatie rondom Time_and_Values niet goed afhandelt. Als namelijk vóór het keyword Time_and_Values geen keyword Series is opgegeven, worden de waarden

©


pagina 12 van 26



onder Time_and_values niet herkend, en loopt de invoerverwerking dus verkeerd. Deze fout is nu hersteld: er wordt nu altijd gecontroleerd of de combinaties van Series plus Frame en Values, en Series plus Time_and_Values, correct en volledig zijn.

2.2.8

W04003 - User’s Guide Couple in User’s Guide Waqua

De “User’s Guide for Parallel Waqua/Triwaq and for Domain Decomposition”, in de wandeling ook wel “User’s Guide Couple” genoemd, was een apart, zelfstandig document. Om de verbondenheid met het Waqua pakket te benadrukken, is dit document nu opgenomen in de bundel “User’s Guide Waqua”.

2.2.9

W05006 - Indeling van R-code onder Roughcombination

In de Waqpre invoerstruktuur onder het keyword Friction/Roughcombination (zie paragraaf 2.8.1.5 van de User’s Guide Waqpre) staat het keyword R-code. Dit keyword kent een bepaald bereik, met een indeling naar betekenis. In deze indeling is een kleine wijziging aangebracht: de range 1201-1300 is veranderd naar 1201-1400. De range 1401-1500 blijft ‘not defined’.

©


pagina 13 van 26



2.3 Intro2005-04, GSC: Generalized Spherical Coordinates 2.3.1

M05069 - Opname Generalized Spherical Coordinates

In Waqua is naast de bestaande roostertypen (Rectilinear, Curvilinear en Spherical) een nieuw roostertype geïmplementeerd: Generalized Spherical Coordinates (GSC); zie paragraaf 2.6.1 van de User’s Guide Waqpre. Met GSC kunnen de voordelen van een kromlijnig en een sferisch rooster gecombineerd worden: het is een kromlijnig rooster op de bolvorm van de aarde. Dit nieuwe roostertype kan gebruikt worden in combinatie met parallel rekenen, space varying wind and pressure, en met horizontale en verticale domeindecompositie. NB De implementatie van deze nieuwe functionaliteit heeft een gering effect op bestaande modellen: er kunnen kleine verschillen in resultaten optreden. NB Wanneer dit type grid wordt opgegeven middels een Rgf-bestand moet de gebruiker bedacht zijn op afrondfouten. Bij fijnmazige modellen, waarbij de celgrootte kleiner is dan enkele kilometers is het raadzaam om de coördinaten in dubbele precisie te zetten. Zie ook de ‘Note’ in paragraaf 2.6.1.5 van de User’s Guide Waqpre. Voor achtergrondinformatie wordt verwezen naar de volgende rapporten van Vortech: • Memo CvV 04.049, Generalized Spherical Coordinates in Waqua/Triwaq. • Memo CvV/M04.077, Verslag volledig doorvoeren GSC in Waqpre/Waqpro.

©


pagina 14 van 26



2.4 Intro2005-05, Onderhoudsintroductie 2.4.1

M05045 - Backtracking optie in Simpar

Simpar is uitgebreid met de ‘backtracking’ optie; zie keyword Ibck in paragraaf 3.4.1.1 van de User’s Guide Simpar. Met deze optie is het mogelijk om particle banen te berekenen in omgekeerde tijdsvolgorde. De gebruiker geeft begin- en eindtijd van de simulatie en het programma rekent vervolgens terug van eindtijd naar begintijd. Dus met het eindtijdstip als startpunt worden de particle verplaatsingen berekend met “omgekeerde” snelheidsvectoren, terug tot aan het begin van de simulatie. NB Backtracking is alleen mogelijk onder de volgende voorwaarden: (1) alleen advectie, geen diffusie; (2) alleen floating particles; (3) geen cyclische snelheidsvelden; (4) geen restart-Sds-bestanden.

2.4.2

M05074 - Langere enclosures mogelijk

Met het oog op de grotere parallelle of domein-decompositie modellen is het maximum aantal elementen voor enclosures verhoogd van 20.000 naar 100.000.

2.4.3

M05083 - PVM vervangen door MPI; nieuw stopcriterium Waqua

Voor de verbetering van de performance van parallel Waqua/Triwaq is het pakket voor de ondersteuning van parallel rekenen, PVM (Parallel Virtual Machine), vervangen door het modernere pakket MPI (Message Passing Interface). Om precies te zijn is het pakket MPICH2 versie 1.0.2p1 geïmplementeerd, voor de Linux en Hp versies van Waqua/Triwaq. N.B. Helaas is door deze wijziging parallel draaien op HP niet meer mogelijk ! Daarnaast is een performanceverbetering aangebracht voor parallel Waqua. Dit betreft het zogeheten “nieuwe stopcriterium”. In grote lijnen komt het erop neer dat per deelrij (of deelkolom) niet meer iteraties worden uitgevoerd dan voor die deelrij noodzakelijk is. NB Deze verbetering heeft effect op de rekenresultaten van Waqua; voor sequentiële berekeningen en voor parallelle berekeningen zonder droogval zijn deze effecten in het algemeen klein; voor parallelle runs met droogval kunnen de verschillen met oude runs wat groter zijn. Voor achtergrondinformatie wordt verwezen naar de volgende rapporten van Vortech: • Memo EV/M05.040, Ontwikkeling van een Mpi-versie van parallel Waqua/Triwaq voor Myrinet. • Technisch rapport TR04-03, Verbetering van de performance van parallel Waqua op het Linux rekencluster van Riza-Arnhem.

2.4.4

M05091 - Nieuwe condities dynamische barriersturing

Ten behoeve van het simuleren van de balgstuw bij Ramspol is een uitbreiding van de dynamische aansturing van barriers gecreëerd. Hiertoe is een extra optie toegevoegd aan de rij mogelijkheden voor het definiëren van een Simple Condition, zie paragraaf 2.8.2.9 Barriers, onder Condition. Aan de rij “computed entities” zijn toegevoegd: Sill_depth, Gate_height, en Barrier_width. Wanneer nu een barrier wordt aangestuurd met een (relatieve) tijdreeks, kan met behulp van de nieuwe condities ervoor worden gezorgd dat de tijdreeks geheel wordt doorlopen zonder “onderbroken” te worden door een andere conditie. Met name voor de balgstuw is dit noodzakelijk.

©


pagina 15 van 26


2.4.5


P05017 - Verdubbelde waarden in Rrsbar

Bij parallelle sommen met modellen met lijnbarriers ging er iets fout bij restart. De waarden van het array Rrsbar, getoond in de Waqpro report print file (onderaan), bleken na restart verdubbeld te zijn. Dit euvel is nu opgelost. Voor achtergrondinformatie wordt verwezen naar het volgende rapport van Vortech: • Memo CvV/M05.015, Foute waardes Rrsbar in restart som.

2.4.6

W05001 - Templatebeschrijving en overbodige return Waqview

Voor Waqview zijn twee problemen opgelost: 1) Aanmaken van een template De beschrijving in paragraaf 6.3 van de User’s Guide Waqview van het aanmaken van een template (bestandstype 5) was niet volledig. Met name was niet duidelijk (a) dat het selectievakje ‘Write a template’ aangeklikt moest worden, (b) dat in het tekstveld ‘Set Version number for template’ indien gewenst de waarde 200 ingevuld moet worden, en (c) dat tenslotte de knop ‘Create file’ moet worden ingedrukt. Dit alles staat wél in paragraaf 7.2 (Waqview Faq’s); het probleem is nu eenvoudig opgelost door in paragraaf 6.3 te verwijzen naar paragraaf 7.2. 2) Return bij invulvelden Voor gebruikers is het lastig en niet intuïtief dat na het invullen van een tekstveld als afsluiting een <enter> gegeven moet worden. Dit probleem is programmatisch opgelost: een <enter> is nu niet meer nodig.

2.4.7

W05003 - Kilometreringen bij paskruizen Waqview

In sommige gevallen beeldde Waqview de kilometeraanduidingen bij paskruizen (Knop Pk, paragraaf 6.2.1 User’s Guide Waqview) aan de rechterzijde (y-as) en de bovenzijde (x-as) niet af. Dit probleem is nu opgelost: in alle gevallen worden de kilometreringen afgebeeld.

2.4.8

W05007 - Nieuwe optie minimumwaarden

Naast de berekening van maximumwaarden is het nu ook mogelijk om minimumwaarden te laten berekenen. Hiertoe is het nieuwe keyword Calcminvalues in het leven geroepen. Hierbij is de gelegenheid aangegrepen om de keywords betreffende de maximum- en minimumberekeningen en bepaling van klassewaarden op een logischer plaats te positioneren, namelijk onder het keyword Sdsoutput. Zie paragraaf 2.13.6, 2.13.7 en 2.13.8 van de User’s Guide Waqpre. Het oude keyword Postproc is verwijderd. NB Gebruikers met modelinvoer waarin de keywords nog onder Postproc staan, moeten hun invoer aanpassen. Mochten ze toch per ongeluk een Waqpre run starten met de oude invoer, dan verschijnt er een melding dat het keyword niet herkend wordt. NB De onderhavige berekeningen kunnen nu niet alleen met Waqua, maar ook bij Triwaq met één laag worden uitgevoerd. NB Om een en ander te realiseren, is de structuur van de uitvoer-arrays voor maximum- en minimumwaarden in de Sds-file gewijzigd. De resultaten komen nu in arrays in het compound array Solution_derived. Zie verder de beschrijving van dit array in het document Lds_flow. Dit houdt wel in dat postprocessingprogrammatuur die gebruikmaakt van deze uitkomsten, aangepast zal moeten worden.

©


pagina 16 van 26


3


Informatie voor programmeurs

In dit hoofdstuk worden de voor de programmeur belangrijkste wijzigingen opgesomd.

3.1 Intro2005-03 3.1.1

P04007 - Getting Started en User’s Guide Dirsys

Dirsys kende de twee aparte documenten “Getting Started Dirsys” en de “Programmer’s Guide Dirsys”. Het eerste document is nu vertaald naar het Engels en daarna zijn de twee documenten samengevoegd tot één document. De naam van het nieuwe document is “Dirsys Getting Started and Programmer’s Guide”.

3.1.2

P05013 - Foutje in Sireft

Het hulpprogramma Sireft liep fout doordat een bepaalde variabele (‘sekese’) niet werd geïnitialiseerd. Deze fout is nu hersteld.

3.2 Intro2005-05 3.2.1

W05007 - Nieuwe optie minimumwaarden

Dit nummer is beschreven in paragraaf 2.4.8, voor gebruikers. In die paragraaf staat ook informatie die van belang is voor programmeurs van postprocessingprogrammatuur.

©


pagina 17 van 26


4


Simona onder Windows XP

4.1 Het Simona pakket Informeel bestaat er reeds geruime tijd een Windows XP-versie van een geselecteerd aantal deelsystemen van Simona, t.w. waqpre, waqpro, waqbhd, waqpan en waqview. Vooralsnog betreft dit een sequentiële versie, waarbij de aansturing van de applicaties verloopt via bat-files. Voor het rekenen met grote(re) modellen worden executables met verschillende (discrete) buffergroottes meegeleverd. Eerder dit jaar is een eerste, op MPI gebaseerde implementatie gerealiseerd van een parallelle versie van het Simona-pakket. De aanroep van de verschillende deelsystemen gebeurt hierbij met behulp van Perl-scripts. Momenteel bevindt deze versie zich nog in het stadium van de acceptatie tests. Om die reden zal in de onderhavige export nog ‘slechts’ een sequentiële versie worden uitgeleverd. In deze release zijn ook de deelsystemen SIMPAR en SLIB3D opgenomen.

4.2 De Ipw Zoals bekend wordt de IPW sinds de vorige export alleen nog maar uitgeleverd voor het Windows platform (Windows2000/XP). Voor het draaien van de IPW moet de gebruiker beschikken over een X-emulator (Exceed v.7.1 of hoger resp. Reflection X v.12.04 of hoger). Uit de dagelijkse praktijk is inmiddels bekend dat de applicatie ook werkt met ‘oudere’ versies van Reflection X (vanaf v.10.x); essentieel is wel dat binnen Reflection de kleurinstelling ingesteld staat op ‘Best Visual’ of ‘True Color’. Een en ander staat ook beschreven in de installatie-handleiding van de IPW. Wat is nieuw in IPW v. 3.08.00 Naast kleine onderhoudsverbeteringen bevat de nieuwste versie van de IPW twee belangrijke aanpassingen: 1. Ondersteuning Baseline opslagstructuur. Bij het bewaren van een met de IPW ingelezen model wordt de modelinvoer nu weggeschreven volgens de directory-structuur van Baseline. 2. Aansluiting bij nieuwe droogvalopties. Onder het menu-item Mesh → Bathymetry.... biedt de IPW nu de mogelijkheid om de verschillende opties samenhangend met de nieuwe droogvalprocedure op te geven resp. te wijzigen. Dit betreft de key-words DPD_GIVEN resp. DPS_GIVEN, Method_DPS en Method_DPUV.

©


pagina 18 van 26


5


Sidonia

In dit hoofdstuk worden de verbeteringen in de Sidonia beschreven. Dit betreft de aanpassingen in de huidige exportversie (V0.96) ten opzichte van de vorige exportversie (V0.93).

5.1 Onderhoud 5.1.1

Sid014 - Realisatie van documentatie

In het kader van dit nummer is, gebruikmakende van reeds beschikbare documenten, documentatie voor Sidonia gerealiseerd. Dit betreft de volgende drie documenten: 1. Gebruikershandleiding Sidonia 2. Programmeerhandleiding Sidonia 3. Installatiegids Sidonia

5.1.2

Sid020 - Variabele root-path

Middels het zetten van de variabele root_path in het bestand properties.txt, wat staat in de directory sidonia-client-xp/conf c.q. sidonia-clientlinux/conf, kan een standaard leespad, bijv. voor de Gui, worden ingesteld. Op XP gaat het instellen van het pad bijvoorbeeld als volgt: root_path=c:\\sidonia_v094\\test\\basis_eems

(let op de dubbele scheidingstekens)

5.1.3

Sid026 - Koppen boven datakolommen min/max waarden

De koppen boven de datakolommen van de min/max waarden worden nu netjes uitgelijnd bij uitvoer van stationsdata in ASCII.

5.1.4

Sid027 - Stationsnaam bij min/max waarden

De stationsnaam wordt nu netjes op dezelfde regel als de data weergegeven bij uitvoer van stationsdata in ASCII.

5.1.5

Sid030 - Padnamen Gui en Getdata

De padnamen worden nu zowel binnen de Gui als binnen Getdata op gelijke, correcte wijze geïnterpreteerd. Er is nog wel een probleem met lange padnamen voor ‘results’ in de GUI.

5.1.6

Sid031 - Inlezen van een Ascii uitvoerreeks in Excel

Het inlezen van een Ascii reeks uit Getdata in een Excel bestand verloopt nu zonder fouten. Aangeraden wordt om hiervoor de volgende procedure te gebruiken: Het inlezen van de data in Excel • Start Excel op, selecteer via !Open het invoerbestand (een import-dialoog verschijnt) • druk op ! volgende (step 2) • druk op ! volgende (step 3) • druk op Geavanceerd ! scheidingsteken duizendtallen moet zijn ‘ ‘ (spatie), ! decimaalteken aanpassen in ‘.’ (punt). ! OK

©


pagina 19 van 26


•


! Voltooien (nu is het bestand, als het goed is correct ingelezen)

Het maken van een afbeelding (tijdreeks) in Excel • selecteer Y-kolommen door te klikken op de “kolomkoppen” • activeer grafiek middels icoon (dialoog verschijnt) • ! kies TAB Reeks en selecteer X-as gegevens door te klikken op eindknop naast invulveld (een twee-regelige wizard verschijnt) • selecteer in Excel-sheet X-kolom door weer op de kop te klikken en activeer deze setting door de wizard te beëindigen (druk op kruis rechtsboven). • toon tijdreeks door ! volgende, ! volgende, ! voltooien, waarna afbeelding door gebruiker verder kan worden afgemaakt. Verder zijn in de ASCII-reeks uitvoer extra (variabele) koppen boven de betreffende kolommen gerealiseerd.

5.1.7

Sid032 - Parameter -d gis in commandoregel

Als naast de optie ‘-d gis’ in de commandoregel van Getdata een Sds-bestand stond met een ‘-d’ in de naam, bijvoorbeeld ‘SDS-d01’, dan werd de commandoregel niet goed door Getdata verwerkt. Dit euvel is nu verholpen.

5.1.8

Sid034 - Runscript getdata.bat

Er is een runscript getdata.bat gemaakt, waardoor het aanroepen van Getdata nu eenvoudiger geworden is.

5.1.9

Sid037 - Verschillen tussen Arcview en Waqview

Er traden verschillen op tussen de waterstandvelden van de postprocessingsprogramma’s Waqview en Arcview. Na onderzoek is gebleken dat de Simarc-uitgang binnen Getdata anders omgaat met deze velden dan Waqview. Door een aanpassing aan de Simarc-uitgang van Getdata is deze ongelijkheid ongedaan gemaakt, waardoor beide programma’s nu dezelfde resultaten te zien geven.

5.2 Nieuwe ontwikkelingen 5.2.1

Ontsluiten informatie m.b.t. dynamische barriers

De informatie met betrekking tot het gedrag van dynamische barriers tijdens een simulatie is beschikbaar als zogenaamde “history-data” op de Waqua SDS-file. De naam van de betreffende array is RRSBAH, dimensie (7,nsluv). “nsluv” is het aantal barriers. Per barrier zijn de volgende 7 gegevens beschikbaar: 1. tijdstap 2. geprefereerde drempelhoogte 3. actuele drempelhoogte 4. geprefereerde schuifhoogte 5. actuele schuifhoogte 6. geprefereerde doorstroombreedte 7. actuele doorstroombreedte Het gedrag van de dynamische barriers kan via getdata als volgt worden opgevraagd: getdata -f <SDS-file> -s -v RRSBAH

©


pagina 20 van 26



De uitvoer bestaat per barrier uit een standaard header met vervolgens per regel de gelezen tijdstap van de SDS-file, gevolgd door de 7 bovenstaande gegevens. Via de SIDONIA gui kan de barrier informatie (RRSBAH) worden opgevraagd via het “Barriers” profiel. Als kandidaat stations worden daarin de namen van de beschikbare barriers getoond.

5.2.2

Het kunnen ophalen van de initiële tijdstap voor veldwaarden

Het opvragen van de initiële tijdstap van velddata was met SIDONIA niet mogelijk. Dit werd veroorzaakt door de wijze waarop de –t optie werd geïnterpreteerd in getdata. Er kon alleen een reguliere tijdreeks worden opgevraagd vanaf tfmapf t/m tlmapf met een opgegeven tijdstap (default timapf). De starttijd van de simulatie, tstart, kan hier echter vóór liggen. Ook is het mogelijk dat er restart-tijdstippen zijn die buiten de reguliere map-tijdreeks liggen (op restarttijdstippen wordt ook map-data geschreven). Tenslotte is het nog mogelijk dat (een gedeelte van) de map tijdstippen niet op de SDS staan, omdat bijvoorbeeld de simulatie nog loopt. De gehele tijdstippen afhandeling is gewijzigd in getdata. Er wordt niet meer gekeken naar de waarden voor tfmapf e.d., maar er wordt uitgegaan van de tijdstippen die impliciet in de SIMONA administratie voor tijdafhankelijke arrays is opgeslagen. Getdata gaat nu als volgt om met de invoer via de –t optie. •

Geen –t optie of geen waarden opgegeven bij de –t optie : Alle beschikbare tijdstippen worden uitgevoerd. Getdata vraagt de beschikbare tijdstippen op aan Dirsys (SIMONA).

•

-t : Alle beschikbare tijdstippen vanaf worden uitgevoerd.

•

-t , : Alle beschikbare tijdstippen vanaf t/m worden uitgevoerd.

•

-t ,, : Alle beschikbare tijdstippen vanaf t/m met interval worden uitgevoerd. Het gaat hier niet om de dichtstbijzijnde tijdstippen, maar om exact die tijdstippen die voldoen aan de gevraagde reguliere reeks (exact met een epsilon verschil).

Indien er om verschil-data wordt gevraagd door een tweede SDS-file te specificeren dan worden de tijdstippen bepaald uit de doorsnede van de beschikbare tijdstippen in beide SDS-files. Hierbij wordt geen rekening gehouden met de referentiedatum (itdate in geval van Waqua, niet beschikbaar i.g.v. Slib3D).

5.2.3

Ontsluiten van Slib3D SDS-files

SIDONIA is zodanig opgezet dat het in principe mogelijk is om ieder type SDS file te benaderen. Aanpassen van de dlds, gegevenswoordenboek en de gui views zou voldoende moeten zijn. Dit bleek echter niet voldoende: bij NetCDF en Ascii uitvoer wordt er een header gevuld met daarin o.a. de referentiedatum. Deze is niet beschikbaar in Slib3D SDSfiles en getdata liep daarop stuk. Nadat dit probleem was verholpen (default

©


pagina 21 van 26



referentiedatum ingevuld) was het mogelijk om ook Slib3D SDS-files via SIDONIA te benaderen.

5.2.4

Ontsluiten van wind SDS-files

Net zoals bij het ontsluiten van Slib3D SDS-files geldt ook hier dat aanpassen van de dlds, gegevenswoordenboek en de gui views voldoende zou moeten zijn. Ook hier bleek dit niet voldoende vanwege de referentiedatum. Deze is wel beschikbaar in wind SDS-files, maar heeft een andere naam (refdat i.p.v. itdate). Nadat dit probleem was verholpen was het mogelijk om ook wind SDS-files via SIDONIA te benaderen.

5.2.5

Implementatie van de verbeterde afvoerpotentiaal

De code voor de verbeterde afvoerpotentiaal is aangeleverd door RIKZ (Robert Vos, Dirsys routine dic004). Deze routine heeft een gewijzigde parameterlijst t.o.v. de voorgaande versie. Dit is aangepast in Dirsys en de dlds file. De belangrijkste wijziging ten aanzien van SIDONIA betrof het kunnen invoeren van gebruikerswaarden voor parameters t.b.v. Dirsys “compute” variabelen. Dirsys zelf biedt die mogelijkheid, maar de parameter waarden moeten via getdata op gebruiksvriendelijke manier gezet kunnen worden. In getdata is het definiëren van parameterwaarden voor Dirsys “compute” variabelen als volgt mogelijk gemaakt: getdata –f <SDS-file> -v ‘psi(ioptsc=3,mg=100,ng=50)’

Ofwel: Per opgevraagde variabele kunnen parameterwaarden worden gespecificeerd door deze tussen ronde haken achter de variabele naam te zetten. Een parameterwaarde wordt opgegeven als <parameternaam>=<parameterwaarde> en de verschillende parameterwaarde specificaties moeten worden gescheiden door een komma. De consequentie van het gebruik van ronde haken is dat minimaal dat onderdeel tussen enkele quotes moet staan bij gebruik op Linux. Anders wil de Linux shell er iets mee doen en resulteert dat in een fout. In getdata worden de parameterwaarden in de Dirsys administratie gezet direct na het parsen van de opdrachtregel. Hierdoor is het ook beschikbaar voor de GIS uitvoer en voor de verschilfuncties.

5.2.6

Het kunnen ophalen van verschillen tussen twee SDS-files

De nieuwe versie van getdata kan voor alle gewenste uitvoerformaten nu ook de verschillen uitvoeren tussen twee SDS files. De aansturing hiervoor is als volgt: getdata –f <SDS 1>,<SDS 2> <etc>

Ofwel : de optie –f gevolgd door twee SDS-filenamen gescheiden door een komma. Indien getdata gevraagd wordt om verschillen te berekenen wordt het volgende gedaan: • Verschillen worden alleen berekend voor variabelen die voorkomen in SIMONA characteristics beginnend met SOLUTION, RESTART of HISTORY en voor alle “compute” variabelen. Dit is eenvoudig uit te breiden. In geval van GIS uitvoer worden verschillen berekend voor VANGLE, VSPEED, WLEVEL, PSI en CONCENTRATIONS. ©


pagina 22 van 26


•

• • • •

•


Voor alle andere variabelen en voor alle character variabelen worden de waarden uit de eerste SDS-file teruggegeven. Het verschil tussen VANGLE en VSPEED waarden wordt bepaald door de hoek resp. snelheid eerst uit te rekenen en vervolgens af te trekken. Bij VANGLE vindt geen correctie plaats om het resultaat in een bepaalde range te brengen, zoals 0-360 graden. Van de beschikbare tijdstappen op de twee SDS-files wordt een doorsnede bepaald, waarna de gebruikersinvoer m.b.t. de tijden wordt verwerkt. Er wordt geen rekening gehouden met een referentiedatum. Bij het berekenen van de verschillen worden de waarden op dezelfde tijdstippen (in minuten) op de twee SDS-files van elkaar afgetrokken. Bij het opvragen van meta-data worden alleen de gegevens van de eerste SDSfile uitgevoerd. In geval van mapdata wordt droogval in de verschillen verwerkt door in de eerste SDS-file alle permanent droge punten op een dummy waarde te zetten en alle tijdelijk droge punten op 0.0. In de tweede SDS-file worden zowel alle permanente als tijdelijk droge punten op 0.0 gezet. In het verschil krijgen daarmee alle permanent droge punten in de eerste SDS-file een dummy waarde. Bij tijdreeksen wordt de positie van een gevraagd station in het gevraagde variabele-array (bijv. ZWL) bepaald aan de hand van de eerste SDS-file. Dezelfde positie wordt gehanteerd in de tweede SDS-file. De variabele-array zelf moet in beide SDS-files dezelfde dimensie hebben. Met andere woorden: om verschillen tussen gelijke stations te kunnen bepalen moet in beide SDSfiles dezelfde reeks stations staan in dezelfde volgorde.

Voor ontwikkelaars: De basis voor de implementatie van verschilberekeningen vormt een uitbreiding van Dirsys. Er zijn een 7-tal nieuwe functies aan Dirsys toegevoegd, gebaseerd op bestaande functies, die een verschil kunnen uitrekenen tussen twee gelijksoortige variabelen op twee verschillende SDS-files. Deze variabelen dienen gelijke dimensie te hebben, gelijk type, en moeten op dezelfde tijdstap beschikbaar zijn in beide SDSfiles. Daar wordt in de nieuwe functies op getest. Het verschil tussen de nieuwe routines en de routines waarop ze zijn gebaseerd is in de aanroep alleen twee extra parameters: de naam van de tweede SDS-file en het experiment. De nieuwe routines zijn zodanig opgezet dat ze de bestaande routines kunnen vervangen. Als namelijk geen naam wordt opgegeven voor de tweede SDSfile, dan retourneren de routines de gelezen waarden van de eerste SDS-file. De nieuwe routines zijn: Nieuwe routine didiar.f didipa.f didrar.f didrpa.f ididsi.f rdidsi.f didasc.f

Gebaseerd op digiar.f digipa.f digrar.f digrpa.f idigsi.f rdigsi.f dipasc.f

De naamgeving is in het algemeen gewijzigd van “dig” (Dirsys Get) in “did” (Dirsys Difference).

©


pagina 23 van 26


6


Configuratie

6.1 Platforms Simona wordt op drie platforms ondersteund: • HP9000/700 onder HP-UX 11 • PC (cluster) onder Linux • PC onder Windows XP of Windows 2000

6.2 Configuratie-eisen Voor informatie over de configuratie-eisen, zie de Installatiegids SIMONA. Enkele hoofdpunten worden hieronder vermeld. Het Simona pakket is gegenereerd en getest in de volgende omgevingen: Platform Compiler Bijzonderheden Linux - Suse Gnu vs 3.2 Linux - Suse Intel vs 9.0 Alleen voor Rikz. Linux - Slackware 9 Intel vs 8.0 HP9000 onder HP-UX 11 HP compiler vs B.11.11 Draait niet meer parallel. Windows XP Visual Studio vs ++6.0 Draait nog niet parallel. Het Sidonia pakket is gegenereerd en getest in de volgende omgevingen: Platform Compiler Bijzonderheden Linux - Suse Lahey Fortran vs 6.1 Java JRE vs 1.5.0 Nodig voor het draaien. Windows XP Visual Studio ++6.0 Nodig voor generatie. Oracle J-Developer Nodig voor generatie. vs 10.1.2.0.0 (1811) Java JRE vs 1.5.0 Nodig voor het draaien.

6.3 Producten De standaard Simona levering bestaat uit de onderstaande onderdelen: • Simona productie-versie • Documentatie in PDF • Release Notes

©


pagina 24 van 26



6.4 Bestanden Hieronder volgt een opsomming en uitleg van de belangrijkste bestanden uit de Simona-omgeving. Met $SIMONADIR wordt de productie-directory van Simona op het lokale systeem bedoeld. Bestandsnaam

Omschrijving

$SIMONADIR/bin/*.exe $SIMONADIR/bin/*.run $SIMONADIR/bin/*.hlp $SIMONADIR/bin/*.arr $SIMONADIR/bin/simona.env* $SIMONADIR/bin/SIMETF $SIMONADIR/util/sicheck* $SIMONADIR/util/simkf* $SIMONADIR/util/keywords* $SIMONADIR/util/sireft* referentiebestanden $SIMONADIR/util/sirecovr* $SIMONADIR/util/sidsview* $SIMONADIR/lib/*.a $SIMONADIR/lib/*.o $SIMONADIR/lib/*.f $SIMONADIR/lib/*.i $SIMONADIR/usedoc/readme.txt $SIMONADIR/usedoc/waqua $SIMONADIR/usedoc/donsim $SIMONADIR/usedoc/tools

programma-executables runprocedures helpinformatie bestanden referentiebestanden waqriv/waqua simona environmentfile (ascii) simona foutmeldingenbestand (binair) utility sicheck (source controle) utility simkf (aanmaken makefile) hulpbestanden keyword controle utility voor het aanmaken van

$SIMONADIR/sysdoc/waqua $SIMONADIR/sysdoc/donsim $SIMONADIR/test/* .* $SIMONADIR/testbank

©


utility herstellen sds-bestanden utility bekijken sds-bestanden (statische) bibliotheken objects hoofdprogramma's source hoofdprogramma's include-bestanden fortran beschrijving afdrukken documentatie WAQUA User's Guide DONSIM User's Guide SIMONA Programmer’s Guide, systeemdocumentatie van de tools LDS WAQUA, sysdoc WAQUA LDS DONSIM, sysdoc DONSIM invoerbestanden testmodellen Beschrijving testbankmodellen

pagina 25 van 26


7


Versienummers

Na Export2005-05 gelden de volgende versienummers; COPPRE COEXEC COPPOS DIFTOSDS SDSTOINC SIDSVIEW SIERWR SIMPAR SIRECOVR SIREFT SLIB3D WAQBHD WAQCLV IPW WAQPAN WAQPRE WAQPRO WAQRIV WAQVIEW WAQWND

1.16.00 1.13.00 1.08.00 1.05.03 1.04.10 3.07.05 3.02.10 1.10.00 3.03.10 3.05.00 1.07.03 1.10.03 1.00.02 3.08.00 8.36.05 9.24.00 9.28.00 8.31.05 3.04.00 9.02.03

Opmerkingen: • De versienummers zijn terug te vinden in de message files van de deelsystemen. • Waqview kent geen message file, dus is de versie alleen te zien op het display. • De handleiding van de Ipw is nog niet geheel bijgewerkt, en beschrijft daarom nog de Ipw versie 3.06.

©


pagina 26 van 26

SIMONA B&O. Release Notes Export

Recommend Documents