### Navigation

- index
- modules |
- modules |
- homepage|
- search|
- documentation »
- Module code »

# Source code for spacepy.irbempy

    # -*- coding: utf-8 -*-

    """
    module wrapper for irbem_lib
    Reference for this library
    https://sourceforge.net/projects/irbem/
    D. Boscher, S. Bourdarie, P. O'Brien, T. Guild, IRBEM library V4.3, 2004-2008

    Most functions in this module use an options list to define the models used and the settings
    that define the quality level of the result. The options list is a 5-element list and is defined
    as follows.

    Options
    -------
       - 1st element: 0 - don't compute L* or phi ;  1 - compute L*; 2- compute phi
       - 2nd element: 0 - initialize IGRF field once per year (year.5);
           n - n is the  frequency (in days) starting on January 1st of each year
           (i.e. if options(2nd element)=15 then IGRF will be updated on the following days of the
           year: 1, 15, 30, 45 ...)
       - 3rd element: resolution to compute L* (0 to 9) where 0 is the recomended value to ensure a
           good ratio precision/computation time
           (i.e. an error of ~2% at L=6). The higher the value the better will be the precision, the
           longer will be the computing time. Generally there is not much improvement for values
           larger than 4. Note that this parameter defines the integration step (theta)
           along the field line such as dtheta=(2pi)/(720*[options(3rd element)+1])
       - 4th element: resolution to compute L* (0 to 9). The higher the value the better will be
           the precision, the longer will be
           the computing time. It is recommended to use 0 (usually sufficient) unless L* is not
           computed on a LEO orbit. For LEO orbit higher values are recommended. Note that this
           parameter defines the integration step (phi) along the drift shell such as
           dphi=(2pi)/(25*[options(4th element)+1])
       - 5th element: allows to select an internal magnetic field model (default is set to IGRF)
           - 0 = IGRF
           - 1 = Eccentric tilted dipole
           - 2 = Jensen&Cain 1960
           - 3 = GSFC 12/66 updated to 1970
           - 4 = User-defined model (Default: Centred dipole + uniform [Dungey open model] )
           - 5 = Centred dipole

    The routines also require specification of the external magnetic field model. The default is the
    Tsyganenko 2001 storm-time model. The external model is always specified using the extMag keyword
    and the following options exist.

    extMag
    ------
        - '0' = No external field model
        - 'MEAD' = Mead and Fairfield
        - 'T87SHORT' = Tsyganenko 1987 short (inner magnetosphere)
        - 'T87LONG' = Tsyganenko 1987 long (valid in extended tail region)
        - 'T89' = Tsyganenko 1989
        - 'OPQUIET' = Olsen-Pfitzer static model for quiet conditions
        - 'OPDYN' = Olsen-Pfitzer static model for active conditions
        - 'T96' = Tsyganenko 1996
        - 'OSTA' = Ostapenko and Maltsev
        - 'T01QUIET' = Tsyganenko 2001 model for quiet conditions
        - 'T01STORM' = Tsyganenko 2001 model for active conditions
        - 'T05' = Tsyganenko and Sitnov 2005 model
        - 'ALEX' = Alexeev model

    Authors
    -------
    Josef Koller, Steve Morley

    Copyright 2010 Los Alamos National Security, LLC.
    """

    from .irbempy import *

### Navigation

- index
- modules |
- modules |
- homepage|
- search|
- documentation »
- Module code »