WallReaction *wall_f = new WallDynamics();

    mesh.ReactDiffuse_New(transport_f, cellreaction_f, wall_f, par.rd_dt);*/

				mesh.ReactDiffuse(par.rd_dt);

    Plot(par.resize_stride);

  }

  i++;

  return mesh.getTime();

}

/* Called if a cell is clicked */

void Cell::OnClick(QMouseEvent *e) {

  /* #ifdef HAVE_QWT

	// Launch DataPlot window

	QStringList curvenames;

    for (int i=0;i<NChem();i++) {

		curvenames += QString("Chem #%1").arg(i);

	}

    PlotDialog *plot = new PlotDialog((Main *)main_window, QString("Monitor for Cell %1").arg(Index()), curvenames);

    QObject::connect(this, SIGNAL(ChemMonValue(double, double *)),

                     plot, SLOT(AddValue(double,double *)));

#endif

	*/

	//getMesh().plugin->OnClick(*this);

}

/* Custom message handler - Default appends a newline character to the end of each line. */ 

 void vlMessageOutput(QtMsgType type, const char *msg)

 {

     switch (type) {

     case QtDebugMsg:

       //fprintf(stderr, "Debug: %s\n", msg);

       cerr << msg << flush;

       break;

     case QtWarningMsg:

       //fprintf(stderr, "Warning: %s\n", msg);

       cerr << "Warning: " << msg << flush;

       break;

     case QtCriticalMsg:

       fprintf(stderr, "Critical: %s\n", msg);

       cerr << "Critical: " << msg << flush;

       break;

     case QtFatalMsg:

       //fprintf(stderr, "Fatal: %s\n", msg);

       cerr << "Fatal: " << msg << flush;

       abort();

     }

 }

Parameter par;

int main(int argc,char **argv) {

  try {

    int c;

    char *leaffile=0;

    char *modelfile=0;

    while (1) {

      //int this_option_optind = optind ? optind : 1;

      int option_index = 0;

      static struct option long_options[] = {

	{"batch", no_argument, NULL, 'b'},

	{"leaffile", required_argument, NULL, 'l'},

	{"model", required_argument, NULL, 'm'} 

      };

      // short option 'p' creates trouble for non-commandline usage on MacOSX. Option -p changed to -P (capital)

      static char *short_options = "blm";

      c = getopt_long (argc, argv, "bl:m:",

		       long_options, &option_index);

      if (c == -1)

	break;

      if (c==0) {

	printf ("option %s", long_options[option_index].name);

	if (optarg)

	  printf (" with arg %s", optarg);

	printf ("\n");

	c = short_options[option_index];

      }

      switch (c) {

      case 'b':

	cerr << "Running in batch mode\n";

	batch=true;

	break;

      case 'l':

	leaffile=strdup(optarg);

	if (!leaffile) {

	  throw("Out of memory");

	}

	printf("Reading leaf state file '%s'\n", leaffile);

	break;

      case 'm':

	modelfile=strdup(optarg);

	if (!modelfile) {

	  throw("Out of memory");

	}

	break;

      case '?':

	break;

      default:

	printf ("?? getopt returned character code 0%o ??\n", c);

      }

    }

    if (optind < argc) {

      printf ("non-option ARGV-elements: ");

      while (optind < argc)

	printf ("%s ", argv[optind++]);

      printf ("\n");

    }

    bool useGUI = !batch;

    qInstallMsgHandler(vlMessageOutput); // custom message handler

    QApplication app(argc,argv,useGUI);

    QPalette tooltippalette = QToolTip::palette();

    QColor transparentcolor = QColor(tooltippalette.brush(QPalette::Window).color());

    tooltippalette.setBrush (QPalette::Window, QBrush (transparentcolor) );

    QToolTip::setPalette( tooltippalette );

    QGraphicsScene canvas(0,0,8000,6000);

    if (useGUI) {

         main_window=new Main(canvas, mesh);

      if ( QApplication::desktop()->width() > ((Main *)main_window)->width() + 10

	   && QApplication::desktop()->height() > ((Main *)main_window)->height() +30 ) {

	((Main *)main_window)->show();

	((Main *)main_window)->resize( ((Main *)main_window)->sizeHint());

    } else {

        ((Main *)main_window)->showMaximized();

    }

  } else {

        main_window=new MainBase(canvas, mesh);

 }

    canvas.setSceneRect(QRectF());

    if (!batch) {

      QObject::connect( qApp, SIGNAL(lastWindowClosed()), qApp, SLOT(quit()) );

    }

	  // Load plugins

	  /*QVector<SimPluginInterface *> plugins = LoadPlugins();

	  InstallPlugin(plugins[0], main_window);

	  cerr << "List of models:" << endl;

	  foreach (SimPluginInterface *p, plugins) {

		  cerr << p->ModelID().toStdString() << endl;

	  }

	  */

    // Install model or read catalogue of models

    ModelCatalogue model_catalogue(&mesh, useGUI?(Main *)main_window:0,modelfile);

    if (useGUI)

      model_catalogue.PopulateModelMenu();

    model_catalogue.InstallFirstModel();

    main_window->Init(leaffile);

    Cell::SetMagnification(1);

    Cell::setOffset(0,0);

    main_window->FitLeafToCanvas();

    main_window->Plot();

    if (batch) {

      double t=0.;

      do {

	t = main_window->TimeStep();

      } while (t < par.maxt);

    } else

      return app.exec();

  } catch (const char *message) {

    if (batch) { 

      cerr << "Exception caught:" << endl;

      cerr << message << endl;

      abort();

    } else {

      QString qmess=QString("Exception caught: %1").arg(message);

      QMessageBox::critical(0, "Critical Error", qmess);

      abort();

    }

  } catch (ios_base::failure) {

    stringstream error_message;

    error_message << "I/O failure: " << strerror(errno);

    if (batch) {

      cerr << error_message.str() <<endl;

      abort();

    } else {

      QString qmess(error_message.str().c_str());

      QMessageBox::critical(0, "I/O Error", qmess );

      abort();

    }

  }

}