adj  =adjustment to parameter in units of parameter
sig  =standard deviation of parameter estimate in units of paramtr
nwv  =new value of parameter in units of parameter
fract=adj/sig (large number to printout **** if sig=0.0_10)
ntype=1 if abs(adj).ge.eps(9)*sig
ntype=2 if abs(adj).ge.eps(10)*sig
ntype=3 if neither of the above
nkind(ntype)=running total of three types of adjustments
n    =index for solution vector of normal equations
line =line counter
il=1 blank line not yet written as separator
il=2 blank line written as separator
fradj= fraction (between zero and one) by which adjustment
is to be multiplied (input as eps(14))
statf= statistics for fract
nf   = number of points included in fract statistics