 }
 
 /**
-* \brief Write numerical data at a given iteration number.
-*
-* It creates a directory whose name is the iteration number
-* and stores the fluxes, the pressures (and potentially the body forces)
-* into the files \a q, \a lambda, respectively.
-*
-* \param caseFolder the path of the folder case
-* \param n the iteration number
-* \param q Cusp array that contains the fluxes
-* \param lambda Cusp array that contains the pressures (and potentially the body forces)
-* \param D information about the computational grid
-*/
+ * \brief Writes numerical data at a given iteration number.
+ *
+ * It creates a directory whose name is the iteration number
+ * and writes the flux, the pressure (and potentially the body forces)
+ * into the files \a q, \a lambda, respectively.
+ *
+ * \param caseFolder directory of the simulation
+ * \param n the iteration number
+ * \param q Cusp array that contains the fluxes
+ * \param lambda Cusp array that contains the pressures (and potentially the body forces)
+ * \param D information about the computational grid
+ */
 template <>
 void writeData<vecD>(std::string &caseFolder, int n, vecD &q, vecD &lambda, domain &D)//, bodies &B)
 {

//##############################################################################

/**
* \brief Read inputs by parsing the command-line and simulation files
*
* \param DB database that contains all the simulation parameters
* \param D contains information about the computational grid
*/
void readInputs(int argc, char **argv, parameterDB &DB, domain &D)
{
	// get a default database

}

/**
 * \brief Writes the mesh points into a file named \a grid.
 *
 * \param caseFolder the path of the case folder
 * \param D information about the computational grid
 */
void writeGrid(std::string &caseFolder, domain &D)
{
	std::stringstream out;
