free_radiation_state Subroutine

   subroutine free_radiation_state(state)
      type(radiation_state_t), intent(inout) :: state
      if (allocated(state%temperature)) deallocate(state%temperature)
      if (allocated(state%pressure)) deallocate(state%pressure)
      if (allocated(state%Y)) deallocate(state%Y)
      if (allocated(state%scalars)) deallocate(state%scalars)
      if (allocated(state%phase)) deallocate(state%phase)
      state%step = 0
      state%time = zero
      state%dt = zero
      state%ncells = 0
      state%n_species = 0
      state%n_scalars = 0
   end subroutine free_radiation_state