2010-11-30    <guravage@petitdru.sen.cwi.nl>

	* canvas.cpp (snapshot): Preface each image format extension with

	an asterisk so files with those extensions to appear in the file

	dialogue.

	* VirtualLeaf-install-windows.nsi: Include imageformats folder

	alongside executable.

	* canvas.cpp (snapshot): PNG is the default image format.

2010-11-29    <guravage@petitdru.sen.cwi.nl>

	* canvas.cpp (snapshot): Query and display supported snapshot file

	formats.

	* VirtualLeaf-install-windows.nsi: To facilitate debugging I added

	the debug versions of all the windows DLLs.

2010-10-19    <guravage@caterpie.sen.cwi.nl>

	* mainbase.cpp (Save): Use format specified in function prototype.

	* VirtualLeaf.cpp (Plot): Render PNG instead of JPEG. Write image

	and XML files dependent only on getTime().

2010-10-18    <guravage@caterpie.sen.cwi.nl>

	* mainbase.cpp (Save): QDir::toNativeSeparators(fname). Invoke

	save with format argument set to zero to force QImage to guess the

	format by looking at fileName's suffix.

	* VirtualLeaf-install.nsi: Put uninstaller in top directory.

	* modelcatalogue.cpp (InstallModel): For all OS-es. Move from "bin" directory to root application folder.

	* VirtualLeaf.pro: MAXOSX, look for gpl3.txt in the ../doc directory.

	* Makefile (Makefile.libplugin): Removed -makefile qmake option.

	* VirtualLeaf.cpp (TIMESTEP): Removed getIterations().

	* mesh.h (Mesh): Removed {increment,get}Iterations().

	* canvas.cpp (TimeStepWrap): Replaced getIterations() with getTime().

2010-10-15    <guravage@caterpie.sen.cwi.nl>

	* mesh.h (Mesh): Added iterations. incrementIterations() and

	getIterations().

	* VirtualLeaf.cpp (Plot): Replaced local frame counter with

	mesh.incrementIterations().

	* canvas.cpp (TimeStepWrap): Replaced local counter with

	mesh.incrementIterations().

2010-10-14    <guravage@caterpie.sen.cwi.nl>

	* VirtualLeaf.pro: Turned debug off to make all profiles must be

	consistant.

	* VirtualLeaf-install.nsi: Tweaked paths to coincide with UNIX

	distribution, i.e. VirtualLeaf.exe, its libraries, models, and the

	uninstaller all go in the bin directory. And leaves directory

	placed under data directory.

	* canvas.cpp (exportCellData): Added a check to inquire before

	overwritting an existing file.

2010-10-14    <merks@cwi.nl>

	* mesh.cpp: In response to referees' comments, added new parameter

	"yield_threshold" instead of fixed parameter '4' for yield

	threshold.

	* wallitem.cpp: when clicking a wall, both the wall type was

	cycled and the transporterdialog popped up. Corrected this - for

	wall type cycling, hold the Control button while left

	clicking. TransporterDialog only pops up for left click. Also made

	sure the wall is redrawn after changing the transporter values.

	* cell.cpp (DivideWalls): accomodated for rename of Circumference -> WallCircumference

	* hull.h: added an operator< to sort Points

	* hull.cpp: added an operator< to sort Points

	* cellbase.cpp (ExactCircumference): I added a new function

	ExactCircumference, yielding the circumference of the cell along

	its wall_elements

	* VirtualLeaf.cpp: adjust info_string to accomodate for new name

	of function CellBase::Circumference -> CellBase::WallCircumference

	* mesh.cpp: corrected Mesh::Compactness, the boundary coordinates

	need to be sorted in x,y order for the convex hull algorithm

	(thanks Margriet!). I updated CSVExportCellData so it exports the

	circumferences of hull and boundary_polygon.

2010-10-14    <guravage@caterpie.sen.cwi.nl>

	* VirtualLeafpar.tmpl (export_fn_prefix): changed to 'cell.'

	* canvas.cpp (TimeStepWrap): Removed TimeStamp(). Cell data

	filename now incorporates iteration number.

	* canvas.h (MainBase): Removed TimeStamp().

