ULYSSES In Space

 

Sun Banner
ULYSSES

Ulysses HISCALE Data Analysis Handbook

 

Appendix 10. Effect of Backscattered Electrons on the Geometric Factors of the LEMS30 Telescope (Hong MS Thesis)

 

 

Most of the code in this project is adapted from Buckley's thesis. Traj4pt is the main program reproduced from Buckley's Traj6pt to trace the trajectories of the electrons. It needs the following programs to make it complete: Chechitsen4.for, Refsub4.for, Trajrt4bdet.for, Trajrt4.for,Tracksub4.for, Timer.for, and Fdmod1.for. Most of the program names are retained except the variation of the number. For example, Chechitsen4.for is the modification of Chechitsen1.for. Trajrt4bdet.for and Trajrt4.for are the same programs except that Trajrt4bdet.for is used to handle the backscattering situation and Trajrt4.for is for non-scattering. Refsub4.for is a new program which includes all the subroutines for the backscattering simulation. There is a common file named Proba.cmn to support Refsub4.for, which calculates the backscattering probability according to the backscattering coefficient formula.

 

To run Traj4pt one must compile and link all programs mentioned above. Shown here is a brief description of what is needed for the input data once Traj4pt is run. For detailed descriptions please refer to Appendix 9.

 

ENTER THE CHOICE OF MODEL OF THE SCATTERING 
NON SCATTERING: 0
SPECULAR SCATTERING: 1
DIFFUSE SCATTERING: 2
   Input one of the numbers to select a particular backscattering mode. 
ENTER THE NAME OF THE FILE OF COORDINATES
   Input coord.dat
ENTER THE NAME OF THE FILE OF COEFFICIENTS
   Input coeff.dat
ENTER THE INITIAL COORDINATES OF THE e IN INCHES
   Input the starting coordinates in inches/100
ENTER THE LOWER BOUND ON TIME 
   Input 0
ENTER THE UPPER BOUND ON TIME 
   Input a long time (1 or 2 seconds would be good) 
ENTER THE TIME STEP 
   Input the time step in which an electron moves according to the differential equation solver. The time step must be chosen in such a way that it won't take too much time but will get an accurate trajectory.
ENTER THE ERROR BOUND
   The selection of error bound is about .05 of the time step.
ENTER THE ENERGY OF PARTICLE IN MEV
   Input the value as instruction.
ENTER MID. THETA, THETAMAX, DEL THETA 
   Input the values of the range of the polar angle.
ENTER MIN. PHI, NO. OF STEPS, DEL PHI
   Input the values of the range of the azimuthal angle.

 

The system which creates the name of the output files is retained. Those file names are of the following form, (coordinate).E(energy) IN. An output file records the escaped electron's starting polar and azimuthal angles for the particular coordinate and energy. One can use the files to calculate the geometric factor for a particular energy after the program has run through all small areas. These file names are long and not easy to type, so a new program named Cname.pas was written in PASCAL to generate a whole list of the file names in a file named filename.dat for the inut energy. FactorD.for is written to calculate the geometric factor for a particular energy after filename.dat is generated.

 

Traj3pt.for with Chechitsen4.for, Refsub4.for, Timer.for, Fdmod1.for, Tracksub3, Trajrt3bdet.for, Trajrt3.for would record the coordinates of the trajectories. It is modified from Buckley's Traj2pt. And the difference between Traj4pt and Traj3pt is the output files. Traj3pt records the coordinates of the trajectories which produce the graphs in Section A10.10. Another new program named Geom.for combined with Cname.pas is written to read (coordinate).E(energy)IN files and to generate a com file to run Traj3pt; it produces only those trajectories which escape the open aperture.

 

For detailed code see Section A10.11.

 

 

Next: A10.9 Appendix - Number of Escapes

 

Return to the Table of Contents for Hong's MS Thesis


Updated 8/8/19, Cameron Crane

QUICK FACTS

Manufacturer: ESA provided the Ulysses spacecraft, NASA provided the power supply, and various others provided its instruments.

Mission End Date: June 30, 2009

Destination: The inner heliosphere of the sun away from the ecliptic plane

Orbit:  Elliptical orbit transversing the polar regions of the sun outside of the ecliptic plane