#include <cassert>
#include <vector>

extern "C" void create_z_grid( float* za,    int* nl,     int* ml,
	float* ztop,  float*  zbot,   float* deltaz )
{
	assert( *ztop > *zbot );
	
	std::vector<float> Z;
	Z.push_back( *ztop );
	while( Z.back()- *deltaz >= *zbot)  Z.push_back( Z.back()- *deltaz );
	*nl = (int)Z.size();
	assert( *nl <= *ml);

	std::copy(Z.begin(),Z.end(), za);

}