#include <background.h>
temporary parameters and workspace passed to phase space distribution function