Actions

Difference between revisions of "Elle basic plus"

From Lrose Wiki

Line 1: Line 1:
  
== Basic+ elle tutorial ==
+
=== '''Basic+ elle tutorial''' ===
  
 
Like the basic elle tutorial, this workflow will convert a single raw radar file to the CfRadial format and estimate the surface rain rate, with the added tasks of estimating the sounding at the radar from GFS analysis and running RadxBeamBlock. Note, the current parameter files are currently optimized to work on Unix systems. This workflow is one way to test that the applications have been installed properly.
 
Like the basic elle tutorial, this workflow will convert a single raw radar file to the CfRadial format and estimate the surface rain rate, with the added tasks of estimating the sounding at the radar from GFS analysis and running RadxBeamBlock. Note, the current parameter files are currently optimized to work on Unix systems. This workflow is one way to test that the applications have been installed properly.
  
=== Run each program one by one on the command line ===
+
=== '''Run each program one by one on the command line''' ====
 
<ol start="1" style="list-style-type: decimal;">
 
<ol start="1" style="list-style-type: decimal;">
 
<li>Download the [https://drive.google.com/file/d/1FnxZMv1cmr6enP9S9nU-bqYLynq05OtU/view?usp=sharing elle_sound_bb.tar.gz] file</li></ol>
 
<li>Download the [https://drive.google.com/file/d/1FnxZMv1cmr6enP9S9nU-bqYLynq05OtU/view?usp=sharing elle_sound_bb.tar.gz] file</li></ol>
Line 51: Line 51:
 
<li>In the elle_sound_bb directory, you should now see 3 new directories: convert, rate, and qpe. You can look through the the files in those directories to inspect the new files and check the quality of the data.</li></ol>
 
<li>In the elle_sound_bb directory, you should now see 3 new directories: convert, rate, and qpe. You can look through the the files in those directories to inspect the new files and check the quality of the data.</li></ol>
  
=== Run programs from a sample script ===
+
==== '''Run programs from a sample script''' ====
 
<ol start="1" style="list-style-type: decimal;">
 
<ol start="1" style="list-style-type: decimal;">
 
<li>Follow steps 1-3 and 5-7 from the section above.</li></ol>
 
<li>Follow steps 1-3 and 5-7 from the section above.</li></ol>
Line 67: Line 67:
 
<li>In the elle_sound_bb directory, you should now see 3 new directories: convert, rate, and qpe. You can look through the the files in those directories to inspect the new files and check the quality of the data.</li></ol>
 
<li>In the elle_sound_bb directory, you should now see 3 new directories: convert, rate, and qpe. You can look through the the files in those directories to inspect the new files and check the quality of the data.</li></ol>
  
=== Example output ===
+
==== '''Example output''' ====
 
Shown below are some HawkEye snapshots of data from the files created during each step.
 
Shown below are some HawkEye snapshots of data from the files created during each step.
  

Revision as of 21:18, 11 January 2021

Basic+ elle tutorial

Like the basic elle tutorial, this workflow will convert a single raw radar file to the CfRadial format and estimate the surface rain rate, with the added tasks of estimating the sounding at the radar from GFS analysis and running RadxBeamBlock. Note, the current parameter files are currently optimized to work on Unix systems. This workflow is one way to test that the applications have been installed properly.

Run each program one by one on the command line =

  1. Download the elle_sound_bb.tar.gz file
  1. Extract contents into the desired directory
tar -zxvf elle_sound_bb.tar.gz

You'll end up with a directory structure that includes files to set environment variables, the raw radar files, and a sample script that will run all the programs in one go.

Elle sound bb dirs.png

  1. In /elle_sound_bb/env_dirs, edit the directories (note: there is purposely no slash at the beginning of PROJ). Together, HOMED/PROJ should direct you to the elle_sound_bb directory. RAW is the full path to the raw radar data in /elle_sound_bb.
export HOMED="/path/to/home/directory"
export PROJ=“short/path/to/elle_sound_bb”
export RAW="/path/to/elle_sound_bb/raw"
  1. Set the variables in the terminal and set the date. Print out the directory structure to make sure it's pointing to the right directory.
source env_dirs
source env_vars
export days=“20170602”
export daysbb="20170602"
echo $HOMED/$PROJ
  1. Download the GFS analysis from this website. Select the AIRS option for the 0.5º domain, which is the first GFS option available. Since the tutorial includes two different files, select all UTC cycles and enter the start and end dates (02 Jun 2017). Select the batch option and enter your email. Download the data when your order is ready and put it in a folder named gfs within the elle_sound_bb directory.
  1. Download the topography data using the following command for this particular case. FYI: the folder is approximately 6 GB. Unzip the contents and make note of the resulting directory.
wget -r https://dds.cr.usgs.gov/srtm/version2_1/SRTM3/Eurasia/
  1. In the params folder, modify line 53 of RadxBeamBlock.RCWF.new to indicate the directory of the unzipped topographic data.
  1. Convert the raw file to a CfRadial file
/path/to/RadxConvert -params ./params/RadxConvert.* -f $RAW/$RADAR_NAME/$days/*.raw*
  1. Run RadxBeamBlock for the RCWF radar location
/path/to/RadxBeamBlock -params ./params/RadxBeamBlock.$RADAR_NAME.*
  1. Convert the GFS Grib2 files to Mdv (Grib2toMdv) and then extract the sounding closest to the RCWF radar (Mdv2SoundingSpdb)
/path/to/Grib2toMdv -params ./params/Grib2toMdv.* -f ./gfs/*.grb2
/path/to/Mdv2SoundingSpdb -params ./params/Mdv2SoundingSpdb.* -f ./gfs/mdv/$i/*.mdv
  1. Run RadxRate to calculate Kdp, the PID, and various rain rates
/path/to/RadxRate -params ./params/RadxRate.* -f ./convert/$RADAR_NAME/$days/*.nc
  1. Run RadxQpe to estimate the near-surface rainfall using the hybrid method
/path/to/RadxQpe -params ./params/RadxQpe.* -f ./rate/$RADAR_NAME/$days/*.nc
  1. In the elle_sound_bb directory, you should now see 3 new directories: convert, rate, and qpe. You can look through the the files in those directories to inspect the new files and check the quality of the data.

Run programs from a sample script

  1. Follow steps 1-3 and 5-7 from the section above.
  1. Edit directories in run_radx.sh script (lines 3 and 12)
cd /PATH/TO/elle_sound_bb
rdir="/PATH/TO/lrose/bin"
  1. Make the file executable, if it isn't already
chmod u+x run_radx.sh
  1. Run shell script
./run_radx.sh
  1. In the elle_sound_bb directory, you should now see 3 new directories: convert, rate, and qpe. You can look through the the files in those directories to inspect the new files and check the quality of the data.

Example output

Shown below are some HawkEye snapshots of data from the files created during each step.

Example of the beam blockage fraction (0-1) from RadxBeamBlock. In this tutorial, output from RadxQpe will look for the lowest elevation angle with < 0.25 blockage.

Radxbeamblock beame.png

Example of the reflectivity field displayed in HawkEye created by RadxConvert.

Radxconvert dbz.png

Example of the Kdp field displayed in HawkEye created by RadxRate.

Radxrate kdp.png

Example of the hybrid rain rate field displayed in HawkEye created by RadxRate.

Radxrate hybrid low.png

Example of the estimated surface rain rate field displayed in HawkEye created by RadxQpe.

Radxqpe hybrid.png