ther      = output thrust level and attitude errors and offsets
ther(1)   = signed thrust level (+ in positive z2 roll direction,
- in negative z2 roll direction) (micropounds)
ther(2)   = attitude error about z1 pitch axis (radians)
ther(3)   = attitude error about z2 roll axis (radians)
ther(4)   = attitude error about z3 yaw axis (radians)
offset(1) = offset angle about z1 pitch axis (radians)
offset(2) = offset angle about z2 roll axis (radians)

at tabular point i (i=1,2,3)
jdles(i)  = julian day number
frles(i)  = ephemeris time fraction of day
ithrst(i) = indicate if tangential thrusters are firing
0   no firing
+1   start or middle of firing (thrust is interpolated
between tabular points with +1, or last point for
interpolation can be -1)
-1   end of firing (can only occur after a +1 has occured,
thrust is nonzero before the -1 and zero after)

thr(1,i)  = thrust (micropounds)
if thrust positive, thrusters on -z2 roll face are firing
with dv in +z2 roll direction
(west face on les9, dv with orbital velocity)
(east face on les8, dv against orbital velocity)
if thrust negative, thrusters on +z2 roll face are firing
with dv in -z2 roll direction
(east face on les9, dv against orbital velocity)
(west face on les8, dv with orbital velocity)
note: les8 is turned upside down from les9
thr(2,i)  = attitude error about z1 pitch axis (radians)
thr(3,i)  = attitude error about z2 roll axis (radians)
thr(4,i)  = attitude error about z3 yaw axis (radians)
off(1,i)  = offset angle about z1 pitch axis (radians)
off(2,i)  = offset angle about z2 roll axis (radians)
les       = data set reference number for thrust and attitude
information that was written in subroutine lesin of
input link from input cards after units were changed
nrecls      = record number that has been read on data set les
0   nothing read yet, get three records into storage
-1   end of file encountered, no more to be gotten
from this tape
>0   we are in middle of data set, be careful not to
back into load point