AstroBEAR subroutine in scrambler.f90

1. Call AMRInit
   1. read in global data
2. Call PhysicsInit
   1. read in physics data
3. Call LevelsInit
   1. Call SinkParticleInit
      1. initializes variables for sink module
4. Call ModuleProblemInit
   1. Sets refinement criteria
   2. Call ProblemModuleInit
5. Call PhysicsFinalizeInit
   1. Uses EOS and finalizes q array
6. Begin time-stepping

Amr Subroutine in Amr_Control.f90:

1. Call ParticlePreUpdate
   1. Call ClearParticleMomementa
   2. Call CalcMoments
   3. Call SynchronizeMoments
   4. Call CheckForNewParticles
      1. If find a new particle: 1. Call CreateNewParticle, 2. Call AddNewSinkParticle
   5. Call CollectNewParticles
      1. Collects particles across processors (using mpi_allgather)
      2. Call CreateParticle: initializes variables for the new particle
      3. Call AddSinkParticle: adds to global sink list
      4. Call CreatePointGravityObject and sets these parameters
   6. Call DoAccretions
   7. Call SynchronizeAccretions
   8. Call FinalizeAccretions
2. Call ParticlePostUpdate
   1. Call SynchronizeGasForces
   2. Call AdvanceParticles
   3. Call UpdateParticles (only called when on max level)

Question, where does particle steps fit in with hydro, elliptic, radiative?

Elliptic → Hydro → Particle → Radiative?