aboutsummaryrefslogtreecommitdiff
path: root/test/cross.py
diff options
context:
space:
mode:
Diffstat (limited to 'test/cross.py')
-rwxr-xr-xtest/cross.py81
1 files changed, 0 insertions, 81 deletions
diff --git a/test/cross.py b/test/cross.py
deleted file mode 100755
index 5f79a39..0000000
--- a/test/cross.py
+++ /dev/null
@@ -1,81 +0,0 @@
-#!/usr/bin/env python3
-
-import matplotlib.pyplot as plt
-import matplotlib as mpl
-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, 0.9, 0.9),
- (1.0, 0.3, 0.3),
- ),
- 'green': (
- (0.0, 0.5, 0.5),
- (0.1, 0.9, 0.9),
- (1.0, 0.3, 0.3),
- ),
- 'blue': (
- (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
-ax = fig.add_subplot(5,1,(1,2))
-
-ax.set_ylabel('Pressure level [hPa]')
-
-#clc = ax.plot(data.step.values.astype('float64'), data.isobaricInhPa, data.ccl.transpose())
-#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, 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', pad=0.0, fraction=0.015)
-
-
-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, 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()
-