#include <iostream>

extern "C" void dump_initial( int* nl, int* ngas, int* igas, float* zt, 
	float* pt,  float* tt, float** qmix,int *iblk )
{
if( *iblk != 7) return;

	for(int i = 0; i< *nl; ++i)   {

		std::cout << // << pt[i] // qmix[0][i] 
                  tt[i] <<"   " << zt[i] ; // "  " <<   pt[i] << "  " << tt[i] ;
/*
		for( int j =0; j < *ngas ; ++j ) {
			std::cout << "  " << igas[j] << " " << qmix[j][i] ;
		}
*/
		std::cout << std::endl;

	}
//exit(23);
		std::cout << "&" << std::endl;
}