netcdf util uses data type specific interface

Author: mjreno

flopy/mf6/data/mfdataarray.py

@@ -1602,10 +1602,10 @@ def _set_storage_netcdf(self, nc_dataset, create=False):
         input_tag = f"{m}/{p}/{n}".lower()
 
         if create:
-            # cache data and update file before resetting storage
+            # add array to netcdf dataset
             nc_varname = f"{p}_{n}".lower()
             data = self._get_data()
-            nc_dataset.create_var(nc_varname, input_tag, 0, data, self.structure.shape)
+            nc_dataset.create_array(nc_varname, input_tag, 0, data, self.structure.shape)
 
         storage._set_storage_netcdf(
             nc_dataset, input_tag, self.structure.layered, storage.layer_storage.get_total_size() - 1

flopy/mf6/data/mffileaccess.py

@@ -473,7 +473,7 @@ def load_netcdf_array(
         nc_tag,
         layer,
     ):
-        return nc_dataset.data(nc_tag, layer + 1)
+        return nc_dataset.array(nc_tag, layer + 1)
 
     def get_data_string(self, data, data_type, data_indent=""):
         layer_data_string = [str(data_indent)]

flopy/mf6/data/mfstructure.py

@@ -1425,7 +1425,7 @@ def __init__(self, data_item, model_data, package_type, dfn_list):
             or "nodes" in data_item.shape
             or len(data_item.layer_dims) > 1
         )
-        # TODO: only gwf, gwt, gwe models
+        # TODO: only gwf, gwt, gwe models NETCDF-DEV
         self.netcdf = (
             ("ncol" in data_item.shape
             or "nrow" in data_item.shape

flopy/mf6/mfmodel.py

@@ -937,19 +937,26 @@ def load_base(
                 dis_str = {
                     "dis6": "structured",
                     "disv6": "vertex",
-                    "disu6": "unstructured",
+                    #"disu6": "unstructured",
                 }
                 dis_type = None
                 for t in instance.name_file.packages.get_data():
                     if t[0].lower().startswith("dis"):
                         dis_type = t[0].lower()
                         break
-                if dis_type:
+                if dis_type and dis_type in dis_str:
                     nc_fpth = os.path.join(instance.model_ws, nc_filerecord[0][0])
                     instance._nc_dataset = open_netcdf_dataset(nc_fpth, dis_type=dis_str[dis_type])
                 else:
-                    pass
-                    # TODO: set error "invalid dis for netcdf input"
+                    message = (
+                        "Invalid discretization type "
+                        f"provided for model {modelname} "
+                        "NetCDF input"
+                    )
+                    raise MFDataException(
+                        model=modelname,
+                        message=message,
+                    )
 
         # order packages
         vnum = mfstructure.MFStructure().get_version_string()