Esempio n. 1
0
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)
Esempio n. 2
0
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
Esempio n. 3
0
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
ages = np.linspace(1.4, 1.7, 20)
tau = 10**badsimha['logtau'][i]
delt_trunc = badsimha['delt_trunc'][i]
sf_slope = np.tan(badsimha['sf_tanslope'][i])

# filters to project onto
Esempio n. 4
0
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