type :: radiation_source_t real(rk), allocatable :: qrad(:) ! [W/m^3], positive adds energy to gas end type radiation_source_t