@@ -82,6 +82,14 @@ def layered(self):
82
82
return res
83
83
84
84
def open (self , nc_fpth : str ) -> None :
85
+ """
86
+ Open an existing dataset. Assumes the dataset has been annotated
87
+ with the necessary attributes to read and update, including global
88
+ attributes modflow_model and modflow_grid.
89
+
90
+ Args:
91
+ nc_fpth (str): Path to an existing NetCDF file.
92
+ """
85
93
fpth = Path (nc_fpth ).resolve ()
86
94
self ._fname = fpth .name
87
95
@@ -95,6 +103,16 @@ def open(self, nc_fpth: str) -> None:
95
103
def create (
96
104
self , modeltype : str , modelname : str , nc_type : str , fname : str , modelgrid : Grid
97
105
) -> None :
106
+ """
107
+ Create a new dataset.
108
+
109
+ Args:
110
+ modeltype (str): A model type, e.g. GWF6.
111
+ modelname (str): The model name.
112
+ nc_type (str): A supported NetCDF file type: mesh2d or structured.
113
+ fname (str): The generated NetCDF file name.
114
+ modelgrid (Grid): A FloPy derived discretization object.
115
+ """
98
116
self ._modelname = modelname .lower ()
99
117
self ._modeltype = modeltype .lower ()
100
118
self ._nc_type = nc_type .lower ()
@@ -971,10 +989,10 @@ def open_dataset(nc_fpth: str, grid_type: str) -> ModelNetCDFDataset:
971
989
ModelNetCDFDataset: A dataset derived from the base class.
972
990
"""
973
991
nc_dataset = None
974
- dis_str = grid_type .lower ()
992
+ grid_t = grid_type .lower ()
975
993
976
994
# grid_type corresponds to a flopy.discretization type
977
- if dis_str != "vertex" and dis_str != "structured" :
995
+ if grid_t != "vertex" and grid_t != "structured" :
978
996
raise Exception (
979
997
"Supported NetCDF discretication types "
980
998
'are "vertex" (DISV) and "structured" '
@@ -990,10 +1008,10 @@ def open_dataset(nc_fpth: str, grid_type: str) -> ModelNetCDFDataset:
990
1008
else :
991
1009
modelname = dataset .attrs ["modflow_model" ].split (":" )[0 ].lower ()
992
1010
gridtype = dataset .attrs ["modflow_grid" ].lower ()
993
- if dis_str == "vertex" :
1011
+ if grid_t == "vertex" :
994
1012
if gridtype == "layered mesh" :
995
1013
nc_dataset = DisvNetCDFMesh2d ()
996
- elif dis_str == "structured" :
1014
+ elif grid_t == "structured" :
997
1015
if gridtype == "layered mesh" :
998
1016
nc_dataset = DisNetCDFMesh2d ()
999
1017
elif gridtype == "structured" :
@@ -1007,7 +1025,7 @@ def open_dataset(nc_fpth: str, grid_type: str) -> ModelNetCDFDataset:
1007
1025
raise Exception (
1008
1026
f"Unable to load netcdf dataset for file grid "
1009
1027
f'type "{ gridtype } " and discretization grid '
1010
- f'type "{ dis_str } "'
1028
+ f'type "{ grid_t } "'
1011
1029
)
1012
1030
1013
1031
return nc_dataset
0 commit comments