IPDS.DAT

Basisinformationen

Datei-Typ

ipds.dat

Datei-Form

FORMATTED

Version

1.3

Beschreibung

Juni 2007

Bedeutung der Datei

Eine Datei des Typs ipds.dat enthält Initialwerte für physikalischeDatensätze Die Möglichkeiten, die damit verbunden sind (z.B. zur Definition von Anfangssituationen), mögen folgende Beispiele verdeutlichen. Neuerungen mit Version 1.3 sind wie folgt markiert: Neu in V1.3

Datei-Inhalt

1. Block date_and_time
   • Datumsangabe
   • Es ist genau ein Block für die Datumsangabe erlaubt
   • Ein-Zeilen-Angabe des Datums mit Uhrzeit mit folgenden erlaubten Formaten:

   • DD.MM.YYYY-hh:mm:ss
   • DD.MM.YYYY-hh:mm:ss.nnnnnnnnn
   • DD.MM.YYYY-hh:mm:ss.nnnnnnnnn-ZONE

   • Legende:

   • DD = Tag des Monats (01-31)
   • MM = Monat des Jahrs (01-12)
   • YYYY = Jahr (vierstellig)
   • hh = Stunde des Tages (00-23)
   • mm = Minute der Stunde (00-59)
   • ss = Sekunde der Minute (00-59)
   • nnnnnnnnn = Nanosekundenanteil (9stellig)
   • ZONE = Zeitzone (UTC/MEZ/MESZ/+3H/..)

   • Keyzeile lautet datetime = DD.MM.YYYY-hh:mm:ss.nnnnnnnnn-ZONE
2. Block global_constant_values
   • Definition von physikalischen Konstantwerten für das Gesamtgebiet
   • Genau ein Block mit Default-Einstellungen für verschiedene physikalische Größen ist erlaubt
   • Die Default-Einstellungen können von den nachfolgenden regional gültigen Einstellungen (Block regional values) überschrieben werden
   • Implementierte physikalische Größen:

   • Wasserstand

   • Keyzeile: waterlevel = ..wert..
   • Genau eine Wertangabe möglich

   • Salzgehalt

   • Keyzeile: salinity = ..wert..
   • Genau eine Wertangabe möglich

   • Temperatur

   • Keyzeile: temperature = ..wert..
   • Genau eine Wertangabe möglich

   • Schwebstoffgehalt

   • Keyzeile für eine Schwebstoff-Fraktion: suspendedload = ..sediment_name.. ..gehalt..
   • Es können beliebig viele Schwebstoff-Fraktionen angegeben werden.

   • Sediment-Fraktionen und Porosität

   • Keyzeile für eine Sediment-Fraktion: sediment_fraction = ..sediment_name.. ..anteil..
   • Es können beliebig viele Sediment-Fraktionen beschrieben werden
   • Keyzeile für die Porosität: porosity = ..anteil..
   • Für die Porosität ist genau eine Wertangabe möglich
   • Die hier verwendeten Sedimentbezeichnungen müssen mit denen aus der verwendeten Sediment-Klassifikations-Datei (Typ soil.dat) übereinstimmen
   • Die Summe aller relativen Vorkommen der Sedimentklassen muss 1 sein

   • Riffelwellenzahl

   • Die Wellenzahl k bestimmt Richtung und Wellenlänge (λ) der Riffel (k = 2 π / λ)
   • Zur Definition einer Wellenzahl werden seine x- und y-Komponente angegeben
   • Keyzeile: ripple_wavenumber = ..x-komponente.. ..y-komponente..
   • Genau eine Angabe zur Riffelwellenzahl möglich

   • Riffelhöhe

   • Keyzeile: ripple_height = ..wert..
   • Genau eine Wertangabe möglich

   • Dünenwellenzahl

   • Die Wellenzahl k bestimmt Richtung und Wellenlänge (λ) der Dünen (k = 2 π / λ)
   • Zur Definition einer Wellenzahl werden seine x- und y-Komponente angegeben
   • Keyzeile: dune_wavenumber = ..x-komponente.. ..y-komponente..
   • Genau eine Angabe zur Dünenwellenzahl möglich

   • Dünenhöhe

   • Keyzeile: dune_height = ..wert..
   • Genau eine Wertangabe möglich

   • Tiefenlage des unerodierbaren Bodens

   • Die Angaben zur Tiefenlage des unerodierbaren Boden können absolut ("absolute") oder relativ ("relative") vorgenommen werden
   • Absolute Angaben werden dann gemacht, wenn sich der unerodierbare Boden in der angegebenen (absoluten) geodätischen Höhe befinden soll
   • Die Angabe kann auch relativ zur ungestörten Wassertiefe erfolgen, der unerodierbare Boden befindet sich dann um den angegebenen Betrag unterhalb (negative Werte) des Wasserkörper-Bodens
   • Keyzeile: rigid_layer_depth = relative|absolute ..wert..
   • Genau eine Angabe zur Tiefenlage des unerodierbaren Bodens möglich

   • Rauheit des unerodierbaren Bodens

   • Die Rauheit des unerodierbaren Bodens kommt immer dann zum Tragen, wenn keine weiteren Sedimente auf dem unerodierbaren Boden liegen
   • Zwei Angaben zur Definition der Rauheit sind erforderlich:

   • Angabe zum verwendeten Bodenschubspannungsgesetz (z.Z. ist nur "nikuradse" implementiert)
   • Angabe des Rauheitswertes

   • Keyzeile: rigid_layer_roughness = nikuradse ..wert..
   • Genau eine Angabe zur Rauheit des unerodierbaren Bodens möglich

   • Rauheit des Bodens

   • Keyzeile: z0_roughnesslength = ..wert..

