import flare
import flare.filters
import synthobs
from synthobs.sed import models

import synthobs.morph.images
import synthobs.morph.PSF



cosmo = flare.default_cosmo()
z = 8.


test = synthobs.test_data() # --- read in some test data

# --- calculate V-band (550nm) optical depth for each star particle
A = 5.2
test.tauVs = (10**A) * test.MetSurfaceDensities



model = models.define_model('BPASSv2.2.1.binary/ModSalpeter_300') # DEFINE SED GRID -
model.dust = {'slope': -1.0} # define dust curve





Exemplo n.º 2
0
import sys
import os

import synthobs
from synthobs.sed import models

import flare

# --- initialise model with SPS model and IMF. Set verbose = True to see a list of available lines.

# m = models.EmissionLines('BPASSv2.2.1.binary/ModSalpeter_300', verbose = False)
m = models.EmissionLines('BPASSv2.2.1.binary/ModSalpeter_300', verbose=False)

# --- read in test data based on BLUETIDES

test = synthobs.test_data()

# UNITS:
# masses = star particle mass in M_sol
# ages = star particle age in Myr
# Z = star mass fraction in metals
# tauV = V-band (550nm) optical depth for each star particle

# --- calculate intrinsic quantities

o = m.get_line_luminosity('HI6563',
                          test.Masses,
                          test.Ages,
                          test.Metallicities,
                          verbose=True)  # intrinsic line luminosities