Procedure: tt01 Purpose: tplot wrapper for the functional interface to the IDL Geopack implementation of the Tsyganenko 2001 and IGRF model. Input: pos_gsm_tvar: the tplot variable storing the position in gsm coordinates Keywords: pdyn(optional): Solar wind pressure(nanoPascals) should either be a string naming a tplot variable or an array or a single value. If a tplot input is used it will be interpolated to match the time inputs from the position var. Non-tplot array values must match the number of times in the tplot input for pos_gsm_tvar dsti(optional): DST index(nanoTeslas) should either be a string naming a tplot variable or an array or a single value. If a tplot input is used it will be interpolated to match the time inputs from the position var. Non-tplot array values must match the number of times in the tplot input for pos_gsm_tvar yimf(optional): y component of the interplanetary magnetic field should either be a string naming a tplot variable or an array or a single value. If a tplot input is used it will be interpolated to match the time inputs from the position var. Non-tplot array values must match the number of times in the tplot input for pos_gsm_tvar zimf(optional): z component of the interplanetary magnetic field should either be a string naming a tplot variable or an array or a single value. If a tplot input is used it will be interpolated to match the time inputs from the position var. Non-tplot array values must match the number of times in the tplot input for pos_gsm_tvar g1(optional): index describes solar wind conditions in the previous hour, should either be a string naming a tplot variable or an array or a single value. If a tplot input is used it will be interpolated to match the time inputs from the position var. Non-tplot array values must match the number of times in the tplot input for pos_gsm_tvar g2(optional): index describes solar wind conditions in the previous hour should either be a string naming a tplot variable or an array or a single value. If a tplot input is used it will be interpolated to match the time inputs from the position var. Non-tplot array values must match the number of times in the tplot input for pos_gsm_tvar parmod(optional): can input the Nx10 parmod array used by the Fortran Tsyganenko model. If passed as a raw array it will not be modified or interpolated so be sure its has the correct number of entries. It can also be passed as a tplot variable name in which case it will be interpolated. If values are passed individually and as par the par values will be overwritten. period(optional): the amount of time between recalculations of geodipole tilt in seconds(default: 60) increase this value to decrease run time get_nperiod(optional): Return the number of periods used in the time interval newname(optional):the name of the output variable. (default: pos_gsm_tvar+'_bt01') This option is ignored if globbing is used. error(optional): named variable in which to return the error state of this procedure call. 1 = success, 0 = failure get_tilt(optional): Set this value to a tplot variable name in which the geodipole tilt for each period will be returned. One sample will be returned for each period with time at the center of the period. set_tilt(optional): Set this to a tplot variable name or an array of values containing the dipole tilt that should be used. If a tplot input is used it will be interpolated to match the time inputs from the position var. Non-tplot array values must match the number of times in the tplot input for pos_gsm_tvar add_tilt(optional): Set this to a tplot variable name or an array of values containing the values to be added to the dipole tilt that should be used for each period. If a tplot input is used it will be interpolated to match the time inputs from the position storm(optional): Set this to use the storm-time version of the T01 model geopack_2008 (optional): Set this keyword to use the latest version (2008) of the Geopack library. Version 9.2 of the IDL Geopack DLM is required for this keyword to work. Output: Stores the result of the field model calculations in tplot variables Notes: 1. converts from normal gsm to rgsm by dividing vectors by earth's radius(6374 km) ie inputs should be in km 2. Input must be in GSM coordinates 3. Haje Korth's IDL/Geopack DLM must be installed for this procedure to work 4. either the variables setting parmod or the variables setting the individual parameter arrays should be set because the defaults aren't scientifically accurate 5. model parameters that are input as tplot variables they will be interpolated to match the time values on the input position 6.Find more documentation on the inner workings of the model, any gotchas, and the meaning of the arguments at: http://geo.phys.spbu.ru/~tsyganenko/modeling.html -or- http://ampere.jhuapl.edu/code/idl_geopack.html 7. Definition of G1 and G2 can be found at: http://modelweb.gsfc.nasa.gov/magnetos/data-based/Paper220.pdf http://modelweb.gsfc.nasa.gov/magnetos/data-based/Paper219.pdf $LastChangedBy: egrimes $ $LastChangedDate: 2014-03-18 13:03:40 -0700 (Tue, 18 Mar 2014) $ $LastChangedRevision: 14577 $ $URL: svn+ssh://thmsvn@ambrosia.ssl.berkeley.edu/repos/spdsoft/trunk/external/IDL_GEOPACK/t01/tt01.pro $