Description of pyrad config files#

0 : Plot a volume scan (pyrad_config/00_plot_vol)#


This pyrad configuration file does the following:

  • Read polar data from JMA

  • Perform a series of plots from the reflectivity data

    • PPI plot

    • PseudoRHI plot

    • Bscope plot

    • Histogram plot

    • Plot along specific polar coordinates

How to run

cd pyrad_config/00_plot_vol; main_process_data.py --cfgpath ./ --starttime 20230801000000 --endtime 20230802000000 jma_plot_vol.txt

1 : Create a Cartesian grid and plot it (pyrad_config/01_grid)#


This pyrad configuration file does the following:

  • Read polar reflectivity and RhoHV data from JMA

  • Convert polar data to a Cartesian grid

  • Perform a series of plots from the Cartesian data

    • Histogram plot

    • Surfaces at specified heights

    • Latitude slice

    • Longitude slice

    • Cross-section between two points

How to run

cd pyrad_config/01_grid; main_process_data.py --cfgpath ./ --starttime 20230801000000 --endtime 20230802000000 jma_grid.txt

2 : Ground-clutter removal (pyrad_config/02_clutter_filter)#


This pyrad configuration file does the following:

  • Read polar data from MeteoFrance radar (Zh, RhoHV, ZDR, PhiDP)

  • Use this data to detect ground-clutter

  • Remove ground-clutter from the radar data

  • Make PPI plots of filtered Zh, RhoHV, ZDR and PhiDP

How to run

 cd pyrad_config/02_clutter_filter; main_process_data.py --cfgpath ./ --starttime 20210114101500 --endtime 20210114120000 mf_pam_mtcy_filter.txt

3 : Compute KDP (pyrad_config/03_phidp_kdp)#


This pyrad configuration file does the following:

  • Read polar data from MeteoFrance radar (Zh, RhoHV, ZDR, PhiDP)

  • Read a GRIB iso0 height file from MeteoFrance

  • Correct ground clutter as in 02_clutter_filter

  • Perform smoothing of raw PhiDP

  • Estimate KDP with different methods:

    • Vulpiani filtering method

    • Maesaka variational method

    • Linear programming method

    • Kalman filter method

    • Moving least-squares filter

  • Perform attenuation correction

    • ZPHI method

    • Philinear method

  • Perform PPI plots of KDP, attenuation corrected Zh, specific attenuation, PIA and PIDA

How to run

cd pyrad_config/03_phidp_kdp; main_process_data.py --cfgpath ./ --starttime 20210114101500 --endtime 20210114120000 mf_pam_mtcy_kdp.txt

4 : Classifiy hydrometeors (pyrad_config/04_hydroclass)#


This pyrad configuration file does the following:

  • Read polar data from MeteoFrance radar (Zh, RhoHV, ZDR, PhiDP)

  • Read a GRIB iso0 height file from MeteoFrance

  • Compute hydrometeor classification

    • With semi-supervised method (Besic, 2016)

    • By mapping MeteoFrance operational classification

  • Plot PPIs of hydrometeor classification

How to run

cd pyrad_config/04_hydroclass; main_process_data.py --cfgpath ./ --starttime 20200516200500 --endtime 20200516220000 mf_odim_coll_hydro.txt

5 : Melting layer detection (pyrad_config/05_ml)#


This pyrad configuration file does the following:

  • Read polar data from MeteoFrance radar (Zh, RhoHV, ZDR, PhiDP)

  • Read a GRIB iso0 height file from MeteoFrance

  • Compute melting layer (ML) extent

    • With Giangrande method

    • With MeteoFrance method

  • Estimate VPR (vertical profile of reflectivity) with MeteoFrance algorithm

  • Correct reflectivity with VPR

  • Create relevant products

    • Timeseries of ML and VPR

    • PPIs of corrected reflectivity and ML extent overlaid on radar moments

    • Altitude plots of apparent VPR and theoretical VPR

How to run

cd pyrad_config/05_ml; main_process_data.py --cfgpath ./ --starttime 20200516200500 --endtime 20200516220000 mf_odim_coll_ml.txt

6 : Quantitative precipitation estimation (pyrad_config/06_rr)#


This pyrad configuration file does the following:

  • Read polar data from MeteoFrance radar (Zh, RhoHV, ZDR, PhiDP)

  • Read a GRIB iso0 height file from MeteoFrance

  • Compute attenuation with ZPHI method

  • Estimate precipitation with 4 methods

    • From Zh with R-Z relationship

    • From KDP with R-KDP relationship

    • From Ah (specific attenuation) with R-A relationship

    • From hydrometeor classification using different estimators for every class

  • Compute rainfall accumulation over one hour

How to run

cd pyrad_config/06_rr; main_process_data.py --cfgpath ./ --starttime 20200516200500 --endtime 20200516220000 mf_odim_coll_rr.txt

7 : Compute velocity azimuth display (pyrad_config/07_VAD)#


This pyrad configuration file does the following:

  • Read Doppler data from MeteoFrance radar

  • Compute velocity azimuth display (VAD) from Doppler velocity

  • Plot and save (csv) vertical profiles of wind:

    • U,V,W wind components

    • Wind speed and direction

How to run

cd pyrad_config/07_VAD; main_process_data.py --cfgpath ./ --starttime 20210114101500 --endtime 20210114120000 mf_pag_mtcy_doppler.txt

8 : Compute dual-Doppler estimation (pyrad_config/08_DDA)#


This pyrad configuration file does the following:

  • Read polar data from MeteoSwiss C-band and X-band radar

  • Perform ground clutter removal for X-band data

  • Load C-band operational dealiased Doppler velocity

  • Perform dual-Doppler estimation using the pyDDA wrapper with sounding data from wmo:06610 for initial conditions

  • Generate plots:

    • 2D wind arrow maps at 1000, 2000 and 3000 m.a.s.l

    • Cross-sections of U,V and W wind components along runway 32 (RW32) of Zurich Airport

How to run

cd pyrad_config/08_DDA; main_process_data.py --cfgpath ./ --starttime 20231127174000 --endtime 20231127175000 dda.txt

9 : Quasi-vertical profiles (pyrad_config/09_QVP)#


This pyrad configuration file does the following:

How to run

cd pyrad_config/09_QVP; main_process_data.py --cfgpath ./ --starttime 20200516200500 --endtime 20200516220000 mf_odim_coll_qvp.txt

10 : Radar visibility from a DEM with the GECSX tool (pyrad_config/10_visibility/)#


This pyrad configuration file does the following:

  • Read a DEM of the Swiss topography in geotiff format

  • Estimate the static radar visibility and clutter power/reflectivity using either data from an X-band mobile radar (visibility.txt) or using user specified radar parameters and no actual radar data (visibility_norad.txt).

  • Create PPI plots of the clutter and the visibility

  • Create Grid map plots of several other products generated by GECSX (e.g. slope, aspect, min. visible elevation angle,…)

How to run

cd pyrad_config/10_visibility/; main_process_gecsx.py --cfgpath ./ visibility_norad.txt # no radar data
cd pyrad_config/10_visibility/; main_process_gecsx.py --cfgpath ./ --starttime 20231127170000 --endtime 20231127180000 visibility.txt # with radar data