2010-10-08    <guravage@caterpie.sen.cwi.nl>

	* pardialog.h:

	* pardialog.cpp:

	* parameter.h:

	* parameter.cpp: Regenerated to include export_interval and export_fn_prefix.

	* VirtualLeafpar.tmpl: Appended export_interval and export_fn_prefix.

	* canvas.h (MainBase): Declared polymorphic exportCellData() functions.

	* canvas.cpp:

	(TimeStamp): New private TimeStamp() function.

	(TimeStepWrap): Added invocation of exportCellData().

	(exportCellData): Created two polymorphic functions: one with a

	single QString argument, the other with no argument. The former is

	called from TimeStepWrap() while the latter is called from the

	"Export cell areas" item in the file menu.

2010-10-07    <guravage@caterpie.sen.cwi.nl>

	* canvas.cpp (exportCellData): Added a Q3FileDialog to inquire

	where to write the exportCellData.

2010-06-28    <guravage@caterpie.sen.cwi.nl>

	* VirtualLeaf-install.nsi: Grab gpl3.txt from doc directory.

	* canvas.cpp (gpl): gpl3.txt can be either in an ancestor doc

	directory (Linux) or a decedent doc directory (Windows, via the

	binary installer).

	* VirtualLeaf-install.nsi: Add VirtualLeaf doc directory.

2010-06-25    <guravage@caterpie.sen.cwi.nl>

	* gpl3.txt: Moved gpl3.txt from doc to src directory.

	* VirtualLeaf.pro: Added -Wno-write-strings and -Wno-unused-parameter to QMAKE_CXXFLAGS.

	* libplugin.pro: Ditto.

	* parameter.cpp: Result of adding datadir changes to make_parameter_source.pl.

	* parameter.h: Ditto.

	* output.h: Declared new function (AppendHomeDirIfPathRelative).

	* output.cpp (AppendHomeDirIfPathRelative): Added new function.

	* canvas.cpp (gpl): Moving gpl3.txt from doc to src obviates the need to docDir.cd("../doc").

	* VirtualLeaf-install.nsi: Grab gpl3.txt from src directory.

	Add missing libiconv/bin, libxml2bin and libz/bin directories.

	Copy libiconv-2.dll, libxml2.dll and zlib1.dll from relative paths.

	* VirtualLeaf.pro: copy gpl3.txt as part of QMAKE_POST_LINK.

2010-06-24    <guravage@caterpie.sen.cwi.nl>

	* libplugin.pro: Use correct library path.

	* VirtualLeaf.pro: Ditto.

	* VirtualLeaf.cpp (DrawCell): Iterate over NChem to construct info_string.

2010-06-23    <guravage@caterpie.sen.cwi.nl>

	* simitembase.cpp: Removed NULL assignments to unused variables.

	* VirtualLeaf.cpp: Ditto.

	* apoplastitem.cpp: Ditto.

	* canvas.cpp: Ditto.

	* cell.cpp: Ditto.

	* cellbase.h: Ditto.

	* forwardeuler.cpp: Ditto.

	* mainbase.h: Ditto.

	* nodeitem.cpp: Ditto.

	* qcanvasarrow.h: Ditto.

	* simitembase.cpp: Ditto.

	* Makefile (clean): Add -f Makefile argument to each make invocation.

	* VirtualLeaf-install.nsi: New gpl license text.

	* VirtualLeaf.pro: Disabled console mode.

	* mesh.cpp (Clear): Added parentheses to qDebug statments.

	(TestIllegalWalls): Replaced qDebug().

	* canvas.cpp (mouseReleaseEvent): Replaced qDebug() with cerr since qDebug complains about *node_set.

	* wall.cpp (CorrectWall): Rplaced gDebug() with cerr in transform call and when printing *this.

2010-06-22    <guravage@caterpie.sen.cwi.nl>

	* Makefile (tutorials): Add tutorials target.

2010-06-21    <guravage@caterpie.sen.cwi.nl>

	* parameter.cpp: Added particular reassignment of datadir.

	* canvas.cpp (gpl): Added GPL3 License text. Display detail text only if the source text file exists.

2010-06-18    <guravage@caterpie.sen.cwi.nl>

	* canvas.cpp (gpl): Added gpl slot to display GPL license.

	* VirtualLeaf.pro: Changed default LIBXML2DIR, LIBICONVDIR and LIBZDIR to corresponding distribution lib directories.

	* libplugin.pro: Ditto.

	* Makefile (clean): add if stmt not to `touch` on windows.

