For MYSTIC with polarization the noise for Q in clear-sky atmosphere is large, because photons reflected at the surface only contribute to I, not to Q –> include some importance sampling method to enhance Rayleigh scattering?
Optical properties (aerosol_files can be used for scalar calculations) can not be specified in the required way (spectral ssa, spectral altitude dependent tau, constant phase matrix)
best solution now: optprop_file including spectral ssa, const phamat, ext=1, use wc.dat to specify tau (but this could only be done sequentially for all wavelengths, much overhead).
may be better: make a aerosol_files option for phase matrix (but this requires too much memory).