anchor.diff:       real(dp),allocatable :: j_nu_var(:)          ! independent emissivity variable
        real(dp),allocatable :: log10_j_nu_var(:)    ! independent emissivity variable [Log10]
        type(pdf_dp),allocatable :: j_nu(:)          ! emissivity
  -     type(pdf_dp),allocatable :: e_nu(:)          ! emissivity for diffusion (includes scattering)
   
        ! integer :: beta ! power of the photon energy sampling
        ! real(dp),allocatable  :: a(:) ! Energy of the emitted photon = a*nu^beta * incoming energy
  
positive.code:   public :: dust_scatter_peeloff
    public :: dust_sample_emit_probability
    public :: dust_sample_emit_frequency
    public :: dust_sample_emit_frequency_eta_nu
    public :: dust_jnu_var_pos_frac
  
    type dust
  
negative.code:     nu = 10._dp**nu
  
    end subroutine dust_sample_emit_frequency
  
    subroutine dust_scatter_peeloff(d,nu,a,s,a_req)
      implicit none
  