2010-06-17    <guravage@caterpie.sen.cwi.nl>

	* VirtualLeaf.pro: Removed perl references.

	* libplugin.pro: Ditto.

2010-06-15    <guravage@caterpie.sen.cwi.nl>

	* VirtualLeaf.pro: Removed xmlwritecode.cpp from SOURCES list.

	* xmlwrite.cpp (XMLSave): Removed references to XMLWriteLeafSourceCode and XMLWriteReactionsCode.

	* xmlwrite.h (XMLIO): Ditto!

	* mesh.cpp (findNextBoundaryNode): Initialize Node *next_boundary_node = NULL;

	* xmlwrite.cpp (XMLReadSimtime): Removed unused variable cur

	(XMLReadWalls): viz_flux need not be declared twice; default value of 0.0.

	(XMLReadCells): Removed unused count variable.

	(XMLReadSimtime): Removed unused cur variable.

	(XMLRead): Removed unused v_str variable.

	* simitembase.cpp (userMove): Use assignment merely to obviate compilation warning.

	(SimItemBase) Ditto.

	* qcanvasarrow.h (QGraphicsLineItem): Use assignment merely to obviate compilation warning.

	* output.cpp (OpenWriteFile): Removed unused par variable.

	* nodeitem.cpp (paint): Use assignment merely to obviate compilation warning.

	* forwardeuler.cpp (odeint): Use assignment merely to obviate compilation warning.

	* cell.cpp (DivideOverGivenLine): Use assignment merely to obviate compilation warning.

	* canvas.cpp (FigureEditor): Use assignments merely to obviate compilation errors.

	(mousePressEvent): Removed unused item variable.

	* apoplastitem.cpp

	(ApoplastItem): Removed unused par variable.

	(OnClick): Use NULL assignment merely to obviate compilation warning.

	* mainbase.h (MainBase): Use assignment merely to obviate compilation warning.

	* cellbase.h (CellsStaticDatamembers): Use assignment merely to obviate compilation warning.

	* cell.cpp: Wrapped diagnostic output in QDEBUG blocks.

	* VirtualLeaf.cpp ditto.

	* canvas.cpp ditto.

	* cell.cpp ditto.

	* data_plot.cpp ditto.

	* forwardeuler.cpp ditto.

	* mesh.cpp ditto.

	* mesh.h

	* random.cpp ditto.

	* wall.cpp ditto.

	* wallbase.cpp ditto.

	* wallitem.cpp ditto.

2010-06-07    <guravage@caterpie.sen.cwi.nl>

	* VirtualLeaf.pro: Removed explicit perl invocation to regerenerate parameter files.

	* libplugin.pro: ditto.

2010-06-03    <guravage@caterpie.sen.cwi.nl>

	* pardialog.h: Added default versions of this automatically generated file.

	* pardialog.cpp: ditto.

	* parameter.h: ditto.

	* parameter.cpp: ditto.

	* VirtualLeaf.pro: delete/generate  parameter.{h,cpp}and pardialog.{h,cpp} only if perl is installed.

 	* libplugin.pro: dito.

	* Makefile: Added top-level Makefile

2010-05-10    <guravage@caterpie.sen.cwi.nl>

	* VirtualLeaf.pro: Added -fPIC option to QMAKE_CXXFLAGS.

	mesh.TestIllegalWalls();

	mesh.RepairBoundaryPolygon();

#ifdef QDEBUG

	qDebug() << "Done RepairBoundaryPolygon" << endl;

#endif

	mesh.TestIllegalWalls();

	mesh.CleanUpWalls();

#ifdef QDEBUG

	qDebug() << "Done CleanUpWalls" << endl;

#endif

	mesh.TestIllegalWalls();

      }

      dynamic_cast<Main *>(parent())->Plot();

#ifdef QDEBUG

      cerr << "NodeSet of cutting line: " << *node_set << endl;

#endif

    }

  } else 

    if (e->button()==Qt::RightButton) {

      if (intersection_line /* && !angle_line*/) {

	QPointF p = mapToScene(e->pos());

	QPointF sp = intersection_line->line().p1();

	viewport()->setMouseTracking( TRUE );

      } 

    }

}

