old_contents: #Import statements so that you will have access to the necessary methods
  import numpy
  from landlab import RasterModelGrid
  from landlab.plot.imshow import imshow_node_grid
  from pylab import show, figure
  
  #Create a raster grid with 25 rows, 40 columns, and cell spacing of 10 m
  mg = RasterModelGrid(25, 40, 10.0)
  
  #Create a field of node data (an array) on the grid called elevation.  
  #Initially populate this array with zero values.
  z = mg.add_zeros('node', 'elevation')
  
  #Check the size of the array
  len(z)
  
  #Create a diagonal fault across the grid
  fault_y = 50.0 + 0.25*mg.node_x
  upthrown_nodes = numpy.where(mg.node_y>fault_y)
  z[upthrown_nodes] += 10.0 + 0.01*mg.node_x[upthrown_nodes]
  
  #Illustrate the grid
  imshow_node_grid(mg, z, cmap='jet', grid_units=['m','m'])
  show()
  
  #Define paramaters
  kd = 0.01   # 0.01 m2 per year
  dt = 0.2*mg.dx*mg.dx/kd   # CFL condition
  
  #Set boundary conditions
  mg.set_closed_boundaries_at_grid_edges(False, True, False, True)
  
  #Get id values of the cord nodes on which you will operate
  interior_nodes = mg.get_core_nodes()
  
  #Evolve landscape
  for i in range(25):
   	g = mg.calculate_gradients_at_active_links(z)
   	qs = -kd*g
   	dqsdx = mg.calculate_flux_divergence_at_nodes(qs)
   	dzdt = -dqsdx
   	z[interior_nodes] += dzdt[interior_nodes]*dt
  
  #Plot new landscape
  figure()
  imshow_node_grid(mg, z, cmap='jet', grid_units=['m','m'])
  show()
  
new_contents: #Import statements so that you will have access to the necessary functions
  import numpy
  from landlab import RasterModelGrid
  from landlab.plot.imshow import imshow_node_grid
  from pylab import show, figure
  
  #Create a raster grid with 25 rows, 40 columns, and cell spacing of 10 m
  mg = RasterModelGrid(25, 40, 10.0)
  
  #Create a field of node data (an array) on the grid called elevation.  
  #Initially populate this array with zero values.
  z = mg.add_zeros('node', 'elevation')
  
  #Check the size of the array
  len(z)
  
  #Create a diagonal fault across the grid
  fault_y = 50.0 + 0.25*mg.node_x
  upthrown_nodes = numpy.where(mg.node_y>fault_y)
  z[upthrown_nodes] += 10.0 + 0.01*mg.node_x[upthrown_nodes]
  
  #Illustrate the grid
  imshow_node_grid(mg, z, cmap='jet', grid_units=['m','m'])
  show()
  
  #Define paramaters
  kd = 0.01   # 0.01 m2 per year
  dt = 0.2*mg.dx*mg.dx/kd   # CFL condition
  
  #Set boundary conditions
  mg.set_closed_boundaries_at_grid_edges(False, True, False, True)
  
  #Get id values of the cord nodes on which you will operate
  interior_nodes = mg.get_core_nodes()
  
  #Evolve landscape
  for i in range(25):
   	g = mg.calculate_gradients_at_active_links(z)
   	qs = -kd*g
   	dqsdx = mg.calculate_flux_divergence_at_nodes(qs)
   	dzdt = -dqsdx
   	z[interior_nodes] += dzdt[interior_nodes]*dt
  
  #Plot new landscape
  figure()
  imshow_node_grid(mg, z, cmap='jet', grid_units=['m','m'])
  show()
  
subject: Update scarp_diffusion_no_component.py