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