anchor.diff:       real(dp), allocatable :: log10_emissivity(:,:)
     end type dust_helper
   
  +  type filter
  +
  +     real(dp) :: n_wav
  +     real(dp), allocatable :: nu(:), tr(:)
  +
  +  end type filter
  +
     type image
   
        ! IMAGE PARAMETERS
  
positive.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
               if(img%uncertainties) then
                  img%sed2(inu,ir,im,io,:) = img%sed2(inu,ir,im,io,:) + img%tmp_stokes**2._dp * p%energy**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
  
    subroutine image_bin_raytraced(img,p,x_image,y_image,im,spectrum)
  
  
negative.code:     real(dp),intent(in) :: x_image, y_image
      integer,intent(in) :: im ! sub-image to bin into
      real(dp) :: log10_nu_image
      integer :: inu,io ! Bins
      integer :: iorig, ifilt
      real(dp) :: transmission
  
      if(img%compute_image.and..not.allocated(img%img)) call error('bin_photon','Image not allocated')
      if(img%compute_sed.and..not.allocated(img%sed)) call error('bin_photon','SED not allocated')
  