
Physical Sciences Lab - PSLC (VAX 8600)

Username: SANCHEZ
Password: 
            Welcome to VAX/VMS Version 5.1-1 (DECwindows)
                  University of Wisconsin - Madison
                     Physical Sciences Laboratory

    Last interactive login on Friday, 28-JUL-1989 14:34
    Last non-interactive login on Friday, 28-JUL-1989 10:47

    The time is Friday 28-JUL-1989 14:52:18.98 CST +0/-3 Seconds
    Node PSLC, Port _RTA1:, Default directory D5:[SANCHEZ]
$ SET DEF [.DEMO]
$ SHADOW
 
 
 
 
 
Defining Logical Names and Environmental Parameters
Graphic device set at TT:. Enter $SETPLOT to modify.


All SHADOW logical names and command defined.
Type $NEWS for more info on SHADOW.


$ BRAGG
All crystal structures are refered to a cubic unit cell.
Bravais lattice type : 
0	for 	ZincBlende
1	for 	Rocksalt
2	for	simple FCC
3       for     CsCl structure
Then ? 
0
Lattice constant (Angs) ? 
5.43
Index of crystal plane of reflection H,K,L : 
1
1
1
***********************************************************
The ZINCBLENDE structure is defined by atom A located at 
(0,0,0) and atom B at (1/4,1/4,1/4) of the fcc lattice.
***********************************************************
Enter 2-letters (capitalized) atomic symbol for atom A : 
SI
Enter 2-letters (capitalized) atomic symbol for atom B : 
SI
***********************************************************
Atomic scattering factor is defined by fo + f' + if", where
    fo = fo(SIN(theta)/Lambda) is the non-dispersive part
    f', f" (Lambda) are the dispersive part.
***********************************************************
We need fo at 3 different values of SIN(theta)/Lambda, which
should cover the range of interest and center around :
SIN(theta)/Lambda =   0.1594890246380182     ratio.
***********************************************************
Please enter 1) SIN(theta)/Lambda, 2) fo ,
For atom A, first set : 
.1,12.16
          , second set : 
.2,9.67
          , third set : 
.3,8.22
For atom B, first set : 
.1,12.16
          , second set : 
.2,9.67
          , third set : 
.3,8.22
***********************************************************
f', f" is furnished from optical constant library within ...
minimum photon energy (eV) : 
7700
maximum photon energy (eV) : 
7800
energy step (eV) : 
10
Do you want to include crystal absorption [1/0] ? 
1
Temperature (Debye-Waller) factor : 
1.0
Output file name (for SHADOW) : 
SI.PAR
Do you want to generate a rocking curve [1/0] ? 
1
... at what energy (eV) ? 
7750
 
So far, we are working with:
Lambda 		 =    1.599809032258065     Angstroms
Theta (graz) 	 =    14.78259099147192     degrees
Bragg angle      =    14.78449376730854     degrees
Structure factor  =  (61.29837293588675,1.976536272774998)
Refraction index = 1 - delta - i*beta :
           delta =   8.1919893381993125E-06
            beta =   2.0040488614427989E-07
Absorption coeff =    157.4164179247638     cm-1
 
Do you want an asymmetrical cut [1/0] ? 
1
Angle (deg) between face and bragg planes (CW)= 
7
The width of the Rock Curve in function of incident angle is
1/2 width for s-pol  =  :    10.76098572573549     microradians
1/2 width for p-pol  =  :    9.359850847544060     microradians
The width of the Rock Curve in function of reflected angle is
1/2 width for s-pol  =  :    29.48906494122322     microradians
1/2 width for p-pol  =  :    25.64943923522630     microradians
Asymmetric factor b=  :    2.740368372638810    
    
+/- how many microradians : 
80
How many points : 
150
Do you want to use sec [1/0] ? 
1
Do you want the R.C. not centered [1/0] ? 
1
Do you want to try another energy [1/0] ? 
0
$
\
\
\
\  Notice that the file SI.PAR used for SHADOW do NOT
\  depend if asymmetrical or symmetrical diffraction is chosen.
\  You do not have to run BRAGG again for change the asymmetrical
\  cut in TRACE procedures.
\
\
\
\
$ DIR

Directory D5:[SANCHEZ.DEMO]

