def main():
"""Run main function."""
# Runtime: 8 sec
# get file path
config_path = 'physics_settings/configb.sh'
config = updated_config(config_path)
bd_flux_rat = config['bd_flux_rat']
print("bd_flux_rat = {}".format(bd_flux_rat))
# get two scalar values: scaled_bulge_ratio and scaled_disk_ratio
frb, frd = np.genfromtxt(bd_flux_rat, unpack=True, dtype=float)
print("frb = {}, frd = {}".format(frb, frd))
# get bulge and disk weights
b, d = get_bulge_disk_weights(config_path)
# get psf path names for weighted psf files to be written.
psfb, psfd, psfm = config['psfb'], config['psfd'], config['psfm']
# create bulge and disk psf
create_psfbd(b, psfb)
create_psfbd(d, psfd)
# create monochromatic psf
create_psfmono(frb, frd, psfb, psfd, psfm)
def debug():
config = updated_config('physics_settings/configd.sh')
rescaled_outfile = config['rescaled_outfiled']
print("Debugging:")
print(config['catalog_file']) # jedisim_out/out0/scaled_disk/trial1_catalog.txt
print(config['dislist_file']) # jedisim_out/out0/scaled_disk/trial1_dislist.txt
run_process("jeditransform", ['./executables/jeditransform',
config['catalog_file'],
config['dislist_file']
])
def main():
"""Run main function."""
# Config dictionaries
configb = updated_config('physics_settings/configb.sh')
configd = updated_config('physics_settings/configd.sh')
configm = updated_config('physics_settings/configm.sh')
# print("configb['90_output_folder'] = {}".format(configb['90_output_folder']))
# Get convolved scaled files g_csb, g_csd , and g_csm
lsst_TDCR(configb, configb['psfb'], configb['90_rescaled_outfileb'])
lsst_TDCR(configd, configd['psfd'], configd['90_rescaled_outfiled'])
lsst_TDCR(configm, configm['psfm'], configm['90_rescaled_outfilem'])
# get final monochromatic image
# jedisim_out/out0/scaled_bulge_disk/trial1_lsst_mono.fits
lsst_monochromatic(configm)
# get final chromatic image
# jedisim_out/out0/scaled_bulge_disk/trial1_lsst.fits
lsst_chromatic(configb, configd, configm)
def main():
"""Run main function."""
# Runtime: 1 min 30 sec
config_path = 'physics_settings/configb.sh'
config = updated_config(config_path)
# Create output dirs
replace_outfolder('simdatabase/scaled_bulge_f8')
replace_outfolder('simdatabase/scaled_disk_f8')
replace_outfolder('simdatabase/scaled_bulge_disk_f8')
# Now write files into them.
scale_galaxies(NUM_GALS, config)
def main():
"""Run main function."""
# Get file name to write output
# ofile = 'physics_settings/bd_flux_rat_z1.5.txt'
config_path = 'physics_settings/configb.sh'
config = updated_config(config_path)
ofile = config['bd_flux_rat']
print("Creating: {}".format(ofile))
# Now run the main program
fr, frb, frd = sbd_flux_ratio(ofile)
# Print output values
print("fr = {:.4f} frb = {:.4f} frd = {:.8f}".format(fr, frb, frd))
def jedisim_3cats():
"""Run main function."""
# Create catalogs
config_paths = [
'physics_settings/config{}.sh'.format(i) for i in list('bdm')
]
for config_path in config_paths:
config = updated_config(config_path)
run_jedicatalog(config_path)
# Update 3 catalogs for rotated case.
rotate_outfolder_angle_catalog(config)
rotate_outfolder_convolvedlist(config)
rotate_outfolder_distortedlist(config)
def main():
# Usage: python a3_bd_weights.py -z 1.5
# Runtime: 1.6 sec
# Config parameters
config = updated_config('physics_settings/configb.sh')
z = float(config['fixed_redshift'])
# Variables
infileb = 'sed/ssp_pf_interpolated_z{}.csv'.format(z)
infiled = 'sed/exp9_pf_interpolated_z{}.csv'.format(z)
outfile = 'physics_settings/bd_weights_z{}.txt'.format(z)
# Run the main program
write_bd_weights(config, infileb, infiled, outfile)
def main():
"""Run main function."""
# Config dictionaries
configb = updated_config('physics_settings/configb.sh')
configd = updated_config('physics_settings/configd.sh')
configm = updated_config('physics_settings/configm.sh')
# NOTE: we should first create output dirs and catalogs before running
# this jedisimulate program.
# Get convolved scaled files g_csb, g_csd , and g_csm
lsst_TDCR(configb, configb['psfb'], configb['rescaled_outfileb']) # out 3
lsst_TDCR(configd, configd['psfd'], configd['rescaled_outfiled']) # out 4
lsst_TDCR(configm, configm['psfm'], configm['rescaled_outfilem']) # out 5
# get final monochromatic image
# jedisim_out/out0/scaled_bulge_disk/trial1_lsst_mono.fits # main out 2
lsst_monochromatic(configm)
# get final chromatic image
# jedisim_out/out0/scaled_bulge_disk/trial1_lsst.fits # main out 1
lsst_chromatic(configb,configd,configm)
def main():
# Usage: python a3_bd_weights.py -z 1.5
# Runtime: 1.6 sec
# Config parameters
config = updated_config('physics_settings/configb.sh')
z = float(config['fixed_redshift'])
# Variables
sed_b = 'sed/ssp_pf_interpolated_z{}.csv'.format(z)
sed_d = 'sed/exp9_pf_interpolated_z{}.csv'.format(z)
file_bulge = 'sed/ssp_pf_interpolated_z{}.csv'.format(z)
file_disk = 'sed/exp9_pf_interpolated_z{}.csv'.format(z)
# Run the main program
bf_df(sed_b, sed_d, config)
def d90_replace_outfolders_jedi(config_path):
"""Replace old output folders before running new simulation."""
config = updated_config(config_path)
outfolders = [
config["90_output_folder"], config['90_output_folder'] + 'convolved/'
]
for outfolder in outfolders:
replace_outfolder(outfolder)
# Make stamps and distorted folders
# e.g. jedisim_out/out90/stamp_0 to stamp_12
# e.g. jedisim_out/out90/distorted_0 to stamp_12
for x in range(0, int(math.ceil(float(config['num_galaxies']) / 1000))):
postage_path = config['90_output_folder'] + "stamp_" + str(x)
distorted_path = config['90_output_folder'] + "distorted_" + str(x)
if not os.path.exists(postage_path):
os.makedirs(postage_path)
if not os.path.exists(distorted_path):
os.makedirs(distorted_path)
def jedisim_outfolders():
config = updated_config(config_path)
jouts = {}
z = config['fixed_redshift']
# five galaxies, three psfs, and two texts.
keys = [
'lsst', 'mono', 'gcsb', 'gcsd', 'gcsm', 'psfb', 'psfd', 'psfm',
'catalog', 'dislist'
]
odirs = [
'jedisim_output/jout_z{}_2017'.format(z) +
time.strftime("_%b_%d_%H_%M/z{}/{}/".format(z, key)) for key in keys
]
for i, odir in enumerate(odirs):
# replace output dirs
replace_outfolder(odir)
replace_outfolder(odir[0:-1] + '90/')
# also create dictionary
jouts[keys[i]] = odirs[i]
jouts[keys[i] + '90'] = odirs[i][0:-1] + '90/'
return jouts
def main():
# Global Variables
config_path = 'physics_settings/configb.sh'
config = updated_config(config_path)
lambda1 = 1000
lambda2 = 12000
infileb = 'sed/ssp_pf.cat'
infiled = 'sed/exp9_pf.cat'
# Derived names
z = float(config['fixed_redshift'])
outfileb = 'sed/ssp_pf_interpolated_z{}.csv'.format(z)
outfiled = 'sed/exp9_pf_interpolated_z{}.csv'.format(z)
# Parameters
age = util.age(4.0) - util.age(z)
age = int(round(age)) # Nearest Gyr
age = age - 1 # column 0 is 1 Gyr
age = abs(age)
# Interpolate
interpolate_flux(config, infileb, outfileb, lambda1, lambda2)
interpolate_flux(config, infiled, outfiled, lambda1, lambda2)
