ONERA-DESP LIBRARY
Modification history


V2.4: April,12 2005:
trace_field_line: A new routine has been added. It allows to trace magnetic field lines. Full description is available in the user guide.

make_lstar and make_lstar_shell_splitting: a minor bug has been fixed -> when any spacecraft was at the magnetic equator within double real precision the two routines were returning wrong number (sometimes NaN or L* between ~0.1 to ~0.8). Note that make_lstar_shell_splitting was returning bad value only for 90 degree pitch angle.

make_lstar: another minor bug has been fixed -> when a serie of input postion was submitted to the routine, the first value of L* returned after a data gap could have been set to bad data flag.

Download onera_desp_lib_2.4.tar.gz


V2.5: June,7 2005:
IGRF: IGRF 2005 has been implemented.

Tsyganenko 2004 storm model has been added

Alexeev (paraboloid model - Submitted to ISO) has been added.

Download onera_desp_lib_2.5.tar.gz


V3.0: June,1st 2006:
Geographic to heliospheric coordinate transformations have been implemented.

Heliospheric coordinate transformations have been added

NASA AE8/AP8 and AFRL CRRESELE/CRRESPRO have been added.

Download onera_desp_lib_3.0.tar.gz


V3.1: July,12 2006:
Geographic to heliospheric coordinate transformations have been implemented.

Bug fix in Heliospheric coordinate transformations

Easier to use AFRL CRRESELE/CRRESPRO function.

MATLAB wrappers have been added (see new User's guide).

Download onera_desp_lib_3.1.tar.gz

V3.2: December,15 2006:
In this new release bugs have been fixed in the Matlab wrappers only. This upgrade resolves fatal errors experienced only by some users. To resolve incompatibilities with older versions of Matlab, we have made a minor change to how onera_desp_lib_rotate calls anonymous functions. Also, to accommodate Matlab's unreliable handling of spaces in path names, we have introduced a new "prototype" function onera_desp_lib_proto.m that provides Matlab the same information that used to be provided by the C-style header file onera_desp_lib.h. This approach avoids the need for Matlab to invoke a C compiler to read the .h file. The code for this prototype function was automatically built by Matlab's loadlibrary routine, so it is not documented. The new function, onera_desp_lib_proto, is used by onera_desp_lib_load for loading the onear_desp_lib shared library file. Most users will never explicitly use either onera_desp_lib_load or onera_desp_lib_proto, since these functions are called automatically by the more usual library routines.
Bugs in the matlab wrappers have been fixed (this affected all matlab distributions). This latter problem caused calls to onera_desp_lib_rotate to not work properly when called with a time series of points to be rotated.
None of these changes alters the calling syntax of the main library routines. The modified and new .m files should be installed in the same directory as the library's other .m files. No change to the FORTRAN or shared library files was needed.

Download onera_desp_lib_3.2.tar.gz

V4.1: April,16 2007:
In this new release bugs have been fixed, the list is given below:

- In GET_MLT, xGEO was corrupted when it was returned by the library. This is fixed.

- In FLY_IN_AFRL_CRRES the returned electron fluxes were wrong by only a factor 1000! (was per keV rather than per MeV as it was set in the manuel). Now the returned numbers are correct and the flux are per MeV.

- In FLY_IN_AFRL_CRRES various coordinates system were allowed as input but because coordinate rotations are time dependent when GSM,GSE,SM or GEI are involved this was not handle propoerly by the routine. This is now fixed but time when input position is provided is now required. Look carefully at the new manual to track the new call sequence.

- In FLY_IN_AEAP various coordinates system were allowed as input but because coordinate rotations are time dependent when GSM,GSE,SM or GEI are involved this was not handle propoerly by the routine. This is now fixed but time when input position is provided is now required. Look carefully at the new manual to track the new call sequence.

New features:

- In FLY_IN_AFRL_CRRES it is now possible to compute fluxes for several energies at once (up to 25). This new feature impose a different calling sequence. Check with the new manual to track new calling sequence.

- In FLY_IN_AEAP it is now possible to compute fluxes for several energies at once (up to 25). This new feature impose a different calling sequence. Check with the new manual to track new calling sequence.

- COORD_TRANS : this is a new generic and easy to use routine to perform any coordinate rotations. More details are available in the manual.

- COORD_TRANS_VEC : same as COORD_TRANS but handle up to 100000 rotations in one call. More details are available in the manual.

- Two orbit propagators have been added:

SGP4_TLE allows to propagate NORAD two lines elements.

SGP4_ELE allows to propagate any orbit from different type of elements.

RV2COE allows to find the classical orbital elements given the Geocentric Equatorial Position and Velocity vectors.

- A Makefile is now provided to build the library on various platforms. It is highly recommended to use the make command now because in the future the scripts provided in the compile directory will not be maintained any more.


Download onera_desp_lib_4.1.tar.gz



V4.2: May,28 2008:
From this new version, the library is distributed under LGPL licence (see http://www.gnu.org/licenses/)
In this new release bugs have been fixed in the matlab wrappers and in the fortran source codes. There were no major bugs only minor ones, so they are not listed here.

New features:

- In MAKE_LSTAR and MAKE_LSTAR_SHELL_SPLITTING, there is a switch to compute Φ rather than L*.

- GET_HEMI : this routine allows to find in which magnetic hemisphere a location is.

- LSTAR_PHI : this routine allows one to convert from L* to Φ and vice versa.

- GET_AE8_AP8_FLUX : This function allows one to compute NASA AE8 min/max and AP8 min/max flux for any B/Bo, L position. The output can be differential flux or flux within an energy range or integral flux.

- GET_CRRES_FLUX : This routine allows one to compute AFRL CRRESPRO and CRRESELE flux for any B/Bo, L position. The output can be differential flux or flux within an energy range or integral flux.

- FLY_IN_IGE : This routine allows one to fly any geostationary spacecraft in IGE (International Geostationary Electron) models. The use of the model is limited to geostationary altitude as it is based on LANL-GEO bird series (1976 to 2006) and JAXA-DRTS spacecraft (added in IGE2006). Three versions of the model are provided.

- Atmospheric models have been added: MSIS86, MSISE90, NRLMSIS00



Download onera_desp_lib_4.2.tar.gz