anchor.diff:      real(dp) :: nu
   
       type(angle3d_dp) :: a_scat
  +    type(angle3d_dp) :: a_final
   
       real(dp) :: P1,P2,P3,P4,norm
   
  
positive.code:     S%U = - sin_2_beta * RLS2 + cos_2_beta * RLS3
      S%V = RLS4
  
      a = a_final
  
    end subroutine scatter_stokes
  
    elemental real(dp) function B_nu(nu,T)
  
negative.code:        P2 = interp2d(d%mu,d%nu,d%P2,a_scat%cost,nu)
         P3 = interp2d(d%mu,d%nu,d%P3,a_scat%cost,nu)
         P4 = interp2d(d%mu,d%nu,d%P4,a_scat%cost,nu)
         call scatter_stokes(s,a,a_scat,a_req,P1,P2,P3,P4)
      end if
      a = a_req
    end subroutine dust_scatter_peeloff
  