Inclined_sweeps_HAARP

Get source data

Zenodo record

Shindin, Alexey. (2019). Properties of Stimulated Electromagnetic Emissions during Inclined HF Pumping of the Ionosphere near the 4th Electron Gyroharmonic at HAARP - data set (Version 1.0.0) [Data set]. Zenodo. http://doi.org/10.5281/zenodo.2600237

File "raw_data_archive.urls" contains all urls for the data set above.

"Run all" bash script with all commands below for Linux OS to do all data processing and make all figures

Scipt below tested in Ubuntu 18.04 with Anaconda 2018.12 python3 distribution on PC (CPU: AMD Ryzen 7 1800X, RAM: 32 GB). On such PC about 30 GB of free disc space is needed and about 48 hours to do all processing.

bash run_all.sh

Download the data

i.e. using wget:

wget -i raw_data_archive.urls

Check the data

In Linux:

sed 's/\r$//' checksums.sha | sed 's/\r/\n'/ | sha1sum -c

In Windows one can use Double commander to do this.

Unpack the data

Third-party software needed: 7zip for '*.7z' files, wavpack for '*.wv' files.

wvunpack -r *.wv
7z -y x "*.7z"

Reproducing the figures

Setup environment

All code tested with Anaconda 2018.12 python distribution:

Python 3.7.1
numpy-1.15.4
pandas-0.23.4
matplotlib-3.0.2
scipy-1.1.0

To extract RF data from '*.DRXS12' files GNU Octave (or Matlab) is requied. Octave code tested with GNU Octave 5.1.0.

Extract RF Data from '*.DRXS12'

octave-cli extract_spew.m

Get PSD spectrograms

python get_spectrograms_A.py
python get_spectrograms_B.py
python get_spectrograms_C.py

Concatenate spectrograms to get one spectrogram per site

python concat_spectrograms.py

Process spectrograms to get f0-shifted spectrograms

python inline_sweeps.py

Create database of spectrograms

python fill_database.py

Process DM SEE feature

python dm_proc.py

Process BUM SEE feature

python bum_proc.py

Figure 1

python figure1.py

Figure 2

python figure2.py

Figure 3

python figure3.py

ashindin/Inclined_sweeps_HAARP