anchor.diff:       integer :: n_dust
   
        integer :: group_id = 999
  +
  +     logical :: use_filters
  +     type(filter), allocatable :: filters(:)
   
     end type image
   
  
positive.code:        io = 1
      end if
  
      if(img%n_stokes == 4) then
         img%tmp_stokes = [p%s%i,p%s%q,p%s%u,p%s%v]
      else
  
negative.code:        if(img%compute_sed) then
            call find_sed_bin(img,x_image,y_image,ir)
            if(ir >= 1 .and. ir <= img%n_ap) then
               img%sed(inu,ir,im,io,:) = img%sed(inu,ir,im,io,:) + img%tmp_stokes * p%energy * transmission
               if(img%uncertainties) then
                  img%sed2(inu,ir,im,io,:) = img%sed2(inu,ir,im,io,:) + (img%tmp_stokes * p%energy * transmission) ** 2._dp
                  img%sedn(inu,ir,im,io,:) = img%sedn(inu,ir,im,io,:) + 1._dp
               end if
            end if
         end if
      end if
  
    end subroutine image_bin_single
  
    subroutine image_bin_raytraced(img,p,x_image,y_image,im,spectrum)
  
  