[GitHub] [incubator-sdap-nexus] fgreg commented on a change in pull request #68: SDAP-192 Created DOMS netCDF reader tool.

[GitBox <gi...@apache.org>]

fgreg commented on a change in pull request #68: SDAP-192 Created DOMS netCDF reader tool.
URL: https://github.com/apache/incubator-sdap-nexus/pull/68#discussion_r269743469
 
 

 ##########
 File path: tools/doms/doms_reader.py
 ##########
 @@ -0,0 +1,133 @@
+# Licensed to the Apache Software Foundation (ASF) under one or more
+# contributor license agreements.  See the NOTICE file distributed with
+# this work for additional information regarding copyright ownership.
+# The ASF licenses this file to You under the Apache License, Version 2.0
+# (the "License"); you may not use this file except in compliance with
+# the License.  You may obtain a copy of the License at
+#
+#   http://www.apache.org/licenses/LICENSE-2.0
+#
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+
+import argparse
+from netCDF4 import Dataset
+import sys
+import datetime
+import csv
+from collections import OrderedDict
+
+def assemble_matches(filename):
+    """
+    Read a DOMS netCDF file and return a list of matches.
+    
+    Arguments: 
+        filename (string): the DOMS netCDF file name.
+    
+    Returns:
+        matches (list): List of matches. Each list element is a dictionary:
+            For netCDF group GROUP (SatelliteData or InsituData) and group
+                variable VARIABLE:
+            matches[GROUP]['matchID']: MatchedRecords dimension ID for the match
+            matches[GROUP]['GROUPID']: GROUP dim dimension ID for the record
+            matches[GROUP][VARIABLE]: variable value 
+    """
+    # Open the netCDF file.
+    try:
+        doms_nc = Dataset(filename, 'r')
+    except OSError as err:
+        sys.exit("Error reading netCDF file " + filename + ": " + str(err))
+    
+    # Check that the number of groups is consistent w/ MatchedGroups dimension
+    matched_groups = doms_nc.dimensions['MatchedGroups'].size
+    groups = []
+    for group in doms_nc.groups:
+        groups.append(group)
+    
+    assert len(groups) == matched_groups, \
+        ("Number of groups isn't the same as MatchedGroups dimension.")
+    
+    matches = []
+    matched_records = doms_nc.dimensions['MatchedRecords'].size
+    
+    # Check that the time units are 'seconds since 1970-01-01 00:00:00 0:00'.
+    # Necessary to create the datetime field using utcfromtimestamp function.
+    for group in groups:
+        time_units = "seconds since 1970-01-01 00:00:00 0:00"
+        assert doms_nc.groups[group]['time'].units == time_units, \
 
 Review comment:
   Don't need this if using `num2date` see comment below.

----------------------------------------------------------------
This is an automated message from the Apache Git Service.
To respond to the message, please log on to GitHub and use the
URL above to go to the specific comment.
 
For queries about this service, please contact Infrastructure at:
users@infra.apache.org


With regards,
Apache Git Services