To unclutter PPIs:

First run the script that creates tmp files (there is a bug in code and this is the only way to create a tmp directory):
createTmpFiles.pl inputdir outputdir
(the full path is /home/disk/blitz/bin/createTmpFiles.pl)

Then on a 32-bit machine: rainband, clutter,dbz,eyewall,etc (except beam for some reason) run:

rmClutter_ppi_cluttermap_32
(The full path is /home/disk/blitz/bin or /home/disk/blitz/src/swp-src/)

This code uses a clutter map and the clutter product of the particle identification (PID) algorithm to remove ground clutter. If the PID variables are not available, as it is often the case for PPI, then only the clutter map is used.

Usage: rmClutter_ppi_cluttermap_32 [swpFlist] [reflName] [velName] [pdName] [wpdName] [clrAirFlist] [clrReflName] [reflThres] [elevThres] [outputdir]

swpFlist = file that contains a list of the sweep files to unclutter

reflName = name of the reflectivity variable, normally DZ (use 'readswp_headers' to get info on the headers of a sweepfile)

velName = name of the radial velocity variable, normally VE or VR (use 'readswp_headers' to get info on the headers of a sweepfile)

pdName and wpdName = names of particle type fields, which are often not available for PPI, so set them to '-'

clrAirFlist = file that contains the name or names of the sweepfile that will be use as a clutter map. The first line in this file needs to have the total number of clutter map files (normally one if all the PPI correspond to the same elevation angle)

clrReflName = name of the reflectivity variable in the clutter map file (normally DZ, use 'readswp_headers' to get info on the headers of a sweepfile

reflThres = This is the most important papameter of the algorithm: radar gates that have relfectivity values below this threshold in the clutter map file will be removed in the uncluttered fields.

elevThres = This is a parameter that is useful for removing clutter in RHI, not in PPI. It indicates that above this elevation, no clutter will be removed. For PPI, it should be set to zero.

outputdir = Output directory

The output files has both original and uncluttered reflectivity and radial velocity. The the name of the uncluttered variables strat with a 'C' for 'corrected'.

The only files that get uncluttered are the ones that match the clutter map in number of range bins, if they are different, the files don't get uncluttered (i.e., will need a different clutter map for long range and short range swps. Again 'readswp_headers' is useful to know the number of gates (numbins) in each file).