compute_RPS function

Compute RPS for a single forecast