// returns a vector of pointer to cells colliding with intersection line

vector <CellItem *> FigureEditor::getIntersectedCells(void)

{ 

  vector <CellItem *> colliding_cells;

  QList<QGraphicsItem *> l = intersection_line->collidingItems( );

#ifdef QDEBUG

  qDebug() <<  "l.size() = " << l.size() << endl;

#endif

  for (QList<QGraphicsItem *>::Iterator it=l.begin(); it!=l.end(); ++it) {

#ifdef QDEBUG

    qDebug() << typeid(**it).name() << endl;

#endif

    if ( !strcmp(typeid(**it).name(),"8CellItem") ) {

      colliding_cells.push_back(dynamic_cast<CellItem *>(*it));

    }

  }

  delete intersection_line;

  intersection_line = 0;

  return colliding_cells;

}

void FigureEditor::FullRedraw(void)

{

  QList<QRectF> rl;

  rl << sceneRect();

  updateScene(rl);   

}

NodeItem *FigureEditor::selectedNodeItem(QList<QGraphicsItem *> graphicItems) const

{

  NodeItem *nodeItem = 0;

  // graphicItems is a list of the QgraphicItems under the mouse click event 

  QList<QGraphicsItem *>::Iterator it = graphicItems.begin();

  for (; it != graphicItems.end(); ++it) {

    if ( !strcmp(typeid(**it).name(),"8NodeItem")) {

      // the object under the mouse click is a Nodeitem

      nodeItem = dynamic_cast<NodeItem *>(*it);

      // indicate we've selected it

      nodeItem->setBrush(dark_red);

      break;

    }

  }

  return nodeItem;

}

void FigureEditor::insertNode(QPointF p)

{

  Node *node = new Node(p.x(), p.y(), 0);

  mesh.AddNode(node);

  scene()->clearSelection();

  dynamic_cast<Main *>(parent())->Plot();

  FullRedraw();

#ifdef QDEBUG  

  qDebug() << "Node: " << p << endl;

#endif

}

static uint mainCount = 0;

Main::Main(QGraphicsScene& c, Mesh &m, QWidget* parent, const char* name, Qt::WindowFlags f) :

  Q3MainWindow(parent,name,f),

  MainBase(c,m),

  mesh(m)

