Function to compute the CRPS

This function was developed by Luis to compute the Continuous Ranked Probability Score for ensemble forecasts.