SuperDARN / pyDARNio

Python Library for read in SuperDARN data
GNU Lesser General Public License v3.0
8 stars 2 forks source link

Multi-beam Record Conversion #34

Closed RemingtonRohel closed 2 years ago

RemingtonRohel commented 2 years ago

Extended conversion to dmap dictionaries to allow for multiple beams in a record.

Scope

Converting HDF5 files to DMAP format.

Issue: Borealis files with multiple beams in a record only have the last beam per record converted.

Approval

Number of approvals: 1

Test

import pydarnio

infile = "path to borealis file with multi-beam records"
file_type = "rawacf" or "bfiq"
file_structure = "site" or "array"
sdarn_file = "path to output sdarn file"

converter = pydarnio.BorealisConvert(infile, file_type, sdarn_file, borealis_file_structure=file_structure)

After this, you can try plotting the data for all beams to verify that the data has been correctly converted to DMAP. See "Plotting Lag-0" in the pyDARN documentation for how to plot (if your data is rawacf): https://pydarn.readthedocs.io/en/main/user/range_time/

RemingtonRohel commented 2 years ago

Tested and was able to plot all beams from iqdat and rawacf dmap files.

kevinkrieger commented 2 years ago

Tested with /home/dataman/widebeam_3tx/20211213/multi_record/20211213.2034.52.pgr.0.rawacf.hdf5.site.avg16 on maxwell. Used the script shown above then plotted all beams. Here's three of the plots:

bm0_20211213 2034 52 pgr 0 rawacf avg16 bm7_20211213 2034 52 pgr 0 rawacf avg16 bm15_20211213 2034 52 pgr 0 rawacf avg16