function constant, expr
    c = 3d10         ; cm per second
    h = 6.626d-27    ; ergs seconds
    kB = 1.381d-16   ; ergs per kelvin
    pi = 3.14159265d
    void = execute('value='+expr)
    return, value
end