Beispiel #1
0
        from radis.levels.dunham import EvJ
        dunham_coeffs = get_dunham_coefficients('CO', 1, 'X1SIG+')
        
        # Now calculate energy
        EvJ(v=0, J=0, **dunham_coeffs)
    
    See Also
    --------
    
    Get Dunham coefficients for a molecule: :py:func:`~radis.db.utils.get_dunham_coefficients`

    Use Herzberg coefficients instead: :py:func:`~radis.levels.dunham.Gv`, 
    :py:func:`~radis.levels.dunham.Fv`
        
    '''

    E = 0

    for ykl_name, Ykl in Ykl_dict.items():
        k = int(ykl_name[1])
        l = int(ykl_name[2])
        E += Ykl * (v + 0.5)**k * (J * (J + 1))**l

    return E


if __name__ == '__main__':

    from radis.test.phys.test_dunham import _run_all_tests
    print('Testing Dunham.py: ', _run_all_tests())
Beispiel #2
0
        from radis.levels.dunham import EvJ
        dunham_coeffs = get_dunham_coefficients('CO', 1, 'X1SIG+')

        # Now calculate energy
        EvJ(v=0, J=0, **dunham_coeffs)

    See Also
    --------

    Get Dunham coefficients for a molecule: :py:func:`~radis.db.utils.get_dunham_coefficients`

    Use Herzberg coefficients instead: :py:func:`~radis.levels.dunham.Gv`,
    :py:func:`~radis.levels.dunham.Fv`
    """

    E = 0

    for ykl_name, Ykl in Ykl_dict.items():
        k = int(ykl_name[1])
        l = int(ykl_name[2])
        E += Ykl * (v + 0.5) ** k * (J * (J + 1)) ** l

    return E


if __name__ == "__main__":

    from radis.test.phys.test_dunham import _run_all_tests

    print("Testing Dunham.py: ", _run_all_tests())