# Input Parameters: radius in cells nrad (<= 8) tstop, modulate input energy?; no maps?; 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, ratio of ordered helical to total mag. field, pitch angle of helical field
  6 2      0 1   # nrad =radius in cells nrad (<= 8), tstop, modulate input energy?; 1=no maps (0=make maps)
0.030     # zred = redshift
0.13      # dgpc = distance in Gpc
0.4       # alpha = optically thin power-law slope of synchrotron spectrum
0.2       # pn = normalized dispersion of density about mean value
0.2       # pb = normalized dispersion of magnetic field about mean value
1.0       # bave = mean turbulent magnetic field upstream of shock, in Gauss
1.0       # -PSD slope if input energy is modulated according to a PSD
0.25      # E(electrons)/E(B)
0.0007    # rsize = radius of cells in pc
2.0e5     # gmaxmin = minimum value of maximum electron energy after plasma crosses shock
20.0      # gmrat = (maximum/minimum) values of maximum electron energy after plasma crosses shock
2000.0    # gmin = minimum electron energy after plasma crosses shock
0.99875d0 # betaup = laminar flow velocity in units of c upstream of shock
0.1       # betat = turbulent vortex rotation velocity in units of c upstream of shock
6.0d0     # zeta = angle between recollimation shock and jet axis
2.0d0     # thlos = angle between jet axis and line of sight
0.5d0     # opang = jet opening half-angle
1200.0    # tdust = blackbody temperature of dust in the torus
0.001     # 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.5       # bratio = ratio of ordered helical to total (helical+turbulent) magnetic field
90.0      # bpitch = pitch angle of helical magnetic field (90 = toroidal, 0 = longitudinal)