Actions

Difference between revisions of "Elle grid"

From Lrose Wiki

 
(13 intermediate revisions by the same user not shown)
Line 1: Line 1:
=== '''elle grid tutorial''' ===
+
=== '''Elle grid tutorial''' ===
  
 
This workflow will convert several raw radar files to the CfRadial format, interpolate the data to a cartesian grid, and apply a convective stratiform separation algorithm.
 
This workflow will convert several raw radar files to the CfRadial format, interpolate the data to a cartesian grid, and apply a convective stratiform separation algorithm.
  
==== '''Run each program one by one on the command line''' ====
+
 
 +
=== '''Run program 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/1Jd4LIZoy2CvZ2t3fsgLmVd1QFlIPcLhk/view?usp=sharing elle_grid.tar.gz] file</li></ol>
 
<li>Download the [https://drive.google.com/file/d/1Jd4LIZoy2CvZ2t3fsgLmVd1QFlIPcLhk/view?usp=sharing elle_grid.tar.gz] file</li></ol>
Line 19: Line 20:
 
<li>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.</li></ol>
 
<li>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.</li></ol>
 
  export RADAR_NAME="KHGX"
 
  export RADAR_NAME="KHGX"
  export days=“20170827”
+
  export days=20170827
 
  export rdir="/path/to/lrose/bin"
 
  export rdir="/path/to/lrose/bin"
 
  echo $HOMED/$PROJ
 
  echo $HOMED/$PROJ
Line 27: Line 28:
 
<ol start="6" style="list-style-type: decimal;">
 
<ol start="6" style="list-style-type: decimal;">
 
<li>Regrid the data to a cartesian grid and apply the convective stratiform algorithm</li></ol>
 
<li>Regrid the data to a cartesian grid and apply the convective stratiform algorithm</li></ol>
  $rdir/Radx2Grid -params ./params/RadxGrid.* -f convert/$RADAR_NAME/$days/*.nc.
+
  $rdir/Radx2Grid -params ./params/RadxGrid.* -f convert/$RADAR_NAME/$days/*.nc
  
  
 +
=== '''Expected output''' ===
 
After running these two programs, you should end up with NetCDF files that show the following output.   
 
After running these two programs, you should end up with NetCDF files that show the following output.   
  
[[File:Radx2grid_dbz.png|250px]]         [[File:Radxgrid_convstrat.png|250px]]
+
[[File:Grid_tutorial_ref_leg.png|275px]][[File:Grid_tutorial_ref_field.png|275px]]  [[File:Grid_tutorial_cs_leg.png|275px]][[File:Grid_tutorial_cs_field.png|275px]]
  
==== '''Make modifications to the gridding procedure in RadxGrid.new''' ====
+
 
 +
=== '''Modify the gridding procedure in RadxGrid.new''' ===
 
<ol start="1" style="list-style-type: decimal;">
 
<ol start="1" style="list-style-type: decimal;">
 
<li>Select the interpolation mode in interp_mode on line 225. The current setting is a cartesian grid (INTERP_MODE_CART).</li></ol>
 
<li>Select the interpolation mode in interp_mode on line 225. The current setting is a cartesian grid (INTERP_MODE_CART).</li></ol>
Line 46: Line 49:
 
<li>Set the grid_projection on line 498 (the current setting is PROJ_FLAT).</li></ol>
 
<li>Set the grid_projection on line 498 (the current setting is PROJ_FLAT).</li></ol>
 
<ol start="6" style="list-style-type: decimal;">
 
<ol start="6" style="list-style-type: decimal;">
<li>Choose whether to censor output based on polarimetric variables (e.g., PID category or Rhohv values) on line 1278 (current setting is TRUE and currently censors data where Rhohv is outside [0.8, 1.05] deg/km and PID is outside [1, 14]).</li></ol>
+
<li>Choose whether to censor output based on polarimetric variables (e.g., PID category or Rhohv values) on line 1278 (current setting is TRUE and currently censors data where Rhohv is outside [0.8, 1.05] deg/km and REF is outside [-10, 100]).</li></ol>
 
<ol start="7" style="list-style-type: decimal;">
 
<ol start="7" style="list-style-type: decimal;">
 
<li>Choose whether to remove long range rays in NEXRAD data on line 1468 (current setting is TRUE).</li></ol>
 
<li>Choose whether to remove long range rays in NEXRAD data on line 1468 (current setting is TRUE).</li></ol>
 
<ol start="8" style="list-style-type: decimal;">
 
<ol start="8" style="list-style-type: decimal;">
 
<li>Choose whether to identify convective and stratiform precipitation and set the thresholds in lines 2416-2597 (current setting is TRUE). The threshold for definitive convection is on line 2491, the radius of convection is on line 2504, and choosing to write the partition to the file is on line 2559. However, I encourage suggest reading through the entire section and determining whether the values are appropriate for the precipitation type.</li></ol>
 
<li>Choose whether to identify convective and stratiform precipitation and set the thresholds in lines 2416-2597 (current setting is TRUE). The threshold for definitive convection is on line 2491, the radius of convection is on line 2504, and choosing to write the partition to the file is on line 2559. However, I encourage suggest reading through the entire section and determining whether the values are appropriate for the precipitation type.</li></ol>
 +
 +
[[File:Radx2Grid_censor.png|550px]]

Latest revision as of 17:13, 3 February 2023

Elle grid tutorial

This workflow will convert several raw radar files to the CfRadial format, interpolate the data to a cartesian grid, and apply a convective stratiform separation algorithm.


Run program on the command line

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

You'll end up with a directory that includes the raw radar files and parameter files needed to run RadxConvert and Radx2Grid.

  1. Set the directories (note: there is purposely no slash at the beginning of PROJ). Together, HOMED/PROJ should direct you to the elle_grid directory. RAW is the full path to the raw radar data in /elle_grid.
export HOMED="/path/to/home/directory"
export PROJ=“short/path/to/elle_grid”
export RAW="/path/to/elle_basic/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.
export RADAR_NAME="KHGX"
export days=20170827
export rdir="/path/to/lrose/bin"
echo $HOMED/$PROJ
  1. Convert the raw files to CfRadial files
$rdir/RadxConvert -params ./params/RadxConvert.* -f $RAW/$RADAR_NAME/$days/*.ar2v
  1. Regrid the data to a cartesian grid and apply the convective stratiform algorithm
$rdir/Radx2Grid -params ./params/RadxGrid.* -f convert/$RADAR_NAME/$days/*.nc


Expected output

After running these two programs, you should end up with NetCDF files that show the following output.

Grid tutorial ref leg.pngGrid tutorial ref field.png Grid tutorial cs leg.pngGrid tutorial cs field.png


Modify the gridding procedure in RadxGrid.new

  1. Select the interpolation mode in interp_mode on line 225. The current setting is a cartesian grid (INTERP_MODE_CART).
  1. In lines 227-312, there are various parameters for setting the interpolation the user might be interested in (e.g., whether to used the measured or fixed angles, setting the beamwidth fraction).
  1. Choose the vertical level spacing, either using evenly spaced levels (lines 354-358) or a custom array (386-400).
  1. Set the domain size and horizontal grid spacing in lines 439-446.
  1. Set the grid_projection on line 498 (the current setting is PROJ_FLAT).
  1. Choose whether to censor output based on polarimetric variables (e.g., PID category or Rhohv values) on line 1278 (current setting is TRUE and currently censors data where Rhohv is outside [0.8, 1.05] deg/km and REF is outside [-10, 100]).
  1. Choose whether to remove long range rays in NEXRAD data on line 1468 (current setting is TRUE).
  1. Choose whether to identify convective and stratiform precipitation and set the thresholds in lines 2416-2597 (current setting is TRUE). The threshold for definitive convection is on line 2491, the radius of convection is on line 2504, and choosing to write the partition to the file is on line 2559. However, I encourage suggest reading through the entire section and determining whether the values are appropriate for the precipitation type.

Radx2Grid censor.png