From 6e683bfbd5d58d2e9ec53ee37f1472b8dac4a7bb Mon Sep 17 00:00:00 2001
From: Jonas Gunz <himself@jonasgunz.de>
Date: Sun, 15 Oct 2023 23:16:59 +0200
Subject: test/cross

---
 test/cross.py | 46 +++++++++++++++++++++++++++++++++-------------
 1 file changed, 33 insertions(+), 13 deletions(-)

(limited to 'test')

diff --git a/test/cross.py b/test/cross.py
index 6194630..5f79a39 100755
--- a/test/cross.py
+++ b/test/cross.py
@@ -6,35 +6,46 @@ from matplotlib.colors import LinearSegmentedColormap
 
 import xarray as xr
 from metpy.interpolate import cross_section
+import metpy.calc as mpcalc
+from metpy.units import units
 
 clcov_cmap = {
     'red': (
         (0.0, 0.0, 0.0),
-        (0.1, 1.0, 1.0),
-        (1.0, 0.2, 0.2),
+        (0.1, 0.9, 0.9),
+        (1.0, 0.3, 0.3),
     ),
     'green': (
-        (0.0, 0.0, 0.0),
-        (0.1, 1.0, 1.0),
-        (1.0, 0.2, 0.2),
+        (0.0, 0.5, 0.5),
+        (0.1, 0.9, 0.9),
+        (1.0, 0.3, 0.3),
     ),
     'blue': (
-        (0.0, 1.0, 1.0),
-        (0.1, 1.0, 1.0),
-        (1.0, 0.2, 0.2),
+        (0.0, 0.9, 0.9),
+        (0.1, 0.9, 0.9),
+        (1.0, 0.3, 0.3),
     ),
 }
 
 mpl.colormaps.register(LinearSegmentedColormap('clcov', clcov_cmap))
 
+# backend_kwargs={'filter_by_keys':{'typeOfLevel': 'heightAboveGround','level':2}}
 data = xr.load_dataset('dwd_icon-eu/combined.grib2', engine='cfgrib')
+
+
 lat, lon = (47.96, 11.99)
 data = data.sel(latitude=lat, longitude = lon, method='nearest')
+
+#data = data.assign_coords(step=(data.step / (10**9 * 3600)))
+data = data.assign_coords(step=(data.step.values.astype(float) * units('ns')).to('hour'))
+
 print(data)
 
+fig = plt.figure(figsize=(5, 5), layout="constrained")
+
 
 # start figure and set axis
-fig, ax = plt.subplots(figsize=(5, 5))
+ax = fig.add_subplot(5,1,(1,2))
 
 ax.set_ylabel('Pressure level [hPa]')
 
@@ -42,20 +53,29 @@ ax.set_ylabel('Pressure level [hPa]')
 #clc = ax.imshow(data.ccl.transpose(), extent=(data.step.values.astype(float).min(), data.step.values.astype(float).max(), data.isobaricInhPa.min(), data.isobaricInhPa.max()), aspect='auto', cmap='Blues_r', vmin=0, vmax=100)
 #plt.colorbar(clc, label='clcov')
 # Blues_r
-clc = ax.contourf(data.step.values.astype('float64'), data.isobaricInhPa, data.ccl.transpose(), cmap='clcov', vmin=0, vmax=100, levels=9)
+clc = ax.contourf(data.step, data.isobaricInhPa, data.ccl.transpose(), cmap='clcov', vmin=0, vmax=100, levels=9)
 
 # use Format parameter for n/8
-plt.colorbar(clc, label='cloudcov', extendfrac=None, ticks=[100*n/8 for n in range(9)], format=lambda x,_: f'{int(x/12.5)}/8')
+plt.colorbar(clc, label='cloudcov', extendfrac=None, ticks=[100*n/8 for n in range(9)], format=lambda x,_: f'{int(x/12.5)}/8', pad=0.0, fraction=0.015)
 
 
-cf = ax.contour(data.step.values.astype('float64'), data.isobaricInhPa, data.t.metpy.convert_units('degC').transpose())
+cf = ax.contour(data.step, data.isobaricInhPa, data.t.metpy.convert_units('degC').transpose())
 ax.clabel(cf, inline=True, fontsize=10)
 #plt.colorbar(cf, pad=0, aspect=50)
 #plt.colorbar(cf)
 
-barbs = ax.barbs(data.step.values.astype('float64'), data.isobaricInhPa, data.u.transpose(), data.v.transpose())
+barbs = ax.barbs(data.step, data.isobaricInhPa, data.u.transpose(), data.v.transpose())
 #ax.barbs(data.u, data.v, color='black', length=5, alpha=0.5)
 
 ax.invert_yaxis()
+
+### Second plot
+
+ax2 = fig.add_subplot(5,1,3,sharex=ax)
+ax2.plot(data.step, data.t2m.metpy.convert_units('degC').transpose(), color='red', label='Temperature (2m)')
+ax2.plot(data.step, mpcalc.dewpoint_from_relative_humidity(data.t2m, data.r2).transpose(), color='blue', label='Dewpoint (2m)')
+ax2.set_ylabel('Temperature [degC]')
+ax2.legend(loc='lower right')
+
 plt.show()
 
-- 
cgit v1.2.3