Added transmittance computation

Pretty confined it works, although I parameterise and store it not like
Bruneton and Neyret