From 27b78272512cf229fa5dfa52a376a04877cc7dfd Mon Sep 17 00:00:00 2001 From: Jonas Gunz Date: Fri, 25 Aug 2023 12:48:51 +0200 Subject: move stuff --- grib.py | 56 ------------------------------------------------------- grib_info.py | 7 ------- test/grib.py | 56 +++++++++++++++++++++++++++++++++++++++++++++++++++++++ test/grib_info.py | 7 +++++++ 4 files changed, 63 insertions(+), 63 deletions(-) delete mode 100755 grib.py delete mode 100755 grib_info.py create mode 100755 test/grib.py create mode 100755 test/grib_info.py diff --git a/grib.py b/grib.py deleted file mode 100755 index 767f3d7..0000000 --- a/grib.py +++ /dev/null @@ -1,56 +0,0 @@ -#!/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') - -# GRIB-Objekt wichtige attribute: -# analDate -# validDate -# units -# level -# typeOfLevel -# name -# shortName - -for grb in grib: - #print(grb) - vals = grb.values - - lats, lons = grb.latlons() - - for k in grb.keys(): - print(k) - - print(grb.analDate) - print(grb.validDate) - print(grb.units) - print(grb.typeOfLevel) - break - - # TODO Data in xarray? - - 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()) - diff --git a/grib_info.py b/grib_info.py deleted file mode 100755 index 1f3b5c7..0000000 --- a/grib_info.py +++ /dev/null @@ -1,7 +0,0 @@ -#!/usr/bin/env python3 -import pygrib -import sys - -grib = pygrib.open(sys.argv[1]) -for grb in grib: - print(sys.argv[1], grb) diff --git a/test/grib.py b/test/grib.py new file mode 100755 index 0000000..767f3d7 --- /dev/null +++ b/test/grib.py @@ -0,0 +1,56 @@ +#!/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') + +# GRIB-Objekt wichtige attribute: +# analDate +# validDate +# units +# level +# typeOfLevel +# name +# shortName + +for grb in grib: + #print(grb) + vals = grb.values + + lats, lons = grb.latlons() + + for k in grb.keys(): + print(k) + + print(grb.analDate) + print(grb.validDate) + print(grb.units) + print(grb.typeOfLevel) + break + + # TODO Data in xarray? + + 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()) + diff --git a/test/grib_info.py b/test/grib_info.py new file mode 100755 index 0000000..1f3b5c7 --- /dev/null +++ b/test/grib_info.py @@ -0,0 +1,7 @@ +#!/usr/bin/env python3 +import pygrib +import sys + +grib = pygrib.open(sys.argv[1]) +for grb in grib: + print(sys.argv[1], grb) -- cgit v1.2.3