##of.flow_at_one_node(rg, z, study_node)
##of.update_at_one_point(rg, rainfall_duration=900, model_duration=900,rainfall_intensity=0.0000193333)
##of.update_at_one_point(rg, rainfall_duration=480, model_duration=480,rainfall_intensity=0.0000183333)
# Once run, we can plot the output
##of.plot_at_one_node()
# Now the flow_across_grid() method will be discussed. The commands needed to
# run this method are triple-commented ('###') for convenience. All flow_at_one_node()
# methods calls *should* be commented out to run this driver quicker.
# Again, this function reads in data using the ModelParameterDictionary and the default
# input file. The only required arguments for this method are the RasterModelGrid instance
# and the the initial elevations.
# This is the standard call to the flow_across_grid() method
of.flow_across_grid(rg,z)
of.update_across_grid(rg, rainfall_duration=900, model_duration=900,rainfall_intensity=0.0000193333)
of.update_across_grid(rg, rainfall_duration=480, model_duration=900,rainfall_intensity=0.0000183333)
plt.figure('Total Erosion, m')
imshow_grid(rg, of.total_dzdt, show=True)
endtime = time.time()
print(endtime - start_time, "seconds")
##of.update_at_one_point(rg, rainfall_duration=480, model_duration=480,rainfall_intensity=0.0000183333)
# Once run, we can plot the output
##of.plot_at_one_node()
# Now the flow_across_grid() method will be discussed. The commands needed to
# run this method are triple-commented ('###') for convenience. All flow_at_one_node()
# methods calls *should* be commented out to run this driver quicker.
# Again, this function reads in data using the ModelParameterDictionary and the default
# input file. The only required arguments for this method are the RasterModelGrid instance
# and the the initial elevations.
# This is the standard call to the flow_across_grid() method
of.flow_across_grid(rg, z)
of.update_across_grid(rg,
rainfall_duration=900,
model_duration=900,
rainfall_intensity=0.0000193333)
of.update_across_grid(rg,
rainfall_duration=480,
model_duration=900,
rainfall_intensity=0.0000183333)
plt.figure('Total Erosion, m')
imshow_grid(rg, of.total_dzdt, show=True)
endtime = time.time()
print endtime - start_time, "seconds"
# This is the standard way to call the flow_at_one_node method using the input file.
##of.flow_at_one_node(rg, z, study_node)
##of.update_at_one_point(rg, rainfall_duration=900, model_duration=3000,rainfall_intensity=0.0000093133)
# Once run, we can plot the output
##of.plot_at_one_node()
# Now the flow_across_grid() method will be discussed. The commands needed to
# run this method are triple-commented ('###') for convenience. All flow_at_one_node()
# methods calls *should* be commented out to run this driver quicker.
# Again, this function reads in data using the ModelParameterDictionary and the default
# input file. The only required arguments for this method are the RasterModelGrid instance
# and the the initial elevations.
# This is the standard call to the flow_across_grid() method
of.flow_across_grid(rg, z)
of.update_across_grid(rg,
rainfall_duration=900,
model_duration=3000,
rainfall_intensity=0.0000093133)
plt.figure('Total Erosion, m')
imshow_grid(rg, of.total_dzdt, show=True)
endtime = time.time()
print endtime - start_time, "seconds"
