#include "stdlib.h" #include "amath.h" int* sort_m(float *values, int npts, float missing) { int *ordero, *orderi, *order, nptsi, i; float *valuei; orderi = calloc(sizeof(int),npts); order = calloc(sizeof(int),npts); valuei = calloc(sizeof(float),npts); nptsi = 0; for (i=0; i 0) { ordero = sort(valuei, nptsi); for (i=0; i < nptsi; i++) order[i] = orderi[ordero[i]]; } for (i=nptsi; i 0) free(ordero); return order; }