///
/// This file is part of Shalava.
///
/// Copyright (C) 2019 Noe Bernabeu <Noe.Bernabeu@gmail.com>
///
/// Shalava is free software; you can redistribute it and/or modify
/// it under the terms of the GNU General Public License as published by
/// the Free Software Foundation; either version 2 of the License, or
/// (at your option) any later version.
///
/// Shalava is distributed in the hope that it will be useful,
/// but WITHOUT ANY WARRANTY; without even the implied warranty of
/// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
/// GNU General Public License for more details.
///
/// You should have received a copy of the GNU General Public License
/// along with Shalava; if not, write to the Free Software
/// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
/// 
/// =========================================================================
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;