Exemple #1
0
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]),