| https:///weradocs/wera16_beam.shtml |
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 .
| As of V2.0.0 (November 2009), the time series read from the .SORT file can be splitted into shorter overlapping fractions, which can be configured similar to wera16_grid_ol . Check the values in params.cfg as described below. Two different integration times can be handeled simultaneously: CUR and TSU. |
wera16_beam 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 params.cfg:
C-------------------------------------------------------------------------------
C Enable .SPEC/.rad files with different integration times
C-------------------------------------------------------------------------------
CUR_SPEC = .TRUE. !enable writing of current files (*.CUR.SPEC)
c CUR_SPEC = .FALSE. !enable writing of current files (*.CUR.SPEC)
TSU_SPEC = .TRUE. !enable writing of tsunami detection (*.TSU.SPEC)
c TSU_SPEC = .FALSE. !enable writing of tsunami detection (*.TSU.SPEC)
C-------------------------------------------------------------------------------
C Define integration times
C-------------------------------------------------------------------------------
! the following parameters must be multiples of 128
NSER_AVG = 512 ! length of combined .USORT fractions
NCUR_AVG = 2048 ! integration time for current measurements, minimum 512
NTSU_AVG = 512 ! Integration time for Tsunami detection, minimum 256
C-------------------------------------------------------------------------------
C Result file time steps within a coherent .USORT stream
C-------------------------------------------------------------------------------
! the following parameters must be multiples of 128
NCUR_STEP = 512 ! Samples between .CUR.SPEC / .crad result files
NTSU_STEP = 128 ! Samples between .TSU.SPEC / .crad Tsunami result files
C-------------------------------------------------------------------------------
C Name of WERA Gain File
C-------------------------------------------------------------------------------
GAIN_NAME = '/home/wera/etc/wera_gain.dat'
C-------------------------------------------------------------------------------
C The following 7 parameters allow simultaneous use of a linear and
C a square array within the same WERA system. They provide a mapping
C of receiver channels to antenna numbers. A value of '0' indicates that
C this antenna configuration does not exist in the actual WERA system.
C-------------------------------------------------------------------------------
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]
!Name of file for FullCalibration
FCAL_001 = '/home/wera/etc/LaRevellata.fullcal'
!Name of file for linear array spacing
DIST_001 = '/home/wera/etc/antdist_LaRevellata.asc'
!Name of file for antenna locations
APOS_001 = '/home/wera/etc/antpos_LaRevellata.asc'
ADJ_RAN_001 = 200.0 !Range to add to WERA measurenent [m]
ADJ_ANG_001 = 2.3 !Angle to add to WERA measurenent [degrees]
|
Look here for all options and their parameters to define the linear receive antenna array and beam forming.
When CUR_SPEC or TSU_SPEC are defined to be .TRUE., the resulting file name has to
include the seconds and integration time. The time code is generated from the starting time of the fraction
of the time series and the <path_to_.rad_file> is modified from e.g.
/home/wera/data/Hawaii/2009274/20092710115_kok.rad to
/home/wera/data/Hawaii/2009274/2009271011545-0512_kok.rad
i.e. the target directory, site code, and extension are kept.
Program expects 'FMSOCO' header code and writes 'FMRADB' header code.
Look here for a complete list of header codes.
The following box shows a typical output generated by wera16_beam:
$ wera16_beam /home/wera/data/Hawaii/2009274/20092710000_kok.SORT /home/wera/data/Hawaii/2009274/20092710000_kok.rad params_hawaii_hioos.cfg WERA16_Beam V2.0.0 Linux: *** Reading configuration from /home/wera/etc/params_hawaii_hioos.cfg *** The time is 27-NOV-2009 11:11:10 *** Now processing /home/wera/data/Hawaii/2009274/20092741500_kok.SORT nant_sort = 16 antennas Frequency = 16.0499992 MHz NSER = 2048 Samples T_chirp = 0.352083 s F_bragg = 0.408730 Hz Last RC = 100 range cells RANGEresol = 1.5000 km RAN_OFF = 0.5000 RC RAN_OFF = 750 m TRUE_NORTH = 310 degrees DDS_Clock = 181.48430769 MHz MT = 2080 samples SSR = 0.000000000000 seconds --- NSER_AVG = 512 --- NCUR_AVG = 2048 --- NTSU_AVG = 512 *** MAX_RAD_CUR is set to 3.50 m/s. *** VARX is set to ....... 1.50 *** Number of averages for current processing: 13 in 1 file. *** Number of averages for current processing: 1 in 13 files. *** Reading antenna locations from /home/wera/etc/antpos_kok.asc W16beam: 2048 SAMPLES 01-OCT-09 15:00 UTC kok FMSOCO FREQUENZ 16.050 *** FullCalib does not exist ....: /home/wera/etc/kok.fullcal *** Reading WERA_Gain ...........: /home/wera/etc/wera_gain.dat *** Calculating complex window function for beam forming. *** Calculating S/N of each antenna. *** Normalizing antenna gains. *** Forming beams: 0%...10%...20%...30%...40%...50%...60%...70%...80%...90%...100% *** Result in /home/wera/data/Hawaii/2009274/2009274150000-2048_kok.rad *** Result in /home/wera/data/Hawaii/2009274/2009274150000-0512_kok.rad *** Result in /home/wera/data/Hawaii/2009274/2009274150045-0512_kok.rad *** Result in /home/wera/data/Hawaii/2009274/2009274150130-0512_kok.rad *** Result in /home/wera/data/Hawaii/2009274/2009274150215-0512_kok.rad *** Result in /home/wera/data/Hawaii/2009274/2009274150300-0512_kok.rad *** Result in /home/wera/data/Hawaii/2009274/2009274150345-0512_kok.rad *** Result in /home/wera/data/Hawaii/2009274/2009274150430-0512_kok.rad *** Result in /home/wera/data/Hawaii/2009274/2009274150515-0512_kok.rad *** Result in /home/wera/data/Hawaii/2009274/2009274150600-0512_kok.rad *** Result in /home/wera/data/Hawaii/2009274/2009274150645-0512_kok.rad *** Result in /home/wera/data/Hawaii/2009274/2009274150730-0512_kok.rad *** Result in /home/wera/data/Hawaii/2009274/2009274150815-0512_kok.rad *** Result in /home/wera/data/Hawaii/2009274/2009274150900-0512_kok.rad *** The time is 27-NOV-2009 11:12:31 |
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_beam.F Compilation and linking with '. do_wera16_beam.sh' Call with > wera16_beam <path_to_.sort_file> <path_to_.rad_beam_file> [params_*.cfg] e.g. > wera16_beam /directory/20003341355_TOR.SORT /directory/20003341355_TOR.rad_beam params_hawaii.cfg You should use the same time code in the .SORT and .rad_beam filenames. |
VMS:
Source is in wera16_beam.F
Compilation and linking with @wera16_beam
Program needs to be defined as an external command to take parameters
from the command line, e.g.
$ WERA16_BEAM :== $$USER1:[WERA.Fortran]WERA16_BEAM.EXE
Call with
$ wera16_beam <path_to_.sort_file> <path_to_.rad_beam_file> [param_uv*.cfg] e.g.
$ WERA16_BEAM -
$USER1:[WERA.Fortran]20003341355_TOR.SORT -
$USER1:[WERA.Fortran]20003341355_TOR.rad_beam -
param_uv_hawaii.cfg
|
