Ali A., Harries T. J., Douglas T. A., 2018, MNRAS, 477, 5422 (doi:10.1093/mnras/sty1001)

ionrp (ionization + radiation pressure)

notes
- ???? indicates dump number (0076-0392)
- radiation initiated from 0076 
- dumps separated by 0.5e4 yr
- ### indicates source number (001-029, descending in initial mass)

----------------

grids????-????.zip (grouped in sets of 5)
- dump_????.grid (torus grid file)

dumpvtu????-????.zip (grouped in sets of 10)
- dump_????.vtu (vtk files for grid, e.g. density, temperature, etc.)

nbody.zip (all)
- nbody_????.vtk (vtk files for sources)

source.zip (all)
- source_????.dat (torus source file)

logs.zip
- parameters.dat (torus parameters file)
- log (torus print statements)

radial.zip
nb: columns are (1) r (2) rho (3) u (4) rhoe (5) p (6) phi_stars (7) phi_gas (8) kappa*Flux (9) radmom (10) radpressure (11) temperature (12) tdust
- radial+X????.dat (cell values from grid centre to boundary, advancing in +x direction)  
- radial+Y????.dat (+y direction)
- radial+Z????.dat (+z direction)
- radial+XYZ????.dat (+xyz direction)

columnz.zip
- columnz_????.fits (fits files showing column density in g/cm2)

sourceinfo.zip
- sourceinfo_###.fits (info for star number ###, e.g. mass, radius, luminosity, position)

synthetic observations (t = colatitude theta, p = azimuth phi in degrees)
- halpha.zip 
- sii6731.zip 
- oiii5007.zip
- freefree.zip 
- dustimages.zip - as well as observations, also includes...
--- ????_avgtdust_*.fits (mean dust temp weighted by mass, or no weighting)
--- ????_normvarianceTdust_*.fits (variance normalised by mean value, for different weights)


analysis.zip
- plotting scripts
   - starburst.py (star positions and column density at starburst)
   - opacity.py & albedo01.dat (opacity vs wavelength)
   - makeplot.py (bulk grid properties)
   - multidensityfits.py (column density snapshots)
   - multirgb.py (HST plot)
   - densitypdfSinglePlot.py (column density PDFs)
   - globaltempplot.py (T vs time)
   - globalneplot.py (ne vs time)
   - habingplot.py (G0 around star 1 vs time)
   - habingplot2.py (G0 around star 25 vs time)
   - luminosityhalpha.py (Halpha luminosity/extinction vs time)
   - nlyman.py (NLy vs time)
   - bbratio.py (450 and 850 um plots)
- tables of data
   - tableoftotals.dat (total of variables on grid, eg mass)
   - valuesOnGrid.dat (more totals)
   - tableofnlyman.dat (NLy from integrated stellar spectra)
   - estimatedLymanFlux.dat (NLy from recombination rate - RHS of eq. 30, assuming case A, case B, Rubin1968 coefficients)
   - avgtempGlobalRms*.dat (average/stddev of T, ne, for different weights)
   - averageHabingFlux.dat (G0 around star 1 in Habing units)
   - averageHabingFlux_allStars.dat (G0 around each star in Habing units)