HIP
Go to https://hip.dataforsyningen.dk/
Images illustrating how to download potential and topography from HIP
Go to https://hip.dataforsyningen.dk/
Images illustrating how to find ETRS89/UTM32N coordinates
Alternatively use https://koordinattransformation.dk/Denmark to find the coordinate or convert from another coordinate representation. With this approach you need to truncate the coordinates because we need integer values.
Install daisy tools
git clone git@github.com:daisy-model/tools.git
cd tools
pip install .
Use prepare_hip_data_for_daisy to extract data. There are various options that can be printed with --help
prepare_hip_data_for_daisy --help
Assume we have downloaded topography data to topology/DK1_2020_100m_layers.nc and pressure data to head-elevation/dk1_2020_100m_head_10km_613_71_potentiale.nc. We can then extract data from coordinate 714450, 6135250 and display it in centimeters with the following
prepare_hip_data_for_daisy topology/DK1_2020_100m_layers.nc head-elevation/dk1_2020_100m_head_10km_613_71_potentiale.nc --x 714450 --y 6135250 --unit cm
which should output
============================= Soil column =============================
dk_model X Y terrain_height layer aquifer elevation thickness unit
0 DK1 714450 6135250 2671.237061 CompLayer_1 False 2471.237061 200.000000 cm
1 DK1 714450 6135250 2671.237061 CompLayer_2 False 2248.258545 222.978516 cm
2 DK1 714450 6135250 2671.237061 CompLayer_11 True -2751.741455 5000.000000 cm
============================= Terrain height =============================
2671.237 cm
============================= Top aquifer =============================
elevation layer name: CompLayer_11
head_elevation layer name: 0
dk2019 layer name: kalk
aquifer layer name: glw6
============================= Head elevation =============================
time pressure unit
0 1989-12-28 06:00:00 143.741943 cm
1 1990-01-12 06:00:00 143.800049 cm
2 1990-01-27 06:00:00 132.112305 cm
3 1990-02-11 06:00:00 125.428711 cm
4 1990-02-26 06:00:00 124.750732 cm
.. ... ... ...
725 2019-10-07 06:00:00 115.349609 cm
726 2019-10-22 06:00:00 107.164307 cm
727 2019-11-06 06:00:00 103.935791 cm
728 2019-11-21 06:00:00 107.864502 cm
729 2019-12-06 06:00:00 105.703125 cm
[730 rows x 3 columns]
We can also store the data in csv files by specifying and output directory
prepare_hip_data_for_daisy topology/DK1_2020_100m_layers.nc head-elevation/dk1_2020_100m_head_10km_613_71_potentiale.nc --x 714450 --y 6135250 --unit cm --outdir data-for-daisy/714450x-6135250y
Which should produce two csv files
data-for-daisy/714450x-6135250y/pressure.csv
data-for-daisy/714450x-6135250y/soil_column.csv