# Input Parameters: radius in cells nrad (<= 12), no. time steps, modulate input energy?; zred, dgpc, alpha, pn,pb, bavg, -PSD slope, E(e)/E(B), rsize, gmaxmn, gmrat, gmin, betaup(laminar), betaup(turb), zeta (shock angle), thlos, opang, dust temp., obs. dust luminosity in units of 1e45 erg/s, dust torus distance, torus x-section rad., zdist0, vdm (vol. of Mach disk/cell vol.), ratio of ordered helical to total mag. field, pitch angle of helical field
  8 5000    1   # nrad =radius in cells nrad (<= 8), tstop, modulate input energy?; 1=no maps (0=make maps)
0.069     # zred = redshift
0.307     # dgpc = distance in Gpc
0.65      # alpha = optically thin power-law slope of synchrotron spectrum
0.3       # pn = normalized dispersion of density about mean value
0.3       # pb = normalized dispersion of magnetic field about mean value
0.15      # bave = mean turbulent magnetic field upstream of shock, in Gauss
2.0       # -PSD slope if input energy is modulated according to a PSD
1.0e0     # E(electrons)/E(B)
0.003     # rsize = radius of cells in pc
15000.0   # gmaxmin = minimum value of maximum electron energy after plasma crosses shock
10.0      # gmrat = (maximum/minimum) values of maximum electron energy after plasma crosses shock
1000.0    # gmin = minimum electron energy after plasma crosses shock
0.986d0   # betaup = laminar flow velocity in units of c upstream of shock
0.577     # betat = turbulent vortex rotation velocity in units of c upstream of shock
8.0d0     # zeta = angle between recollimation shock and jet axis
3.0d0     # thlos = angle between jet axis and line of sight
1.9d0     # opang = jet opening half-angle
1200.0    # tdust = blackbody temperature of dust in the torus
0.1       # ldust = luminosity of dust IR emission in units of 1E45 erg/s
1.5       # dtdist = radius (distance from black hole) of central circle of dust torus in pc
0.3       # dtrad = cross-sectional radius of dust torus in pc
1.0       # zdist0 = distance from jet vertex where standing shock starts
0.1       # vdm (vol. of Mach disk divided by the volume of "normal" cells)
0.4       # bratio = ratio of ordered helical to total (helical+turbulent) magnetic field
90.0      # bpitch = pitch angle of helical magnetic field (90 = toroidal, 0 = longitudinal)