last = current[:]
current = []
for c in hydroseqs:
if flowlines[hydroseqs[c]].down in last:
comids.append(c)
current.append(hydroseqs[c])
print(('found {} flowlines including:\n'.format(len(comids))))
for c in comids:
print(c)
print('')
# make an instance of the FtableCalculator to use for the data from the file
calculator = FtableCalculator(gageid)
# calculate log-log regressions of flow and width vs depth
calculator.calculate_regressions()
# get the length and average flow from the flowline VAAs
length = flowlines[hydroseqs[comid]].length
qref = flowlines[hydroseqs[comid]].flow * 0.3048**3
# calculate the FTABLE for the gage station using the length of the reach
ftable = calculator.create_ftable(length, units=units)
print(('FTABLE for gage station {}:\n'.format(comid)))
# find the record indices of the comid and reach length in km in the file
i = [f[0] for f in reader.fields].index('LENGTHKM') - 1
j = [f[0] for f in reader.fields].index('COMID') - 1
# get the reach length and common identifier
length = record[i]
comid = record[j]
it = comid, length
print('comid {} is closest to the gage and has a length of {} km\n'.format(*it))
# make an instance of the FtableCalculator to use for the data from the file
calculator = FtableCalculator(gageid)
# calculate log-log regressions of flow and width vs depth using the observed
# values from the NWIS database
calculator.calculate_regressions()
# calculate the FTABLE using the length of the reach and the regressions for:
#
# the width (surface area = width * length),
# the volume (volume = width * depth * length)
# the discharge
#
# this approach utilizes all available data to infer discharge rate for a given
# depth of flow needed for continuous simulation
i = [f[0] for f in reader.fields].index('LENGTHKM') - 1
j = [f[0] for f in reader.fields].index('COMID') - 1
# get the reach length and common identifier
length = record[i]
comid = record[j]
it = comid, length
print(
'comid {} is closest to the gage and has a length of {} km\n'.format(*it))
# make an instance of the FtableCalculator to use for the data from the file
calculator = FtableCalculator(gageidpath)
# calculate log-log regressions of flow and width vs depth using the observed
# values from the NWIS database
calculator.calculate_regressions()
# calculate the FTABLE using the length of the reach and the regressions for:
#
# the width (surface area = width * length),
# the volume (volume = width * depth * length)
# the discharge
#
# this approach utilizes all available data to infer discharge rate for a given
# depth of flow needed for continuous simulation
while len(current) > 0:
last = current[:]
current = []
for c in hydroseqs:
if flowlines[hydroseqs[c]].down in last:
comids.append(c)
current.append(hydroseqs[c])
print('found {} flowlines including:\n'.format(len(comids)))
for c in comids: print(c)
print('')
# make an instance of the FtableCalculator to use for the data from the file
calculator = FtableCalculator(gageid)
# calculate log-log regressions of flow and width vs depth
calculator.calculate_regressions()
# get the length and average flow from the flowline VAAs
length = flowlines[hydroseqs[comid]].length
qref = flowlines[hydroseqs[comid]].flow * 0.3048**3
# calculate the FTABLE for the gage station using the length of the reach
ftable = calculator.create_ftable(length, units = units)
print('FTABLE for gage station {}:\n'.format(comid))
