text: import openmoc
  import openmoc.plotter as plotter
  from openmoc.options import Options
  from lattices import universes
  
  ###############################################################################
  #                          Main Simulation Parameters
  ###############################################################################
  
  opts = Options()
  
  openmoc.log.set_log_level('NORMAL')
  
  openmoc.log.py_printf('TITLE', 'Simulating the LRA Benchmark Problem...')
  
  ###############################################################################
  ##########################     Creating Cmfd mesh    ##########################
  ###############################################################################
  
  cmfd = openmoc.Cmfd()
  cmfd.setSORRelaxationFactor(1.5)
  cmfd.setLatticeStructure(110, 110)
  cmfd.setCentroidUpdateOn(False)
  
  
  ###############################################################################
  ##########################   Creating the Geometry   ##########################
  ###############################################################################
  
  openmoc.log.py_printf('NORMAL', 'Creating geometry...')
  
  geometry = openmoc.Geometry()
  geometry.setCmfd(cmfd)
  geometry.setRootUniverse(universes['Root'])
  
  
  ###############################################################################
  #                          Creating the TrackGenerator
  ###############################################################################
  
  openmoc.log.py_printf('NORMAL', 'Initializing the track generator...')
  
  track_generator = openmoc.TrackGenerator(geometry, opts.num_azim,
                                           opts.track_spacing)
  track_generator.setNumThreads(opts.num_omp_threads)
  track_generator.generateTracks()
  
  
  ###############################################################################
  #                            Running a Simulation
  ###############################################################################
  
  solver = openmoc.CPUSolver(track_generator)
  solver.setConvergenceThreshold(opts.tolerance)
  solver.setNumThreads(opts.num_omp_threads)
  solver.computeEigenvalue(opts.max_iters)
  solver.printTimerReport()
  
  
  ###############################################################################
  #                             Generating Plots
  ###############################################################################
  
  openmoc.log.py_printf('NORMAL', 'Plotting data...')
  
  openmoc.plotter.plot_materials(geometry, gridsize=500)
  openmoc.plotter.plot_cells(geometry, gridsize=500)
  openmoc.plotter.plot_flat_source_regions(geometry, gridsize=500)
  openmoc.plotter.plot_spatial_fluxes(solver, energy_groups=[1,2])
  
  openmoc.log.py_printf('TITLE', 'Finished')
  