1 #include "global.h"
2 
get_aver_elev(struct Ortho_Control_Points * cpz,double * aver_z)3 int get_aver_elev(struct Ortho_Control_Points *cpz, double *aver_z)
4 {
5     double meanz;
6     double *cp = cpz->z2;
7     int i, n;
8 
9     /*  Need 1 control points */
10     if (cpz->count <= 0) {
11 	return (-1);
12     }
13 
14     /* set average elevation from mean values of CONZ points */
15     meanz = 0;
16     n = 0;
17     for (i = 0; i < cpz->count; i++) {
18 	if (cpz->status[i] <= 0)
19 	    continue;
20 
21 	n++;
22 	meanz += *(cp++);
23 	G_debug(3, "In ortho meanz = %f", meanz);
24     }
25 
26     *aver_z = meanz / n;
27 
28     G_debug(1, "In ortho aver_z = %f", *aver_z);
29 
30     return 0;
31 }
32