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)