- Save R_sys_ns and R_sys to file
- Remove R_sys_ns and R_sys and replace with refl(...)
- Return D_sys for ground measurement
- Allow calculation of absolute values of Raman cross sections, not finalized yet
- replace FORTRAN77 atm.f with FORTRAN90 atm.f90
- make ss_inelm.f90
- Replace 1+d0m with Fourier expansion quantities (e.q. Bplus, etc.) in all subroutines
- Remove layer dependancy of various optical properties
- Use plm_0, plm_p, plm_m of gsf.f90 instead of setplm.f90 for multiple scatterigng calc, and remove setplm.f90
- double/add from top to bottom instead of bottom to top as in current implementaion
- Check for empty wavenr bins (when dnu is very fine ~1cm-1).
- Use single precision "real(kind=sp)" for wavenr/wavelength arrays
- Start with second order Raman scattering instead of single scattering in basiclayer.f90 (as in elastic doubling adding code of De Haan, 1984) with corresponding 2nd order wavenrgrid
- Extend to include polarization, make DAS a vector-code
- Include vibrational Raman scattering by air (involves much larger wavenr shifts than +-800 cm-1, and involves many empty wavenrbins)
- Include aerosol scattering. For Fourier moments m=3 and higher: make sure to use only elastic adding and doubling.
- Include vibrational Raman scattering
- Alternative radiative transfer model approach for dnu < 1cm-1 (e.g needed to simulate O2A band) use large (1:nwave_fine,1:nwave_fine) matrix and fill it instead of calculating multiple times (1:nwave_fine,1:nwave_fine) for every n=1,nwave_coarse as is done in the current implementation.
- Add example of bidrectional reflection function (BDRF)
- Ocean interface (medium with different refraction index and rough surface) with oceanic Raman scattering