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().c_str() << 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",

		      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);

  }

}

void Main::readLastStateXML()

{

  // 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 );

    next_file = xml_files.back();

    next_file = currentFile_path+"/"+next_file;

    readStateXML((const char*)next_file);

  }

}

void Main::readFirstStateXML()

{

  // 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 );

    next_file = xml_files.front();

    next_file = currentFile_path+"/"+next_file;

    readStateXML((const char*)next_file);

  }

}

QDir Main::GetLeafDir(void) {

  QDir LeafDir(QApplication::applicationDirPath()); 

  QStringList plugin_filters; // filter for plugins, i.e "*.dll", "*.dylib"

#if defined(Q_OS_WIN) 

  if (LeafDir.dirName().toLower() =="debug" 

      ||LeafDir.dirName().toLower() =="release") 

    LeafDir.cdUp(); 

  //plugin_filters << "*.dll";

#elif defined(Q_OS_MAC) 

  if (LeafDir.dirName() =="MacOS"){ 

    LeafDir.cdUp(); 

    LeafDir.cdUp(); 

    LeafDir.cdUp(); 

  }

#endif

  // for all OS-es. Move from "bin" directory to root application folder.

  if (LeafDir.dirName() == "bin") {

    LeafDir.cdUp();

  }

  LeafDir.cd("data/leaves");

  /* if (!LeafDir.cd("data/leaves")) {

    QString status_message = QString("No directory data/leaves");

    statusBar()->showMessage(status_message);

    return LeafDir;

  } 

  */

  return LeafDir;

}

void Main::readStateXML()

{

  //  extern Mesh mesh;

  stopSimulation();

#ifdef QDEBUG

  qDebug() << "Trying to open an OptionFileDialog" << endl;

#endif

  OptionFileDialog *fd = new OptionFileDialog( this, "read dialog", TRUE );

  fd->setMode( OptionFileDialog::ExistingFile );

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

  if (working_dir) {

    fd->setDir(*working_dir);

  } else {

    fd->setDir(par.datadir);

  }

  QString fileName;

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

    fileName = fd->selectedFile();

    if (working_dir) {

      delete working_dir;

    }

    working_dir = fd->dir();

    if (readStateXML((const char *)fileName,fd->readGeometryP(), fd->readParametersP()) )

      return readStateXML(); // user can try again

  }

}

void Main::clear()

{

  editor->clear();

}

void Main::about()

{

  static QMessageBox* about = new QMessageBox

       <p>\

         An Open Source framework for cell-based modeling of plant\

         tissue growth and development.\

       </p>\

       <p>(c) 2005-2008, Roeland Merks <i>et al.</i>\

         <a href=\"http://www.psb.vib-ugent.be\">VIB Department Plant Systems Biology</a>,\

         Ghent, Belgium.\

       </p>\

         <a href=\"http://www.cwi.nl/~merks\">Roeland Merks <i>et al.</i></a>\

         <a href=\"http://www.cwi.nl\">Centrum Wiskunde & Informatica (CWI)</a> and\

         <a href=\"http://www.ncsb.nl\"> Netherlands Consortium for Systems Biology (NCSB)</a>,\

         Amsterdam, the Netherlands.\

       </p>\

       <p>\

         VirtualLeaf 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 3 of the\

         License, or (at your option) any later version.\

       </p>\

       <p>\

         If you use this code for your projects, please cite our paper in\

         <a href=\"http://www.plantphysiol.org\">Plant Physiology</a>, &lsquo;\

         Roeland M. H. Merks, Michael Guravage, Dirk Inzé, and Gerrit T.S. Beemster,\

         <a href=\"http://www.plantphysiol.org/cgi/reprint/155/2/656.pdf\">\

         VirtualLeaf: an Open Source framework for cell-based modeling of plant tissue growth and development</a>,<br>\

         Plant Physiology 2011: 155:656-666.\

       </p>\

       <p>\

         Please share your model plugins and extensions at\

         <a href=\"http://virtualleaf.googlecode.com\">http://virtualleaf.googlecode.com</a>.\

       </p>",

      QMessageBox::Information, 1, 0, 0, this, 0, FALSE );

  about->setButtonText( 1, "Dismiss" );

  about->setIconPixmap(QPixmap( leaficon_small ));

  about->show();

}