3. Block region
   • Definition von Regionen
   • Es können beliebig viele Regionen definiert werden
   • Es muss mindestens eine Region definiert werden, wenn abweichend von den Konstantwerten für das Gesamtgebiet regional andere Werte definiert werden sollen
   • Es müssen keine Regionen definiert werden
   • Jede Region braucht einen Namen

   • Dieser Name wird im Folgenden als Identifikationsschlüssel verwendet
   • Der Name sollte sich an geographischen Gegebenheiten orientieren
   • Groß- und Kleinschreibung beachten
   • Die Keyzeile lautet: region_name = ..name..

   • Beschreibung des diese Region umschließenden Randpolygons

   • Das Randpolygon muss aus mindestens drei Punkten bestehen
   • Die maximale Anzahl Randpolygonpunkte ist nicht beschränkt
   • Das Randpolygon wird automatisch geschlossen
   • Die Keyzeile für einen Randpolygonpunkt lautet: border_point = ..x-koordinate.. ..y-koordinate..

4. Block sampling_point
   • Definition von (Mess-)stationen mit Messwerten
   • Es sind beliebig viele Messstationen erlaubt
   • Es muss mindestens eine Messstation definiert werden, wenn abweichend von den Konstantwerten für das Gesamtgebiet regional andere Werte definiert werden sollen
   • Es müssen keine Messstationen definiert werden
   • Jede Messstationen braucht einen Namen

   • Dieser Name wird im Folgenden als Identifikationsschlüssel verwendet
   • Der Name sollte sich an geographischen Gegebenheiten orientieren
   • Groß- und Kleinschreibung beachten
   • Die Keyzeile lautet: sampling_point_name = ..name..

   • Positionsangabe der Messstation

   • Angabe der x- und y-Koordinate
   • Nur eine Positionsangabe verwenden
   • Die Keyzeile lautet: sampling_point_xy = ..x-koordinate.. ..y-koordinate..

   • Definition der an der Messstation gültigen Werte der unterschiedlichen physikalischen Größen

   • Es können die gleichen physikalischen Größen wie im Block global constant values beschrieben verwendet werden

