def plot_tardis(): radial1d_mdl = model_radial_oned.Radial1DModel(config_object, atom_data) radial1d_mdl.create_packets() out_nu, out_energy, jsestimator, nubarestimator = montecarlo_multizone.montecarlo_radial1d(radial1d_mdl) spectrum = synspec.get_lambda_spec(out_nu, out_energy, 500*1e-8, 20000*1e-8, samples=1000) ax.plot(np.log10(spectrum.wave/1e4), np.log10(spectrum.flux), label=config_object.plasma_type) ax.set_xlim(-1,0) ax.set_xlabel(r'$\log10{\lambda} [\log10{\mu m}]$') ax.set_ylabel(r'$\log10{Flux}$') return out_nu, out_energy, jsestimator, nubarestimator
def plot_tardis(): radial1d_mdl = model_radial_oned.Radial1DModel(config_object, atom_data) radial1d_mdl.create_packets() out_nu, out_energy, jsestimator, nubarestimator, line_id_in, line_id_out = montecarlo_multizone.montecarlo_radial1d( radial1d_mdl) spectrum = synspec.get_lambda_spec(out_nu, out_energy, 500 * 1e-8, 20000 * 1e-8, samples=1000) return out_nu, out_energy, jsestimator, nubarestimator, line_id_in, line_id_out, radial1d_mdl
def plot_tardis(): radial1d_mdl = model_radial_oned.Radial1DModel(config_object, atom_data) radial1d_mdl.create_packets() out_nu, out_energy, jsestimator, nubarestimator = montecarlo_multizone.montecarlo_radial1d( radial1d_mdl) spectrum = synspec.get_lambda_spec(out_nu, out_energy, 500 * 1e-8, 20000 * 1e-8, samples=1000) ax.plot(np.log10(spectrum.wave / 1e4), np.log10(spectrum.flux), label=config_object.plasma_type) ax.set_xlim(-1, 0) ax.set_xlabel(r'$\log10{\lambda} [\log10{\mu m}]$') ax.set_ylabel(r'$\log10{Flux}$') return out_nu, out_energy, jsestimator, nubarestimator
def plot_tardis(): radial1d_mdl = model_radial_oned.Radial1DModel(config_object, atom_data) radial1d_mdl.create_packets() out_nu, out_energy, jsestimator, nubarestimator, line_id_in, line_id_out = montecarlo_multizone.montecarlo_radial1d(radial1d_mdl) spectrum = synspec.get_lambda_spec(out_nu, out_energy, 500*1e-8, 20000*1e-8, samples=1000) return out_nu, out_energy, jsestimator, nubarestimator, line_id_in, line_id_out, radial1d_mdl