New routine for molecular calculation from first principles.
Now the read_sonde_* should only read the sounding data and make not calculation.
Inverted order of beta_mol. Now it is (levels, wavelength).
Inverted alpha_mol which is now (bins, wlen). In the raman analysis,
stopped doing the substituition ray_ext=alpha_mol and simply uses
alpha_mol everywhere.
Raman code was using zet_0 as the first bin from ground to start doing raman
analysis. This was very confusing because in Ansmann papers zet_0 is the
reference height, ie, the heighest bin. Changed zet_0 for bin1st.
Verified the start and end of loops and transformed some of them to simple matrix
multiplication. Raman_Manaus now uses my own function for the derivative... but it has
the same problem as the one before: fails at the start and end points. We should think
of introducing a polynomial fit and then taking the derivative.