DEMO.LOG;1                0  28-JUL-1989 14:52:10.57  (RWED,RWED,RE,)
PLOT.TOP;1                1  28-JUL-1989 14:51:53.90  (RWED,RWED,RE,)
ROCK_CURVE.PAR;1          3  28-JUL-1989 14:54:22.92  (RWED,RWED,RE,)
ROCK_CURVE_INC.P;1       15  28-JUL-1989 14:54:01.80  (RWED,RWED,RE,)
ROCK_CURVE_INC.S;1       15  28-JUL-1989 14:54:01.64  (RWED,RWED,RE,)
ROCK_CURVE_REF.P;1       15  28-JUL-1989 14:54:02.08  (RWED,RWED,RE,)
ROCK_CURVE_REF.S;1       15  28-JUL-1989 14:54:01.93  (RWED,RWED,RE,)
SI.PAR;1                  4  28-JUL-1989 14:53:50.49  (RWED,RWED,RE,)

Total of 8 files, 68 blocks.
$ TY PLOT.TOP
MOD NOECHO
@ROCK_CURVE_INC.S
JOIN
@ROCK_CURVE_INC.P
PLOT
@ROCK_CURVE_REF.S
JOIN
@ROCK_CURVE_INC.P
PLOT
$ TDFILE
TOP DRAWER (10/10/80) NOW EXECUTING

@PLOT.TOP
READING FROM INDIRECT COMMAND FILE: PLOT.TOP 
           1  MOD NOECHO                     
                                                    
READING FROM INDIRECT COMMAND FILE: ROCK_CURVE_INC.S
           2    -11.82631107269572      4.5692957405253544E-03  
                                                                           
END OF INDIRECT COMMAND FILE
         152  JOIN                                              
                                                                
READING FROM INDIRECT COMMAND FILE: ROCK_CURVE_INC.P            
         153    -11.82631107269572      3.4491750826317568E-03  
                                                                           
 PRINTING SUPPRESSED AFTER     0 POINTS.  (PLOT CONTINUES.)
END OF INDIRECT COMMAND FILE
         303  PLOT                                              
                                                                
READING FROM INDIRECT COMMAND FILE: ROCK_CURVE_REF.S            
         304    -3.690421467485779      3.6519038661560184E-02  
                                                                           
 PRINTING SUPPRESSED AFTER     0 POINTS.  (PLOT CONTINUES.)
END OF INDIRECT COMMAND FILE
         454  JOIN                                         
                                                                           
READING FROM INDIRECT COMMAND FILE: ROCK_CURVE_INC.P       
         455    -11.82631107269572      3.4491750826317568E-03
                                                                           
 PRINTING SUPPRESSED AFTER     0 POINTS.  (PLOT CONTINUES.)
END OF INDIRECT COMMAND FILE
         605  PLOT                                            
                                                                           
END OF INDIRECT COMMAND FILE

FORTRAN STOP
$ DTEKON
$ UGTEK

  UG element file: 
SAVE
  Plot which pictures (default=all)? 

  Enter picture scale factor (default=1.0): 

$ DTEKOFF
$ go trace


  :::::::  :::    ::            :::    :::::::        :::::     :::  ::  ::
 ::::::::  :::    ::           ::::    ::::::::      :::::::    :::  ::  ::
:::        :::    ::          ::: :    :::   :::    :::    ::   :::  ::  ::
:::        :::    ::         ::: ::    :::    ::    :::    ::   :::  ::  ::
 :::::::   :::::::::        :::  ::    :::    ::    :::    ::   :::  ::  ::
  :::::::  :::::::::       :::   ::    :::    ::    :::    ::   :::  ::  ::
       ::  :::    ::      :::    ::    :::    ::    :::    ::   :::  ::  ::
      :::  :::    ::     :::     ::    :::   :::    :::   :::   ::: ::: :::
::::::::   :::    ::    :::      ::    ::::::::      :::::::     :::::::::
:::::::    :::    ::   :::       ::    :::::::        :::::       ::: :::




VERSION     88.2                                               1 June  1988




Hit <RETURN> when ready: 

Ray Tracing Selected. Begin procedure.

$ go source

SOURCE selected. Begin procedure.