{

  editor = new FigureEditor(canvas,mesh, this);

#ifdef QDEBUG

  qDebug() << "Interactive = " << editor->isEnabled();

#endif

  working_dir = 0;

  QObject::connect( editor, SIGNAL(MousePressed()), this, SLOT(PauseIfRunning()));

  QObject::connect( editor, SIGNAL(MouseReleased()), this, SLOT(ContIfRunning()));

  QMenuBar* menu = menuBar();

  Q3PopupMenu* file = new Q3PopupMenu( menu );

  file->insertItem("&Read leaf", this, SLOT(readStateXML()));

  file->insertItem("&Save leaf", this, SLOT(saveStateXML()));

  file->insertItem("Snapshot", this, SLOT(snapshot()), Qt::CTRL+Qt::SHIFT+Qt::Key_S);

  file->insertSeparator();

  file->insertItem("Read next leaf", this, SLOT(readNextStateXML()), Qt::Key_PageDown);

  file->insertItem("Read previous leaf", this, SLOT(readPrevStateXML()), Qt::Key_PageUp);

  file->insertItem("Read last leaf", this, SLOT(readLastStateXML()), Qt::Key_End);

  file->insertItem("Read first leaf", this, SLOT(readFirstStateXML()), Qt::Key_Home);

  file->insertItem("Export cell data", this, SLOT(exportCellData()));

  file->insertSeparator();

  file->insertItem("&Print...", this, SLOT(print()), Qt::CTRL+Qt::Key_P);

  file->insertSeparator();

  file->insertItem("E&xit", qApp, SLOT(quit()), Qt::CTRL+Qt::Key_Q);

  menu->insertItem("&File", file);

  Q3PopupMenu* edit = new Q3PopupMenu( menu );

  edit->insertItem("Reset Chemicals and Transporters", this, SLOT( CleanMesh()), Qt::CTRL+Qt::Key_R );

  edit->insertItem("Reset Chemicals", this, SLOT( CleanMeshChemicals()) );

  edit->insertItem("Reset Transporters", this, SLOT( CleanMeshTransporters()) );

  edit->insertItem("Randomize PIN1 Transporters", this, SLOT( RandomizeMesh()) );

  edit->insertItem("Cut away SAM", this, SLOT( CutSAM() ));

  menu->insertItem("&Edit", edit);

  run = new Q3PopupMenu( menu );

  running = false;

  paused_id = run->insertItem("&Simulation paused", this, SLOT(togglePaused()), Qt::Key_S);

  run->setItemChecked(paused_id, FALSE);

  menu->insertItem("&Run", run);

  view = new Q3PopupMenu( menu );

  view->insertItem("&Zoom in", this, SLOT(zoomIn()), Qt::CTRL+Qt::Key_Equal);

  view->insertItem("Zoom &out", this, SLOT(zoomOut()), Qt::CTRL+Qt::Key_Minus);

  view->insertSeparator();

  com_id = view->insertItem("Show cell &centers", this, SLOT(toggleShowCellCenters()));

  view->setItemChecked(com_id, FALSE);

  mesh_id = view->insertItem("Show &nodes", this, SLOT(toggleShowNodes()), Qt::CTRL+Qt::SHIFT+Qt::Key_N);

  view->setItemChecked(mesh_id, TRUE);

  node_number_id = view->insertItem("Show node numbers", this, SLOT(toggleNodeNumbers()), Qt::CTRL+Qt::SHIFT+Qt::Key_M);

  view->setItemChecked(node_number_id, FALSE);

  cell_number_id = view->insertItem("Show cell numbers", this, SLOT(toggleCellNumbers()));

  view->setItemChecked(cell_number_id, FALSE);

  hide_cells_id = view->insertItem("Hide cells", this, SLOT(toggleHideCells()));

  view->setItemChecked(hide_cells_id, FALSE);

  border_id = view->insertItem("Show &border cells", this, SLOT(toggleShowBorderCells()));

  view->setItemChecked(border_id, FALSE);

  cell_axes_id = view->insertItem("Show cell &axes", this, SLOT(toggleCellAxes()));

  cell_strain_id = view->insertItem("Show cell &strain", this, SLOT(toggleCellStrain()));

  view->setItemChecked(cell_axes_id, FALSE);

  fluxes_id = view->insertItem("Show &fluxes", this, SLOT(toggleShowFluxes()));

  view->setItemChecked(fluxes_id, FALSE);

  cell_walls_id = view->insertItem("Show transporters", this, SLOT(toggleShowWalls()));

  view->setItemChecked(cell_walls_id, FALSE);

 // apoplasts_id = view->insertItem("Show apoplasts", this, SLOT(toggleShowApoplasts()));

 // view->setItemChecked(apoplasts_id, FALSE);

  view->insertSeparator();

  only_boundary_id = view->insertItem("Show only leaf &boundary", this, SLOT(toggleLeafBoundary()));

  view->insertSeparator();

  movie_frames_id = view->insertItem("Start saving movie &frames", this, SLOT(toggleMovieFrames()));

  view->setItemChecked(movie_frames_id, par.movie);

  view->setItemChecked(only_boundary_id, FALSE);

  menu->insertItem("&View", view);

  options = new Q3PopupMenu( menu );

  dyn_cells_id = options->insertItem("Cell growth", this, SLOT(toggleDynCells()));

  options->setItemChecked(dyn_cells_id, true);

  options->insertItem("Edit &parameters", this, SLOT(EditParameters()), Qt::CTRL+Qt::Key_E);

  rotation_mode_id = options->insertItem("Rotate leaf", this, SLOT(EnterRotationMode()), Qt::CTRL + Qt::SHIFT + Qt::Key_R);

  options->setItemChecked(rotation_mode_id, false);

  menu->insertItem("&Options",options);

  // Menu of models

  modelmenu = new QMenu( menu );

  menu->insertItem("&Models", modelmenu);

  menu->insertSeparator();

  helpmenu = new Q3PopupMenu( menu );

  tooltips_id = helpmenu->insertItem("Show Cell&Info", this, SLOT(Refresh()));

  helpmenu->setItemChecked(tooltips_id, true);

  helpmenu->insertSeparator();

    //helpmenu->insertSeparator();

  helpmenu->insertItem("&LICENSE", this, SLOT(gpl()) );

  helpmenu->insertItem("About", this, SLOT(about()) ); //, Key_F1);

	menu->insertItem("&Help",helpmenu);

  statusBar();

  setCentralWidget(editor);

  printer = 0;

  init();

  // Start timer which repetitively invokes

  // a simulation time step

  timer = new QTimer( this );

  connect( timer, SIGNAL(timeout()), SLOT(TimeStepWrap()) );

  stopSimulation();

  statusBar()->addWidget(new QLabel("Ready."));

  setCaption(caption);

  gifanim = 0;

  infobar = new InfoBar();

  addDockWindow(infobar);

}

