import sys, os, time
import numpy as np
import matplotlib.pyplot as pl
import fsps
from prospect.sources import CompositeSFH
sfhtype = {1: 'tau', 4: 'delaytau', 5: 'simha'}
compute_vega_mags = False
zcontinuous = 1
sps = fsps.StellarPopulation(compute_vega_mags=compute_vega_mags,
zcontinuous=zcontinuous)
mysps = CompositeSFH(sfh_type='tau',
interp_type='logarithmic',
flux_interp='linear',
compute_vega_mags=compute_vega_mags,
zcontinuous=zcontinuous)
mysps.configure()
sspages = np.insert(mysps.logage, 0, 0)
mysps.mint_log = -3
wlo = 1e3
whi = 1.2e4
def main():
tage = 1.4
# vary sf_trunc
sf_trunc = tage * np.linspace(0.90, 1.02, 9)
import sys, os, time
import numpy as np
import matplotlib.pyplot as pl
import fsps
from prospect.sources import CompositeSFH
from sedpy import observate
sfhtype = {1: 'tau', 4: 'delaytau', 5: 'simha'}
# build FSPS and Prospector sps objects
zcontinuous = 1
sps = fsps.StellarPopulation(zcontinuous=zcontinuous)
tres = np.round(len(sps.ssp_ages) / 94.)
mysps = CompositeSFH(sfh_type='tau',
interp_type='logarithmic',
mint_log=5.45,
flux_interp='linear',
zcontinuous=zcontinuous)
mysps.configure()
# Save the Prospector SSP time axis
sspages = np.insert(mysps.logage, 0, mysps.mint_log)
# Set up some parameters that cause trouble in FSPS
pname = 'tage' # the parameter that will vary
badsimha = {
'logtau': [1.34, 0.62],
'delt_trunc': [0.91, 0.98],
'sf_tanslope': [1.24, -1.54],
}
i, sfh = 1, 5
