Aktionen

NetCDF Unstrukturiertes Gitter mit SubGrid

Aus BAWiki

Version vom 14. September 2010, 07:03 Uhr von imported>Lang Guenther (→‎SubPolygone)

Kurze Beschreibung

Ein aus beliebigen Polygonen (Dreiecke und Vierecke) aufgebautes Gitter wird beschrieben. Zusätzlich sind typische Daten des SubGrid vorhanden. Der Dateityp hierfür ist utrsub_grid.dat. Die Speicherung von Daten auf diesem Gitter, z. B. der Tiefe (Bathymetrie), wird an anderer Stelle beschrieben (siehe NetCDF Synoptische Daten im unstrukturierten Gitter mit SubGrid). 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 vollständig für die staggered data Positionen node, edge, poly, subedge und subpoly gestaltet.

Inhalt der Datei

Dimensionen

  • Siehe NetCDF Unstrukturiertes Gitter, Abschnitt "Dimensionen".
  • Zusätzlich zu den dort angegebenen Dimensionen werden noch benötigt:
    dimensions:
    nMesh2_subnode = number of sub grid nodes ;
    nMesh2_subedge = number of sub edges ;
    nMesh2_subpoly = number of sub polygons ;
    nMaxMesh2_edge_subedges = maximum number of sub edges along an edge ;
    nMaxMesh2_poly_subpolys = maximum number of sub polygons within a polygon ;
    nMaxMesh2_subpoly_subnodes = maximum number of nodes within a subpolygon ;

Lokale Koordinaten

Berechnungsgitter

SubGrid

(optional) SubKnoten

double Mesh2_subnode_x(nMesh2_subnode) ;
Mesh2_subnode_x:standard_name = "projection_x_coordinate" ;
Mesh2_subnode_x:long_name = "x-coordinate of 2D mesh sub grid nodes" ;
Mesh2_subnode_x:units = "m" ;
double Mesh2_subnode_y(nMesh2_subnode) ;
Mesh2_subnode_y:standard_name = "projection_y_coordinate" ;
Mesh2_subnode_y:long_name = "y-coordinate of 2D mesh sub grid nodes" ;
Mesh2_subnode_y:units = "m" ;
Hinweis: Prüfen, ob diese Daten überhaupt erforderlich sind. Den SubKnoten werden keine Daten zugeordnet. SubPolygone und SubKanten werden als Zellen durch ihre Abmessungen (Bounds) beschrieben. Falls kein Attribut FillValue in den Bounds-Variablen zulässig ist, so werden die SubKnoten in jedem Fall benötigt.

SubKanten

double Mesh2_subedge_x(nMesh2_subedge) ;
Mesh2_subedge_x:standard_name = "projection_x_coordinate" ;
Mesh2_subedge_x:long_name = "x-coordinate of 2D mesh sub grid edges, e.g. center of sub edge" ;
Mesh2_subedge_x:units = "m" ;
Mesh2_subedge_x:bounds = "Mesh2_subedge_x_bnd"
double Mesh2_subedge_y(nMesh2_subedge) ;
Mesh2_subedge_y:standard_name = "projection_y_coordinate" ;
Mesh2_subedge_y:long_name = "y-coordinate of 2D mesh sub grid edges, e.g. center of sub edge" ;
Mesh2_subedge_y:units = "m" ;
Mesh2_subedge_y:bounds = "Mesh2_subedge_y_bnd"

SubPolygone

double Mesh2_subpoly_x(nMesh2_subpoly) ;
Mesh2_subpoly_x:standard_name = "projection_x_coordinate" ;
Mesh2_subpoly_x:long_name = "x-coordinate of 2D mesh sub polygons, e.g. any position within sub polygon" ;
Mesh2_subpoly_x:units = "m" ;
Mesh2_subpoly_x:bounds = "Mesh2_subpoly_x_bnd"
double Mesh2_subpoly_y(nMesh2_subpoly) ;
Mesh2_subpoly_y:standard_name = "projection_y_coordinate" ;
Mesh2_subpoly_y:long_name = "y-coordinate of 2D mesh sub polygons, e.g. any position within sub polygon" ;
Mesh2_subpoly_y:units = "m" ;
Mesh2_subpoly_y:bounds = "Mesh2_subpoly_y_bnd"

Begrenzung für SubKanten

double Mesh2_subedge_x_bnd(nMesh2_subedge,Two) ; \\ no metadata required
double Mesh2_subedge_y_bnd(nMesh2_subedge,Two) ; \\ no metadata required

Begrenzung für SubPolygone

double Mesh2_subpoly_x_bnd(nMesh2_subpoly,nMaxMesh2_subpoly_subnodes) ; \\ no metadata required
Mesh2_subpoly_x_bnd:_Fillvalue = fillvalue ; \\ for missing subnodes
double Mesh2_subpoly_y_bnd(nMesh2_subpoly,nMaxMesh2_subpoly_subnodes) ; \\ no metadata required
Mesh2_subpoly_y_bnd:_Fillvalue = fillvalue ; \\ for missing subnodes

Geografische Koordinaten

Berechnungsgitter

SubGrid

(optional) SubKnoten

