def main():
# INITIALIZE
# Name of parameter input file
input_file_name = 'test_inputs_for_diffusion_model.txt'
# Open input file and read run-control parameters
mpd = ModelParameterDictionary(input_file_name)
run_duration = mpd.get('RUN_DURATION', ptype=float)
# Create and initialize a grid
mg = create_and_initialize_grid(mpd)
# Create and initialize a diffusion component
dc = LinearDiffuser(mg)
dc.initialize(mpd)
# RUN
# Run the diffusion component until it's time for the next output
dc.run_until(run_duration)
# FINALIZE
# Display results to screen
mg.imshow('node', 'landscape_surface__elevation')
import pylab
pylab.show()
from landlab.io.netcdf import write_netcdf
write_netcdf('diffusion_example.nc', mg)
def main():
# INITIALIZE
# Name of parameter input file
input_file_name = 'test_inputs_for_diffusion_model.txt'
# Open input file and read run-control parameters
mpd = ModelParameterDictionary(input_file_name)
run_duration = mpd.get('RUN_DURATION', ptype=float)
# Create and initialize a grid
mg = create_and_initialize_grid(mpd)
# Create and initialize a diffusion component
dc = LinearDiffuser(mg)
dc.initialize(mpd)
# RUN
# Run the diffusion component until it's time for the next output
dc.run_until(run_duration)
# FINALIZE
# Display results to screen
mg.imshow('node', 'landscape_surface__elevation')
import pylab
pylab.show()
from landlab.io.netcdf import write_netcdf
write_netcdf('diffusion_example.nc', mg)
