setup and read quantities into temporary storage , then
rearrange into labelled common, such that corrections belonging
to a given planet are grouped together.  the order of the
planets is as follows

for j=1,...,numphs we have
phsit(j)= first four characters of eight character optical
observing site name
phser(j)= four characters giving observation series name
nplphs(j)= planet number of observed body
0= sun              1= mercury          2= venus
4= mars             5= jupiter          6= saturn
7= uranus           8= neptune          9= pluto
10= moon
aphase(i,j)= dimensionless phase coefficients
ncphs(j)= number of coefficients in phase model
lphs(i,j) =0 aphase(i,j) not adjusted
lphs(i,j) =1 aphase(i,j) adjusted in least squares analysis

viii. planetary phase corrections  (phsred)
for each planet and observation series
card 1
columns
1- 4  observing site name  first four characters        (1a4)
6- 9  observation series name  first four characters    (1a4)
10-12  planet number                                     (i3)
13-15  number of coefficients in phase correction mode   (i3)
limitation is 9 (if less than or equal 3, only one card read)
16-17  l1=1,0  first   coefficient is adjusted or not    (i2)
18-19  l2=1,0  second  coefficient is adjusted or not    (i2)
20-21  l3=1,0  third   coefficient is adjusted or not    (i2)
22-23  l4=1,0  fourth  coefficient is adjusted or not    (i2)
24-25  l5=1,0  fifth   coefficient is adjusted or not    (i2)
26-27  l6=1,0  sixth   coefficient is adjusted or not    (i2)
28-29  l7=1,0  seventh coefficient is adjusted or not    (i2)
30-31  l8=1,0  eigth   coefficient is adjusted or not    (i2)
32-33  l9=1,0  nineth  coefficient is adjusted or not    (i2)
37-48  first   phase coefficient                         (e12.5)
49-60  second  phase coefficient                         (e12.5)
61-72  third   phase coefficient                         (e12.5)
card 2
1-12  fourth  phase coefficient                         (e12.5)
13-24  fifth   phase coefficient                         (e12.5)
25-36  sixth   phase coefficient                         (e12.5)
37-48  seventh phase coefficient                         (e12.5)
49-60  eigth   phase coefficient                         (e12.5)
61-72  nineth  phase coefficient                         (e12.5)
limitation  no more than u_mxphs sites allowed