removed debug variable countabit

removed check for ( (derx_los[j * nx + i] + dery_los[j * nx + i]) == 0) in computeLOSderivative

Performed a check on count_mask before assigning a value to ERRMSHA to avoid a divide by zero error.
If count_mask is indeed zero then ERRMSHA is now NaN.

reverted to 1.35 and added comment to SAVNCPP header

out of bounds problem still remained in Function 8; fixed it

changed the numerical derivative. here is a description of the problem that i fixed:

In computeJz(), the newly malloc'd jz array has its interior
points filled with values but its boundary points are left
alone:

    for (i = 1; i <= nx-2; i++)
    {
        for (j = 1; j <= ny-2; j++)
        {
            // calculate jz at all points
            ...

In all subsequent functions the array jz is used whole, including
the uninitialized border.

I ran your test case in the debugger.  On n02, jz was 0 on the
border; on solar3, I got wild values

(gdb) p jz[0]
$16 = 1.37765386e+10
(gdb) p jz[1]
$17 = 1.35754556e-19
(gdb) p jz[2]
$18 = 1.35631939e-19

If the code relies on the border being 0, then it was pure luck
that it ran correctly on x86_64.

added lorentz force function to calculate keys TOTFX, TOTFY, TOTFZ, TOTBSQ, EPSX, EPSY, EPSZ

No code change. Added documentation to the header of the R function.

reverted back to v1.25

remedied an incorrect check

added R parameter (adapted from Schrijver 2007)

: 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