void Main::RefreshInfoBar(void)

{

  infobar->SetText(mesh.ModelID());

}

void Main::UserMessage(QString message, int timeout)

{

  statusBar()->showMessage(message, timeout);

}

void Main::init()

{

  clear();

  static int r=24;

  srand(++r);

  mainCount++;

}

Main::~Main()

{

  delete printer;

  if ( !--mainCount ) {

  }

  //EndGifAnim();

}

void Main::newView()

{

  // Open a new view... have it delete when closed.

  Main *m = new Main(canvas, mesh, 0, 0, Qt::WDestructiveClose);

  qApp->setMainWidget(m);

  m->show();

  qApp->setMainWidget(0);

}

void Main::EditParameters()

{

  ParameterDialog *pardial = new ParameterDialog(this, "stridediag");

  // Make sure the values in the parameter dialog are updated after a file is read 

  // each method changing the parameters (reading XML or PAR files) should

  // emit this signal

  QObject::connect(   this, SIGNAL( ParsChanged() ), pardial, SLOT( Reset() ) );

}

void Main::savePars()

{

  stopSimulation();

  Q3FileDialog *fd = new Q3FileDialog( this, "file dialog", TRUE );

  fd->setMode( Q3FileDialog::AnyFile );

  fd->setFilter( "Parameter files (*.par)");

  fd->setDir(par.datadir);

  QString fileName;

  if ( fd->exec() == QDialog::Accepted ) {

    fileName = fd->selectedFile();

    ofstream parfile((const char *)fileName);

    par.Write(parfile);

  }

  startSimulation();

}

void Main::readPars()

{

  stopSimulation();

  Q3FileDialog *fd = new Q3FileDialog( this, "file dialog", TRUE );

  fd->setMode( Q3FileDialog::ExistingFile );

  fd->setFilter( "Parameter files (*.par)");

  fd->setDir(par.datadir);

  QString fileName;

  if ( fd->exec() == QDialog::Accepted ) {

    fileName = fd->selectedFile();

    par.Read((const char *)fileName);

  }

  emit ParsChanged();

}

void Main::saveStateXML()

{

  stopSimulation();

  Q3FileDialog *fd = new Q3FileDialog( this, "file dialog", TRUE );

  fd->setMode( Q3FileDialog::AnyFile );

  fd->setFilter( "LeafML files (*.xml)");

  fd->setDir(par.datadir);

  QString fileName;

  if ( fd->exec() == QDialog::Accepted ) {

    fileName = fd->selectedFile();

    // extract extension from filename

    QFileInfo fi(fileName);

    QString extension = fi.suffix();

    if (extension.isEmpty()) {

      extension = "xml";

      fileName += ".";

      fileName += extension;

    }

    if ( QFile::exists( fileName ) &&

	 QMessageBox::question(

			       this,

			       tr("Overwrite File? -- Leaf Growth"),

			       tr("A file called %1 already exists."

				  " Do you want to overwrite it?")

			       .arg( fileName ),

			       tr("&Yes"), tr("&No"),

			       QString::null, 1, 1 ) ) {

      return saveStateXML();

    } else {

      mesh.XMLSave((const char *)fileName, XMLSettingsTree());

      QString status_message;

      status_message = QString("Wrote LeafML to %1").arg(fileName);

      cerr << status_message.toStdString() << endl;

      statusBar()->showMessage(status_message);

    }

  }

}

void Main::snapshot()

