6. Other executables in DD system

6.1 crlparam, crlparamF

One can calculate the CORALL plasma parameters (in case there is the CORALL data in DDBASE) using crlparam procedure:
crlparam YYYYDDDHHMMSSMLS TimeInterval
Here DDD is Day of Year -1 and MLS is milliseconds. TimeInterval is in seconds.

Output:

YYYDDDHHMMSSMLS T[eV] Vx[km/s] Vy[km/s] Vz[km/s] N[1/cm^3]

Parameters are calculated in the spacecraft' frame at the moment YYYDDDHHMMSSMLS (as if the spacecraft stops its rotation at YYYDDDHHMMSSMLS).
To obtain parameters in GSE, GSM or SM frames one can use crlparamF procedure:
crlparamF YYYYDDDHHMMSS Frame (sc,gse,gsm,sm)

Output:

YYYDDDHHMMSSMLS Vx[km/s] Vy[km/s] Vz[km/s] T[eV] N[1/cm^3] flag

If the "flag" is greater than 0 the parameters are meaningfull. The real density value is N*9.4

6.2 GetGSE

Spacecraft attitude is returned by GetGSE procedure:
GetGSE YYYYDDDHHMMSSMLS
Example:

 
GetGSE 1995266000000000

returns

 
    Xgse   Ygse    Zgse   Ygsm   Zgsm 
Xsc 0.995  -0.100  0.012  -0.098  0.026 
Ysc -0.076  -0.662  0.745  -0.548  0.833 
Zsc -0.067  -0.742  -0.667  -0.831  -0.552 
Rotation time = 117.9

6.3 delay

To get time lag between WIND and INTERBALL (TAIL and AURORAL) or Magnetopause Subsolar Point one can use delay procedure. Time lag is calculated under the assumption of propagation with solar wind speed.

delay YYYYDDDHHMM Here DDD is Day of Year -1.

Example:

  delay 19960001000

returns

 
Time's delay from WIND to Magnetopause : 1182.698364.  
Time's delay from WIND to Auroral Sat. : 1386.966064.   
Time's delay from WIND to Tail Sat.    : 1602.507324.

6.4 DDday

Procedure returns Day of Year minus 1 (DD day).
DDday YYYY MM DD
Here YYYY is Year, MM - month (1 - 12), DD - day (1 - 31)

6.5 fm3text, miftext

Procedures return the values of magnetic field components from FM3 or MIF magnetometers.
fm3text YYYYDDDHHMMSS DHHMMSS frame miftext YYYYDDDHHMMSS DHHMMSS frame
Here YYYYDDDHHMMSS is the start time , DHHMMSS - time interval and frame - the frame of reference=[sc,gse,gsm,sm].

6.6 GetVBKP

Program returns the frame needed to plot CORALL's ion distribution function in frame {B,V}. B is put along Z axis of new frame and Y axis lies along velocity of convection.
Program GetVBKP uses 2 minutes averaged magnetic field data, Corall data and attitude data of INTERBALL-TAIL.
Program has interactive input. Firstly you have to input the start time in usual format YYYYDDDHHMMSSMLS exactly as in the distribution function request. Then you should input the kind of frame. "0" means real {B, V} , but "1" means frame of {Byz, Vyz}, where Byz is the projection of B onto YZ of GSE reference frame, and Vyz is projection of V onto YZ of GSE.
The output of the program is:

Spin period of s/c, sec
X, Y and Z components of s/c frame in new {B,V}. frame at the time of closest sunpulse as it is necessary for distribution function objects.

Example

GetVBKP

  Type Time in YYYYDDDHHMMSSMLS format:

1996048030000000

 Type 1 if you would like to use YZ plane or 0 otherwise:

1

output

 
Time of sunpuls is 1996048030055713 
 
Spin period is       118.000 sec 
 
Spacecraft axis in VB frame: 
Xn     -0.993000   -0.0213120     0.112151 
Yn    -0.0890000    -0.474073    -0.876048 
Zn     0.0720000    -0.879952     0.468998

6.7 Tsyganenko96 model magnetic field

Program calculates and draws Tsyganenko96 model magnetic field (GSM) lines projected onto XgsmYgsm,XgsmZgsm and YgsmZgsm planes for given time. Lines are drawn in the following manner: XZ projection - field lines which intersect line parallel to Zgsm and pass through the point (XsatYsat). XY and YZ projections - field lines which intersect line parallel to Ygsm and pass through the point (XsatZsat). Step of field lines drawing (in Re) is defined by user.
User can choose different satellites - Interball-Tail, Interball-Auroral or arbitrary. In case of arbitrary satellite user should define satellite position (in GSM) for given time in resource file.
Also drawn is the Magnetopause location predicted by Shue (1997) model. Actually it is the cross-section of magnetopause by plane containing satellite and parallel to XY, YZ or XZ (GSM) planes.
The location of satellite in the graph is marked by red label.

Program is run by entering T96. It uses resource file T96.res.

An example of T96.res file:

 
1996295221000        
1300.  -40. 
-40. 10. -25. 25. -25. 25.                
0. 22. -22. 22. -22. 22. 2. 
0                                        
T               
-20 5.5 6

Here

1 string - Year, DDday, Hour, Minute, Second - moment of time for magnetic field calculation
2 string - Time delay WIND-Subsolar magnetopause (sec) and Dst value (nT)
3 string - Range for Graphs axes (Re) - Xmin, Xmax, Ymin, Ymax, Zmin, Zmax (GSM)
4 string - Step and Ranges for field lines calculations (Re). Z_start_point, Z_end_point (for XZ projection), Y_start_point, Y_end_point (for YZ projection), Y_start_point, Y_end_point (for XY projection), Step for magnetic field lines calculation is the same for all three graphs.
5 string - Hardcopy parameter: 0 - only screen output, 1 - screen and GIF output, 2 - only PostScript output.
6 string - What satellite: T - Interball-Tail, A - Interball-Auroral, any other letter - arbitrary satellite. If user chooses T or A then coordinates for the given moment for chosen satellite are taken automatically.
7 string - If user has chosen neither Tail nor Auroral Satellite then coordinates in GSM (Re) for the moment of magnetic field calculation should be defined.

If 'hardcopy' parameter (5th string) in T96.res file is nonequal to 0 then the output will be 3 files 'xz.gif'(ps),'yz.gif'(ps),'xy.gif'(ps) in the current working directory.

An output for above T96.res file: