1. Compile and run scalarhb.f for data production.

1. Compile and run anascalar.f for the integrated autocorrelation
   time of the field: iobs=2,ntau=2**6,nskip=1,nstep=1,lmean=.false..
   Plot with a1tau.plt. Histogram with h01.plt or i01.plt.

3. Energy estimates from correlation functions:

3.1 Compile and run ana_scorj.f. Plot with cor0.plt and cor1.plt.

3.2 Fits first with massfit.f. Copy c*000.d data sets to fort.10.
    Plot with gfit.plt and gfitln.plt. Everything together: cor.plt.

3.3 Jackknife analysis with massfitj.f.
    CJ="c0j": XMM  mass estimate (momentum k=0). Run first for xmj.d file.
    CJ="c1j": EML  energy estimate from correlation function.
              EL   energy estimate from m and lattice momentum 2*sin(k1/2).
              E1   energy estimate from m and k1=2pi/N momentum.
              XMM  mass used.
              DEL  deviation (EL-EML)/EL in %.
              XK1M momentum extrapolation k1=sqrt(EML**2-m**2).
              XKL  lattice momentum l2=2*sin(k1/2).
              XK1  k1=2pi/N momentum.
              DEL  deviation (XKL-XK1m)/XKL in %.