(filterIR, 0.02, 0, 0),
(filterUV, 0, 0, 4)],
postfix="_augmented",
fmin=fmin,
fmax=fmax)
move_visualization_files(skirtrun, filenames)
# -----------------------------------------------------------------
# build SED plots for each SKIRT-run
if vistype == 'seds':
from pts.core.plot.seds import plotseds
print "Building SED plots for {} SKIRT-runs".format(len(skirtruns))
for skirtrun in skirtruns:
print "Building SED plot for SKIRT-run {}...".format(skirtrun.runid())
plotseds(skirtrun.simulation())
move_visualization_files(skirtrun, filenames)
# -----------------------------------------------------------------
# build temperature histograms
if vistype == 'temperature':
from pts.eagle.plottemperature import plottemperature
print "Building temperature histograms for {} SKIRT-runs".format(
len(skirtruns))
for skirtrun in skirtruns:
print "Building temperature histogram for SKIRT-run {}...".format(
skirtrun.runid())
plottemperature(skirtrun)
# -----------------------------------------------------------------
# Import standard modules
import sys
# Import the relevant PTS classes and modules
from pts.core.simulation.simulation import createsimulations
from pts.core.plot.seds import plotseds
# -----------------------------------------------------------------
# a value of None means that the axis is autoscaled;
# alternatively specify a range through a tuple with min and max values
xlim = None
ylim = None
#xlim = ( 5e-2, 1e3 )
#ylim = ( 1e-13, 1e-9 )
# -----------------------------------------------------------------
print "Starting plotseds..."
# get the command-line argument specifying the simulation(s)
argument = sys.argv[1] if len(sys.argv) > 1 else ""
# construct the list of simulation objects and make the plots
for simulation in createsimulations(argument):
plotseds(simulation, xlim=xlim, ylim=ylim)
print "Finished plotseds"
# -----------------------------------------------------------------
# convert to W / [m2 * Hz]
fnus = []
for i in range(len(flambda)):
fnu = flambda[i] / spectral_factor_hz_to_micron(
lambdas[i] * u("micron")) * u("W / (m2 * Hz)")
fnu = fnu.to("Jy").value
fnus.append(fnu)
# Create the SED
sed = SED(photometry_unit="Jy")
for i in range(len(lambdas)):
sed.add_point(lambdas[i] * u("micron"), fnus[i] * u("Jy"))
plotter = SEDPlotter()
plotter.add_sed(sed, "...", residuals=False)
#min_wavelength = None
#max_wavelength = None
#min_flux = None
#max_flux = None
#plotter.run(min_wavelength=min_wavelength, max_wavelength=max_wavelength, min_flux=min_flux, max_flux=max_flux)
plotseds(simulation)
# -----------------------------------------------------------------
print "Building gas particle plots for {} SKIRT-runs".format(
len(skirtruns))
for skirtrun in skirtruns:
print "Building gas particle plot for SKIRT-run {}...".format(
skirtrun.runid())
plotgasparticles(skirtrun)
# -----------------------------------------------------------------
# build SED plots for each SKIRT-run
if vistype == 'seds':
from pts.core.plot.seds import plotseds
print "Building SED plots for {} SKIRT-runs".format(len(skirtruns))
for skirtrun in skirtruns:
print "Building SED plot for SKIRT-run {}...".format(skirtrun.runid())
plotseds(skirtrun.simulation(), output_path=skirtrun.vispath())
# -----------------------------------------------------------------
# build temperature histograms
if vistype == 'temperature':
from pts.eagle.plottemperature import plottemperature
print "Building temperature histograms for {} SKIRT-runs".format(
len(skirtruns))
for skirtrun in skirtruns:
print "Building temperature histogram for SKIRT-run {}...".format(
skirtrun.runid())
plottemperature(skirtrun)
# -----------------------------------------------------------------
fmin,fmax = makeintegratedrgbimages(skirtrun.simulation(),
[ (filterR, 1,0,0), (filterG, 0,1,0), (filterB, 0,0,1) ], postfix="_optical")
makeintegratedrgbimages(skirtrun.simulation(),
[ (filterR, 1,0,0), (filterG, 0,1,0), (filterB, 0,0,1), (filterIR, 0.02,0,0), (filterUV, 0,0,4) ],
postfix="_augmented", fmin=fmin, fmax=fmax)
move_visualization_files(skirtrun, filenames)
# -----------------------------------------------------------------
# build SED plots for each SKIRT-run
if vistype=='seds':
from pts.core.plot.seds import plotseds
print "Building SED plots for {} SKIRT-runs".format(len(skirtruns))
for skirtrun in skirtruns:
print "Building SED plot for SKIRT-run {}...".format(skirtrun.runid())
plotseds(skirtrun.simulation())
move_visualization_files(skirtrun, filenames)
# -----------------------------------------------------------------
# build temperature histograms
if vistype=='temperature':
from pts.eagle.plottemperature import plottemperature
print "Building temperature histograms for {} SKIRT-runs".format(len(skirtruns))
for skirtrun in skirtruns:
print "Building temperature histogram for SKIRT-run {}...".format(skirtrun.runid())
plottemperature(skirtrun)
# -----------------------------------------------------------------
# build wavelength movies for each SKIRT-run
if vistype=='particles':
from pts.eagle.plotgasparticles import plotgasparticles
print "Building gas particle plots for {} SKIRT-runs".format(len(skirtruns))
for skirtrun in skirtruns:
print "Building gas particle plot for SKIRT-run {}...".format(skirtrun.runid())
plotgasparticles(skirtrun)
# -----------------------------------------------------------------
# build SED plots for each SKIRT-run
if vistype=='seds':
from pts.core.plot.seds import plotseds
print "Building SED plots for {} SKIRT-runs".format(len(skirtruns))
for skirtrun in skirtruns:
print "Building SED plot for SKIRT-run {}...".format(skirtrun.runid())
plotseds(skirtrun.simulation(), output_path=skirtrun.vispath())
# -----------------------------------------------------------------
# build temperature histograms
if vistype=='temperature':
from pts.eagle.plottemperature import plottemperature
print "Building temperature histograms for {} SKIRT-runs".format(len(skirtruns))
for skirtrun in skirtruns:
print "Building temperature histogram for SKIRT-run {}...".format(skirtrun.runid())
plottemperature(skirtrun)
# -----------------------------------------------------------------
# build wavelength movies for each SKIRT-run
if vistype=='wavemovie':
