• UpdateRadiation() | Radiation
• ir(ncells3, ncells2, ncells1, n_ang*nfreq): intensity
• ir_avg(nfreq, ncells3, ncells2, ncells1): average intensity over all angles | —radiation={const,six_ray} | radiation_task_list.hpp
• RadiationIntegratorTaskList | pmb→prad pmb→prad→pradintegrator | ir_avg is used as an parameter for chemistry (heating and chemical reaction rates) ir variable coordinates ordering will have problems with AMR functions (see rad_newtonian implementation) | | rad_newtonian | radiation/ radiation/integrators/ | RadIntegrator
• CalculateFluxes(AthenaArray<Real> &ir, const int order) | Radiation
• ir(ncells3, ncells2, ncells1, nangnoctnfreq): specific intensity
• ir_gray(ncells3, ncells2, ncells1, nang): integrated over all frequencies | -radiation -implicit_radiation | **im_rad_task_list.hpp
• IMRadTaskList** | pmb→prad pmb→prad→pradintegrator | | | gr_rad | radiation/ | Radiation
• CalculateFluxes(AthenaArray<Real> &prim_rad, const AthenaArray<Real> &prim_hydro, int order, Real dt) | Radiation
• prim(nang, ncells3, ncells2, ncells1): primitive intensity
• cons(nang, ncells3, ncells2, ncells1): conserved intensity | -r | time_integrator.cpp
• if (RADIATION_ENABLED) | pmb→prad | | | ART (Athena4.2) | radps/ | radps_transfer(&Mesh) | Erad[ncells3][ncells2][ncells1][nfrqe]: intensity | -enable-radps | N/A | N/A | Will be moved to Athena++. Need to be coupled with chemistry Merge six-ray with ART? | | MonteCarlo? | | | | | | | Shane Davis Coordinate with ART? |

Things need to be coordinated:

• File structure
• Radiation Variable Class
  • Class naming
  • Variable naming
  • Pointer
• Radiation Integrator Class
  • Class naming
  • Pointer
• Tasklist
  • Filenames
  • Class naming
  • Pointers (in main.cpp)
• Configure options
• Macros
• Logistics and order of merging to the master branch
• Wiki documentation

Option 1: treat them as completely different modules - The PLAN

• File structure
  • chem_rad/, newtonian_rad/, gr_rad/
  • Also for input files and regression tests, the folder naming should follow above
• Radiation Variable Class
  • Class naming: ChemRadiation, NewtonianRadiation, GRRadiation
  • Variable naming: up to the developer
  • Pointer: pmb→pchem_rad, pmb→pnewtonian_rad, pmb→pgr_rad
• Radiation Integrator Class
  • Class naming: ChemRadIntegrator, NewtonianRadIntegrator
  • Pointer: pmb→pchem_rad→pradintegrator, pmb→pnewtonian_rad→pradintegrator
• Tasklist
  • Filenames: chem_rad_task_list.hpp/cpp, newtonian_rad_task_llist.cpp
    • Can gr_rad be in a separate task list or does it have to be in time_integrator.cpp?
  • Class naming: ChemRadTaskList, NewtonianRadTaskList
  • Pointers: pchem_rad_list, pnewtonian_rad_list
• Configure options
  • —chem_radiation=[choice]
  • -newtonian_radiation, -newtonian_radiation_implicit
  • -gr_radiation
• Macros
  • CHEM_RAD_ENABLED, NEWTONIAN_RAD_ENABLED, GR_RAD_ENABLED
• Logistics and order of merging to the master branch
  • Since they are separate modules, this doesn’t matter so much
• Wiki documentation
  • “Radiation” page with subpages
    • Radiation for Chemistry
    • Newtonian Radiation
    • GR Radiation