Examples

Reflectivity quicklook

The quickest way to inspect a processed product is a time-range view of Ze or another 2D field. 

fig, ax = mrr.quicklook(variable="Ze", source="raprompro", vmin=0, vmax=40)

Doppler spectrum and spectrogram

The package can extract a single-gate Doppler spectrum or a range-by-velocity spectrogram for a selected time. 

import datetime

fig, path = mrr.plot_spectrum(
    datetime.datetime(2025, 3, 8, 12, 5, 0),
    target_range=2880.0,
    spectrum_var="spectrum_raw",
)

fig, path = mrr.plot_spectrogram(
    datetime.datetime(2025, 3, 8, 12, 5, 0),
    spectrum_var="spectrum_raw",
)

Microphysical process analysis

Layer-scale trends in Dm, Nw and LWC are analysed non-parametrically by default and then projected to RGB for hexagram-based classification. 

analysis = mrr.rain_process_analyze(
    period=(datetime.datetime(2025, 3, 8, 12, 0, 0), datetime.datetime(2025, 3, 8, 12, 10, 0)),
    layer=(1000.0, 2000.0),
    k=11,
    ze_th=-5.0,
    trend_method="kendall_theilsen",
    tau_zero_tol=0.05,
    min_points_trend=10,
)

classified = mrr.classify_rain_process(
    analysis=analysis,
    min_tau_strength=0.10,
)

fig, path = mrr.plot_classified_processes_on_hexagram(
    classified=classified,
    analysis=analysis,
    savefig=False,
)
Hexagram classification example
Example of classified rain-process samples plotted on the RGB hexagram.

For method comparison or backwards checks, pass trend_method="ols" to rain_process_analyze().