example_model_sunshine: shouldn't the NA filling be performed by era?

[luca-s]

I am referring to this line. I believe the median value of the features should be computed on a per-era basis to avoid introducing biases.