Pastas vs mcmc vd Spek paper

[Camu-git]

No description provided.

[martinvonk]

Hi @Camu-git,

Looks very nice!

Maybe we could look into getting the data from the pastas-data repository. That way we could use the pastas load_dataset funciton:

ds = ps.load_dataset("spek_2017")
series = ds["heads"]["B32E0031_PUTTEN"].dropna()
precname = ds["metadata"].at[series.name, "Prec [m/d]"]
evapname = ds["metadata"].at[series.name, "Evap [m/d]"]
precipitation = ds["precipitation"].loc[:, precname].multiply(1e3)
evaporation = ds["evaporation"].loc[:, evapname].multiply(1e3)

This makes it easier to also test other head series because the precipitation and evaporation series are automatically coupled to the head series via the metadata table.

[martinvonk]

@Camu-git send me an email with a notebook that runs using the pastas-data and load_dataset. Maybe add that one to this pr? (and remove the old notebook if everything works.

[mbakker7]

I thought @Camu-git used more recent data. Are you sure the pastas-data datasets are up to date? Otherwise maybe move the updated datasets to pastas-data?