This procedure generates a SOURCE for SHADOW.
Mode selected [ ? <ret> for HELP ] ? : prompt
Call to INPUT
---------------------------------- S H A D O W  ------------------------------- 
Jun 1988                                         F.Cerrina        ECE - UW      

Defining source : 
When prompted for a yes/no answer, you may enter:
for   YES  answer      Y, 1
for   NO   answer      anything else




Do you want a verbose [ 1 ] or terse [ 0 ] output ?
1
------------- SOURCE SPECS ------------------
Options available:
Random in  BOTH REAL and MOMENTUM space		      0
Grid       BOTH					      1
Grid       REAL, random  MOMENTUM	              2
Random     REAL,  grid   MOMENTUM		      3
Ellipses in PHASE space, random around each ellipse   4
Ellipses in PHASE space,   grid around each ellipse   5
Source modelling type [ 0-5 ] ? 
0
How many rays [ 1 - 5 000 ] ? 
1500
Seed [ odd, 1000 - 1 000 000 ] ? 
133587
Do you want to optimize the source ? 
0
Source type : [ 0 ] regular source
              [ 1 ] wiggler
              [ 2 ] undulator
Then ? 
0
 The source is specified in the laboratory reference frame. The program will then rotate the set of 
rays in the mirror frame.

Type of source,now.

         use  ( 0 ) for point source  
              ( 1 ) for rectangular s.
              ( 2 ) for elliptical  s.
              ( 3 ) for gaussian    s.

X-Z plane source type [ 0-3 ] ? 
0
Source depth. The actual source will be centered on the no-depth position. Use 
         	(1) for no depth,
         	(2) for flat depth distribution,
         	(3) for gaussian depth distribution,
                (4) for a synchrotron source depth distr.

                   	Then ? 
Source Depth [ 1-4 ] ? 
1

O.K., got it so far.
         Source distribution now. We may use 
             ( 1 ) for a flat source
             ( 2 )       uniform   s.
             ( 3 )       gaussian  s.
             ( 4 )       synchrotron 
             ( 5 )       conical 
             ( 6 )       exact synchrotron 

Source Angle Distribution [ 1-6 ] ? 
0
Horizontal half-divergence [ (+)x, rads ] ? 
0
                           [ (-)x, rads ] ? 
0
Vertical                   [ (+)z, rads ] ? 
.0002
                           [ (-)z, rads ] ? 
.0002
Do you want a Photon energy [ Y/N ] ? 
1
We have these choices :
Single line ......................... 1
Several lines ....................... 2
Uniform source....................... 3
Energy distribution [ 1-3 ] ? 
1
Photon Energy [ 0 ] or Angstroms [ 1 ] ? 
0
Energy [ eV ] ? 
7750
Do you want to store the optical paths (OPD) [Y/N] ? 
0
Do you want to generate the A vectors (electric field) [Y/N] ?
0
Exit from INPUT_SOURCE

Generated          250 rays out of         1500
                   500
                   750
                  1000
                  1250
                  1500
Exit from SOURCE
SOURCE => Source has been succesfully generated.

SOURCE procedure completed.

Shadow:: exit

Exit to DCL

$ go trace

Ray Tracing Selected. Begin procedure.

Mode selected [ ? <ret> for HELP ] ? : prompt
PROMPT selected.
Call to RESET
Exit from RESET
Mode selected is: 
PROMPT                                                                          
Options: to start anew                 [ 0 ] 
            to restart from a given OE [ 1 ] 
Then ? 
0
Call to INPUT_OE
---------------------------------- S H A D O W  ------------------------------- 
Jun 1988                                         F.Cerrina        ECE - UW      
 
When prompted for a yes/no answer, you may enter:
for   YES  answer      Y, 1
for   NO   answer      anything else
 

Defining Optical Element:  1
Continue ? [ ^Z or %EXIT to terminate OS ] 




Do you want a verbose [ 1 ] or terse [ 0 ] output ?
1
You may save disk space by not writing out the  intermediate STAR or MIRR data f
iles. In general you will not need them unless you have specific needs (footprin
ts, etc.) 
 
Files to write out. Options: 
All............................ [ 0 ] 
Mirror only.....................[ 1 ] 
Image at CP only................[ 2 ] 
None............................[ 3 ] 
 Then ? 
0
	Let's define the optical or central axis of the system
for this particular optical element.
	By this I mean a "virtual" line that runs throughout
