################################################################################
##########
########## Population model setup
##########
################################################################################
logger.setLevel(-1)
args = parse_args()
logger.info("Loading base %s", path.basename(args.baseFile))
t1 = time.clock()
base = StarlightBase(args.baseFile, args.baseDir)
l_ssp = np.arange(3650.0, 6850.0, 2.0)
f_ssp = base.f_sspResam(l_ssp)
logger.info("Took %.2f seconds to read the base (%d files)" % (time.clock() - t1, base.sspfile.size))
wl_norm_window = (l_ssp < 5680.0) & (l_ssp > 5590.0)
################################################################################
##########
########## Morphology model setup
##########
################################################################################
logger.info("Creating original B-D model.")
norm_model = get_model(args.trueModel, with_default=True)
norm_params = np.array(norm_model.getParams(), dtype=norm_model.dtype)
logger.info("Original model at normalization window:\n%s\n" % str(norm_model))
logger.info('Creating PSF (FWHM = %.2f ", beta = %.2f)' % (args.psfFWHM, args.psfBeta))
gs.tight_layout(fig, rect=[0, 0, 1, 0.95])
plt.savefig('plots/tese/simulation_initmodel.pdf')
################################################################################
##########
########## Plot original population model
##########
################################################################################
logger.debug('Plotting original population model.')
#index_norm = find_nearest_index(l_ssp, 5635.0)
logger.info('Loading base %s', path.basename(args.baseFile))
t1 = time.clock()
base = StarlightBase(args.baseFile, args.baseDir)
wl = np.arange(3650.0, 6850.0, 2.0)
f_ssp = base.f_sspResam(wl)
logger.info('Took %.2f seconds to read the base (%d files)' % (time.clock() - t1, base.sspfile.size))
wl_norm_window = (wl < 5680.0) & (wl > 5590.0)
t0 = base.ageBase.max()
r = np.arange(30)
tau = tau_r(args.tau0, args.dtau_dr, r)
print tau
bulge_sfh_tau0 = get_synth_sfh(t0, args.tau0, base.ageBase)
bulge_sfh_tau1 = get_synth_sfh(t0, 4e9, base.ageBase)
bulge_sfh_tau2 = get_synth_sfh(t0, 8e9, base.ageBase)
disk_sfh_tauInf = get_synth_sfh(t0, 1e12, base.ageBase)
bulge_flux_tau0 = (f_ssp * bulge_sfh_tau0[:, np.newaxis]).sum(axis=1).sum(axis=0)
