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

WERA HF Radar Software - .SORT Data to Radial Current Components (square array)

wera4_dir

Call: wera4_dir <path_to_.sort_file> <path_to_.rad_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):

Please check for correct setting of the parameters in params.cfg . 

Values in PARAM_UV.CFG:

	MIN_CONTRIB     =  3    !Direction Finding, filter parameter
	NUM_TRIES       =  4    !Direction Finding, filter parameter
	VARX            =  1.5  !Variance of radial comp must be less [m/s]
	MAX_RAD_CUR     =  2.5  !Maximum 2-dim current allowed [m/s]
	WE4_SN		=  6.0  !min. S/N wera4_dir (Factor, not dB)

C       Name and size (NX, NY) of Grid File
C
        NX              = 20
        NY              = 20
        GRID_NAME       = '$user1:[WERA.Fortran]Maasmond_grid_20_20.asc'

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

C       Name of WERA Phase File (required for square array direction finding)
C
        PHASE_NAME      = '$user1:[WERA.Fortran]PHAS_THEO_360_180deg.asc'

        !The following 6 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
        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        ='$data1:[WERA.Fortran]Ouddorp.fullcal'

Values in params.cfg:

	MIN_CONTRIB     =  3    !Direction Finding, filter parameter
	NUM_TRIES       =  4    !Direction Finding, filter parameter
	VARX            =  1.5  !Variance of radial comp must be less [m/s]
	MAX_RAD_CUR     =  2.5  !Maximum 2-dim current allowed [m/s]
	WE4_SN		=  6.0  !min. S/N wera4_dir (Factor, not dB)

C       Name and size (NX, NY) of Grid File
C
        NX              = 20
        NY              = 20
        GRID_NAME       = '/home/wera/etc/maasmond_grid_20_20.asc'

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

C       Name of WERA Phase File (required for square array direction finding)
C
        PHASE_NAME      = '/home/wera/etc/PHAS_THEO_360_180deg.asc'

        !The following 6 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
        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        ='/home/wera/etc/Ouddorp.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 wera4_dir : 

wera@ifmwera2:~/Fortran> wera4_dir /home/wera/data/examples/20032202220_rav.SORT /home/wera/data/examples/20032202220_rav.RAD_dsq params_Adria.cfg
 WERA4_Dir 1.1.1 Linux:
 *** Now processing /home/wera/data/examples/20032202220_rav.SORT
 NFFT       =   2048 Samples
 F_BRAGG    =  0.4074547 Hz
 Last RC    =  96
 RANGEresol =  1.5000 km
 RAN_OFF    =  0.50 RC
 TRUE_NORTH =  163 degrees
 Square array antennas:   1  2  3  4
 Number of frac-series =  193
 W4dir: 2048 SAMPLES  08-AUG-03 22:20 UTC  Ravenna     FMSOCO FREQUENZ  15.950
 *** FullCalib does not exist ....: /home/wera/etc/Ravenna.fullcal
 *** Reading WERA_Gain ...........: /home/wera/etc/wera_gain.dat
 *** Reading WERA_Phase: /home/wera/etc/PHAS_THEO_360_180deg.asc
 *** MAX_RAD_CUR is set to    1.50 m/s.
 *** WE4_SN is set to .....   6.00
 *** VARX is set to .......   1.50
 *** MIN_CONTRIB is set to    3
 *** NUM_TRIES is set to ..   4
 *** The time is  25-JUN-2004  15:12:35
 *** Result in /home/wera/data/examples/20032202220_rav.RAD_dsq
 *** The time is  25-JUN-2004  15:12:35

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

Compilation from source: 
Linux:

Source is in wera4_dir.F
Compilation and linking with '. do_wera4_dir.sh'

Call with
> wera4_dir <path_to_.sort_file> <path_to_.spec_file> params*.cfg  e.g.
> wera4_dir /directory/20003341355_oud.sort /directory/20003341355_oud.rad_dsqu params_Maasmond.cfg

You should use the same time code in the .sort and .rad_dsqu filenames.


VMS:

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

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

Call with
$ wera4_dir <path_to_.sort_file> <path_to_.spec_file> param_uv*.cfg e.g.
$ WERA4_DIR -
         $USER1:[WERA.Fortran]20003341355_oud.sort -
         $USER1:[WERA.Fortran]20003341355_oud.rad_dsqu -
         param_uv_Maasmond.cfg

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