5. Block regional_values
   • Definition von regionel abweichenden physikalischen Werten
   • Spezifizierung von Regionen, für die die Werte der physikalischen Größen an Messstationen abgelegt sind
   • Für die Zwischenräume der Regionen werden die Werte der physikalischen Größen interpoliert
   • Es sind beliebig viele dieser Blöcke erlaubt
   • Es müssen keine Blöcke definiert werden (die Konstantwerte für das Gesamtgebiet werden dann nicht verändert)
   • Die Blöcke sind so zu sortieren, dass die Blöcke mit der geringsten Priorität zuerst erfolgen (nur relevant bei Gebietsüberschneidungen)
   • Angabe der Region

   • Die Keyzeile lautet: region = inside|outside ..regionname..
   • Der Name der Region (zweite Angabe) muss übereinstimmen mit einem der oben definierten Schlüsselnamen für Regionen
   • Zusätzlich wird kenntlich gemacht, ob die hier gewünschte Region innerhalb ("inside") oder außerhalb ("outside") des definierten Randpolygons liegt
   • Definition der geodätischen Schicht, in der die angegebenen Werte gültig sind

   • Die Keyzeilen lauten: lower_validity_level = ..unterer_Grenzwert..
     upper_validity_level = ..oberer_Grenzwert..
   • Die Angabe beider Werte ist optional

   • Angabe der Messstationen

   • Es können beliebig viele Messstationen angegeben werden
   • Es muss mindestens eine Messstation angegeben werden
   • Bei Angabe von genau einer Messstation erhält man Konstantwerte für die ganze Region
   • Der Name jeder Messtation muss übereinstimmen mit einem der oben definierten Schlüsselnamen für Messstationen
   • Die Keyzeile lautet: sampling_point = ..messstationname..

   • Angabe des max. Abstands der Messstationen von einem Punkt P. Nur diejeniegen Messstationen, die sich innerhalb des Abstandsradius um P befinden, werden für die Interpolation der Messdaten verwendet.

   • Es kann nur ein max. Abstand angegeben werden
   • Der Wert für den max. Abstand muss größer Null sein
   • Die Keyzeile lautet: sampling_point_maxdist = ..maximaler_Abstand_zu_Messstationen..
   • Die Angabe der Zeile ist optional. Wird die Zeile nicht angegeben, wird der Default-Wert (40.E6) verwendet. Damit ist sichergestellt, dass sich mit IPDS-Dateien der Version 1.1 trotz neuer Programmversion identische Resultate erzielen lassen.
   • Die folgenden Beispielgrafiken sollen den Einfluss des max. Abstandes dokumentieren: Dargestellt ist in allen Grafiken die aus 1382 Einzelmessungen berechnete Verteilung des Feinsandes (.063-.177 mm Korndurchmesser) in der Emsmündung, im Dollart und in der Unterems. Die Einzelmessungen sind dem Sedimentatlasdes Rijksinstituut voor Kust en Zee (RIKZ) entnommen. Der Farbverlauf bildet den rel. Anteil des Feinsandes ab: Rot entspricht 100%, Gelb 75%, Grün 50% und Dunkelblau 0% Feinsand-Anteil.

• Der max. Abstand zu den Messstationen (300m) entspricht der internen Gitterauflösung
  

Die Interpolation findet nur in eng begrenzter Nähe zu den Messstationen statt. Für alle weiteren Gitternetz-Punkte werden keine Messstationen gefunden. Die Lage der Messstationen kommt deutlich zum Vorschein. (642x768 Pixel, 21kB).

• Der max. Abstand zu den Messstationen (1000m) beträgt etwa das 3-4 fache der internen Gitterauflösung
  

Die Interpolation findet flächendeckend dort statt, wo Messstationen nicht weit entfernt sind. Das Ergebnis ist dort in erster Näherung naturähnlich. In weiten Bereichen der Emsmündung findet keine Interpolation statt, die Messpunkte sind dort zu weit entfernt. (642x768 Pixel, 30kB).

• Ohne Angabe eines max. Abstandes Die Interpolation findet flächendeckend für jeden Punkt des Gitternetzes statt. Weit entfernt von der Lage der Messstationen ist das Ergebnis stark von künstlichen Strukturen geprägt. (642x768 Pixel, 85kB).

• Angabe des Interpolationsverfahrens

• Die Keyzeile lautet:interpolation_method = nearest_points_in_sectors|triangular_interpolation
• Zur Zeit sind nur diese zwei Verfahren implementiert
• Das Interpolationsverfahren nearest_points_in_sectors ist ein im Haus der BAW entwickeltes Verfahren, welches für einen gewünschten Ort möglichst viele, möglichst dicht benachbarte Messstationen zur Interpolation heranzieht, wenn diese Messstationen in unterschiedlichen Richtungen zu finden sind. Die physikalischen Werte der Messstationen werden gewichtet zur Mittelung heran gezogen. Die Wichtung erfolgt über das Quadrat der Abstände zu den Messstationen.
• Neu in V1.3 Für das Interpolationsverfahren triangular_interpolation werden alle zur Verfügung stehenden Messstationen zu einem Dreiecksgitter vermascht. Die Vermaschung erfolgt nach Delaunay-Kriterien und bildet eine konvexe Hülle um alle Messstationen. Die Vermaschung wird im Haus der BAW mit Hilfe des externen Software-Paketes GEOMPACK (Autor: Barry Joe) generiert und über die Schnittstelle ext_ipds_setup_mespos_tria bereitgestellt.
• Andere Verfahren können integriert werden

Programme, welche diesen Datei-Typ benutzen

TELEMAC-2D [SediMorph], UNS, UNTRIM, UNTRIM2007

Beispiel-Datei

$PROGHOME/examples/ipds/ipds.dat

---

zurück zu: Dateikennblätter

---

Strukturübersicht