** Add changes at the beginning! **

15/11/2002 A. Moralejo

ph2cph.c: Fixed BUG in timing calculation!!! The time to be subtracted 
(mirror till ground) had the wrong sign!!!

reflector.c: in the wavelength range check, changed the range limits to 
the fixed values 290 and 600 nm. The StarFieldAdder cer files could not be 
processed otherwise, because the w range is not written in their headers.

ph2cph.c: changed sign of phi angle to make it positive. This was necessary
since now the vector indicating the reflected photon trajectory has also 
changed sign! (see note below, 11/11).

14/11/2002 A. Moralejo
reflector.c, parms.c Added wobble mode option. Added wobble flag to the
output and also an atmospheric_model flag.

13/11/2002 A. Moralejo
reflector.c, header.h: Introduced 3 counters to keep track of what proportion 
of the C-photons in each event have been produced by electrons, muons or 
other particles.

11/11/2002 A. Moralejo
ph2cph.c: BUG found! and fixed: the direction cosines in the cer files are 
those of the downgoing directions of photons, hence third component is 
negative. Before it was assumed positive. This, together with a wrong 
reflection in the individual mirrors (also changed) resulted in a subtle 
error which produced a quite de-focused reflector!

reflector.c: for the same reason as above, the sign of theta in the calls 
to makeOmega had to be changed (+pi added to the phi argument).

geometry.c, ph2cph.c, reflector.c, init.h: eliminated array ct_Focal[] 
which contained the focal lengths of individual mirrors. It was redundant
since the same information is available in ct_data[].f

Oct/2002 A. Moralejo
added attach.c to attach the files magic.def, axisdev.dat and
reflectivity.dat to the end of the reflector output.

Sept/2002 A. Moralejo
header.h, reflector.c: added run header (=Corsika's), changed event header.

3/7/2002 A. Moralejo
reflector.c: Introduced NaN check in the photon loop. 
If NAns are found in a photon data block, it is not processed.

ph2cph.c: Introduced "check of positiveness" before taking the sqrt 
in the calculation of the photon trajectory intersection with the
global paraboloid.