mpage: 1 for light curves, 2 for radial velocity curves, 5 for maps
nref: number of reflections if mref=2 (1)
mref: specifies whether simplified or detailed reflection is to be applied (1)
ifsmv1: whether primary spots are allowed to move in longitude or not (0)
ifsmv2: whether secondary spots are allowed to move in longitude or not (0)
icor1: computed RV’s are to include proximity and eclipse corrections or not (1)
icor2: computed RV’s are to include proximity and eclipse corrections or not (1)
if3b: is there a light time effect due to a kinematic third body (0)
ld1: primary limb darkening law (linear=1, log=2, or square root=3) (1)
ld2: secondary limb darkening law (linear=1, log=2, or square root=3) (1)
kspev: whether spots are to grow and decay or not (0)
kspot: whether the old simple spot model or a more accurate one applies (0)
nomax: spot growth and decay timewise profile is trapezoidal or triangular (1)
ifcgs: light curves are to be computed in arbitrary units or physical units (0)
ktstep: only for mpage=6

jdphs: whether time (1) or phase (2) is the indep. var. in light and RV curves
hjd0: ephemeris reference time – usually heliocentric Julian Date
pzero: Orbit period in mean solar days at time HJD0
dpdt: Rate of period change, assumed constant (dimensionless)
pshift: constant shift applied to move a light or RV curve in phase
delph: phase spacing of Gauss. quadrature abscissas for phase/time smearing sim.
nga: number of Gaussian quadrature abscissas in phase/time smearing simulation
stdev: Fractional std deviation of simulated scatter in output light or RV (0)
noise: how output scatter (simulated noise) is to vary with light level.
seed: number that determines individual noise excursions of light or RV points.

hjdst: Start time for output points
hjdsp: Stop time for output points
hjdin: Increment for times of output points
phstrt: Starting phase for output points. Applied only if JDPHS=2.
phstop: Stop phase for output points.
phin: Increment for phases of output points
phn: Phase of normalization for columns of normalized light and set magnitude
phobs: Phase at which a spectroscopic temperature was estimated
lsp: Tells whether TOBS and TAVSP are for star 1 or star 2.
tobs: Spectroscopic temperature at phase PHOBS
tavsp: Computed flux-weighted mean surface temperature

MODE: model constraints to be applied (wdtype)
IPB: compute L2 from T2 (0); decouple L from T (1)
IFAT1: blackbody (0) or stellar atmosphere (1)
IFAT2: blackbody (0) or stellar atmosphere (1)
N1: number of latitude rows per hemisphere (10)
N2: number of latitude rows per hemisphere (10)
perr0: (initial) argument of periastron for star 1 in radians
dperdt: The first time derivative of ω.
the: semi-duration of primary eclipse in phase units
VUNIT: radial velocity input unit in kilometers per second

E: eccentricity
A: length of the semi-major axis in solar radii of the relative orbital ellipse
F1: ratio of axial rotation rate to orbital rate (1)
F2: ratio of axial rotation rate to orbital rate (1)
VGA: Gamma (systemic) velocity of barycenter (double or triple)
XINCL: orbital inclination
GR1: exponents in the bolometric gravity brightening
GR2: convective envelopes: 0.3, radiative: 1.0
abunin: metallicity as usually defined logarithmically, relative to the Sun
fspot1: spot angular drift rates in longitude
fspot2: rate 1.0000 means that drift just matches the mean orbital angular rate

tavh: mean surface effective temperature (hot) star
tavc: mean surface effective temperature (cool) star
alb1: bolometric albedo
alb2: convective envelopes: 0.5-1.0, radiative:1.0
pot1: star 1: modified surface ’potentials’ 
pot2: star 2: = semi-major axis / radius
rm: mass ratio of binary components, m2/m1
xbol1: linear coefficient in bolometric limb darkening law
xbol2: linear coefficient in bolometric limb darkening law
ybol1: 2nd coefficient for logarithmic or quadratic law
ybol2: 2nd coefficient for logarithmic or quadratic law
dpclog: logarithm (base 10) of distance (d) in parsecs

With tertiary component only:
a3b: semi-major axis of the relative third body orbit in solar radii.
p3b: orbit period in mean solar days
xinc3b: orbital inclination
e3b: eccentricity
perr3b: Argument of periastron of third body orbit, in radians.
tconj3b: HJD of superior conjunction of the EB center of mass with the tertiary

iband: Photometric band identifier
HL: band luminosity (star 1)
CL: band luminosity (star 1)
X1: Limb darkening coefficients in the linear cosine law, star1
x2: Limb darkening coefficients in the linear cosine law, star2
y1: Coefficients in the non-linear limb darkening terms, star1
y2: Coefficients in the non-linear limb darkening terms, star2
EL3: third light
opsf: opacity in a specific band, related to ISM
ZERO: A reference magnitude that can shift a magnitude light curve
FACTOR: scaling factor for LC output curves in flux
wl: Wavelength of light, in microns, now used only for circumstellar matter
aextinc: Interstellar extinction in the designated photometric band.
calib: Flux calibration constant in erg s−1 cm−3 for a star of magnitude 0

XLAT(KP,I):
XLONG(KP,I):
RADSP(KP,I):
TEMSP(KP,I):
tstart(kp,i): Onset time of spot
tmax1(kp,i): Start time of constant maximum size for spot
tmax2(kp,i): End time of constant maximum size for spot
tfinal(kp,i): Disappearance time of spot

Spectral line profiles only (mpage=3):
binwm1,sc1,sl1,nf1
wll1(iln),ewid1(iln),depth1(iln),kks(1,iln)
binwm2,sc2,sl2,nf2
wll2(iln),ewid2(iln),depth2(iln),kks(2,iln)
xcl(i),ycl(i),zcl(i),rcl(i),op1(i),fcl(i),edens(i),xmue(i),encl(i)

Backup before suffix "15" change: imlc.bck.f
Backup before "cleaning" change: imlc15.bck.f
Tests:
lctest.f: parms for lcout.active
