From ae709fdf1eb0df384835509a6c9500b69301d14c Mon Sep 17 00:00:00 2001 From: Jonas Gunz Date: Sun, 13 Aug 2023 21:02:42 +0200 Subject: asdf --- grib.py | 37 +++++++++++++++++++++++++++++++++++++ 1 file changed, 37 insertions(+) create mode 100755 grib.py (limited to 'grib.py') 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 + +# ///icon-d2_regular-lat-lon_pressure-level____ + +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()) + -- cgit v1.2.3