; Spreadsheet array: (n, i)
; spreadsheet(0,0) is max_n
; i is as follows:
; 1 : File number for crescent image (unique for day)
; 2 : File number for fulldisk image (unique for day)
; 3 : Time crescent fiducial spot read (s)
; 4 : Time fulldisk fiducial spot read (s)
; --- "Moonshine" values are nearest MS points (in time) to ES points
; 5 : Earthshine spot intensity (corrected only for radius)
; 6 : Moonshine spot intensity (corrected only for radius & NG filter)
; 7 : Earthshine I_0 (point corrected for airmass; no distance correction)
; 8 : Moonshine I_0 (point corrected for airmass; no distance correction)
; 9 : Earthshine I_0 (airmass fit intercept exp(beta) from Beer's Law; everything corrected)
; 10 : Moonshine I_0 (airmass fit intercept exp(beta) from Beer's Law; everything corrected)
; 11 : Earthshine alpha (atmospheric absorptivity parameter, fit result)
; 12 : Moonshine alpha (atmospheric absorptivity parameter, fit result)
; 13 : Earthshine chi squared (fit result - 'sigma')
; 14 : Moonshine chi squared (fit result - 'sigma')
; 15 : ES coaltitude parameter	- almanac
; 16 : MS coaltitude parameter	- almanac
; 17 : ES airmass parameter for this time/date	- reference fn
; 18 : MS airmass parameter for this time/date	- reference fn
; 19 : ES image Julian date - reference fn
; 20 : MS image Julian date - reference fn
;; 21 : ES phase angle	- almanac (in radians, topocentric)	
; 21 : ES phase angle	- almanac (in radians, geocentric) ; changed to geoc @QJ	
; 22 : MS phase angle	- almanac (in radians, topocentric)
; 23 : ES alpha libration (topocentric)	- almanac
; 24 : ES beta libration (topocentric)	- almanac
; 25 : MS alpha libration (topocentric)	- almanac
; 26 : MS beta libration (topocentric)	- almanac
; 27 : ES radius earth-moon (topocentric)	- almanac
; 28 : MS radius earth-moon (topocentric)	- almanac
; 29 : Lambert phase function value at ES phase angle
; 30 : Lunar phase function value at MS phase angle
; 31 : radius earth-moon (average of ES/MS from almanac)
;;32 : A* without airmass correction, from nearest MS point
; 32 : New: put asun here, selenographic longitude of the sun.
; 33 : [A*] without airmass correction, from nearest MS point
;;34 : A* with airmass correction, from nearest MS point	@Qj
; 34 : New: put bsun here, selenographic latitude of the sun.	@Qj
; 35 : [A*] with airmass correction, from nearest MS point
; 37 : Lunar phase function value at ES phase angle
; 38 : MS fit point ("airmass corrected"; N.B.: no distance correction 'cause got to be used with es_i [RAW!])
; 39 : A* using ES point & MS fit
; 40 : [A*] using ES point & MS fit
; 41 : year
; 42 : month
; 43 : day
; 44 : CS alpha (airmass absorptivity from fitting). @QJ
; 45 : CS intensity (EXP(cs_beta) with everything corrected). @QJ
; 46 : ES Sun Earth distance  @QJ
; 47 : MS Sun Earth distance  @QJ
; 48 : ES moon sun distance @QJ
; 49 : MS moon sun distance @QJ