print_input.icc

derr << "south = " << south << endl;
derr << "est = " << est << endl;
derr << "north = " << north << endl;
derr << "west = " << west << endl;
derr << "reso = " << reso << endl;
derr << "rho = " << rho << endl;
derr << "mu0 = " << mu0 << endl;
derr << "tau = " << tau << endl;
derr << "nb_vents = " << nb_vents << endl;
for (size_t i=1; i<=nb_vents; i++){
	derr << "flow_rate_filename_"+itos(i)+" = " << flow_rate_filenames[i] << endl;
	derr << "Q_"+itos(i)+" = " << Q[i] << endl;
	derr << "R_"+itos(i)+" = " << R[i] << endl;
	derr << "xe_"+itos(i)+" = (" << point(xe[i])[0] << "," << point(xe[i])[1] << ")" << endl;
	derr << "ti_"+itos(i)+" = " << ti[i] << endl;
	derr << "tf_"+itos(i)+" = " << tf[i] << endl;
}
derr << "presence_of_trees = " << presence_of_trees << endl;
derr << "d = " << d << endl;
derr << "M = " << M << endl;
derr << "arrangement_of_trees = " << arrangement_of_trees << endl;
derr << "nb_regions = " << nb_regions << endl;
for (size_t i=1; i<=nb_regions; i++){
	derr << "nb_vertex_"+itos(i)+" = " << nb_vertex[i] << endl;
	for (size_t j=1; j<=nb_vertex[i]; j++){
	derr << "vertex_"+itos(i)+itos(j)+" = (" << point(vertex[i][j])[0] << "," << point(vertex[i][j])[1] << ")" << endl;
	}
}
derr << "final_shape = " << final_shape << endl;