def run_swmm():
#=======================================
# setup all the nodes before starting
#=======================================
sim = Simulation('./pipe_test.inp')
#=======================================
# Start the simulation
#=======================================
print(sim._isStarted)
sim.start()
print(sim._isStarted)
node_names = ['Inlet', 'Outlet']
link_names = ['Culvert']
nodes = [Nodes(sim)[names] for names in node_names]
links = [Links(sim)[names] for names in link_names]
# type, area, length, orifice_coeff, free_weir_coeff, submerged_weir_coeff
opening0 = nodes[0].create_opening(4, 1.0, 1.0, 0.6, 1.6, 1.0)
opening0 = nodes[1].create_opening(4, 1.0, 1.0, 0.6, 1.6, 1.0)
print("\n")
print("n0 is coupled? ", nodes[0].is_coupled)
print("n1 is coupled? ", nodes[1].is_coupled)
nodes[0].overland_depth = 1.0
nodes[0].coupling_area = 1.0
# This step_advance should be an integer multiple of the routing step
# which is set in the ,inp file. Currently set to 1s.
# Should be able to interrogate sim to find out what the
# routing stepsize is. Maybe should issue a warning if
# step_advance is set lower than the routing step size.
# Indeed maybe step_advance should just allow advance n routing steps?
for i in range(50):
ids = dict()
volumes = sim.coupling_step(1.0)
print("Step:", i)
print("Current time:", sim.current_time)
#print("volumes:",volumes)
#print(volumes)
for flows in volumes:
#print(flows)
vols = flows[1]
for key in vols:
print("Volume total at node", key, "at time", flows[0], ":",
vols[key])
domain.set_name("anuga_swmm")
for t in domain.evolve(yieldstep=1.0, finaltime=60.0):
print("\n")
#print(f"coupling step: {t}")
domain.print_timestepping_statistics()
if t < stop_release_water_time:
# assume we need to release the water into the domain for first two seconds
pass
#op_inlet.set_rate(flow)
else:
# set the overland_depth
# TODO: set up the overland depth, modify this function
print("total volume: ",domain.get_water_volume())
volumes = sim.coupling_step(1.0)
print(volumes)
nodes[0].overland_depth = get_depth(op_inlet)
print("inlet overland depth: ", get_depth(op_inlet))
volumes_in_out = volumes[-1][-1]
print(volumes_in_out)
if t <= stop_release_water_time+1:
# no water exchange as the first two steps from swmm and anuga did not match.
print("Volume total at node Inlet" ":", volumes_in_out["Inlet"])
print("Oulet: ", nodes[1].total_inflow)
op_inlet.set_rate(0)
op_outlet.set_rate(0)
else:
