setup_scat_test2 Function

    function setup_scat_test2(dict) result(array)
    !! set up scattering test scene 2 with user defined tau and hgg

        use opticalProperties
        use sdfs,             only : sdf, box
        use vector_class

        type(toml_table), intent(inout) :: dict
        type(sdf), allocatable :: array(:)

        type(opticalProp_t) :: opt
        real(kind=wp) :: mus, mua, hgg, n, tau

        allocate(array(1))
        call get_value(dict, "tau", tau)
        call get_value(dict, "hgg", hgg)

        n = 1._wp
        hgg = hgg
        mua = 1e-17_wp
        mus = tau

        opt = mono(mus, mua, hgg, n)
        array(1) = box(vector(200._wp, 200._wp, 200._wp), opt, 2)

    end function setup_scat_test2