double Mesh2_subnode_lon(nMesh2_subnode) ;
Mesh2_subnode_lon:standard_name = "longitude" ;
Mesh2_subnode_lon:long_name = "longitude of 2D mesh sub grid nodes" ;
Mesh2_subnode_lon:units = "degrees_east" ;
double Mesh2_subnode_lat(nMesh2_subnode) ;
Mesh2_subnode_lat:standard_name = "latitude" ;
Mesh2_subnode_lat:long_name = "latitude of 2D mesh sub grid nodes" ;
Mesh2_subnode_lat:units = "degrees_north" ;
Hinweis: Prüfen, ob diese Daten überhaupt erforderlich sind. Den SubKnoten werden keine Daten zugeordnet. SubPolygone und SubKanten werden als Zellen durch ihre Abmessungen (Bounds) beschrieben. Falls kein Attribut FillValue in den Bounds-Variablen zulässig ist, so werden die SubKnoten in jedem Fall benötigt.

SubKanten

double Mesh2_subedge_lon(nMesh2_subedge) ;
Mesh2_subedge_lon:standard_name = "longitude" ;
Mesh2_subedge_lon:long_name = "longitude of 2D mesh sub grid edges, e.g. center of sub edge" ;
Mesh2_subedge_lon:units = "degrees_east" ;
Mesh2_subedge_lon:bounds = "Mesh2_subedge_lon_bnd" \\ optional
double Mesh2_subedge_lat(nMesh2_subedge) ;
Mesh2_subedge_lat:standard_name = "latitude" ;
Mesh2_subedge_lat:long_name = "latitude of 2D mesh sub grid edges, e.g. center of sub edge" ;
Mesh2_subedge_lat:units = "degrees_north" ;
Mesh2_subedge_lat:bounds = "Mesh2_subedge_lat_bnd" \\ optional

SubPolygone

double Mesh2_subpoly_lon(nMesh2_subpoly) ;
Mesh2_subpoly_lon:standard_name = "longitude" ;
Mesh2_subpoly_lon:long_name = "longitude of 2D mesh sub polygons, e.g. any position within sub polygon" ;
Mesh2_subpoly_lon:units = "degrees_east" ;
Mesh2_subpoly_lon:bounds = "Mesh2_subpoly_lon_bnd" \\ optional
double Mesh2_subpoly_lat(nMesh2_subpoly) ;
Mesh2_subpoly_lat:standard_name = "latitude" ;
Mesh2_subpoly_lat:long_name = "latitude of 2D mesh sub polygons, e.g. any position within sub polygon" ;
Mesh2_subpoly_lat:units = "degrees_north" ;
Mesh2_subpoly_lat:bounds = "Mesh2_subpoly_lat_bnd" \\ optional

(optional) Begrenzung für SubKanten

double Mesh2_subedge_lon_bnd(nMesh2_subedge,Two) ; \\ no metadata required
double Mesh2_subedge_lat_bnd(nMesh2_subedge,Two) ; \\ no metadata required

(optional) Begrenzung für SubPolygone

double Mesh2_subpoly_lon_bnd(nMesh2_subpoly,nMaxMesh2_subpoly_subnodes) ; \\ no metadata required
Mesh2_subpoly_lon_bnd:_Fillvalue = fillvalue ; \\ for missing subnodes
double Mesh2_subpoly_lat_bnd(nMesh2_subpoly,nMaxMesh2_subpoly_subnodes) ; \\ no metadata required
Mesh2_subpoly_lat_bnd:_Fillvalue = fillvalue ; \\ for missing subnodes

Topologische Daten

Berechnungsgitter

Knotenverzeichnis der Kanten

Knotenverzeichnis der Polygone

Polygonverzeichnis der Kanten

SubGrid

SubKantenverzeichnis der Kanten

integer Mesh2_edge_subedges(nMesh2_edge,nMaxMesh2_edge_subedges) ;
Mesh2_edge_subedges:standard_name = "???" ; \\ yet to be determined
Mesh2_edge_subedges:long_name = "list of sub edges for all edges" ;
Mesh2_edge_subedges:_FillValue = fillvalue ; \\ for missing sub edges
Bemerkung: Indizes starten mit Null!

SubPolygonverzeichnis der Polygone

integer Mesh2_poly_subpolys(nMesh2_poly,nMaxMesh2_poly_subpolys) ;
Mesh2_poly_subpolys:standard_name = "???" ; \\ yet to be determined
Mesh2_poly_subpolys:long_name = "list of sub polygons for all polygons" ;
Mesh2_poly_subpolys:_FillValue = fillvalue ; \\ for missing sub polygons
Bemerkung: Indizes starten mit Null!

Topologievariable

integer Mesh2 ;
Mesh2:standard_name = "???" ; \\ yet to be determined
Mesh2:long_name = "name of grid can be stored here" ;
Mesh2:dimensionality = 2 ; \\ flat, two-dimensional topology
Mesh2:locations = "node edge poly subedge subpoly" ;
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:poly_node_connectivity = "Mesh2_poly_nodes" ;
Mesh2:edge_poly_connectivity = "Mesh2_edge_polys" ;
Mesh2:edge_subedge_connectivity = "Mesh2_edge_subedges" ;
Mesh2:poly_subpoly_connectivity = "Mesh2_poly_subpolys" ;

Koordinatentransformation

Anmerkungen, Fragen

  • Datei nutzt Deltares-CF-Erweiterungen mit den orange markierten (BAW-) Erweiterungen in der Topologievariable.
  • Können auch für Kanten bounds angegeben werden?
  • Müssen die Indizes unbedingt mit Null starten?
  • Die Beschreibung der Koordinatentransformation ist für UTM- und Gauß-Krüger-Koordinaten geeignet.

zurück zu NetCDF


Strukturübersicht