: fixed error in MEANSHR, SHRGT45, ERRBZ wherein all were giving nans at times -- this was due to some weird condition being satisfied, so now i check on that condition

corrected gamma, alpha, and error in alpha
this is an incomplete version -- errors in errors remain
but i am pushing it over now so it is available on the web

corrected indexing in the computation of jz_err in computeJz()

Changed all sums to doubles instead of floats, since there is a significant loss of precision in summing millions of floats:

P% cat prec.c
#include <stdio.h>
#include <stdlib.h>

main()
  {
  float sum = 0.0;
  double dsum = 0.0;
  int n = 13000000;
  int i;

  for (i=0; i<n; i++)
    {
    sum += 1.333;
    dsum += 1.333;
    }

  printf("avg = %f, davg = %f\n", sum/n, dsum/n);
  }

P% cc prec.c -o prec
P% prec
avg = 1.181585, davg = 1.333000
P%

reinstated old version of calculating mean free energy; total free energy remains unchanged

Corrected factor of 8*PI in total and mean free energy calculation

Rewrote module to compute errors using formal error propagation and VFISV error arrays

corrected error in rsun_ref / rsun_obs

Use a smoothed jz array using fresize.h routines.
Modified the methodology for calculating alpha.

Modified code to select pixels in conf_disambig map that are greater than 70 and pixels in bitmap that are greater than 30.

Update functioning code sharp.c per Monica's request