def SO_s0(nside, pixel_indices=None, mpi_comm=None, nside_template=512): T_map = get_data_from_url("synch_T_ns{}.fits".format(nside_template)) Q_map = get_data_from_url("synch_Q_ns{}.fits".format(nside_template)) U_map = get_data_from_url("synch_U_ns{}.fits".format(nside_template)) A_I = read_map(T_map, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm) return [{ "model": "power_law", "nu_0_I": 0.408, "nu_0_P": 23., "A_I": A_I, "A_Q": read_map(Q_map, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), "A_U": read_map(U_map, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), "spectral_index": np.ones(len(A_I)) * -3.1, }]
def SO_d1(nside, pixel_indices=None, mpi_comm=None, nside_template=512): T_map = get_data_from_url("dust_T_ns{}.fits".format(nside_template)) Q_map = get_data_from_url("dust_Q_ns{}.fits".format(nside_template)) U_map = get_data_from_url("dust_U_ns{}.fits".format(nside_template)) beta_map = get_data_from_url( "variable_spectral_index/beta_dust_ns{}_1deg.fits".format( nside_template)) temperature_map = get_data_from_url( "variable_spectral_index/temperature_dust_ns{}_1deg.fits".format( nside_template)) A_I = read_map(T_map, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm) return [{ "model": "modified_black_body", "nu_0_I": 545., "nu_0_P": 353., "A_I": A_I, "A_Q": read_map(Q_map, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), "A_U": read_map(U_map, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), "spectral_index": read_map(beta_map, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), "temp": read_map(temperature_map, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), "add_decorrelation": False, }]
def SO_a0(nside, pixel_indices=None, mpi_comm=None, nside_template=512): T_map1 = get_data_from_url("ame_T_ns{}.fits".format(nside_template)) T_map2 = get_data_from_url("ame2_T_ns{}.fits".format(nside_template)) return [ { "model": "spdust", "nu_0_I": 22.8, "nu_0_P": 22.8, "A_I": read_map(T_map1, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), "nu_peak_0": 30., "emissivity": loadtxt(get_data_from_url("ame_emissivity.txt"), mpi_comm=mpi_comm, unpack=True), "nu_peak": 18.95, }, { "model": "spdust", "nu_0_I": 41.0, "nu_0_P": 41.0, "A_I": read_map(T_map2, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), "nu_peak_0": 30., "emissivity": loadtxt(get_data_from_url("ame_emissivity.txt"), mpi_comm=mpi_comm, unpack=True), "nu_peak": 33.35, }, ]
def SO_s1(nside, pixel_indices=None, mpi_comm=None, nside_template=512): T_map = get_data_from_url("synch_T_ns{}.fits".format(nside_template)) Q_map = get_data_from_url("synch_Q_ns{}.fits".format(nside_template)) U_map = get_data_from_url("synch_U_ns{}.fits".format(nside_template)) beta_map = get_data_from_url( "variable_spectral_index/beta_synch_ns{}_1deg.fits".format( nside_template)) A_I = read_map(T_map, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm) return [{ "model": "power_law", "nu_0_I": 0.408, "nu_0_P": 23., "A_I": A_I, "A_Q": read_map(Q_map, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), "A_U": read_map(U_map, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), "spectral_index": read_map(beta_map, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), }]
def SO_d0(nside, pixel_indices=None, mpi_comm=None, nside_template=512): T_map = get_data_from_url("dust_T_ns{}.fits".format(nside_template)) Q_map = get_data_from_url("dust_Q_ns{}.fits".format(nside_template)) U_map = get_data_from_url("dust_U_ns{}.fits".format(nside_template)) A_I = read_map(T_map, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm) return [{ "model": "modified_black_body", "nu_0_I": 545., "nu_0_P": 353., "A_I": A_I, "A_Q": read_map(Q_map, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), "A_U": read_map(U_map, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), "spectral_index": np.ones(len(A_I)) * 1.53, "temp": np.ones(len(A_I)) * 19.6, "add_decorrelation": False, }]
def SO_f1(nside, pixel_indices=None, mpi_comm=None, nside_template=512): T_map = get_data_from_url("freefree_T_ns{}.fits".format(nside_template)) return [{ "model": "power_law", "nu_0_I": 30., "A_I": read_map(T_map, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), "spectral_index": -2.14, }]
def SO_a1(nside, pixel_indices=None, mpi_comm=None, nside_template=512): T_map1 = get_data_from_url("ame1_T_ns{}.fits".format(nside_template)) T_map2 = get_data_from_url("ame2_T_ns{}.fits".format(nside_template)) Q_dust = get_data_from_url("dust_Q_ns{}.fits".format(nside_template)) U_dust = get_data_from_url("dust_U_ns{}.fits".format(nside_template)) nu_peak_map = get_data_from_url( "variable_spectral_index/ame_nu0_peak_ns{}_1deg.fits".format( nside_template)) return [ { "model": "spdust_pol", "nu_0_I": 22.8, "nu_0_P": 22.8, "A_I": read_map(T_map1, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), "nu_peak_0": 30., "emissivity": loadtxt(get_data_from_url("ame_emissivity.txt"), mpi_comm=mpi_comm, unpack=True), 'nu_peak': read_map(nu_peak_map, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), 'pol_frac': 0.01, 'angle_q': read_map(Q_dust, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), 'angle_u': read_map(U_dust, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm) }, { "model": "spdust_pol", "nu_0_I": 41.0, "nu_0_P": 41.0, "A_I": read_map(T_map2, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), "nu_peak_0": 30., "emissivity": loadtxt(get_data_from_url("ame_emissivity.txt"), mpi_comm=mpi_comm, unpack=True), "nu_peak": 33.35, 'pol_frac': 0.01, 'angle_q': read_map(Q_dust, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm), 'angle_u': read_map(U_dust, nside, field=0, pixel_indices=pixel_indices, mpi_comm=mpi_comm) }, ]