Pull precipitation series in SWMM format from DWD data

dwd2dat.py

import warnings
warnings.filterwarnings("ignore")

from wetterdienst.provider.dwd.observation import DwdObservationRequest, \
    DwdObservationPeriod, DwdObservationResolution, DwdObservationParameter, DwdObservationDataset

# # all
# print("All available parameters")
# print(
#     DwdObservationRequest.discover()
# )

# # selection
# print("Selection of minute data")
# print(
#     DwdObservationRequest.discover(
#         filter_=DwdObservationResolution.MINUTE_1
#     )
# )

# define request filter by WGS84 position
request = DwdObservationRequest(
    parameter=[
        DwdObservationParameter.MINUTE_1.PRECIPITATION_HEIGHT
    ],
    resolution=DwdObservationResolution.MINUTE_1,
    period=DwdObservationPeriod.HISTORICAL,
    start_date="2020-01-01",
    end_date="2020-12-12"
).filter_by_rank(
    49.93121683122822,
    8.679585503284287,
    1
)

# fetch values for request and drop na-values and zero values
station_data = request.values.all().df
station_data.dropna(axis=0).head()
df = station_data[station_data.value != 0]

# save to file
station_data.to_csv("prec.dat", sep="\t", columns = ['station_id', 'date', 'value'], date_format = '%Y %m %d %H %M', index=False, header=False)