split write_cell_quantities_to_vtk function for vector<double>, added output of N_it results to mosfet_hcd