#include "Level1Data.h"

#include <cassert>

extern "C" void get_l1_za( int* ml,   int*  nl,
	int* ileave,  int* nleave,  float*  za, float *ztop)
{
//std::cout << " in getl1za " << std::endl;
	std::vector<double> altgrid =  ::gLevel1Data.GetImpactAlts() ;
	assert(*ileave >=1 && *nleave >= *ileave && *nleave <= (int)altgrid.size() );

	*nl= 1+  ( (int)altgrid.size() - *ileave)/ *nleave ;

	assert( *nl <= *ml );

    *nl = 0;
	for( unsigned il= *ileave-1; il< altgrid.size(); il+= *nleave) {
        if((float)altgrid[il] > *ztop) continue;
		za[ *nl ] =(float)altgrid[il];
		*nl += 1;
    }
//std::cout << " out of getl1za " << std::endl;
}