the optical system. Along this line are located the 
"continuation" planes, where the OS is subdivided in the
individual OE. This line does not have to coincide with the
true optical axis, as it is used mainly for bookeeping the
data, but it helps greatly in the data analysis if this
identity is preserved as everything in the program is referredto it.
	Once established, you still have complete freedom of
"moving" around the mirrors. 
In the case of a grating, you will have several choices. The 
program may override your specifications for the central axis 
and locate the source and image at the "best" position. 
You will be prompted later for this option.
It is recommended to use CM as units. This is not critical for
 most cases, but it is in the case of diffraction elements.
Optical Element definition:
Incidence Angle ? 
10
Source Distance ? 
5000
Reflection Angle? 
10
Image Distance  ? 
5000
Reflector [ 0 ] or refractor [ 1 ] ? 
0
Lets define the mirror. I may compute its parameters, like the
radius or the axes. This will not affect the rest of the calcu
lations; all the geometrical parameters may be modified later.
Or, you may wish to specify the mirror parameters yourself.


What kind of surface are we dealing with ?
spherical    = 1
elliptical   = 2
toroidal     = 3
paraboloid   = 4
plane        = 5
Codling slit = 6
hyperbolical = 7
cone	      = 8
polynomial   = 9

Mirror surface [ 1-9] ? 
5
Is the mirror Cylindrical ? 
0
Is this optical element a Fresnel Zone Plate ? 
0
Are we dealing with a Grating ? 
0
Are we dealing with a crystal [ Y/N ] ? 
1
File containing crystal parameters ?
si.par
Is the crystal asymmetric [ Y/N ] ? 
1
Enter the angle formed by the crystal planes with the crystal surface. Use a pos
itive value if the plane normal lies along the +Y direction (focusing), negative
 otherwise.
Planes angle [ deg ] ? 
7
Automatic Tuning of Crystal [ Y/N ] ?
1
You must then supply the wavelength (or photon energy) at 
which you want the crystal tuned.
Energy, in eV, [ 0 ] or wavelength, in Angs., [ 1 ] ? 
0
Photon Energy ? 
7750
Is the mirror convex [ Y/N ] ? 
0
Reflectivity of Surface. SHADOW may solve the Fresnel equations locally.
Available options:'No reflectivity              .......... 0 
Full polarization dependence .......... 1 
No        "           "      .......... 2 
Reflectivity mode [ 0,1,2 ] ? 
0

Mirror orientation angle. Angles are measured CCW, in deg,
referring to the mirror normal. Alpha=0 is the mirror
sitting at the origin, facing up. Alpha = 90 is the 
mirror standing at your right and facing left when you 
look along the beam STANDING ON THE PREVIOUS MIRROR and so on.
Orientation Angle [ Alpha ] ? 
0
Mirror Dimensions finite [ Y/N ] ?
0
 The Optical Element and the the relative mirror are now fully 
defined. The mirror pole is now located at the "center" of the 
optical element. It is possible to override this situation and 
"move" the mirror without affecting the rest of the system.
 It is also possible to move the "source" without affecting the
rest of the system.
 The movements are expressed in the DEFAULT Mirror Ref. Frame.
 so that if you move BOTH source and mirror the relative
movement is the vector sum of the individual ones.
A word of caution: SOURCE movements and MIRROR movements are
NOT equivalent from the point of view of the whole system.
Do you want to move the Source [ Y/N ] ? 
0
   
--- Mirror rotations and position. ---
We define three angles, as rotations around the three axis.
These rotation are  defined in the program as corrections to
the mirror nominal position; that is, they modify the mirror
position relative to the Default Mirror Reference Frame, where
all the calculations are performed. Remember that rotations
do NOT commute. I apply them in the same order of entry.
CW ROTATIONS are (+) angles. A translation can be also applied to the mirror.
Do you want to move the mirror itself [ Y/N ] ? 
0
Distorted surface [ Y/N ] ? 
0
Any screens in this OE [ Y/N ] ? 
0
Slit at continuation plane [ Y/N ] ? 
0
Extra Image plates [ Y/N ] ? 
0
File containing the source array ? 
begin
Exit from INPUT
Tracing optical element #            1
Call to SETSOUR
Exit from SETSOUR
Call to IMREF
Exit from IMREF
Call to OPTAXIS
Exit from OPTAXIS
Call to MSETUP
Exit from MSETUP
Call to RESTART
Exit from RESTART
Call to MIRROR
Exit from MIRROR
Call to IMAGE
Exit from IMAGE
Call to RESET
Exit from RESET
Do you want to change input mode ? 
0
Call to INPUT_OE
---------------------------------- S H A D O W  ------------------------------- 
Jun 1988                                         F.Cerrina        ECE - UW      

