resolvedvelocities is a Python implementation of the Heinselman and Nicolls Bayesian reconstruction algorithm [1] used to resolved 3D ion drift velocity and electric field vectors from AMISR line-of-sight measurements. There are two flavors of this algorithm, one which bins points by magnetic latitude and is typically used in the F-region to get local plasma convection velocities, and one which bins point by altitude and is typically used in the E-region to get altitude profiles of velocity. This package contains programs to run both.
This package requires numpy and apexpy and can be installed with pip:
pip install git+https://github.com/amisr/resolvedvelocities.git
This code is designed to be run from the command line with a configuration file. The configuration file specifies the input and output file names as well as a variety of other parameters.
Latitude-Binned Resolved Velocities (Vvels-Lat):
resolvedvelocities-lat config.ini
Altitude-Binned Resolved Velocities (Vvels-Alt):
resolvedvelocities-alt config.ini
Note:
The apexpy dependency MUST be installed from commit 94e63c3af524b60dd652eb1c92df786e7ac9076a. This version includes the function bvectors_apex
, which returns the value of Be3, necessary to convert the plasma drift velocity to the electric field.
Full documentation is available at https://resolvedvelocities.readthedocs.io/en/latest/.
Source code is available from the resolvedvelocities GitHub page. Code is open source under GNU General Public License v3.0. Please report bugs by submitting an issue!
[1] | Heinselman, C. J., and Nicolls, M. J. (2008), A Bayesian approach to electric field and E‐region neutral wind estimation with the Poker Flat Advanced Modular Incoherent Scatter Radar, Radio Sci., 43, RS5013, doi:10.1029/2007RS003805. |