anchor.diff:      s = stokes_dp(1._dp,0._dp,0._dp,0._dp)
   
       energy_scaling = 1.
  -    ! energy_scaling = d%a(i)*nu**(d%beta)
   
       call random_sphere_angle3d(a)
   
  
positive.code:      type(pdf_dp),allocatable :: b_nu(:)          ! emissivity divided by opacity (blackbodies for LTE dust)
  
       logical :: is_lte ! Whether the emissivities assume therma emission from LTE dust
  
       ! integer :: beta ! power of the photon energy sampling
       ! real(dp),allocatable  :: a(:) ! Energy of the emitted photon = a*nu^beta * incoming energy
  
    end type dust
  
  
negative.code:        call set_pdf(d%b_nu(i),emiss_nu,emiss_jnu(i,:) / interp1d_loglog(d%nu, d%kappa_nu, emiss_nu),log=.true.)
      end do
  
    end subroutine dust_setup
  
    subroutine dust_jnu_var_pos_frac(d,specific_energy,jnu_var_id,jnu_var_frac)
  