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
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