hspf 2.0.1__py3-none-any.whl → 2.0.2__py3-none-any.whl

hspf/hbn.py

@@ -12,9 +12,8 @@ import math
 from struct import unpack
 from numpy import fromfile
 from pandas import DataFrame
-from datetime import datetime, timedelta, timezone
+from datetime import datetime, timedelta #, timezone
 from collections import defaultdict
-import numpy as np
 #from pathlib import Path
 
 
@@ -350,11 +349,13 @@ class hbnClass:
         rows = []
         times = []
         nvals = len(self.mapn[operation, id, activity]) # number constituent timeseries
-        utc_offset = timezone(timedelta(hours=-6)) #UTC is 6hours ahead of CST
+        #utc_offset = timezone(timedelta(hours=-6)) #UTC is 6hours ahead of CST
         for (index, reclen) in self.mapd[operation, id, activity, tcode]:
             yr, mo, dy, hr, mn = unpack('5I', self.data[index + 36: index + 56])
             hr = hr-1
-            dt = datetime(yr, mo, dy, 0, mn ,tzinfo=utc_offset) + timedelta(hours=hr)
+            #dt = datetime(yr, mo, dy, 0, mn ,tzinfo=utc_offset) + timedelta(hours=hr)
+            dt = datetime(yr, mo, dy, 0, mn ) + timedelta(hours=hr)
+
             times.append(dt)
 
             index += 56

hspf/parser/graph.py

