diff options
| -rw-r--r-- | config.yaml | 44 | ||||
| -rwxr-xr-x | plotter/meteogram.py | 26 |
2 files changed, 63 insertions, 7 deletions
diff --git a/config.yaml b/config.yaml index d9e5d45..74a9377 100644 --- a/config.yaml +++ b/config.yaml @@ -21,8 +21,10 @@ aggregator: - hbas_con - htop_con - htop_dc + - tot_prec pressure_levels: [1000, 950, 925, 900, 875, 850, 825, 800, 775, 700, 600, 500, 400, 300, 250, 200, 150, 100] - steps: [0, 3, 6, 9, 12, 15, 18, 21, 24, 27, 30, 33, 36, 39, 42, 45, 48] + steps: [0, 3, 6, 9, 12, 15, 18, 21, 24, 27, 30, 33, 36, 39, 42, 45, 48, 51, 54, 57, 60, 63, 66, 69, 72] + description: FORECAST DWD ICON-EU plotter: - module: plotter.meteogram aggregator: icon_eu @@ -31,6 +33,9 @@ plotter: - name: meteogram_antersberg lat: 47.96 lon: 11.99 + - name: meteogram_munich + lat: 48.16 + lon: 11.57 - module: 'plotter.vertical_from_grib' aggregator: icon_eu @@ -40,10 +45,19 @@ plotter: lon: 11.99 name: skewt_antersberg analysis: lcl + - lat: 11.57 + lon: 48.16 + name: skewt_munich + analysis: lcl - module: 'plotter.horizontal' aggregator: icon_eu output: web/data plots: + - name: dry_con_top + area: null + layers: + - layertype: raster + field: HTOP_DC - name: t_fi_850 area: null layers: @@ -58,6 +72,34 @@ plotter: plot_units: '_gpdm' contours: 5 clabels: true + - name: t_fi_500 + area: null + layers: + - layertype: raster + field: t + level: 500 + plot_units: degC + colormap: jet + - layertype: contour + field: z + level: 500 + plot_units: '_gpdm' + contours: 5 + clabels: true + - name: t_fi_200 + area: null + layers: + - layertype: raster + field: t + level: 200 + plot_units: degC + colormap: jet + - layertype: contour + field: z + level: 200 + plot_units: '_gpdm' + contours: 5 + clabels: true - name: pmsl_t850 area: null layers: diff --git a/plotter/meteogram.py b/plotter/meteogram.py index a360fef..a04fe34 100755 --- a/plotter/meteogram.py +++ b/plotter/meteogram.py @@ -13,6 +13,8 @@ from metpy.units import units import misc +HEIGHT = 11 + def run(data, plots, output='.'): misc.create_output_dir(output) index = [] @@ -32,11 +34,11 @@ def _plot(data, output, name, lat, lon): init_str = init.strftime('%d %b %Y - %HUTC') init_for_filename = init.strftime('%Y-%m-%d-%HUTC') - fig = plt.figure(figsize=(10, 10), layout="constrained") + fig = plt.figure(figsize=(12, 10), layout="constrained") # start figure and set axis - ax = fig.add_subplot(5,1,(1,2)) + ax = fig.add_subplot(HEIGHT,1,(1,4)) ax.set_ylabel('Pressure level [hPa]') @@ -54,15 +56,16 @@ def _plot(data, output, name, lat, lon): ax.invert_yaxis() ### Temp + Dewpoint - ax2 = fig.add_subplot(5,1,3,sharex=ax) + ax2 = fig.add_subplot(HEIGHT,1,(5,6),sharex=ax) ax2.plot(data.valid_time, data.t2m.metpy.convert_units('degC').transpose(), color='red', label='Temperature (2m)') ax2.plot(data.valid_time, mpcalc.dewpoint_from_relative_humidity(data.t2m, data.r2).transpose(), color='blue', label='Dewpoint (2m)') ax2.plot(data.valid_time, data.sel(isobaricInhPa=850.0).t.metpy.convert_units('degC').transpose(), color='grey', label='Tempreature (850hPa)') ax2.set_ylabel('Temperature [degC]') ax2.legend(loc='lower right') + ## MSLP - ax3 = fig.add_subplot(5,1,4,sharex=ax) + ax3 = fig.add_subplot(HEIGHT,1,(7,8),sharex=ax) ax3.plot(data.valid_time, data.prmsl.metpy.convert_units('hPa').transpose(), color='black', label='Temperature (2m)') ax3.set_ylabel('Mean Sea Level Pressure [hPa]') #ax3.legend(loc='lower right') @@ -77,8 +80,19 @@ def _plot(data, output, name, lat, lon): height=(data.hcct.metpy.convert_units('km')-data.HBAS_CON.metpy.convert_units('km')).transpose(), align='edge', width=np.timedelta64(3, 'h')) - ### Info Lines + # Precip + ax5 = fig.add_subplot(HEIGHT,1,(9,10),sharex=ax) + ax5.set_ylabel('Total precipitation [mm]') + ax5.set_ylim(0, 30) + ax5.bar(data.valid_time[:-1], data.tp.diff('step').transpose(), width=np.timedelta64(3, 'h'), + align='edge', alpha=0.7, color='green') + ax6 = ax5.twinx() + ax6.set_ylabel('Snow depth [m]') + ax6.set_ylim(bottom=0) + ax6.plot(data.valid_time, data.sde.transpose(), color='blue') + + ### Info Lines info_lines = [] init = misc.np_time_convert(data.time.values) init_str = init.strftime('%d %b %Y - %HUTC') @@ -89,7 +103,7 @@ def _plot(data, output, name, lat, lon): if '_description' in data.attrs: info_lines.append(data.attrs['_description']) - ax_text = fig.add_subplot(5, 1, 5) + ax_text = fig.add_subplot(HEIGHT, 1, 11) ax_text.text(0, 0, '\n'.join(info_lines), ha='left', va='center', size=10, fontfamily='monospace') ax_text.axis("off") |