Defining Optical Element:  2
Continue ? [ ^Z or %EXIT to terminate OS ] 

End of session
 
Procedure completed. Return to COMMAND level
 
Procedure completed.

Shadow:: exit

Exit to DCL

$ preplot
PREPLOT> Input file ? 
star.01
Read                 1500 rays.
Each ray has           12 entries.
PREPLOT> How many columns to write out ? 
2
Row [1-12]  : the individual column
Row   [20]  : R = SQRT(X**2 + Y**2 + Z**2)
Row   [21]  : angle from the Y-axis
Row   [22]  : the magnitude of A vector
Row   [23]  : A**2
PREPLOT> Row # 1 : 
6
PREPLOT>       2 : 
23
Options --- Enter
0 for excluding the losses
1 for including losses at a particular O.E.
2 for plotting all the rays .
3 for plotting ONLY the losses (all of them)
4 for plotting ONLY the losses at a given O.E.
PREPLOT> Then ? 
0
***********
Found         1500 good points out of        1500
Output options :
  [ 0 ] store rays in a file
  [ 1 ] plot directly on screen
  [ 2 ] both
Then ? 
0
PREPLOT> Output file ? 
star.top
$ dir

Directory D5:[SANCHEZ.DEMO]

BEGIN.DAT;1             282  31-JUL-1989 09:58:10.25  (RWED,RWED,RE,)
EFFIC.01;1                1  31-JUL-1989 10:00:06.40  (RWED,RWED,RE,)
END.00;1                  7  31-JUL-1989 09:58:12.07  (RWED,RWED,RE,)
END.01;1                 11  31-JUL-1989 10:00:09.44  (RWED,RWED,RE,)
MIRR.01;1               282  31-JUL-1989 10:00:05.16  (RWED,RWED,RE,)
OPTAX.01;1                2  31-JUL-1989 09:59:44.44  (RWED,RWED,RE,)
PLOT.TOP;1                1  28-JUL-1989 15:09:16.57  (RWED,RWED,RE,)
ROCK_CURVE.PAR;1          3  28-JUL-1989 14:54:22.92  (RWED,RWED,RE,)
ROCK_CURVE_INC.P;1       15  28-JUL-1989 14:54:01.80  (RWED,RWED,RE,)
ROCK_CURVE_INC.S;1       15  28-JUL-1989 14:54:01.64  (RWED,RWED,RE,)
ROCK_CURVE_REF.P;1       15  28-JUL-1989 14:54:02.08  (RWED,RWED,RE,)
ROCK_CURVE_REF.S;1       15  28-JUL-1989 14:54:01.93  (RWED,RWED,RE,)
SAVE.DAT;1               24  28-JUL-1989 15:09:25.95  (RWED,RWED,RE,)
SI.PAR;1                  4  28-JUL-1989 14:53:50.49  (RWED,RWED,RE,)
STAR.01;1               282  31-JUL-1989 10:00:07.72  (RWED,RWED,RE,)
STAR.TOP;1               94  31-JUL-1989 10:01:46.62  (RWED,RWED,RE,)
START.00;1                7  31-JUL-1989 09:58:02.68  (RWED,RWED,RE,)
START.01;1               11  31-JUL-1989 09:59:43.28  (RWED,RWED,RE,)

Total of 18 files, 1071 blocks.
$ lo
   Resource             Amount   Dollars      
CPU time (Minutes)         .30       .77    2.550/Minute             
Page swaps               3,977       .17     .042/Thou               
Direct I/O                 299       .04     .136/Thou               
Buffered I/O               782       .03     .042/Thou               
Connect time (Hours)       .11       .02     .210/Hour               
 
TOTAL:                              1.03    Interactive    
  SANCHEZ      logged out at 31-JUL-1989 10:01:59.12
