From b398bb68c6c442f6f23d6caf2854fe00e6834c01 Mon Sep 17 00:00:00 2001
From: Jonas Gunz <himself@jonasgunz.de>
Date: Fri, 25 Aug 2023 12:50:50 +0200
Subject: move xar

---
 vertical_from_grib.py | 71 +++++++++++++++++++++++++++++++++++++++++++++++++++
 xar.py                | 71 ---------------------------------------------------
 2 files changed, 71 insertions(+), 71 deletions(-)
 create mode 100755 vertical_from_grib.py
 delete mode 100755 xar.py

diff --git a/vertical_from_grib.py b/vertical_from_grib.py
new file mode 100755
index 0000000..759bca4
--- /dev/null
+++ b/vertical_from_grib.py
@@ -0,0 +1,71 @@
+#!/usr/bin/env python3
+
+import datetime
+
+import xarray as xr
+from metpy.units import units
+import metpy.calc as mpcalc
+import numpy as np
+
+import skewt
+
+def np_time_convert(dt64):
+    unix_epoch = np.datetime64(0, 's')
+    one_second = np.timedelta64(1, 's')
+    seconds_since_epoch = (dt64 - unix_epoch) / one_second
+
+    return datetime.datetime.utcfromtimestamp(seconds_since_epoch)
+
+config = {
+    'source':'dwd_icon-d2/combined.grib2',
+    'plots':[
+        {
+            'lat':47.9626,
+            'lon':11.9964,
+            'name':'Antersberg'
+        },
+    ]
+}
+
+def run(config):
+    data = xr.load_dataset(config['source'], engine='cfgrib')
+
+    for plot in config['plots']:
+        _plot(data, **plot)
+
+def _plot(data, lat, lon, name):
+    for_temp = data.sel(latitude=lat, longitude = lon, method='nearest')
+    for_temp = for_temp[['r', 't', 'u', 'v']]
+
+    for step in for_temp.coords['step']:
+        this_step = for_temp.sel(step=step)
+
+        p = this_step.coords['isobaricInhPa'].values * units.hPa
+        T = this_step.t.values * units.K
+        relHum = this_step.r.values * units.percent
+        Td = mpcalc.dewpoint_from_relative_humidity(T, relHum)
+        u = this_step.u.values * (units.m / units.s)
+        v = this_step.v.values * (units.m / units.s)
+
+        valid = np_time_convert(step.valid_time.values)
+        init = np_time_convert(step.time.values)
+
+        valid_str = valid.strftime('%d %b %Y - %HUTC')
+        init_str = init.strftime('%d %b %Y - %HUTC')
+        hours_since_init_str = str(int(this_step.step.values / np.timedelta64(1,'h'))).zfill(2)
+
+        skt = skewt.Skewt(p=p, T=T, Td=Td)
+        skt.addWindUV(u, v)
+        skt.addInfo(f'{name} INIT+' + hours_since_init_str)
+        skt.addInfo(f"VALID: {valid_str}")
+        skt.addInfo(f"INIT : {init_str}")
+        skt.addInfo(f"LAT {lat} LON {lon}")
+        skt.addAnalysis(shade=True, analysis='lcl')
+        skt.addInfo("FORECAST DWD ICON-D2")
+
+        init_for_filename = init.strftime('%Y-%m-%d-%HUTC')
+
+        skt.plot(filename=f'skewt_{name}_{init_for_filename}+{hours_since_init_str}.png')
+
+if __name__ == '__main__':
+    run(config)
diff --git a/xar.py b/xar.py
deleted file mode 100755
index 759bca4..0000000
--- a/xar.py
+++ /dev/null
@@ -1,71 +0,0 @@
-#!/usr/bin/env python3
-
-import datetime
-
-import xarray as xr
-from metpy.units import units
-import metpy.calc as mpcalc
-import numpy as np
-
-import skewt
-
-def np_time_convert(dt64):
-    unix_epoch = np.datetime64(0, 's')
-    one_second = np.timedelta64(1, 's')
-    seconds_since_epoch = (dt64 - unix_epoch) / one_second
-
-    return datetime.datetime.utcfromtimestamp(seconds_since_epoch)
-
-config = {
-    'source':'dwd_icon-d2/combined.grib2',
-    'plots':[
-        {
-            'lat':47.9626,
-            'lon':11.9964,
-            'name':'Antersberg'
-        },
-    ]
-}
-
-def run(config):
-    data = xr.load_dataset(config['source'], engine='cfgrib')
-
-    for plot in config['plots']:
-        _plot(data, **plot)
-
-def _plot(data, lat, lon, name):
-    for_temp = data.sel(latitude=lat, longitude = lon, method='nearest')
-    for_temp = for_temp[['r', 't', 'u', 'v']]
-
-    for step in for_temp.coords['step']:
-        this_step = for_temp.sel(step=step)
-
-        p = this_step.coords['isobaricInhPa'].values * units.hPa
-        T = this_step.t.values * units.K
-        relHum = this_step.r.values * units.percent
-        Td = mpcalc.dewpoint_from_relative_humidity(T, relHum)
-        u = this_step.u.values * (units.m / units.s)
-        v = this_step.v.values * (units.m / units.s)
-
-        valid = np_time_convert(step.valid_time.values)
-        init = np_time_convert(step.time.values)
-
-        valid_str = valid.strftime('%d %b %Y - %HUTC')
-        init_str = init.strftime('%d %b %Y - %HUTC')
-        hours_since_init_str = str(int(this_step.step.values / np.timedelta64(1,'h'))).zfill(2)
-
-        skt = skewt.Skewt(p=p, T=T, Td=Td)
-        skt.addWindUV(u, v)
-        skt.addInfo(f'{name} INIT+' + hours_since_init_str)
-        skt.addInfo(f"VALID: {valid_str}")
-        skt.addInfo(f"INIT : {init_str}")
-        skt.addInfo(f"LAT {lat} LON {lon}")
-        skt.addAnalysis(shade=True, analysis='lcl')
-        skt.addInfo("FORECAST DWD ICON-D2")
-
-        init_for_filename = init.strftime('%Y-%m-%d-%HUTC')
-
-        skt.plot(filename=f'skewt_{name}_{init_for_filename}+{hours_since_init_str}.png')
-
-if __name__ == '__main__':
-    run(config)
-- 
cgit v1.2.3