Numerical codes for the time-dependent radiative transfer equation by FDM-DOM

The numerical codes provide the calculation of the time-dependent radiative transfer equation (RTE) using the Galerkin quadrature method with the 6th order level-symmetric even quadrature set based on the finite difference and discrete ordinates methods for highly forward scattering media under the non-reentry boundary condition. 【Numerical schemes for the RTE-calculations】 ・3rd order weighted essentially non oscillatory scheme (spatial discretization) ・3rd order total variation diminishing-Runge-Kutta method (temporal discretization) ・Galerkin quadrature method (angular discretization) 【The medium geometry and optical properties (default)】 ・Rectangular medium with the source-detector distance of 0.4 cm ・Spatial grid size: 0.02 cm and timestep: 0.5 ps ・Absorption coefficient: 0.2 cm-1 ・Scattering coefficient: 100 cm-1 ・Anisotropic factor: 0.9 ・Refractive index: 1.4