plasmaprm

#!/usr/bin/env python

from irfpy.util.formula import *

def main():

    print('wpe [/s] = %.3e * sqrt(n [/cc])'%plasma_afreq_el(1))
    print('fpe [/s] = %.3e * sqrt(n [/cc])'%(plasma_freq_el(1)))

    print('wpi [/s] = %.3e * sqrt(n [/cc])'%plasma_afreq_pr(1))
    print('fpi [/s] = %.3e * sqrt(n [/cc])'%plasma_freq_pr(1))

    print('wge [/s] = %.3e * B [nT]' % gyro_afreq_el(1))
    print('fge [/s] = %.3e * B [nT]' % gyro_freq_el(1))

    print('wgi [/s] = %.3e * B [nT]' % gyro_afreq_pr(1))
    print('fgi [/s] = %.3e * B [nT]' % gyro_freq_pr(1))

    print('Tge [s] = %.3e / B [nT]' % gyro_period_el(1))
    print('Tgi [s] = %.3e / B [nT]' % gyro_period_pr(1))

    print('Rge [km] = %.3e * v [km/s] / B [nT]' % gyro_radius_el(1, 1))
    print('Rgi [km] = %.3e * v [km/s] / B [nT]' % gyro_radius_pr(1, 1))

    print('Lde [km] = %.3e * sqrt(Te [eV] / n [/cc])' % debye_length_el(1,1 ))

if __name__ == '__main__':
    main()