According to http://adsabs.harvard.edu/abs/2015A%26A...576A.105P
equation (8) given by Padovani et al. 2012 ( http://adsabs.harvard.edu/abs/2012A%26A...543A..16P ) is wrong.
So, this line in sourcefunc.c
dSigma2 = maxp_g[pos].mol[ispec].dust[iline]g[pos].mol[ispec].knu[iline](0.5_angle[2]_angle[2]-1./3.);
should be replaced by the following line:
dSigma2 = maxp_g[pos].mol[ispec].dust[iline]g[pos].mol[ispec].knu[iline](angle[2]*angle[2]-2./3.);
According to http://adsabs.harvard.edu/abs/2015A%26A...576A.105P
equation (8) given by Padovani et al. 2012 ( http://adsabs.harvard.edu/abs/2012A%26A...543A..16P ) is wrong.
So, this line in sourcefunc.c
dSigma2 = maxp_g[pos].mol[ispec].dust[iline]g[pos].mol[ispec].knu[iline](0.5_angle[2]_angle[2]-1./3.);
should be replaced by the following line:
dSigma2 = maxp_g[pos].mol[ispec].dust[iline]g[pos].mol[ispec].knu[iline](angle[2]*angle[2]-2./3.);