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
