}
relative_strength_total_degeneracy = {
'10-01': 9.0,
'11-01': 9.0,
'12-01': 9.0,
}
line_names = freq_dict.keys()
ckms = units.speedoflight_ms / 1e3 #2.99792458e5
voff_lines_dict = dict([(k, (v - 88.6318470e9) / 88.6318470e9 * ckms)
for k, v in freq_dict.iteritems()])
hcn_vtau = hyperfine.hyperfinemodel(line_names, voff_lines_dict, freq_dict,
line_strength_dict,
relative_strength_total_degeneracy)
hcn_amp = hcn_vtau.ampfitter
hcn_vtau_fitter = hcn_vtau.fitter
hcn_vtau_vheight_fitter = hcn_vtau.vheight_fitter
hcn_varyhf_fitter = hcn_vtau.varyhf_fitter
hcn_varyhf_amp_fitter = hcn_vtau.varyhf_amp_fitter
hcn_varyhf_amp_width_fitter = hcn_vtau.varyhf_amp_width_fitter
def hcn_radex(xarr,
density=4,
column=13,
xoff_v=0.0,
width=1.0,
grid_vwidth=1.0,
'11-01':3.0,
'12-01':5.0,
}
relative_strength_total_degeneracy = {
'10-01':9.0,
'11-01':9.0,
'12-01':9.0,
}
line_names = freq_dict.keys()
ckms = units.speedoflight_ms / 1e3 #2.99792458e5
voff_lines_dict = dict([(k,(v-83.6318470e9)/83.6318470e9*ckms) for k,v in freq_dict.iteritems()])
hcn_vtau = hyperfine.hyperfinemodel(line_names, voff_lines_dict, freq_dict, line_strength_dict, relative_strength_total_degeneracy)
hcn_amp = hcn_vtau.ampfitter
hcn_vtau_fitter = hcn_vtau.fitter
hcn_vtau_vheight_fitter = hcn_vtau.vheight_fitter
def hcn_radex(xarr, density=4, column=13, xoff_v=0.0, width=1.0,
grid_vwidth=1.0,
grid_vwidth_scale=False,
texgrid=None,
taugrid=None,
hdr=None,
path_to_texgrid='',
path_to_taugrid='',
temperature_gridnumber=3,
debug=False,
verbose=False,