WERA HF Radar Software - Spectra to Radial Currents / Radial Waves

wera16_spec_cwrad

A third optional parameter can be used to pass the name of a configuration file. This file must be located in ~/etc e.g. ~/etc/params_Hawaii.cfg . If no third parameter is supplied, the configuration is read from ~/etc/params.cfg .

The program works like wera16_spec_rad, except that it will create .crad and .wrad files, if requested. The extension on the result file (second parameter) must be .crad .

As of December 2007, this program can write NOAA .ruv compatible files, if in params_xxx.cfg NOAA_File = .TRUE.

Reads all information needed on the measurement from the header (.SPEC) or (.SPEC2) and from

Values in PARAM_UV.CFG / params.cfg: VARX = 1.5 !Variance of radial comp must be less [m/s] MAX_RAD_CUR = 1.0 !Maximum radial current allowed [m/s] STEER_001 = 60 !+/- max steering angle [degrees] ! C Parameters for WERA16_SPEC_CWrad / wrad2uv C WERA_Currents = .TRUE. !.true. = create .crad WERA_Waves = .TRUE. !.true. = create .wrad WERA_Wind = .TRUE. !.true. = request wrad2uv to calculate wind direction NOAA_File = .TRUE. !.true. = write NOAA .ruv file

If the water depth is known at the grid points (see data format of the topography file), the reduced phase speed of the Bragg-scattering waves can be taken care of when calculating the surface currend velocity.

MAX_RAD_CUR is used by wera16_spec_cwrad to limit the radial component of the current. (Grid points with higher current velocities are marked as 'gap'.)

Program writes binary (unformatted) data which are incompatible between Linux <--> VMS platforms. To transfer the data, first convert them to ASCII (Radial2ASCII), then use FTP in ASCII mode, and convert back to binary (ASCII2Radial).

Program expects  'FMSPEC' or 'FMSPC2' header code
and writes       'FMRADG'  header code.

Look here for a complete list of header codes.

The following box shows a typical output generated by wera16_spec_cwrad :

wera@ifmpc205:~/Fortran> wera16_spec_cwrad ~/data/surlitop/bre/spec/2005/246/20052460600_bre.SPEC ~/data/surlitop/bre/cwrad/2005/246/20052460600_bre.crad params_brest.cfg WERA16_Spec_CWRad V2.0.0 Linux: *** Reading configuration from /home/wera/etc/param_uv_brest.cfg *** The time is 13-DEC-2007 15:33:14 *** Now processing : ~/data/surlitop/bre/spec/2005/246/20052492348_bre.SPEC Frequency = 12.4580002 MHz NSER = 2048 Samples T_chirp = 0.260000 s F_bragg = 0.360100 Hz Last RC = 114 range cells RANGEresol = 1.5000 km RAN_OFF = 0.5000 RC RAN_OFF = 750 m TRUE_NORTH = 20 degrees Header: 2048 SAMPLES 06-SEP-05 23:48 UTC Brezellec FMRADG FREQUENZ 12.458 MAX_RAD_CUR = 1.50 m/s. Max Steer = +/- 60 degrees *** Set up table of grid points. *** Calculating Signal and Noise. *** Searching for 1st order Bragg peaks with sufficient S/N. *** Calculate statistics for Doppler peaks *** Surface current data in ~/data/surlitop/bre/cwrad/2005/246/20052492348_bre.crad *** Ocean wave data in ~/data/surlitop/bre/cwrad/2005/246/20052492348_bre.wrad *** NOAA radial data in ~/data/surlitop/bre/cwrad/2005/246/RDL_bre_2005_09_06_234800.ruv *** The time is 13-DEC-2007 15:33:19

Linux: Source is in wera16_spec_cwrad.F Compilation and linking with '. do_wera16_spec_cwrad.sh' Call with > wera16_spec_cwrad <path_to_.spec_file> <path_to_.rad_file> e.g. > wera16_spec_cwrad /directory/20003341355_TOR.SPEC /directory/20003341355_TOR.RAD_Beam You should use the same time code in the .SPEC and .RAD filenames and indicate that the result has been processed by Beam Forming by including '_Beam' in the file name.

Here are the data formats of the
input (.SPEC)
input (.SPEC2)
and
output (.RAD) files.

Compilation from source:

VMS: Source is in wera16_spec_cwrad.F Compilation and linking with @wera16_spec_cwrad Program needs to be defined as an external command to take parameters from the command line, e.g. $ wera16_spec_cwrad :== $$USER1:[WERA.Fortran]wera16_spec_cwrad.EXE Call with $ wera16_spec_cwrad <path_to_.spec_file> <path_to_.rad_file> e.g. e.g. $ wera16_spec_cwrad - $USER1:[WERA.Fortran]20003341355_TOR.SPEC - $USER1:[WERA.Fortran]20003341355_TOR.RAD_Beam You should use the same time code in the .SPEC and .RAD filenames and indicate that the result has been processed by Beam Forming by including '_Beam' in the file name.