void Main::gpl()

{

  static QMessageBox* gpl = new QMessageBox ( "GPL License", "", 

      QMessageBox::Information, 1, 0, 0, this, 0, FALSE );

  QDir docDir(QApplication::applicationDirPath());

  docDir.cd("../doc"); // Where Linux expects gpl3.txt

  QString path = docDir.filePath("gpl3.txt");

  if (!docDir.exists("gpl3.txt")){

    docDir = QApplication::applicationDirPath();

    docDir.cd("doc"); // Where Windows expects gpl3.txt

    path = docDir.filePath("gpl3.txt");

  }

  // At this point path points either to the linux variant, which

  // exists, or the windows variant, which may exist. Testing the

  // ifstream object will determine whether we've found gpl3.txt.

  std::ifstream file(path.toStdString().c_str());

  std::string str;

  if (file) {

    file.seekg(0, std::ios::end);   

    str.reserve(file.tellg());

    file.seekg(0, std::ios::beg);

    str.assign((std::istreambuf_iterator<char>(file)),

	       std::istreambuf_iterator<char>());

    gpl->setDetailedText(QString(str.c_str()));

  }

  gpl->setText(QString( "<h3>GNU GENERAL PUBLIC LICENSE</h3>"

			"<p>Version 3, 29 June 2007</p>"

			"<p>Copyright &copy; 2007 Free Software Foundation, Inc. "

			"&lt;<a href=\"http://fsf.org/\">http://fsf.org/</a>&gt;</p><p>"

			"Everyone is permitted to copy and distribute verbatim copies "

			"of this license document, but changing it is not allowed.</p>"

			"<h2>GNU GENERAL PUBLIC LICENSE</h2>"));

  gpl->setButtonText( 1, "Dismiss" );

  gpl->show();

}

void Main::aboutQt(){

  QMessageBox::aboutQt( this, "Virtual Leaf" );

}

void Main::toggleShowCellCenters()

{

  Plot();

}

void Main::toggleShowWalls()

{

  Plot();

}

void Main::toggleShowApoplasts()

{

  Plot();

}

void Main::toggleShowNodes()

{

  Plot();

}

void Main::toggleNodeNumbers(void)

{

  Plot();

}

void Main::toggleCellNumbers(void)

{

  Plot();

}

void Main::toggleCellAxes(void)

{

  Plot();

}

void Main::toggleCellStrain(void)

{

  Plot();

}

void Main::toggleShowFluxes(void)

{

  Plot();

}

void Main::toggleShowBorderCells()

{

  Plot();

}

void Main::toggleHideCells(void)

{

  Plot();

  editor->FullRedraw();

}

void Main::toggleMovieFrames(){}

void Main::toggleLeafBoundary(){}

void Main::toggleDynCells() {}

void Main::startSimulation(void)

{

  timer->start( 0 );

  statusBar()->message("Simulation started");

  running = true;

}

void Main::stopSimulation(void)

{

  timer->stop();

  cerr << "Stopping simulation" << endl;

  statusBar()->message("Simulation paused");

  running = false;

}

void Main::togglePaused()

{

  bool s = run->isItemChecked(paused_id);

  if (s) {

    cerr << "Calling start simulation" << endl;

    startSimulation();

  } else {

    cerr << "Calling stop simulation" << endl;

    stopSimulation();

  }

}

void Main::setFluxArrowSize(int size)

{

  flux_arrow_size = size/100.;

}

void Main::enlarge()

{

  canvas.setSceneRect( QRectF( 0,0, canvas.width()*4./3., canvas.height()*4./3.) );

}

void Main::shrink()

{

  canvas.setSceneRect( QRectF( 0,0, canvas.width()*3/4, canvas.height()*3/4) );

}

void Main::scale(double factor)

{

  QMatrix m = editor->matrix();

  m.scale(factor, factor);

  editor->setMatrix( m );

}

void Main::zoomIn()