NAME

qc_mom_ww - do quality control on moments via the ETL version, which uses the Weber/Wuertz algorithm.

SYNOPSIS

qc_mom_ww -s startdate -e enddate [-t qcwindow] [-T slidetime]

[-q qclimit] [-x extrapolate] [-z filesize] [-i inputdir] [-o outputdir]

[-y verticaldy] [-n verticalsnr] [-m verticalnmin] [-d verticaldxt]

[-h verticaldxh]

[-Y obliquedy] [-N obliquesnr] [-M obliquenmin] [-D obliquedxt]

[-H obliquedxh]

AVAILABLILTY

this program is in /net/sssf2/profiler/analysis/bin

DESCRIPTION

For all data found in input moment netcdf files between startdate and enddate qc_mom_ww computes quality control values for all velocities, interpolates new velocities for velocities marked as bad (qc values above qclimit), and computes an error estimate on the new velocities. It updates the appropriate variables in the moment netcdf file and writes out a new netcdf file. qc_mom_ww looks for the original moment netcdf files in the directory specified in the INPUT_DIR environmental variable, and writes the new moments netcdf files to the directory specified by the OUTPUT_DIR environmental variable. Files are only created for the time period in which data is found. The files are all hourly files.

REQUIRED ARGUMENTS

-s startdate

 

This is not an optional option. It must be provided. It gives the date and time to start collecting data from. It must be given in the following format: "yyyy-mm-dd hh:nn:ss" where yyyy=year, mm=month, dd=day, hh=hour, nn=minutes, ss=seconds.

 

-e enddate

 

This is not an optional option. It must be provided. It gives the date and time to stop collecting data. It must be given in the following format: "yyyy-mm-dd hh:nn:ss" where yyyy=year, mm=month, dd=day, hh=hour, nn=minutes, ss=seconds.

 

OPTIONS

-t qcwindow

 

The size of the window each instance of QC looks at to find patterns, etc. This must be given in the format "hh:mm:ss" where hh=hours, mm=minutes, ss=seconds. The default is: "01:00:00"

 

-T slidetime

 

How far forward to slide the Qc Window specified by QC Time. This must be given in the format "hh:mm:ss" where hh=hours, mm=minutes, ss=seconds. The default is: "01:00:00"

-q qclimit

 

A Qc value above which values are considered "bad" (it will try to interpolate at these values). default: 10

 

-x extrapolate

 

Allow Extrapolation(1) or not(0). default: 0

 

-z filesize

 

NOTE: this option has not been implemented yet. the amount of time to put into each file. The default is 1 hour. This must be entered in the format: hh:mm:ss where hh=hour, mm=minutes, and ss=seconds (Example: -z "01:00:00")

 

-i inputdir

 

specifies the directory in which to find the .mom files to use for calculating winds. If this option is provided it is used. If not, the program looks for the environmental variable INPUT_DIR. If that is not found then it looks in the current working directory.

 

-o outputdir

 

specifies the directory in which to put the newly created .wnd files. If this option is provided it is used. If not, the program looks for the environmental variable OUTPUT_DIR. If that is not found then it puts the new files in the current working directory.

 

The remaining options are all parameters for the Weber/Wuertz QC. Use with caution.

 

Parameters for vertical wind beams only:

-y verticaldy

 

the allowable change over the neighborhood(dx) (m/s) for the (vertical Wind Beams). default: 1.2

 

-n verticalsnr

 

The lower limit for snr. Data below this limit are flagged with a value of 99 in the error_estimate (vertical Wind Beams). default: -30

 

-m verticalnmin

 

Goes into the equation: x=n/verticalnmin; where x = minimum number of required for confidence (vertical Wind Beams). default: 10

 

-d verticaldxt

 

Neighborhood size in the time dimension. Goes into equation:

x=(verticaldxt) * (minimum time diff)

(vertical Wind Beams). default: 2.1

 

-h verticaldxh

 

neighborhood size in the height dimension. Goes into equation:

x = (verticaldxh) * (minimum height diff)

(vertical Wind Beams) default: 2.1

 

Parameters for oblique wind beams only:

-Y obliquedy

 

the allowable change over the neighborhood(dx) (m/s) for the (Oblique Wind Beams). default: 1.2

 

-N obliquesnr

 

The lower limit for snr. Data below this limit are flagged with a value of 99 in the error_estimate (Oblique Wind Beams). default: -30

 

-M obliquenmin

 

Goes into the equation: x=n/obliquenmin; where x = minimum number of required for confidence (Oblique Wind Beams). default: 10

 

-D obliquedxt

 

Neighborhood size in the time dimension. Goes into equation:

x=(obliquedxt) * (minimum time diff)

(Oblique Wind Beams). default: 2.1

 

-H obliquedxh

 

neighborhood size in the height dimension. Goes into equation:

x = (obliquedxh) * (minimum height diff)

(Oblique Wind Beams) default: 2.1

 

ENVIRONMENT

The INPUT_DIR and OUTPUT_DIR variables must be set.

FILES

input moment netcdf files of the correct format must be found in the INPUT_DIR directory or the directory specified by the -i option.

SEE ALSO

profiler_to_netcdf, spc_to_mom_hs, mom_to_wind_etl

NOTES

BUGS

There is currently no error checking. If the wrong format, or the input/output directory do not exist a core dump will occor.