As it is desireable to compute the RPSS of a weighted ensemble against the unweighted reference I've split the weights variable inside the RPSS function into weights_ref and weights_exp.
If you have any suggestions in the variable naming, documentation or methodology, please let me know.