| https:///weradocs/Data_Format_WERA_Spec2.shtml | Last update: |
|
This is the data format of WERA Power Spectra measured on a rectangular grid.
The grid has the size NX * NY (East * South). A list of coordinates is
stored in Gijon_grid_60_60.asc.
This information is also stored in the
header of each data file.
The structure in each data file is as follows: | | Table of pointers | IX=PX(1) | IX=PX(2) |...| IX=PX(NCOV) | | HEAD | |----------|----------|...|-------------| | | PX,PY(1..NCOV) | IY=PY(1) | IY=PY(2) |...| IY=PY(NCOV) |The table of pointers is read from NCOV/64+1 512 byte blocks with the numbers 2 to (NCOV/64 + 2) in the form (PX(i),PY(i)), i=1,NCOV .
integer*4 NCOV
integer*4 list_of_locations(128) ! 512 bytes = 1 block
integer*4 icnt ! counter for blocks
integer*4 PX(NX*NY) ! Table of pointers IX
integer*4 PY(NX*NY) ! Table of pointers IY
do icnt = 1,NCOV/64+1 ! number of 512 byte blocks
read(iun_spec,REC=1 + icnt) (list_of_locations(i), i=1,128) ! read 512 byte block
do i = 1,64 ! 64 pointers (IX,IY) in one block
if(i+64*(icnt-1).le.NCOV) then
PX(i+64*(icnt-1)) = list_of_locations(2*i-1)
PY(i+64*(icnt-1)) = list_of_locations(2*i )
end if
end do
end do
The powerspectrum in each grid point consists of NFTL values from negative
to positive doppler shifts with DC at position 1+NFTL/2
in IEEE REAL*4 floating point format (4 * NFTL bytes).
File_Size [Bytes] = 512 + (NCOV/64+1) * 512 + NCOV * 4 * NFTL Download 2008077133000_syl_CUR.SPEC as an example data set. |