float* average(float *inarr, long npts, float missing) { float ave, max, min, retarr[3]; long i, n; max = -1.0e32; min = 1.0e32; ave = 0.0; n = 0; for (i=0;i max) max = inarr[i]; ave = ave + inarr[i]; n++; } } if (n > 0) { retarr[0] = min; retarr[1] = ave/(float) n; retarr[2] = max; } else { for (i=0;i<3;i++) retarr[i]=missing; } return retarr; }