misspelled function name "thomsonlosfar"? in line:
PRO ThomsonFarY, Y, Elo, I, It, Itr     ; Elo in radians