asdf

1 files changed, 37 insertions, 0 deletions

diff --git a/grib.py b/grib.py
new file mode 100755
index 0000000..c298450
--- /dev/null
+++ b/grib.py
@@ -0,0 +1,37 @@
+#!/usr/bin/env python3
+import pygrib
+
+# http://opendata.dwd.de/weather/nwp/icon-d2/grib/00/t/icon-d2_germany_regular-lat-lon_pressure-level_2023080200_000_1000_t.grib2.bz2
+
+# <BASE>/<RUN>/<PARAMETER>/icon-d2_regular-lat-lon_pressure-level_<INIT>_<OFFSET>_<LEVEL>_<PARAMETER>
+
+GRIBDIR='dwd_icon-d2'
+
+# Flugplatz Saal a. d. Saale
+target_lat='50.3108796'
+target_lon='10.0646952'
+
+# Find nearest coords in model
+# TODO
+
+grib = pygrib.open('dwd_icon-d2/combined.grib2')
+
+
+for grb in grib:
+    print(grb)
+    vals = grb.values
+
+    lats, lons = grb.latlons()
+
+    print(grb.name)
+    #print(grb.level)
+
+    #print('lats min/max: ', lats.shape, lats.max(), lats.min())
+    #print('lons min/max: ', lons.shape, lons.max(), lons.min())
+
+    #print(vals.shape)
+
+    #print(vals[100][::100])
+
+    #print(grb.latlons())
+