Mostly copy-paste from scoringRules:::crps_edf. Adjusted to the data table format, where the observation is a vector of the same length as the ensemble forecast,
but is just repeated (which is why only y[1]) is used.
crps_aux(y, dat)vector of length m with m identical entries, the observation
vector of length m containing the m ensemble forecasts