(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':