New UGRID netcdf format v0.9 - Revision history

Yjzhang at 20:48, 9 February 2017

2017-02-09T20:48:15Z

Eli at 21:23, 30 June 2015

2015-06-30T21:23:33Z

Eli at 21:15, 30 June 2015

2015-06-30T21:15:18Z

Yjzhang at 19:08, 30 June 2015

2015-06-30T19:08:29Z

Eli at 16:16, 30 June 2015

2015-06-30T16:16:29Z

Eli at 16:11, 30 June 2015

2015-06-30T16:11:23Z

Eli at 23:01, 29 June 2015

2015-06-29T23:01:32Z

Show changes

Yjzhang at 17:37, 18 December 2014

2014-12-18T17:37:53Z

Yjzhang at 15:32, 1 December 2013

2013-12-01T15:32:31Z

Yjzhang at 15:32, 1 December 2013

2013-12-01T15:32:11Z

@@ Line 1: / Line 1: @@
-An IOOS working group led by Dr. Rich Signell have been working on a standard for unstructured grids. The latest is v0.9 as found
+An IOOS working group led by Dr. Rich Signell have been working on a standard for unstructured grids. The latest is v1.0 as found
-[https://github.com/ugrid-conventions/ugrid-conventions/blob/v0.9.0/ugrid-conventions.md here].
+[https://github.com/ugrid-conventions/ugrid-conventions here].
 A few SCHISM user groups have either built their own tools that use this new format or are working on this. Below is a summary of progress so far (as well as links if available):

@@ Line 63: / Line 63: @@
      float SCHISM_hgrid_node_x(nSCHISM_hgrid_node) ;
          SCHISM_hgrid_node_x:long_name = "x_coordinates" ;
-	SCHISM_hgrid_node_x:units = "meters" ;
+	SCHISM_hgrid_node_x:units = "meters" ;       # alternatively, degrees lat/lon
 	SCHISM_hgrid_node_x:mesh = "SCHISM_hgrid" ;
 	SCHISM_hgrid_node_x:location = "node" ;
@@ Line 69: / Line 69: @@
      float SCHISM_hgrid_node_y(nSCHISM_hgrid_node) ;
 	SCHISM_hgrid_node_y:long_name = "y_coordinates" ;
-	SCHISM_hgrid_node_y:units = "meters" ;
+	SCHISM_hgrid_node_y:units = "meters" ;       # alternatively, degrees lat/lon
 	SCHISM_hgrid_node_y:mesh = "SCHISM_hgrid" ;
 	SCHISM_hgrid_node_y:location = "node";

@@ Line 29: / Line 29: @@
 	float time(time) ;
 	time:long_name = "Time" ;
 	time:units = "seconds since 2013-28-08  00:00:00 PST" ;
 	time:base_date = "08/28/2013 00:00:00 PST                         " ;
@@ Line 136: / Line 137: @@
      salt:long_name ="water practical salinity";
      salt:missing_value ="-9999999.0";
-     salt:mesh = "Mesh3D";
+     salt:mesh = "SCHISM_hgrid";
      salt:location ="node";
@@ Line 142: / Line 143: @@
      temp:long_name ="water temperature";
      temp:missing_value ="-9999999.0";
-     temp:mesh = "Mesh3D";
+     temp:mesh = "SCHISM_hgrid";
      temp:location ="node";
@@ Line 152: / Line 153: @@
      u:missing_value = "-999999." ;
      u:missing_value = "-999999." ;
-     u:mesh = "Mesh3D" ;
+     u:mesh = "SCHISM_hgrid" ;
      u:location = "node" ;
@@ Line 158: / Line 159: @@
      v:long_name = "Northward Water Velocity" ;
      v:missing_value = "-999999." ;
-     v:mesh = "Mesh3D" ;
+     v:mesh = "SCHISM_hgrid" ;
      v:location = "node" ;
@@ Line 178: / Line 179: @@
      u:long_name = "Eastward Water Velocity" ;
      u:missing_value = "-999999." ;
-     u:mesh = "Mesh3D" ;
+     u:mesh = "SCHISM_hgrid" ;
      u:location = "edge" ;
@@ Line 184: / Line 185: @@
      v:long_name = "Northward Water Velocity" ;
      v:missing_value = "-999999." ;
-     v:mesh = "Mesh3D" ;
+     v:mesh = "SCHISM_hgrid" ;
      v:location = "edge" ;
@@ Line 196: / Line 197: @@
     w:long_name="vertical velocity at face center at all levels";
     w:missing_value ="-999999.0";
-    w:mesh ="Mesh3D";
+    w:mesh ="SCHISM_hgrid";
     w:location ="face";
@@ Line 205: / Line 206: @@
     salt:long_name ="water practical salinity";
     salt:missing_value ="-9999999.0";
-    salt:mesh = "Mesh3D";
+    salt:mesh = "SCHISM_hgrid";
     salt:location ="volume";
@@ Line 211: / Line 212: @@
     temp:long_name ="water temperature";
     temp:missing_value ="-9999999.0";
-    temp:mesh = "Mesh3D";
+    temp:mesh = "SCHISM_hgrid";
     temp:location ="volume";

@@ Line 111: / Line 111: @@
      float elev(time,nSCHISM_hgrid_node);
-     elev:long_name ="water surface height agove reference datum";
+     elev:long_name ="water surface height above reference datum";
      elev:missing_value = "-999999.0";
      elev:mesh ="SCHISM_hgrid";
@@ Line 171: / Line 171: @@
      dtbe:missing_value ="-999999.";
      dtbe:mesh ="SCHISM_hgrid";
-     dtbe:location="node";
+     dtbe:location="node"; <--edge?
   // type 67 vector at 3d mesh edge center

@@ Line 19: / Line 19: @@
 netcdf example_compressed {
 dimensions:
-	nSCHISM_hgrid_node          = 144773 ;
+    nSCHISM_hgrid_node          = 144773 ;
-	nSCHISM_hgrid_face          = 238116 ;
+    nSCHISM_hgrid_face          = 238116 ;
      nSCHISM_hgrid_edge          = 400000 ;
-	nMaxSCHISM_hgrid_face_nodes = 4 ;
+    nMaxSCHISM_hgrid_face_nodes = 4 ;
-	SCHISM_vgrid_layers         = 51;
+    SCHISM_vgrid_layers         = 51;
-	two                  = 2 ;
+    two                  = 2 ;
-	time = UNLIMITED ;
+    time = UNLIMITED ;
 variables:
 	float time(time) ;
-		time:long_name = "Time" ;
+	time:long_name = "Time" ;
-		time:units = "seconds since 2013-28-08  00:00:00 PST" ;
+	time:units = "seconds since 2013-28-08  00:00:00 PST" ;
-		time:base_date = "08/28/2013 00:00:00 PST                         " ;
+	time:base_date = "08/28/2013 00:00:00 PST                         " ;
-		time:standard_name = "time" ;
+	time:standard_name = "time" ;
      int SCHISM_hgrid ;
-		SCHISM_hgrid:cf_role = "mesh_topology" ;
+	SCHISM_hgrid:cf_role = "mesh_topology" ;
-		SCHISM_hgrid:long_name = "Topology data of 2d unstructured mesh" ;
+	SCHISM_hgrid:long_name = "Topology data of 2d unstructured mesh" ;
-		SCHISM_hgrid:topology_dimension = 2 ;
+	SCHISM_hgrid:topology_dimension = 2 ;
-		SCHISM_hgrid:node_coordinates = "SCHISM_hgrid_node_x SCHISM_hgrid_node_y" ;
+	SCHISM_hgrid:node_coordinates = "SCHISM_hgrid_node_x SCHISM_hgrid_node_y" ;
-		SCHISM_hgrid:face_node_connectivity = "SCHISM_hgrid_face_nodes" ;
+	SCHISM_hgrid:face_node_connectivity = "SCHISM_hgrid_face_nodes" ;
          SCHISM_hgrid:edge_coordinates  = "SCHISM_hgrid_edge_x SCHISM_hgrid_edge_y" ;
          SCHISM_hgrid:face_coordinates  = "SCHISM_hgrid_face_x SCHISM_hgrid_face_y" ;
@@ Line 45: / Line 45: @@
      int Mesh3D;
          Mesh3D:cf_role = "mesh_topology" ;
-		Mesh3D:long_name = "Topology data of 3d unstructured mesh" ;
+	Mesh3D:long_name = "Topology data of 3d unstructured mesh" ;
-		Mesh3D:topology_dimension = 3 ;
+	Mesh3D:topology_dimension = 3 ;
      int SCHISM_hgrid_face_nodes(nSCHISM_hgrid_face, nMaxSCHISM_hgrid_face_nodes) ;
@@ Line 73: / Line 73: @@
      float SCHISM_hgrid_face_x(nSCHISM_hgrid_face) ;
-		SCHISM_hgrid_face_x:long_name = "x_coordinates of 2D mesh face" ;
+	SCHISM_hgrid_face_x:long_name = "x_coordinates of 2D mesh face" ;
-		SCHISM_hgrid_face_x:units = "meters" ;
+	SCHISM_hgrid_face_x:units = "meters" ;
-		SCHISM_hgrid_face_x:mesh = "SCHISM_hgrid" ;
+	SCHISM_hgrid_face_x:mesh = "SCHISM_hgrid" ;
-		SCHISM_hgrid_face_x:location = "face" ;
+	SCHISM_hgrid_face_x:location = "face" ;
      float SCHISM_hgrid_face_y(nSCHISM_hgrid_face) ;
-		SCHISM_hgrid_face_y:long_name = "y_coordinates of 2D mesh face" ;
+	SCHISM_hgrid_face_y:long_name = "y_coordinates of 2D mesh face" ;
-		SCHISM_hgrid_face_y:units = "meters" ;
+	SCHISM_hgrid_face_y:units = "meters" ;
-		SCHISM_hgrid_face_y:mesh = "SCHISM_hgrid" ;
+	SCHISM_hgrid_face_y:mesh = "SCHISM_hgrid" ;
-		SCHISM_hgrid_face_y:location = "face" ;
+	SCHISM_hgrid_face_y:location = "face" ;
      float SCHISM_hgrid_edge_x(nSCHISM_hgrid_edge) ;
-		SCHISM_hgrid_edge_x:long_name = "x_coordinates of 2D mesh edge" ;
+	SCHISM_hgrid_edge_x:long_name = "x_coordinates of 2D mesh edge" ;
-		SCHISM_hgrid_edge_x:units = "meters" ;
+	SCHISM_hgrid_edge_x:units = "meters" ;
-		SCHISM_hgrid_edge_x:mesh = "SCHISM_hgrid" ;
+	SCHISM_hgrid_edge_x:mesh = "SCHISM_hgrid" ;
-		SCHISM_hgrid_edge_x:location = "edge" ;
+	SCHISM_hgrid_edge_x:location = "edge" ;
      float SCHISM_hgrid_edge_y(nSCHISM_hgrid_edge) ;
-		SCHISM_hgrid_edge_y:long_name = "y_coordinates of 2D mesh edge" ;
+	SCHISM_hgrid_edge_y:long_name = "y_coordinates of 2D mesh edge" ;
-		SCHISM_hgrid_edge_y:units = "meters" ;
+	SCHISM_hgrid_edge_y:units = "meters" ;
-		SCHISM_hgrid_edge_y:mesh = "SCHISM_hgrid" ;
+	SCHISM_hgrid_edge_y:mesh = "SCHISM_hgrid" ;
-		SCHISM_hgrid_edge_y:location = "edge" ;
+	SCHISM_hgrid_edge_y:location = "edge" ;
      float depth(nSCHISM_hgrid_node) ;
-		depth:long_name = "Bathymetry" ;
+	depth:long_name = "Bathymetry" ;
-		depth:units = "meters" ;
+	depth:units = "meters" ;
-		depth:positive = "down" ;
+	depth:positive = "down" ;
-		depth:mesh = "SCHISM_hgrid" ;
+	depth:mesh = "SCHISM_hgrid" ;
-		depth:location = "node" ;
+	depth:location = "node" ;
@@ Line 162: / Line 162: @@
-  // type 65 missing
+  // type 65 omit
@@ Line 188: / Line 188: @@
-  //type 68 missing
+  //type 68 omit
@@ Line 219: / Line 219: @@
 // global attributes:
-		:Conventions = "CF-1.0" ;
+    :Conventions = "CF-1.0" ;
 data:

@@ Line 5: / Line 5: @@
 <UL>
-   <LI> CAL-DWR (Eli Ateljevich et al): have built VisIT plug-ins for the native binary and v0.9 format;
+   <LI> CAL-DWR (Eli Ateljevich et al): have built VisIT plug-ins for the native binary and our evolving understanding of the v0.9 format;
    <LI> MetOcean, NZ (Brett Beamsley): have made attempt to output into this format. Wish to build 1-way nesting util based on it;
    <LI> Nat. Univ Singapore (Haihua Xu): have built a [http://ccrm.vims.edu/yinglong/wiki_files/Xu-Netcdf-VisIT/ package] (including VisIT plug-in) based on some older UGRID netcdf standard, for ELCIRC and SELFE;
@@ Line 12: / Line 12: @@
-Below is a proposal for how to present SCHISM data in UGRID format. We are still working on names for the grids and variables, so it isn't too late right now to say you hate what we have. The biggest problem at the moment is a seeming lack of guidance as to how to refer to data that is on a layered 2D grid such as temperature in the example below. Also, though we purport to be CF compliant ... we don't really seem to be headed that way for variable names.
+Below is a proposal for how to present SCHISM data in UGRID format. It is meant as a straw man for discussion, but we (California DWR) will act on it because we need make our VisIt visualization plugin a continuous, reasonably well-maintained offering and so far the only stable version has been based on native binary outputs.
   <nowiki>
@@ Line 46: / Line 48: @@
 		Mesh3D:topology_dimension = 3 ;
+    int SCHISM_hgrid_face_nodes(nSCHISM_hgrid_face, nMaxSCHISM_hgrid_face_nodes) ;
-	int SCHISM_hgrid_face_nodes(nSCHISM_hgrid_face, nMaxSCHISM_hgrid_face_nodes) ;
+        SCHISM_hgrid_face_nodes:long_name = "Horizontal Triangular Element Table" ;
-	    SCHISM_hgrid_face_nodes:long_name = "Horizontal Triangular Element Table" ;
+        SCHISM_hgrid_face_nodes:units = "non-dimensional" ;
-		SCHISM_hgrid_face_nodes:units = "non-dimensional" ;
+	SCHISM_hgrid_face_nodes:cf_role = "face_node_connectivity" ;
-		SCHISM_hgrid_face_nodes:cf_role = "face_node_connectivity" ;
+	SCHISM_hgrid_face_nodes:start_index = 1 ;
          SCHISM_hgrid_face_nodes:_FillValue = -999999;
@@ Line 59: / Line 60: @@
          SCHISM_hgrid_edge_nodes:start_index =1;
-	float SCHISM_hgrid_node_x(nSCHISM_hgrid_node) ;
+    float SCHISM_hgrid_node_x(nSCHISM_hgrid_node) ;
-		SCHISM_hgrid_node_x:long_name = "x_coordinates" ;
+        SCHISM_hgrid_node_x:long_name = "x_coordinates" ;
-		SCHISM_hgrid_node_x:units = "meters" ;
+	SCHISM_hgrid_node_x:units = "meters" ;
-		SCHISM_hgrid_node_x:mesh = "SCHISM_hgrid" ;
+	SCHISM_hgrid_node_x:mesh = "SCHISM_hgrid" ;
-		SCHISM_hgrid_node_x:location = "node" ;
+	SCHISM_hgrid_node_x:location = "node" ;
-	float SCHISM_hgrid_node_y(nSCHISM_hgrid_node) ;
+    float SCHISM_hgrid_node_y(nSCHISM_hgrid_node) ;
-		SCHISM_hgrid_node_y:long_name = "y_coordinates" ;
+	SCHISM_hgrid_node_y:long_name = "y_coordinates" ;
-		SCHISM_hgrid_node_y:units = "meters" ;
+	SCHISM_hgrid_node_y:units = "meters" ;
-		SCHISM_hgrid_node_y:mesh = "mesh" ;
+	SCHISM_hgrid_node_y:mesh = "SCHISM_hgrid" ;
-		SCHISM_hgrid_node_y:location = "node";
+	SCHISM_hgrid_node_y:location = "node";
      float SCHISM_hgrid_face_x(nSCHISM_hgrid_face) ;
@@ Line 97: / Line 98: @@
-	float depth(nSCHISM_hgrid_node) ;
+    float depth(nSCHISM_hgrid_node) ;
 		depth:long_name = "Bathymetry" ;
 		depth:units = "meters" ;
@@ Line 151: / Line 152: @@
      u:missing_value = "-999999." ;
      u:missing_value = "-999999." ;
-     u:mesh = "mesh3D" ;
+     u:mesh = "Mesh3D" ;
      u:location = "node" ;
@@ Line 177: / Line 178: @@
      u:long_name = "Eastward Water Velocity" ;
      u:missing_value = "-999999." ;
-     u:mesh = "mesh3D" ;
+     u:mesh = "Mesh3D" ;
      u:location = "edge" ;

@@ Line 8: / Line 8: @@
    <LI> MetOcean, NZ (Brett Beamsley): have made attempt to output into this format. Wish to build 1-way nesting util based on it;
    <LI> Nat. Univ Singapore (Haihua Xu): have built a [http://ccrm.vims.edu/yinglong/wiki_files/Xu-Netcdf-VisIT/ package] (including VisIT plug-in) based on some older UGRID netcdf standard, for ELCIRC and SELFE;
-   <LI> Joseph Zhang (OHSU & VIMS): the combine scripts can optionally output into netcdf format based on an older UGRID standard proposed >8 years ago;
+   <LI> OHSU & VIMS: the combine scripts can optionally output into netcdf format based on an older UGRID standard proposed >8 years ago;
 </UL>