def zonak2004(): '''Zoank, S. et al. 2004, ApJ, 2004, 606, 196 PG1634+706 HST+Keck spectra MgII, SiIV, SiIII from Table 2. Summing Subsystems A (Model 2) and B Errors estimated by JXP (not reported) SiIII in A may be a model SiIV in B may be a model Total NHI from LL. Taken from Fig 3 caption. Error estimated by JXP Not all EWs in Table 1 included Adopting their M/H ''' # Setup radec = xor.stod1('J163428.9897+703132.422') # SIMBAD lls = LLSSystem(name='PG1634+706_z1.041', RA=radec[0], Dec=radec[1], zem=1.337, zabs=1.0414, vlim=[-200., 30.] * u.km / u.s, NHI=17.23, MH=-1.4, sigNHI=np.array([0.15, 0.15])) # SubSystems lls.mk_subsys(2) # Abundances adict = dict(MgII={ 'clm': log_sum([11.45, 11.90, 12.02, 11.68]), 'sig_clm': 0.05, 'flg_clm': 1 }, SiIII={ 'clm': log_sum([12.5, 12.5, 12.8, 12.7]), 'sig_clm': 0.25, 'flg_clm': 1 }, SiIV={ 'clm': log_sum([10.9, 10.8, 11.2, 11.1]), 'sig_clm': 0.15, 'flg_clm': 1 }) lls.subsys['A']._ionclms = IonClms(idict=adict) bdict = dict(SiIII={ 'clm': log_sum([11.8, 12.8, 12.4]), 'sig_clm': 0.15, 'flg_clm': 1 }, SiIV={ 'clm': log_sum([11.2, 12.2, 11.8]), 'sig_clm': 0.15, 'flg_clm': 1 }) lls.subsys['B']._ionclms = IonClms(idict=bdict) # Total lls._ionclms = lls.subsys['A']._ionclms.sum(lls.subsys['B']._ionclms) lls.Refs.append('Zon04') # Return return lls