storm generator

gistfile1.txt

  for (interval_duration, rainfall_rate) in precip.yield_storm_interstorm_duration_intensity():
      if rainfall_rate != 0.:
          mg = fr.route_flow(grid=mg, runoff_rate=rainfall_rate)
          #print 'Area: ', numpy.max(mg.at_node['drainage_area'])
          mg,_,_ = sp.erode(mg, interval_duration, Q_if_used='water_discharges', K_if_used='K_values')
      #add uplift
      mg.at_node['topographic_elevation'][mg.core_nodes] += uplift*interval_duration
      this_trunc = precip.elapsed_time//out_interval
      if this_trunc != last_trunc: #a new loop round
          print 'made it to loop ', out_interval*this_trunc
          last_trunc=this_trunc