hoedanigheid (vertical reference level) is defined in multiple metadata fields

[veenstrajelmer]

When filtering waterlevel extremes, we have to distinguish between NAP/MSL stations in the Groepering column (or include all). The vertical reference level (NAP or MSL) is already defined in the Hoedanigheid column, so it does not make sense to also distinguish between these in the Groepering : GETETM2/GETETBRKD2/GETETMSL2/GETETBRKDMSL2

import ddlpy

locations = ddlpy.locations()

bool_station = locations.index.isin(['HOEKVHLD'])
bool_grootheid_meas = locations['Grootheid.Code'].isin(['WATHTE'])
bool_grootheid_astro = locations['Grootheid.Code'].isin(['WATHTBRKD'])
bool_hoedanigheid = locations['Hoedanigheid.Code'].isin(['NAP'])
bool_groepering_ts = locations['Groepering.Code'].isin(['NVT'])
bool_groepering_ext = locations['Groepering.Code'].isin(['GETETM2','GETETBRKD2','GETETMSL2','GETETBRKDMSL2'])

# get WATHTE locations (ts) >> measured waterlevel
locs_wathte = locations.loc[bool_station & bool_grootheid_meas & bool_hoedanigheid & bool_groepering_ext]
# get WATHTBRKD locations (ts) >> computed astronomical waterlevel
locs_wathtbrkd = locations.loc[bool_station & bool_grootheid_astro & bool_hoedanigheid & bool_groepering_ext]

[TvLoon-RWS]

This will be partly covered with the ProcesType 'meting' 'astro' or 'verwachting' in Wadar WaterWebservices