{

#ifdef QDEBUG

  qDebug() << "Supported Image Formats: " << QImageWriter::supportedImageFormats() << endl;

#endif

  stopSimulation();

  Q3FileDialog *fd = new Q3FileDialog( this, "Save snapshot", TRUE );

  fd->setDir(par.datadir);

  fd->setMode( Q3FileDialog::AnyFile );

  foreach (QString format, QImageWriter::supportedImageFormats()){

    supported_file_formats += (" *." + format);

  }

  fd->setFilter("Image files (" + supported_file_formats + " )");

  QString fileName;

  if ( fd->exec() == QDialog::Accepted ) {

    fileName = fd->selectedFile();

    // extract extension from filename

    QFileInfo fi(fileName);

    QString extension = fi.suffix();

    if (extension.isEmpty()) {

      extension = "png";

      fileName += ".";

      fileName += extension;

    }

    if ( QFile::exists( fileName ) &&

	 QMessageBox::question(

			       this,

			       tr("Overwrite file? -- Cell data export"),

			       tr("A file called %1 already exists."

				  " Do you want to overwrite it?").arg( fileName ),

			       QMessageBox::Yes, QMessageBox::No 

			       ) == QMessageBox::No

	 ) {

      return snapshot();

    } else {

      // Save bitmaps at 1024x768

      QString status_message;

      if (Save((const char *)fileName, extension, 1024, 768)==0) {

	status_message = QString("Wrote snapshot to %1").arg(fileName);

      } else {

	status_message = QString("Error writing snapshot to %1").arg(fileName);

      }

      cerr << status_message.toStdString() << endl;

      statusBar()->showMessage(status_message);

    }

  }

}

void Main::readPrevStateXML()

{

  // if we have already read a file, read the next file

  if (!currentFile.isEmpty() && working_dir) {

    QString next_file;

    QStringList xml_files = working_dir->entryList("*.xml");

    QString currentFile_nopath = currentFile.section( '/', -1 );

    QString currentFile_path = currentFile.section( '/', 0, -2 );

    QList<QString>::iterator f = xml_files.find( currentFile_nopath );

    if (f == xml_files.end()) {

      return;

    }

    if (f==xml_files.begin()) {

      QMessageBox mb( "Read previous leaf",

		      "No more files",

		      QMessageBox::Information,

		      QMessageBox::Ok | QMessageBox::Default,

		      QMessageBox::NoButton,

		      QMessageBox::NoButton);

      mb.exec();

      return;

    }

    next_file = *(--f);

    next_file = currentFile_path+"/"+next_file;

    readStateXML((const char *)next_file);

  }

}

int Main::readStateXML(const char *filename, bool geometry, bool pars, bool simtime)

{

  try {

    xmlNode *settings;

    mesh.XMLRead((const char *)filename, &settings, geometry, pars, simtime);

#ifdef QDEBUG

    qDebug() << "Reading done."<< endl;

#endif

    XMLReadSettings(settings);

    xmlFree(settings);

    Cell::SetMagnification(1);

    Cell::setOffset(0,0);

    FitLeafToCanvas();

    currentFile =  QString(filename);

    Plot();

    QString status_message = QString("Successfully read leaf from file %1. Time is %2 h.").arg(currentFile).arg(mesh.getTimeHours().c_str());

    cerr << status_message.toStdString() << endl;

    setCaption(caption_with_file.arg(filename));

    statusBar()->message(status_message);

    emit ParsChanged();

#ifdef QDEBUG

    qDebug() << "Done. Returning 0." << endl;

#endif

    return 0;

  } catch (const char *error_message) {

    QMessageBox mb( "Read leaf from XML file",

		    QString(error_message),

		    QMessageBox::Critical,

		    QMessageBox::Ok | QMessageBox::Default,

		    Qt::NoButton,

		    Qt::NoButton);

    mb.exec();

    return 1;

  }

}

void Main::readNextStateXML()

{

  // if we have already read a file, read the next file

  if (!currentFile.isEmpty() && working_dir) {

    QString next_file;

    QStringList xml_files = working_dir->entryList("*.xml");

    QString currentFile_nopath = currentFile.section( '/', -1 );

    QString currentFile_path = currentFile.section( '/', 0, -2 );

    QList<QString>::iterator f = xml_files.find( currentFile_nopath );

    if (f == xml_files.end()) {

      return;

    }

    ++f;

    if (f==xml_files.end()) {

      QMessageBox mb( "Read next leaf",

		      "No more files",