MONTE CARLO MODELING OF COMETARY ATMOSPHERES INCLUDING POLARIZATION