Aktionen

NetCDF Dreiecksgitter

Aus BAWiki

Version vom 26. April 2011, 15:39 Uhr von imported>Lang Guenther (→‎Polygone: fehlenden Hinweis ergänzt)

Kurze Beschreibung

Ein aus Dreiecken aufgebautes Gitter wird beschrieben. Typische Dateitypen hierfür sind z. B. gitter05.dat, gitter05.bin sowie selafin. Es könnte sich aber auch um ein spezielles, ausschließlich aus Dreiecken aufgebautes Gitter des Typs untrim_grid.dat handeln. Die Speicherung von Daten auf diesem Gitter, z. B. der Tiefe (Bathymetrie), wird an anderer Stelle beschrieben (siehe NetCDF Synoptische Daten im Dreiecksgitter. Ebenso werden hier keine Informationen über offene und geschlossene Kanten oder Positionen zum Einsteuern von Randwerten aufgeführt. Es wird ausschließlich die Struktur (Topologie) des Gitters beschrieben.

Die Beschreibung wurde möglichst vollständig für die staggered data Postionen node, edge und poly gestaltet, auch wenn für einzelne der o. g. Dateitypen nur ein Teil der nachfolgend genannten Daten benötigt wird.

Inhalt der Datei

Dimensionen

  1. nMesh2_node : Anzahl der Knoten.
  2. nMesh2_edge : Anzahl der Kanten.
  3. nMesh2_poly : Anzahl der Polygone (hier Dreiecke).
  4. two : Konstante ( = 2 ).
  5. three : Konstante ( = 3 ).

Lokale Koordinaten

Knoten

double Mesh2_node_x(nMesh2_node) ;

Mesh2_node_x:long_name = "x-coordinate of 2D mesh nodes" ;
Mesh2_node_x:units = "m" ;
Mesh2_node_x:name_id = 1650 ;
Mesh2_node_x:standard_name = "projection_x_coordinate" ;

double Mesh2_node_y(nMesh2_node) ;

Mesh2_node_y:long_name = "y-coordinate of 2D mesh nodes" ;
Mesh2_node_y:units = "m" ;
Mesh2_node_y:name_id = 1651 ;
Mesh2_node_y:standard_name = "projection_y_coordinate" ;

Hinweise:

  1. Das Attribut "name_id" entspricht der (BAW) PHYDEF-Code-Kennung der Variablen.

Kanten

double Mesh2_edge_x(nMesh2_edge) ;

Mesh2_edge_x:long_name = "x-coordinate of 2D mesh edges" ;
Mesh2_edge_x:units = "m" ;
Mesh2_edge_x:name_id = 1650 ;
Mesh2_edge_x:bounds = "Mesh2_edge_x_bnd" ;
Mesh2_edge_x:standard_name = "projection_x_coordinate" ;

double Mesh2_edge_y(nMesh2_edge) ;

Mesh2_edge_y:long_name = "y-coordinate of 2D mesh edges" ;
Mesh2_edge_y:units = "m" ;
Mesh2_edge_y:name_id = 1651 ;
Mesh2_edge_y:bounds = "Mesh2_edge_y_bnd" ;
Mesh2_edge_y:standard_name = "projection_y_coordinate" ;

Hinweise:

  1. Das Attribut "name_id" entspricht der (BAW) PHYDEF-Code-Kennung der Variablen.

Polygone

double Mesh2_poly_x(nMesh2_poly) ;

Mesh2_poly_x:long_name = "x-coordinate of 2D mesh polygons, center of gravity" ;
Mesh2_poly_x:units = "m" ;
Mesh2_poly_x:name_id = 1650 ;
Mesh2_poly_x:bounds = "Mesh2_poly_x_bnd" ;
Mesh2_poly_x:standard_name = "projection_x_coordinate" ;

double Mesh2_poly_y(nMesh2_poly) ;

Mesh2_poly_y:long_name = "y-coordinate of 2D mesh polygons, center of gravity" ;
Mesh2_poly_y:units = "m" ;
Mesh2_poly_y:name_id = 1651 ;
Mesh2_poly_y:bounds = "Mesh2_poly_y_bnd" ;
Mesh2_poly_y:standard_name = "projection_y_coordinate" ;

Hinweise:

  1. Das Attribut "name_id" entspricht der (BAW) PHYDEF-Code-Kennung der Variablen.
  2. Als Zentrum wird hier der Schwerpunkt des Polygons benutzt, da dieser bei konvexen Polygonen immer innerhalb des Polygons zu liegen kommt. Demgegenüber muss das Zentrum des Umkreises nicht in das Polygon fallen.

(optional) Begrenzung für Kanten

double Mesh2_edge_x_bnd(nMesh2_edge, two) ;

double Mesh2_edge_y_bnd(nMesh2_edge, two) ;

(optional) Begrenzung für Polygone

double Mesh2_poly_x_bnd(nMesh2_poly, three) ;

double Mesh2_poly_y_bnd(nMesh2_poly, three) ;

Hinweise:

  1. Koordinaten der Knoten des Polygons entgegen dem Uhrzeigersinn angeben.

Geografische Koordinaten

Knoten

double Mesh2_node_lon(nMesh2_node) ;

Mesh2_node_lon:long_name = "longitude of 2D mesh nodes" ;
Mesh2_node_lon:units = "degrees_east" ;
Mesh2_node_lon:name_id = 1653 ;
Mesh2_node_lon:standard_name = "longitude" ;

double Mesh2_node_lat(nMesh2_node) ;

Mesh2_node_lat:long_name = "latitude of 2D mesh nodes" ;
Mesh2_node_lat:units = "degrees_north" ;
Mesh2_node_lat:name_id = 1652 ;
Mesh2_node_lat:standard_name = "latitude" ;

Hinweise:

  1. Das Attribut "name_id" entspricht der (BAW) PHYDEF-Code-Kennung der Variablen.

Kanten

double Mesh2_edge_lon(nMesh2_edge) ;

Mesh2_edge_lon:long_name = "longitude of 2D mesh edge centers" ;
Mesh2_edge_lon:units = "degrees_east" ;
Mesh2_edge_lon:name_id = 1653 ;
Mesh2_edge_lon:bounds = "Mesh2_edge_lon_bnd" ;
Mesh2_edge_lon:standard_name = "longitude" ;

double Mesh2_edge_lat(nMesh2_edge) ;

Mesh2_edge_lat:long_name = "latitude of 2D mesh edge centers" ;
Mesh2_edge_lat:units = "degrees_north" ;
Mesh2_edge_lat:name_id = 1652 ;
Mesh2_edge_lat:bounds = "Mesh2_edge_lat_bnd" ;
Mesh2_edge_lat:standard_name = "latitude" ;

Hinweise:

  1. Das Attribut "name_id" entspricht der (BAW) PHYDEF-Code-Kennung der Variablen.

Polygone

double Mesh2_poly_lon(nMesh2_poly) ;

Mesh2_poly_lon:long_name = "longitude of 2D mesh polygons, center of gravity" ;
Mesh2_poly_lon:units = "degrees_east" ;
Mesh2_poly_lon:name_id = 1653 ;
Mesh2_poly_lon:bounds = "Mesh2_poly_lon_bnd" ;
Mesh2_poly_lon:standard_name = "longitude" ;

double Mesh2_poly_lat(nMesh2_poly) ;

Mesh2_poly_lat:long_name = "latitude of 2D mesh polygons, center of gravity" ;
Mesh2_poly_lat:units = "degrees_north" ;
Mesh2_poly_lat:name_id = 1652 ;
Mesh2_poly_lat:bounds = "Mesh2_poly_lat_bnd" ;
Mesh2_poly_lat:standard_name = "latitude" ;

Hinweise:

  1. Das Attribut "name_id" entspricht der (BAW) PHYDEF-Code-Kennung der Variablen.
  2. Als Zentrum wird hier der Schwerpunkt des Polygons benutzt, da dieser bei konvexen Polygonen immer innerhalb des Polygons zu liegen kommt. Demgegenüber muss das Zentrum des Umkreises nicht in das Polygon fallen.

(optional) Begrenzung für Kanten

double Mesh2_edge_lon_bnd(nMesh2_edge,two) ;

double Mesh2_edge_lat_bnd(nMesh2_edge,two) ;

(optional) Begrenzung für Polygone

double Mesh2_poly_lon_bnd(nMesh2_poly,three) ;

double Mesh2_poly_lat_bnd(nMesh2_poly,three) ;

Hinweise:

  1. Koordinaten der Knoten des Polygons entgegen dem Uhrzeigersinn angeben.

Topologische Daten

Von den nachfolgenden topologischen Daten müssen nicht alle vorhanden sein. Es muss alllerdings sichergestellt sein, dass fehlende Verzeichnisse aus den vorhandenen abgeleitet werden können.

Knotenverzeichnis der Kanten

int Mesh2_edge_nodes(nMesh2_edge,two) ;

Mesh2_edge_nodes:long_name = "list of nodes for all edges, start node - end node" ;

Hinweise:

  1. Ein Standardname sollte gelegentlich ergänzt werden.
  2. Der erste Knoten kennzeichnet den Beginn, der zweite das Ende der Kante.
  3. Indizes starten mit Null!

Knotenverzeichnis der Polygone

int Mesh2_poly_nodes(nMesh2_poly, three) ;

Mesh2_poly_nodes:long_name = "list of nodes for all polygons, counterclockwise" ;

Hinweise:

  1. Ein Standardname sollte gelegentlich ergänzt werden.
  2. Knoten des Polygons entgegen dem Uhrzeigersinn angeben.
  3. Indizes starten mit Null!

Kantenverzeichnis der Polygone

int Mesh2_poly_edges(nMesh2_poly,three) ;

Mesh2_poly_edges:long_name = "list of edges for all polygons, counterclockwise" ;

Hinweise:

  1. Ein Standardname sollte gelegentlich ergänzt werden.
  2. Kanten des Polygons entgegen dem Uhrzeigersinn angeben.
  3. Indizes starten mit Null!

Polygonverzeichnis der Kanten

int Mesh2_edge_polys(nMesh2_edge,two) ;

Mesh2_edge_polys:long_name = "list of (adjacent) polygons for all edges - left and right neigbour polygon" ;
Mesh2_edge_polys:_FillValue = fillvalue ;

Hinweise:

  1. Ein Standardname sollte gelegentlich ergänzt werden.
  2. Der erste Wert kennzeichnet das linke, der zweite das rechte Nachbarpolygon.
  3. Mit dem "_FillValue" können Kanten auf dem Rand, die nur ein Nachbarpolygon haben, beschrieben werden.
  4. Indizes starten mit Null!

Topologievariable

int Mesh2 ;

Mesh2:long_name = "name of grid can be stored here" ;
Mesh2:dimensionality = 2 ;
Mesh2:locations = "node edge poly" ;
Mesh2:node_coordinates = "Mesh2_node_lon Mesh2_node_lat" ;
Mesh2:edge_coordinates = "Mesh2_edge_lon Mesh2_edge_lat" ;
Mesh2:poly_coordinates = "Mesh2_poly_lon Mesh2_poly_lat" ;
Mesh2:edge_node_connectivity = "Mesh2_edge_nodes" ;
Mesh2::edge_poly_connectivity = "Mesh2_edge_polys" ;
Mesh2::poly_node_connectivity = "Mesh2_poly_nodes" ;
Mesh2::poly_edge_connectivity = "Mesh2_poly_edges" ;

Hinweise:

  1. Ein Standardname sollte gelegentlich ergänzt werden.
  2. Mit dem Attribut "locations" werden die für das Gitter relevanten Positionen registriert.
  3. In der Praxis müssen nicht alle Verknüpfungslisten (connectivity) angegeben werden, da manche Informationen redundant sind. Sie sollten aber angegeben werden können.
  4. Die orange markierten Bezeichnungen stellen Abweichungen bzw. Ergänzungen zu den Deltares-CF-Erweiterungen dar.

Koordinatentransformation

Siehe unter NetCDF Grid Mapping Variable.

Anmerkungen, Fragen

zurück zu NetCDF

Strukturübersicht

Abgerufen von „https://wiki.baw.de/de/index.php?title=NetCDF_Dreiecksgitter&oldid=5458“