@@ -332,6 +332,14 @@ CREATE TABLE GenInfo (
 
 #%% Methods using universal node id
 
+# def bypass_node(G, node):
+#     preds = list(G.predecessors(node))
+#     succs = list(G.successors(node))
+#     for u in preds:
+#         for v in succs:
+#             if 
+#             G.add_edge(u, v)
+#     G.remove_node(node)
 
 def _add_subgraph_labels(G,G_sub):
     G_sub.labels = {label:node for label, node in G.labels.items() if node in G_sub.nodes}
@@ -519,7 +527,24 @@ def paths(G,reach_id,source_type = 'RCHRES'):
 def count_ancestors(G,node_type,ancestor_node_type):
     return {node['type_id']:len(ancestors(G,node['id'],ancestor_node_type)) for node in get_nodes(G,node_type)}
 
-
+# def catchment_ids(G):
+#     result = []
+#     for node in get_node_ids(G,'RCHRES'):
+#         upstream_nodes = G.predecessors(node)
+#         if any([G.nodes[up]['type'] in ['PERLND','IMPLND'] for up in upstream_nodes]):
+#             result.append(G.nodes[node]['type_id'])
+#     return result
+
+# Very expensive. Should probably standardize it so routing reaches have no implnds/perlnds
+def catchment_ids(G):
+    result = []
+    for node in get_node_ids(G,'RCHRES'):
+        upstream_nodes = G.predecessors(node)
+        if any([G.nodes[up]['type'] in ['PERLND','IMPLND'] for up in upstream_nodes]):
+            cat = make_catchment(G,G.nodes[node]['type_id'])
+            if area(cat) > 0:
+                result.append(G.nodes[node]['type_id'])
+    return result
 
 # Catchment constructor
 def make_catchment(G,reach_id):
@@ -550,8 +575,28 @@ def make_watershed(G,reach_ids):
     return watershed
 
 
-def catcments(G):
-    return None
+# def catcments(G):
+#     cats = [Catchment(graph.make_catchment(G,reach_id) for reach_id in graph.get_node_type_ids(G,'RCHRES'))]
+
+#     return 
+        
+#         for u, v, edge_data in graph.make_catchment(G,reach_id).edges(data=True):
+#             source_node_attributes = G.nodes[u]
+#             # Add or update edge attributes with source node attributes
+#             edge_data["source_type"] = source_node_attributes.get("type")
+#             edge_data["source_name"] = source_node_attributes.get("name")
+#             edge_data["source_type_id"] = source_node_attributes.get("type_id")
+#             cats.append(edge_data)
+        
+#     return pd.DataFrame(cats)
+                     
+                     
+#                      for node in G.nodes:
+#         upstream_nodes = G.predecessors(node)
+#         if any(G.nodes[up]['type'] in ['PELND','IMPLND'] for up in upstream_nodes):
+#             result.append(node)
+
+#     return None
 # Catchment selectors
 
 '''
@@ -637,9 +682,12 @@ class Catchment():
 #%% Legacy Methods for Backwards compatability
 class reachNetwork():
     def __init__(self,uci,reach_id = None):
+        self.uci = uci
         self.G = create_graph(uci)
+        self.catchment_ids = catchment_ids(self.G)
+        self.routing_reaches = self._routing_reaches()
+        self.lakes = self._lakes()
         self.schematic = uci.table('SCHEMATIC').astype({'TVOLNO': int, "SVOLNO": int, 'AFACTR':float})
-        self.uci = uci
     
     def get_node_type_ids(self,node_type):
         return get_node_type_ids(self.G, node_type)
@@ -709,15 +757,18 @@ class reachNetwork():
     
     def subwatersheds(self,reach_ids = None):
         df = subwatersheds(self.uci)
-        if reach_ids is not None:
-            df = df.loc[df.index.intersection(reach_ids)]
-        return df
+        if reach_ids is None:
+            reach_ids = get_node_type_ids(self.G,'RCHRES')
+        return df.loc[df.index.intersection(reach_ids)]
     
     def subwatershed(self,reach_id):
         return subwatershed(self.uci,reach_id) #.loc[reach_id]
     
     def subwatershed_area(self,reach_id):
-        return self.drainage(reach_id).query("source_type in ['PERLND','IMPLND']")['area'].sum()
+        area = self.drainage(reach_id).query("source_type in ['PERLND','IMPLND']")['area'].sum()
+        # if (reach_id in self.lakes()) & (f'FTABLE{reach_id}' in self.uci.table_names('FTABLES')):
+        #     area = area + self.lake_area(reach_id)
+        return area
     
     def reach_contributions(self,operation,opnids):
         return reach_contributions(self.uci,operation,opnids)
@@ -736,6 +787,15 @@ class reachNetwork():
     def outlets(self):
         return [self.G.nodes[node]['type_id'] for node, out_degree in self.G.out_degree() if (out_degree == 0) & (self.G.nodes[node]['type'] == 'RCHRES')]
 
+    def _lakes(self):
+        return list(self.uci.table('RCHRES','GEN-INFO').query('LKFG == 1',engine = 'python').index.astype(int))        
+    
+    def lake_area(self,reach_id):
+        return self.uci.table('FTABLES',f'FTABLE{reach_id}')['Area'].max()
+    
+    def _routing_reaches(self):
+        return [reach_id for reach_id in self.get_node_type_ids('RCHRES') if reach_id not in self.catchment_ids]
+
     def paths(self,reach_id):
         return paths(self.G,reach_id)
     
@@ -799,6 +859,7 @@ def subwatersheds(uci):
     schematic = schematic[schematic['TVOL'] == 'RCHRES'][['SVOLNO','TVOLNO','AFACTR','MLNO']].astype({'SVOLNO':int,'TVOLNO':int,'AFACTR':float,'MLNO':int})
     schematic.reset_index(inplace=True,drop=False)
     schematic.set_index('TVOLNO',inplace=True)
+    schematic = schematic.loc[catchment_ids(uci.network.G)]
     
     dfs = []
     for operation in ['PERLND','IMPLND']:

hspf/uci.py

@@ -78,7 +78,7 @@ class UCI():
                  
     def table(self,block,table_name = 'na',table_id = 0,drop_comments = True):
         # Dynamic parsing of tables when called by user
-        assert block in ['FILES','PERLND','IMPLND','RCHRES','SCHEMATIC','OPN SEQUENCE','MASS-LINK','EXT SOURCES','NETWORK','GENER','MONTH-DATA','EXT TARGETS','COPY']
+        assert block in ['FILES','PERLND','IMPLND','RCHRES','SCHEMATIC','OPN SEQUENCE','MASS-LINK','EXT SOURCES','NETWORK','GENER','MONTH-DATA','EXT TARGETS','COPY','FTABLES']
         
         table = self.uci[(block,table_name,table_id)] #[block][table_name][table_id]
         #TODO move the format_opnids into the Table class?

{hspf-2.0.1.dist-info → hspf-2.0.2.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: hspf
-Version: 2.0.1
+Version: 2.0.2
 Summary: Python package for downloading and running HSPF models
 Project-URL: Homepage, https://github.com/mfratkin1/pyHSPF
 Author-email: Mulu Fratkin <michael.fratkin@state.mn.us>

{hspf-2.0.1.dist-info → hspf-2.0.2.dist-info}/RECORD

@@ -1,9 +1,9 @@
 hspf/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-hspf/hbn.py,sha256=Nsclgm7GFTdeBHgy754OMIY0F4rCh8heCdeG2aGiSaY,19421
+hspf/hbn.py,sha256=SQMxWllZy5OxWGMkhmjiardb8vbSjrmENJrorLBqTDI,19476
 hspf/helpers.py,sha256=djKc12ZZkJmB_cHEbFm-mk8sp4GAbBNfjXxfp7YAELU,3132
 hspf/hspfModel.py,sha256=8XFPd89niSn9bNTjB2UUpoLNAs6wsD6i6Lb9YKoYjUU,8090
 hspf/reports.py,sha256=DfS9DoNwrnD3UvxO879i-bM2gWh5QUMxrV4mdRDgpfE,51878
-hspf/uci.py,sha256=_0nMVOrEn9VVTqMlGNENea22M-T3IIuSwfHXWK-s3Do,30313
+hspf/uci.py,sha256=kb05yc_pmO7w8e8tDOxD7Ln0PxhYALmU8yeoQZHRw9g,30323
 hspf/wdm.py,sha256=q0hNqsMNrTkxHtKEX0q0wWlIZabXv6UX2HjNCF9WEW4,12734
 hspf/wdmReader.py,sha256=-akKWB9SpUzXvXoQMeHLZNi_u584KaeEOyHB-YolTWM,22848
 hspf/data/ParseTable.csv,sha256=ExqUaZg_uUPF5XHGLJEk5_jadnDenKjbwqC4d-iNX_M,193609
@@ -27,8 +27,8 @@ hspf/data/Timeseries Catalog/RCHRES/OXRX.txt,sha256=NWdRFpJ60LsYzCGHjt8Llay3OI8j
 hspf/data/Timeseries Catalog/RCHRES/PLANK.txt,sha256=0MAehIrF8leYQt0Po-9h6IiujzoWOlw-ADCV-bPiqs0,3508
 hspf/data/Timeseries Catalog/RCHRES/SEDTRN.txt,sha256=SiTgD4_YWctTgEfhoMymZfv8ay74xzCRdnI005dXjyE,659
 hspf/parser/__init__.py,sha256=2HvprGVCaJ9L-egvTj1MI-bekq5CNjtSBZfrCtQi3fs,92
-hspf/parser/graph.py,sha256=nBNEmciZ31DX8RzDxq1b_M_pcPNpywEywwF5thoRwBA,34294
+hspf/parser/graph.py,sha256=zjIhr-6y8Iy4jJM-V__aP7Q3zkRzr8PajuHjgqy3cBk,37085
 hspf/parser/parsers.py,sha256=Uk7RD_9aAX4TTjUOyPfwdURPdmE2M2_YItny2PedGi4,19717
-hspf-2.0.1.dist-info/METADATA,sha256=jKf_h1m-EDXasS3HATdTdXj3gujvmB80adyUwtGQU2o,605
-hspf-2.0.1.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
-hspf-2.0.1.dist-info/RECORD,,
+hspf-2.0.2.dist-info/METADATA,sha256=EOwcX8KhToElTZ44rqHbD7mD53d7Nx33ttyHcyZckQc,605
+hspf-2.0.2.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
+hspf-2.0.2.dist-info/RECORD,,