//#include <netcdfcpp.h>
#include <fstream.h>
#include <fcntl.h>

// void write_to_nc ( NcFile &nc, mat &X, std::string Xn, Array<std::string> A ) {
// 	char tmpstr[200];
// 	int Len = X.rows();
// 
// 	sprintf ( tmpstr,"%s.length",Xn.c_str() );
// 	NcDim* lengt = nc.add_dim ( tmpstr, ( long ) Len );
// 	for ( int j=0; j<X.cols(); j++ ) {
// 		if ( j<A.length() )
// 			sprintf ( tmpstr,"%s_%s",Xn.c_str(), ( A ( j ) ).c_str() );
// 		else
// 			sprintf ( tmpstr,"%s_%d",Xn.c_str(),j );
// 		// Create variables and their attributes
// 		NcVar* P = nc.add_var ( tmpstr, ncDouble, lengt );
// 		const double* Dp = X._data();
// 		P->put ( &Dp[j*Len],Len );
// 	}
// }


void dirfile_write (std::ofstream &form, std::string Dir, mat &X, std::string Xn, Array<std::string> A ) {
  char tmpstr[200];
  char file[200];
  const double* Dp=X._data();
  int Len = X.rows();
  
	for ( int j=0; j<X.cols(); j++ ) {
		if ( j<A.length() )
			sprintf ( tmpstr,"%s_%s",Xn.c_str(), ( A ( j ) ).c_str() );
		else
			sprintf ( tmpstr,"%s_%d",Xn.c_str(),j );
		// save files == from kst-doc-datasource
		sprintf(file,"%s/%s",Dir.c_str(),tmpstr);
		
		int fp;
		fp = open(file, O_CREAT | O_WRONLY | O_TRUNC, 00644);
		write(fp,&Dp[j*Len],Len*sizeof(double));
		close(fp);
  
  		form << tmpstr << " RAW  d 1" << endl;
	}
}