pvtend.scan_omega_blowup

pvtend.scan_omega_blowup(era5_w_glob: str, level_pa: float = 30000.0, threshold: float = 25.0, out_csv: str | PathLike | None = None) → DataFrame[source]

Flag ERA5 timestamps whose max absolute ω at level_pa exceeds threshold.

The threshold is a hard physical cutoff in Pa s⁻¹ — the canonical ±5 Pa/s rule applied at 300 hPa for ERA5 ω blowup detection.

Parameters:

  • era5_w_glob – Glob for ERA5 ω NetCDF files.

  • level_pa – Pressure level [Pa]; default 30000 (300 hPa).

  • threshold – Hard cutoff abs(ω) [Pa s⁻¹]; default 25.0 (raw ERA5 300 hPa ω rarely exceeds 2-3 Pa/s, so |ω|>10 Pa/s in derived ω_dry / ω_moist / ω_LHR is essentially always solver blowup).

  • out_csv – Optional output CSV path.

Returns:

DataFrame with columns timestamp, level_pa, max_abs_omega, threshold, n_exceed, exceedance_ratio for every timestamp whose spatial-max abs(ω) > threshold. n_exceed is the number of grid cells exceeding the threshold and exceedance_ratio is max_abs_omega / threshold (>1 → blowup).