https:///weradocs/wera16_dir.shtml Last update:

WERA HF Radar Software - .SORT Data to Radial Component of Current (Direction Finding)

wera16_dir

Call: wera16_dir <path_to_.sort_file> <path_to_.spec_file> [params*.cfg]

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 .

Reads all information needed on number of channels, length of timeseries, etc. from the header (first 512 bytes of data file) and the following values from PARAM_UV.CFG (VMS) / params.cfg (Linux): 

Values in PARAM_UV.CFG:

C       Name of WERA Gain File
C
        GAIN_NAME       = '$user1:[WERA.Fortran]WERA_Gain.dat'

        !The following 7 parameters allow simultaneous use of a linear and
        !a square array within the same WERA system. They provide a mapping
        !of receiver channels to antenna numbers. A value of '0' indicates that
        !this antenna configuration does not exist in the actual WERA system.
        !
        ARR1_001        =  1            !First antenna# in linear array
        ARR2_001        = 12            !Last antenna# in linear array
        ARRD_001        = 9.35          !Distance between elements [meters]
	STEER_001	= 60		!+/- max steering angle [degrees]
        SQA1_001        = 15            !Antenna 1 of square array
        SQA2_001        = 16            !Antenna 2 of square array
        SQA3_001        = 14            !Antenna 3 of square array
        SQA4_001        = 13            !Antenna 4 of square array
                                        !Name of file for FullCalibration
        FCAL_001        = '$user1:[WERA.Fortran]Koolina.fullcal'


Values in params.cfg:

C       Name of WERA Gain File
C
        GAIN_NAME       = '/home/wera/etc/wera_gain.dat'

        !The following 7 parameters allow simultaneous use of a linear and
        !a square array within the same WERA system. They provide a mapping
        !of receiver channels to antenna numbers. A value of '0' indicates that
        !this antenna configuration does not exist in the actual WERA system.
        !
        ARR1_001        =  1            !First antenna# in linear array
        ARR2_001        = 16            !Last antenna# in linear array
        ARRD_001        = 9.35          !Distance between elements [meters]
	STEER_001	= 60		!+/- max steering angle [degrees]
        SQA1_001        =  0            !Antenna 1 of square array
        SQA2_001        =  0            !Antenna 2 of square array
        SQA3_001        =  0            !Antenna 3 of square array
        SQA4_001        =  0            !Antenna 4 of square array

                                        !Name of file for FullCalibration
        FCAL_001        = '/home/wera/etc/Koolina.fullcal'


Program expects  'FMSOCO'  header code
and writes       'FMRADD'  header code.

Look here for a complete list of header codes.

The following box shows a typical output generated by wera16_dir: 

wera@ifmwera2:~/Fortran> wera16_dir ~/data/examples/20031541930_kna.SORT ~/data/examples/20031541930_kna.RAD_beam
 wera16_dir  V1.1.5 Linux:
 *** The time is  25-JUN-2004  13:11:42
 *** MAX_RAD_CUR is set to    2.50 m/s.
 *** VARX is set to .......   1.50
 *** Now processing /home/wera/data/examples/20031541930_kna.SORT
 W16dir: 1024 SAMPLES  03-JUN-03 19:30 UTC  Kaena       FMSOCO FREQUENZ  16.050
 Frequency  = 16.0499992 MHz
 NFFT       =   1024 Samples
 F_BRAGG    =  0.4087299 Hz
 Last RC    = 128
 RANGEresol =  1.5000 km
 RAN_OFF    =  0.50 RC
 TRUE_NORTH =  299 degrees
 Linear array: Antennas  1 to 16
 AntSpacing =  0.5002 lambda
 *** FullCalib does not exist ....: ~/etc/Site.fullcal
 *** Reading WERA_Gain ...........: /home/wera/etc/wera_gain.dat
 *** Reading complex beam window from /home/wera/etc/window_new.asc
 *** Calculating S/N of each antenna.
 *** Normalizing antenna gains.
 *** Forming beams.
 0%...10%...20%...30%...40%...50%...60%...70%...80%...90%...100
 *** Search for 1st order peaks with sufficient S/N.
 *** Result in /home/wera/data/examples/20031541930_kna.RAD_beam
 *** The time is  25-JUN-2004  13:11:57

This is an example of radial components of the surface current field measured at Kaena/Oahu/Hawaii visualized by Plott_UV.

Here are the data formats of the input (.SORT) and output (.RAD_beam) files.

Compilation from source: 
Linux:

Source is in wera16_dir.F
Compilation and linking with '. do_wera16_dir.sh'

Call with
> wera16_dir <path_to_.sort_file> <path_to_.rad_beam_file> [params*.cfg]  e.g.
> wera16_dir /directory/20003341355_TOR.SORT /directory/20003341355_TOR.rad_dir params_Gijon.cfg

You should use the same time code in the .SORT and .rad_beam filenames.


VMS:

Source is in wera16_dir.F
Compilation and linking with @wera16_dir

Program needs to be defined as an external command to take parameters
from the command line, e.g.
$ WERA16_DIR :== $$USER1:[WERA.Fortran]WERA16_DIR.EXE

Call with
$ wera16_dir <path_to_.sort_file> <path_to_.rad_beam_file> param_uv*.cfg e.g.
$ WERA16_DIR -
         $USER1:[WERA.Fortran]20003341355_TOR.SORT -
         $USER1:[WERA.Fortran]20003341355_TOR.rad_beam -
         param_uv_Gijon.cfg

Comments to: gurgel@ifm.uni-hamburg.de (Klaus-Werner Gurgel)