Return the unit provided by the primary log
function tem_logging_unit_primary() result(unit) ! --------------------------------------------------------------------------- !> integer :: unit ! --------------------------------------------------------------------------- unit = primary%funit(0) end function tem_logging_unit_primary