isos = Iso.objects.filter(molecule__molecID=molec_id).all() print len(isos), 'isotopologues for molec', molec_id transitions = Trans.objects.filter(iso__in=isos).all() ntrans = len(transitions) last_pc = 0 for i, trans in enumerate(transitions): prms = trans.prm_set #for prm_name in ('Sw', 'A', 'gamma_air', 'gamma_self', 'n_air', # 'delta_air'): # exec('trans.%s = prms.get(name="%s")' % (prm_name, prm_name)) Ierr = trans.par_line[127:133] Iref = trans.par_line[133:145] this_trans = HITRANTransition() this_trans.par_line = trans.par_line this_trans.molec_id = molec_id this_trans.iso_id = trans.iso.isoID this_trans.nu = HITRANParam(val=prms.get(name="nu").val, ref=int(Iref[:2]), name='nu', ierr=int(Ierr[0])) this_trans.Sw = HITRANParam(val=prms.get(name="Sw").val, ref=int(Iref[2:4]), name='Sw', ierr=int(Ierr[1]), relative=True) this_trans.A = HITRANParam(val=prms.get(name="A").val, ref=int(Iref[2:4]),