# HG changeset patch
# User Michael Guravage <michael.guravage@cwi.nl>
# Date 2010-10-08 16:58:27
# Node ID bb80657b6912e4a8267fba9bf215e5ad41f71d68
# Parent  b3943d0b8126f31b6a660a4e63704ac73eb99573

Export Cell Data at interval specified in the leafML file or explicitly via the file menu.
Added new parameters to VirtualLeafpar.tmpl
Added GPL boilerplate to make_{parameter,pardialog}_source.pl files.
Added new parameters to auxin_growth.xml. See respective ChangeLog files for details.

 --
 user: Michael Guravage <michael.guravage@cwi.nl>
 branch 'default'
 changed data/leaves/ChangeLog
 changed data/leaves/auxin_growth.xml
 changed src/ChangeLog
 changed src/VirtualLeafpar.tmpl
 changed src/canvas.cpp
 changed src/canvas.h
 changed src/parameter.cpp
 changed src/parameter.h
 changed src/pardialog.cpp
 changed src/pardialog.h
 changed src/perl/ChangeLog
 changed src/perl/make_parameter_source.pl
 changed src/perl/make_pardialog_source.pl

diff --git a/data/leaves/ChangeLog b/data/leaves/ChangeLog
--- a/data/leaves/ChangeLog
+++ b/data/leaves/ChangeLog
@@ -1,3 +1,8 @@
+2010-10-08    <guravage@caterpie.sen.cwi.nl>
+
+	* auxin_growth.xml: Added export_interval and export_fn_prefix parameters.
+
+
 2010-06-23    <guravage@caterpie.sen.cwi.nl>
 
 	* Added: tutorial3_init.xml
diff --git a/data/leaves/auxin_growth.xml b/data/leaves/auxin_growth.xml
--- a/data/leaves/auxin_growth.xml
+++ b/data/leaves/auxin_growth.xml
@@ -26,6 +26,8 @@
 
 <leaf name="/Users/roel/VLeafClean_plugin/Snapshot/VLeafMacOSX/auxin_growth.xml" date="Sat Aug 29 18:28:19 2009" simtime="11660">
   <parameter>
+    <par name="export_interval" val="1000"/>
+    <par name="export_fn_prefix" val="CELLDATA_"/>
     <par name="arrowcolor" val="white"/>
     <par name="arrowsize" val="10"/>
     <par name="textcolor" val="red"/>
diff --git a/src/ChangeLog b/src/ChangeLog
--- a/src/ChangeLog
+++ b/src/ChangeLog
@@ -1,3 +1,24 @@
+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
diff --git a/src/VirtualLeafpar.tmpl b/src/VirtualLeafpar.tmpl
--- a/src/VirtualLeafpar.tmpl
+++ b/src/VirtualLeafpar.tmpl
@@ -124,3 +124,5 @@ b3 = false / bool
 b4 = false / bool
 dir1 = . / directory
 dir2 = . / directory
+export_interval = 0 / int
+export_fn_prefix = CellData_ / string
diff --git a/src/canvas.cpp b/src/canvas.cpp
--- a/src/canvas.cpp
+++ b/src/canvas.cpp
@@ -19,6 +19,7 @@
  *
  */
 
+#include <time.h>
 #include <string>
 #include <fstream>
 #include <streambuf>
@@ -1166,6 +1167,11 @@ void Main::TimeStepWrap(void)
   static int t=0;
   TimeStep();
   t++;
+
+  if ((par.export_interval > 0) && ((t % par.export_interval) == 0)){
+    this->exportCellData(QString(par.datadir) + QString('/') + QString(par.export_fn_prefix) + this->TimeStamp());
+  }
+
   // check number of timesteps
   if (t==par.nit) {
     emit SimulationDone();
@@ -1173,6 +1179,18 @@ void Main::TimeStepWrap(void)
 }
 
 
+QString Main::TimeStamp(){
+  time_t rawtime;
+  struct tm * timeinfo;
+  char buffer [15];
+
+  time ( &rawtime );
+  timeinfo = localtime ( &rawtime );
+  strftime (buffer,15,"%Y%m%d%H%M%S",timeinfo);
+  return QString(buffer);
+}
+
+
 void Main::RestartSim(void)
 {
 
@@ -1366,33 +1384,33 @@ xmlNode *Main::XMLSettingsTree(void)
   return settings;
 }
 
-void Main::exportCellData(void) {
-  Q3FileDialog *fd = new Q3FileDialog( this, "file dialog", TRUE );
-  QString fileName;
+void Main::exportCellData(QString fileName) {
   
-  stopSimulation();
-  fd->setMode( Q3FileDialog::AnyFile );
-
-  if ( fd->exec() == QDialog::Accepted ) {
-    fileName = fd->selectedFile();
-
 #ifdef QDEBUG  
-    qDebug() << "exportCellData filename: " << fileName << endl;
+  qDebug() << "exportCellData fileName: " << fileName << endl;
 #endif
 
-  // perhaps make this more general: a popup window were user selects data he wants to export
-  // can go to "settings" section of XML file as well so this can also be measured off-line
-  // mesh.CSVExport would take an QMap or so to record all options
-  // first line gives legenda
-
-    QFile file(fileName);
-    if ( file.open( IO_WriteOnly ) ) {
-      QTextStream stream( &file );
-      mesh.CSVExportCellData(stream);
-      mesh.CSVExportMeshData(stream);
-      file.close();
-    }
+  QFile file(fileName);
+  if ( file.open( IO_WriteOnly ) ) {
+    QTextStream stream( &file );
+    mesh.CSVExportCellData(stream);
+    mesh.CSVExportMeshData(stream);
+    file.close();
   }
 }
 
+
+void Main::exportCellData() {
+  QString fileName;
+  Q3FileDialog *fd = new Q3FileDialog( this, "file dialog", TRUE );
+
+  stopSimulation();
+  fd->setMode( Q3FileDialog::AnyFile );
+  if ( fd->exec() == QDialog::Accepted ) {
+    fileName = fd->selectedFile();
+  }
+  this->exportCellData(fileName);
+}
+
+
 /* finis */
diff --git a/src/canvas.h b/src/canvas.h
--- a/src/canvas.h
+++ b/src/canvas.h
@@ -207,6 +207,7 @@ class Main : public Q3MainWindow, public
   void readFirstStateXML();
   void readLastStateXML();
   void exportCellData();
+  void exportCellData(QString);
   void saveStateXML();
   void snapshot();
   void savePars();
@@ -226,6 +227,8 @@ class Main : public Q3MainWindow, public
 
   void RandomizeMesh();
 
+  QString TimeStamp();
+
  signals:
   void SimulationDone(void);
   void ParsChanged(void);
diff --git a/src/parameter.cpp b/src/parameter.cpp
--- a/src/parameter.cpp
+++ b/src/parameter.cpp
@@ -19,6 +19,9 @@
  *
  */
 
+// Do not edit. All edits will be discarded.
+
+
 #include "parameter.h"
 #include <cstdio>
 #include <cstring>
@@ -183,40 +186,44 @@ Parameter::Parameter() {
   b4 = false;
   dir1 = strdup(".");
   dir2 = strdup(".");
+  export_interval = 100;
+  export_fn_prefix = strdup("CELLDATA");
 }
 
 Parameter::~Parameter() {
     
-  // destruct parameter object
-  // free string parameter
-  CleanUp();
+// destruct parameter object
+// free string parameter
+CleanUp();
 }
 
 void Parameter::CleanUp(void) {
   if (arrowcolor) 
-    free(arrowcolor);
+     free(arrowcolor);
   if (textcolor) 
-    free(textcolor);
+     free(textcolor);
   if (cell_outline_color) 
-    free(cell_outline_color);
+     free(cell_outline_color);
   if (D) 
-    free(D);
+     free(D);
   if (initval) 
-    free(initval);
+     free(initval);
   if (datadir) 
-    free(datadir);
+     free(datadir);
   if (k) 
-    free(k);
+     free(k);
   if (s1) 
-    free(s1);
+     free(s1);
   if (s2) 
-    free(s2);
+     free(s2);
   if (s3) 
-    free(s3);
+     free(s3);
   if (dir1) 
-    free(dir1);
+     free(dir1);
   if (dir2) 
-    free(dir2);
+     free(dir2);
+  if (export_fn_prefix) 
+     free(export_fn_prefix);
 
 }
 
@@ -340,6 +347,8 @@ void Parameter::Read(const char *filenam
   dir2 = sgetpar(fp, "dir2", ".", true);
   if (strcmp(dir2, "."))
     MakeDir(dir2);
+  export_interval = igetpar(fp, "export_interval", 100, true);
+  export_fn_prefix = sgetpar(fp, "export_fn_prefix", "CELLDATA", true);
 }
 
 const char *sbool(const bool &p) {
@@ -356,18 +365,18 @@ void Parameter::Write(ostream &os) const
 
 
   if (arrowcolor) 
-    os << " arrowcolor = " << arrowcolor << endl;
+  os << " arrowcolor = " << arrowcolor << endl;
   os << " arrowsize = " << arrowsize << endl;
 
   if (textcolor) 
-    os << " textcolor = " << textcolor << endl;
+  os << " textcolor = " << textcolor << endl;
   os << " cellnumsize = " << cellnumsize << endl;
   os << " nodenumsize = " << nodenumsize << endl;
   os << " node_mag = " << node_mag << endl;
   os << " outlinewidth = " << outlinewidth << endl;
 
   if (cell_outline_color) 
-    os << " cell_outline_color = " << cell_outline_color << endl;
+  os << " cell_outline_color = " << cell_outline_color << endl;
   os << " resize_stride = " << resize_stride << endl;
   os << " T = " << T << endl;
   os << " lambda_length = " << lambda_length << endl;
@@ -423,12 +432,12 @@ void Parameter::Write(ostream &os) const
   os << " dt = " << dt << endl;
   os << " rd_dt = " << rd_dt << endl;
   if (datadir) {
-    QDir dataDir = QDir::home().relativeFilePath(datadir);
-    os << " datadir = " << dataDir.dirName().toStdString() << endl;
-  }
-  else {
-    os << "datadir = ." << endl;
-  }
+                                     QDir dataDir = QDir::home().relativeFilePath(datadir);
+                                     os << " datadir = " << dataDir.dirName().toStdString() << endl;
+                                 }
+                                 else {
+                                     os << "datadir = ." << endl;
+                                 }
   os << " movie = " << sbool(movie) << endl;
   os << " nit = " << nit << endl;
   os << " maxt = " << maxt << endl;
@@ -457,1403 +466,1431 @@ void Parameter::Write(ostream &os) const
   os << " i5 = " << i5 << endl;
 
   if (s1) 
-    os << " s1 = " << s1 << endl;
+  os << " s1 = " << s1 << endl;
 
   if (s2) 
-    os << " s2 = " << s2 << endl;
+  os << " s2 = " << s2 << endl;
 
   if (s3) 
-    os << " s3 = " << s3 << endl;
+  os << " s3 = " << s3 << endl;
   os << " b1 = " << sbool(b1) << endl;
   os << " b2 = " << sbool(b2) << endl;
   os << " b3 = " << sbool(b3) << endl;
   os << " b4 = " << sbool(b4) << endl;
 
   if (dir1) 
-    os << " dir1 = " << dir1 << endl;
+  os << " dir1 = " << dir1 << endl;
 
   if (dir2) 
-    os << " dir2 = " << dir2 << endl;
+  os << " dir2 = " << dir2 << endl;
+  os << " export_interval = " << export_interval << endl;
+
+  if (export_fn_prefix) 
+  os << " export_fn_prefix = " << export_fn_prefix << endl;
 }
 
 void Parameter::XMLAdd(xmlNode *root) const {
-  xmlNode *xmlparameter = xmlNewChild(root, NULL, BAD_CAST "parameter", NULL);
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "arrowcolor" );
-    ostringstream text;
+    xmlNode *xmlparameter = xmlNewChild(root, NULL, BAD_CAST "parameter", NULL);
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "arrowcolor" );
+  ostringstream text;
 
-    if (arrowcolor) 
-      text << arrowcolor;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "arrowsize" );
-    ostringstream text;
+  if (arrowcolor) 
+    text << arrowcolor;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "arrowsize" );
+  ostringstream text;
     text << arrowsize;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "textcolor" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "textcolor" );
+  ostringstream text;
 
-    if (textcolor) 
-      text << textcolor;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "cellnumsize" );
-    ostringstream text;
+  if (textcolor) 
+    text << textcolor;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "cellnumsize" );
+  ostringstream text;
     text << cellnumsize;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "nodenumsize" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "nodenumsize" );
+  ostringstream text;
     text << nodenumsize;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "node_mag" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "node_mag" );
+  ostringstream text;
     text << node_mag;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "outlinewidth" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "outlinewidth" );
+  ostringstream text;
     text << outlinewidth;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "cell_outline_color" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "cell_outline_color" );
+  ostringstream text;
 
-    if (cell_outline_color) 
-      text << cell_outline_color;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "resize_stride" );
-    ostringstream text;
+  if (cell_outline_color) 
+    text << cell_outline_color;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "resize_stride" );
+  ostringstream text;
     text << resize_stride;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "T" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "T" );
+  ostringstream text;
     text << T;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "lambda_length" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "lambda_length" );
+  ostringstream text;
     text << lambda_length;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "lambda_celllength" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "lambda_celllength" );
+  ostringstream text;
     text << lambda_celllength;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "target_length" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "target_length" );
+  ostringstream text;
     text << target_length;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "cell_expansion_rate" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "cell_expansion_rate" );
+  ostringstream text;
     text << cell_expansion_rate;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "cell_div_expansion_rate" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "cell_div_expansion_rate" );
+  ostringstream text;
     text << cell_div_expansion_rate;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "auxin_dependent_growth" );
-    ostringstream text;
-    text << sbool(auxin_dependent_growth);
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "ode_accuracy" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "auxin_dependent_growth" );
+  ostringstream text;
+text << sbool(auxin_dependent_growth);
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "ode_accuracy" );
+  ostringstream text;
     text << ode_accuracy;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "mc_stepsize" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "mc_stepsize" );
+  ostringstream text;
     text << mc_stepsize;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "mc_cell_stepsize" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "mc_cell_stepsize" );
+  ostringstream text;
     text << mc_cell_stepsize;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "energy_threshold" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "energy_threshold" );
+  ostringstream text;
     text << energy_threshold;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "bend_lambda" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "bend_lambda" );
+  ostringstream text;
     text << bend_lambda;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "alignment_lambda" );
+  ostringstream text;
+    text << alignment_lambda;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "rel_cell_div_threshold" );
+  ostringstream text;
+    text << rel_cell_div_threshold;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "rel_perimeter_stiffness" );
+  ostringstream text;
+    text << rel_perimeter_stiffness;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "collapse_node_threshold" );
+  ostringstream text;
+    text << collapse_node_threshold;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "morphogen_div_threshold" );
+  ostringstream text;
+    text << morphogen_div_threshold;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "morphogen_expansion_threshold" );
+  ostringstream text;
+    text << morphogen_expansion_threshold;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "copy_wall" );
+  ostringstream text;
+text << sbool(copy_wall);
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "source" );
+  ostringstream text;
+    text << source;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "D" );
+  xmlNode *xmlvalarray = xmlNewChild(xmlpar, NULL, BAD_CAST "valarray", NULL);
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "alignment_lambda" );
     ostringstream text;
-    text << alignment_lambda;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "rel_cell_div_threshold" );
-    ostringstream text;
-    text << rel_cell_div_threshold;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "rel_perimeter_stiffness" );
-    ostringstream text;
-    text << rel_perimeter_stiffness;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "collapse_node_threshold" );
-    ostringstream text;
-    text << collapse_node_threshold;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "morphogen_div_threshold" );
-    ostringstream text;
-    text << morphogen_div_threshold;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "morphogen_expansion_threshold" );
-    ostringstream text;
-    text << morphogen_expansion_threshold;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "copy_wall" );
-    ostringstream text;
-    text << sbool(copy_wall);
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "source" );
-    ostringstream text;
-    text << source;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << D[0];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "D" );
-    xmlNode *xmlvalarray = xmlNewChild(xmlpar, NULL, BAD_CAST "valarray", NULL);
-    {
-      ostringstream text;
-      text << D[0];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << D[1];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << D[2];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << D[3];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << D[4];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << D[5];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << D[6];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << D[7];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << D[8];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << D[9];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << D[10];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << D[11];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << D[12];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << D[13];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << D[14];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
+    ostringstream text;
+    text << D[1];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << D[2];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << D[3];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << D[4];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << D[5];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << D[6];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << D[7];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << D[8];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << D[9];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << D[10];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << D[11];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << D[12];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << D[13];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "initval" );
-    xmlNode *xmlvalarray = xmlNewChild(xmlpar, NULL, BAD_CAST "valarray", NULL);
-    {
-      ostringstream text;
-      text << initval[0];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << initval[1];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << initval[2];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << initval[3];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << initval[4];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << initval[5];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << initval[6];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << initval[7];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << initval[8];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << initval[9];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << initval[10];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << initval[11];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << initval[12];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << initval[13];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << initval[14];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
+    ostringstream text;
+    text << D[14];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "initval" );
+  xmlNode *xmlvalarray = xmlNewChild(xmlpar, NULL, BAD_CAST "valarray", NULL);
+  {
+    ostringstream text;
+    text << initval[0];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << initval[1];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << initval[2];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << initval[3];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << initval[4];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << initval[5];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "k1" );
     ostringstream text;
-    text << k1;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "k2" );
-    ostringstream text;
-    text << k2;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << initval[6];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "r" );
     ostringstream text;
-    text << r;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << initval[7];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "kr" );
     ostringstream text;
-    text << kr;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "km" );
-    ostringstream text;
-    text << km;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << initval[8];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "Pi_tot" );
     ostringstream text;
-    text << Pi_tot;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << initval[9];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "transport" );
     ostringstream text;
-    text << transport;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "ka" );
-    ostringstream text;
-    text << ka;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << initval[10];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "pin_prod" );
     ostringstream text;
-    text << pin_prod;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << initval[11];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "pin_prod_in_epidermis" );
     ostringstream text;
-    text << pin_prod_in_epidermis;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << initval[12];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "pin_breakdown" );
     ostringstream text;
-    text << pin_breakdown;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "pin_breakdown_internal" );
-    ostringstream text;
-    text << pin_breakdown_internal;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << initval[13];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "aux1prod" );
     ostringstream text;
-    text << aux1prod;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "aux1prodmeso" );
-    ostringstream text;
-    text << aux1prodmeso;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "aux1decay" );
-    ostringstream text;
-    text << aux1decay;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "aux1decaymeso" );
-    ostringstream text;
-    text << aux1decaymeso;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "aux1transport" );
-    ostringstream text;
-    text << aux1transport;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "aux_cons" );
-    ostringstream text;
-    text << aux_cons;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "aux_breakdown" );
-    ostringstream text;
-    text << aux_breakdown;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "kaux1" );
-    ostringstream text;
-    text << kaux1;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << initval[14];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "kap" );
-    ostringstream text;
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "k1" );
+  ostringstream text;
+    text << k1;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "k2" );
+  ostringstream text;
+    text << k2;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "r" );
+  ostringstream text;
+    text << r;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "kr" );
+  ostringstream text;
+    text << kr;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "km" );
+  ostringstream text;
+    text << km;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "Pi_tot" );
+  ostringstream text;
+    text << Pi_tot;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "transport" );
+  ostringstream text;
+    text << transport;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "ka" );
+  ostringstream text;
+    text << ka;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "pin_prod" );
+  ostringstream text;
+    text << pin_prod;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "pin_prod_in_epidermis" );
+  ostringstream text;
+    text << pin_prod_in_epidermis;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "pin_breakdown" );
+  ostringstream text;
+    text << pin_breakdown;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "pin_breakdown_internal" );
+  ostringstream text;
+    text << pin_breakdown_internal;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "aux1prod" );
+  ostringstream text;
+    text << aux1prod;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "aux1prodmeso" );
+  ostringstream text;
+    text << aux1prodmeso;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "aux1decay" );
+  ostringstream text;
+    text << aux1decay;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "aux1decaymeso" );
+  ostringstream text;
+    text << aux1decaymeso;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "aux1transport" );
+  ostringstream text;
+    text << aux1transport;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "aux_cons" );
+  ostringstream text;
+    text << aux_cons;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "aux_breakdown" );
+  ostringstream text;
+    text << aux_breakdown;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "kaux1" );
+  ostringstream text;
+    text << kaux1;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "kap" );
+  ostringstream text;
     text << kap;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "leaf_tip_source" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "leaf_tip_source" );
+  ostringstream text;
     text << leaf_tip_source;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "sam_efflux" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "sam_efflux" );
+  ostringstream text;
     text << sam_efflux;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "sam_auxin" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "sam_auxin" );
+  ostringstream text;
     text << sam_auxin;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "sam_auxin_breakdown" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "sam_auxin_breakdown" );
+  ostringstream text;
     text << sam_auxin_breakdown;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "van3prod" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "van3prod" );
+  ostringstream text;
     text << van3prod;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "van3autokat" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "van3autokat" );
+  ostringstream text;
     text << van3autokat;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "van3sat" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "van3sat" );
+  ostringstream text;
     text << van3sat;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "k2van3" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "k2van3" );
+  ostringstream text;
     text << k2van3;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "dt" );
+  ostringstream text;
+    text << dt;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "rd_dt" );
+  ostringstream text;
+    text << rd_dt;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "datadir" );
+  ostringstream text;
+
+  if (datadir) {
+                                     QDir dataDir = QDir::home().relativeFilePath(datadir);
+                                     text << dataDir.dirName().toStdString();
+                                   }
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "movie" );
+  ostringstream text;
+text << sbool(movie);
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "nit" );
+  ostringstream text;
+    text << nit;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "maxt" );
+  ostringstream text;
+    text << maxt;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "storage_stride" );
+  ostringstream text;
+    text << storage_stride;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "xml_storage_stride" );
+  ostringstream text;
+    text << xml_storage_stride;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "rseed" );
+  ostringstream text;
+    text << rseed;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "constituous_expansion_limit" );
+  ostringstream text;
+    text << constituous_expansion_limit;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "vessel_inh_level" );
+  ostringstream text;
+    text << vessel_inh_level;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "vessel_expansion_rate" );
+  ostringstream text;
+    text << vessel_expansion_rate;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "d" );
+  ostringstream text;
+    text << d;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "e" );
+  ostringstream text;
+    text << e;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "f" );
+  ostringstream text;
+    text << f;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "c" );
+  ostringstream text;
+    text << c;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "mu" );
+  ostringstream text;
+    text << mu;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "nu" );
+  ostringstream text;
+    text << nu;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "rho0" );
+  ostringstream text;
+    text << rho0;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "rho1" );
+  ostringstream text;
+    text << rho1;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "c0" );
+  ostringstream text;
+    text << c0;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "gamma" );
+  ostringstream text;
+    text << gamma;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "eps" );
+  ostringstream text;
+    text << eps;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "k" );
+  xmlNode *xmlvalarray = xmlNewChild(xmlpar, NULL, BAD_CAST "valarray", NULL);
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "dt" );
     ostringstream text;
-    text << dt;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "rd_dt" );
-    ostringstream text;
-    text << rd_dt;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "datadir" );
-    ostringstream text;
-
-    if (datadir) {
-      QDir dataDir = QDir::home().relativeFilePath(datadir);
-      text << dataDir.dirName().toStdString();
-    }
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "movie" );
-    ostringstream text;
-    text << sbool(movie);
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "nit" );
-    ostringstream text;
-    text << nit;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "maxt" );
-    ostringstream text;
-    text << maxt;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "storage_stride" );
-    ostringstream text;
-    text << storage_stride;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "xml_storage_stride" );
-    ostringstream text;
-    text << xml_storage_stride;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "rseed" );
-    ostringstream text;
-    text << rseed;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "constituous_expansion_limit" );
-    ostringstream text;
-    text << constituous_expansion_limit;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "vessel_inh_level" );
-    ostringstream text;
-    text << vessel_inh_level;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "vessel_expansion_rate" );
-    ostringstream text;
-    text << vessel_expansion_rate;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "d" );
-    ostringstream text;
-    text << d;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "e" );
-    ostringstream text;
-    text << e;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "f" );
-    ostringstream text;
-    text << f;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "c" );
-    ostringstream text;
-    text << c;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "mu" );
-    ostringstream text;
-    text << mu;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << k[0];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "nu" );
     ostringstream text;
-    text << nu;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << k[1];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << k[2];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "rho0" );
+    ostringstream text;
+    text << k[3];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
     ostringstream text;
-    text << rho0;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << k[4];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << k[5];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << k[6];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "rho1" );
+    ostringstream text;
+    text << k[7];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
     ostringstream text;
-    text << rho1;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << k[8];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << k[9];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "c0" );
     ostringstream text;
-    text << c0;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << k[10];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "gamma" );
     ostringstream text;
-    text << gamma;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << k[11];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "eps" );
     ostringstream text;
-    text << eps;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << k[12];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
+  }
+  {
+    ostringstream text;
+    text << k[13];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
   {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "k" );
-    xmlNode *xmlvalarray = xmlNewChild(xmlpar, NULL, BAD_CAST "valarray", NULL);
-    {
-      ostringstream text;
-      text << k[0];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << k[1];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << k[2];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << k[3];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << k[4];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << k[5];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << k[6];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << k[7];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << k[8];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << k[9];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << k[10];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << k[11];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << k[12];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << k[13];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-    {
-      ostringstream text;
-      text << k[14];
-      xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
-      xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
-    }
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "i1" );
     ostringstream text;
-    text << i1;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+    text << k[14];
+    xmlNode *xmlval = xmlNewChild(xmlvalarray, NULL, BAD_CAST "val", NULL);
+    xmlNewProp(xmlval, BAD_CAST "v", BAD_CAST text.str().c_str());
   }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "i2" );
-    ostringstream text;
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "i1" );
+  ostringstream text;
+    text << i1;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "i2" );
+  ostringstream text;
     text << i2;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "i3" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "i3" );
+  ostringstream text;
     text << i3;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "i4" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "i4" );
+  ostringstream text;
     text << i4;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "i5" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "i5" );
+  ostringstream text;
     text << i5;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "s1" );
-    ostringstream text;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "s1" );
+  ostringstream text;
 
-    if (s1) 
-      text << s1;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "s2" );
-    ostringstream text;
+  if (s1) 
+    text << s1;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "s2" );
+  ostringstream text;
 
-    if (s2) 
-      text << s2;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "s3" );
-    ostringstream text;
+  if (s2) 
+    text << s2;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "s3" );
+  ostringstream text;
 
-    if (s3) 
-      text << s3;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "b1" );
-    ostringstream text;
-    text << sbool(b1);
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "b2" );
-    ostringstream text;
-    text << sbool(b2);
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "b3" );
-    ostringstream text;
-    text << sbool(b3);
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "b4" );
-    ostringstream text;
-    text << sbool(b4);
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "dir1" );
-    ostringstream text;
+  if (s3) 
+    text << s3;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "b1" );
+  ostringstream text;
+text << sbool(b1);
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "b2" );
+  ostringstream text;
+text << sbool(b2);
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "b3" );
+  ostringstream text;
+text << sbool(b3);
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "b4" );
+  ostringstream text;
+text << sbool(b4);
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "dir1" );
+  ostringstream text;
 
-    if (dir1) 
-      text << dir1;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
-  {
-    xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
-    xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "dir2" );
-    ostringstream text;
+  if (dir1) 
+    text << dir1;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "dir2" );
+  ostringstream text;
 
-    if (dir2) 
-      text << dir2;
-    xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
-  }
+  if (dir2) 
+    text << dir2;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "export_interval" );
+  ostringstream text;
+    text << export_interval;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
+{
+  xmlNode *xmlpar = xmlNewChild(xmlparameter, NULL, BAD_CAST "par", NULL);
+  xmlNewProp(xmlpar, BAD_CAST "name", BAD_CAST "export_fn_prefix" );
+  ostringstream text;
+
+  if (export_fn_prefix) 
+    text << export_fn_prefix;
+xmlNewProp(xmlpar, BAD_CAST "val", BAD_CAST text.str().c_str());
+}
 }
 void Parameter::AssignValToPar(const char *namec, const char *valc) {
   QLocale standardlocale(QLocale::C);
   bool ok;
-  if (!strcmp(namec, "arrowcolor")) {
-    if (arrowcolor) { free(arrowcolor); }
-    arrowcolor=strdup(valc);
-  }
-  if (!strcmp(namec, "arrowsize")) {
-    arrowsize = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'arrowsize' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "textcolor")) {
-    if (textcolor) { free(textcolor); }
-    textcolor=strdup(valc);
-  }
-  if (!strcmp(namec, "cellnumsize")) {
-    cellnumsize = standardlocale.toInt(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'cellnumsize' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "nodenumsize")) {
-    nodenumsize = standardlocale.toInt(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'nodenumsize' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "node_mag")) {
-    node_mag = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'node_mag' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "outlinewidth")) {
-    outlinewidth = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'outlinewidth' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "cell_outline_color")) {
-    if (cell_outline_color) { free(cell_outline_color); }
-    cell_outline_color=strdup(valc);
-  }
-  if (!strcmp(namec, "resize_stride")) {
-    resize_stride = standardlocale.toInt(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'resize_stride' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "T")) {
-    T = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'T' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "lambda_length")) {
-    lambda_length = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'lambda_length' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "lambda_celllength")) {
-    lambda_celllength = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'lambda_celllength' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "target_length")) {
-    target_length = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'target_length' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "cell_expansion_rate")) {
-    cell_expansion_rate = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'cell_expansion_rate' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "cell_div_expansion_rate")) {
-    cell_div_expansion_rate = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'cell_div_expansion_rate' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "auxin_dependent_growth")) {
-    auxin_dependent_growth = strtobool(valc);
-  }
-  if (!strcmp(namec, "ode_accuracy")) {
-    ode_accuracy = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'ode_accuracy' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "mc_stepsize")) {
-    mc_stepsize = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'mc_stepsize' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "mc_cell_stepsize")) {
-    mc_cell_stepsize = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'mc_cell_stepsize' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "energy_threshold")) {
-    energy_threshold = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'energy_threshold' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "bend_lambda")) {
-    bend_lambda = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'bend_lambda' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "alignment_lambda")) {
-    alignment_lambda = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'alignment_lambda' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "rel_cell_div_threshold")) {
-    rel_cell_div_threshold = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'rel_cell_div_threshold' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "rel_perimeter_stiffness")) {
-    rel_perimeter_stiffness = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'rel_perimeter_stiffness' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "collapse_node_threshold")) {
-    collapse_node_threshold = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'collapse_node_threshold' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "morphogen_div_threshold")) {
-    morphogen_div_threshold = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'morphogen_div_threshold' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "morphogen_expansion_threshold")) {
-    morphogen_expansion_threshold = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'morphogen_expansion_threshold' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "copy_wall")) {
-    copy_wall = strtobool(valc);
-  }
-  if (!strcmp(namec, "source")) {
-    source = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'source' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "k1")) {
-    k1 = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'k1' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "k2")) {
-    k2 = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'k2' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "r")) {
-    r = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'r' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "kr")) {
-    kr = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'kr' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "km")) {
-    km = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'km' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "Pi_tot")) {
-    Pi_tot = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'Pi_tot' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "transport")) {
-    transport = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'transport' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "ka")) {
-    ka = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'ka' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "pin_prod")) {
-    pin_prod = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'pin_prod' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "pin_prod_in_epidermis")) {
-    pin_prod_in_epidermis = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'pin_prod_in_epidermis' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "pin_breakdown")) {
-    pin_breakdown = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'pin_breakdown' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "pin_breakdown_internal")) {
-    pin_breakdown_internal = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'pin_breakdown_internal' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "aux1prod")) {
-    aux1prod = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'aux1prod' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "aux1prodmeso")) {
-    aux1prodmeso = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'aux1prodmeso' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "aux1decay")) {
-    aux1decay = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'aux1decay' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "aux1decaymeso")) {
-    aux1decaymeso = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'aux1decaymeso' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "aux1transport")) {
-    aux1transport = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'aux1transport' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "aux_cons")) {
-    aux_cons = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'aux_cons' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "aux_breakdown")) {
-    aux_breakdown = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'aux_breakdown' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "kaux1")) {
-    kaux1 = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'kaux1' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "kap")) {
-    kap = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'kap' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "leaf_tip_source")) {
-    leaf_tip_source = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'leaf_tip_source' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "sam_efflux")) {
-    sam_efflux = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'sam_efflux' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "sam_auxin")) {
-    sam_auxin = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'sam_auxin' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "sam_auxin_breakdown")) {
-    sam_auxin_breakdown = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'sam_auxin_breakdown' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "van3prod")) {
-    van3prod = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'van3prod' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "van3autokat")) {
-    van3autokat = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'van3autokat' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "van3sat")) {
-    van3sat = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'van3sat' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "k2van3")) {
-    k2van3 = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'k2van3' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "dt")) {
-    dt = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'dt' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "rd_dt")) {
-    rd_dt = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'rd_dt' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "datadir")) {
-    if (datadir) { free(datadir); }
-    datadir=strdup(valc);
-    datadir = AppendHomeDirIfPathRelative(datadir);
-  }
-  if (!strcmp(namec, "movie")) {
-    movie = strtobool(valc);
-  }
-  if (!strcmp(namec, "nit")) {
-    nit = standardlocale.toInt(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'nit' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "maxt")) {
-    maxt = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'maxt' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "storage_stride")) {
-    storage_stride = standardlocale.toInt(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'storage_stride' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "xml_storage_stride")) {
-    xml_storage_stride = standardlocale.toInt(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'xml_storage_stride' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "rseed")) {
-    rseed = standardlocale.toInt(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'rseed' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "constituous_expansion_limit")) {
-    constituous_expansion_limit = standardlocale.toInt(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'constituous_expansion_limit' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "vessel_inh_level")) {
-    vessel_inh_level = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'vessel_inh_level' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "vessel_expansion_rate")) {
-    vessel_expansion_rate = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'vessel_expansion_rate' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "d")) {
-    d = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'd' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "e")) {
-    e = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'e' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "f")) {
-    f = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'f' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "c")) {
-    c = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'c' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "mu")) {
-    mu = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'mu' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "nu")) {
-    nu = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'nu' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "rho0")) {
-    rho0 = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'rho0' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "rho1")) {
-    rho1 = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'rho1' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "c0")) {
-    c0 = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'c0' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "gamma")) {
-    gamma = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'gamma' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "eps")) {
-    eps = standardlocale.toDouble(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'eps' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "i1")) {
-    i1 = standardlocale.toInt(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'i1' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "i2")) {
-    i2 = standardlocale.toInt(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'i2' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "i3")) {
-    i3 = standardlocale.toInt(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'i3' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "i4")) {
-    i4 = standardlocale.toInt(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'i4' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "i5")) {
-    i5 = standardlocale.toInt(valc, &ok);
-    if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'i5' from XML file.",valc); }
-  }
-  if (!strcmp(namec, "s1")) {
-    if (s1) { free(s1); }
-    s1=strdup(valc);
-  }
-  if (!strcmp(namec, "s2")) {
-    if (s2) { free(s2); }
-    s2=strdup(valc);
-  }
-  if (!strcmp(namec, "s3")) {
-    if (s3) { free(s3); }
-    s3=strdup(valc);
-  }
-  if (!strcmp(namec, "b1")) {
-    b1 = strtobool(valc);
-  }
-  if (!strcmp(namec, "b2")) {
-    b2 = strtobool(valc);
-  }
-  if (!strcmp(namec, "b3")) {
-    b3 = strtobool(valc);
-  }
-  if (!strcmp(namec, "b4")) {
-    b4 = strtobool(valc);
-  }
-  if (!strcmp(namec, "dir1")) {
-    if (dir1) { free(dir1); }
-    dir1=strdup(valc);
-  }
-  if (!strcmp(namec, "dir2")) {
-    if (dir2) { free(dir2); }
-    dir2=strdup(valc);
+if (!strcmp(namec, "arrowcolor")) {
+  if (arrowcolor) { free(arrowcolor); }
+  arrowcolor=strdup(valc);
+}
+if (!strcmp(namec, "arrowsize")) {
+  arrowsize = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'arrowsize' from XML file.",valc); }
+}
+if (!strcmp(namec, "textcolor")) {
+  if (textcolor) { free(textcolor); }
+  textcolor=strdup(valc);
+}
+if (!strcmp(namec, "cellnumsize")) {
+  cellnumsize = standardlocale.toInt(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'cellnumsize' from XML file.",valc); }
+}
+if (!strcmp(namec, "nodenumsize")) {
+  nodenumsize = standardlocale.toInt(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'nodenumsize' from XML file.",valc); }
+}
+if (!strcmp(namec, "node_mag")) {
+  node_mag = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'node_mag' from XML file.",valc); }
+}
+if (!strcmp(namec, "outlinewidth")) {
+  outlinewidth = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'outlinewidth' from XML file.",valc); }
+}
+if (!strcmp(namec, "cell_outline_color")) {
+  if (cell_outline_color) { free(cell_outline_color); }
+  cell_outline_color=strdup(valc);
+}
+if (!strcmp(namec, "resize_stride")) {
+  resize_stride = standardlocale.toInt(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'resize_stride' from XML file.",valc); }
+}
+if (!strcmp(namec, "T")) {
+  T = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'T' from XML file.",valc); }
+}
+if (!strcmp(namec, "lambda_length")) {
+  lambda_length = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'lambda_length' from XML file.",valc); }
+}
+if (!strcmp(namec, "lambda_celllength")) {
+  lambda_celllength = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'lambda_celllength' from XML file.",valc); }
+}
+if (!strcmp(namec, "target_length")) {
+  target_length = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'target_length' from XML file.",valc); }
+}
+if (!strcmp(namec, "cell_expansion_rate")) {
+  cell_expansion_rate = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'cell_expansion_rate' from XML file.",valc); }
+}
+if (!strcmp(namec, "cell_div_expansion_rate")) {
+  cell_div_expansion_rate = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'cell_div_expansion_rate' from XML file.",valc); }
+}
+if (!strcmp(namec, "auxin_dependent_growth")) {
+auxin_dependent_growth = strtobool(valc);
+}
+if (!strcmp(namec, "ode_accuracy")) {
+  ode_accuracy = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'ode_accuracy' from XML file.",valc); }
+}
+if (!strcmp(namec, "mc_stepsize")) {
+  mc_stepsize = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'mc_stepsize' from XML file.",valc); }
+}
+if (!strcmp(namec, "mc_cell_stepsize")) {
+  mc_cell_stepsize = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'mc_cell_stepsize' from XML file.",valc); }
+}
+if (!strcmp(namec, "energy_threshold")) {
+  energy_threshold = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'energy_threshold' from XML file.",valc); }
+}
+if (!strcmp(namec, "bend_lambda")) {
+  bend_lambda = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'bend_lambda' from XML file.",valc); }
+}
+if (!strcmp(namec, "alignment_lambda")) {
+  alignment_lambda = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'alignment_lambda' from XML file.",valc); }
+}
+if (!strcmp(namec, "rel_cell_div_threshold")) {
+  rel_cell_div_threshold = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'rel_cell_div_threshold' from XML file.",valc); }
+}
+if (!strcmp(namec, "rel_perimeter_stiffness")) {
+  rel_perimeter_stiffness = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'rel_perimeter_stiffness' from XML file.",valc); }
+}
+if (!strcmp(namec, "collapse_node_threshold")) {
+  collapse_node_threshold = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'collapse_node_threshold' from XML file.",valc); }
+}
+if (!strcmp(namec, "morphogen_div_threshold")) {
+  morphogen_div_threshold = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'morphogen_div_threshold' from XML file.",valc); }
+}
+if (!strcmp(namec, "morphogen_expansion_threshold")) {
+  morphogen_expansion_threshold = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'morphogen_expansion_threshold' from XML file.",valc); }
+}
+if (!strcmp(namec, "copy_wall")) {
+copy_wall = strtobool(valc);
+}
+if (!strcmp(namec, "source")) {
+  source = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'source' from XML file.",valc); }
+}
+if (!strcmp(namec, "k1")) {
+  k1 = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'k1' from XML file.",valc); }
+}
+if (!strcmp(namec, "k2")) {
+  k2 = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'k2' from XML file.",valc); }
+}
+if (!strcmp(namec, "r")) {
+  r = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'r' from XML file.",valc); }
+}
+if (!strcmp(namec, "kr")) {
+  kr = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'kr' from XML file.",valc); }
+}
+if (!strcmp(namec, "km")) {
+  km = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'km' from XML file.",valc); }
+}
+if (!strcmp(namec, "Pi_tot")) {
+  Pi_tot = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'Pi_tot' from XML file.",valc); }
+}
+if (!strcmp(namec, "transport")) {
+  transport = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'transport' from XML file.",valc); }
+}
+if (!strcmp(namec, "ka")) {
+  ka = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'ka' from XML file.",valc); }
+}
+if (!strcmp(namec, "pin_prod")) {
+  pin_prod = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'pin_prod' from XML file.",valc); }
+}
+if (!strcmp(namec, "pin_prod_in_epidermis")) {
+  pin_prod_in_epidermis = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'pin_prod_in_epidermis' from XML file.",valc); }
+}
+if (!strcmp(namec, "pin_breakdown")) {
+  pin_breakdown = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'pin_breakdown' from XML file.",valc); }
+}
+if (!strcmp(namec, "pin_breakdown_internal")) {
+  pin_breakdown_internal = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'pin_breakdown_internal' from XML file.",valc); }
+}
+if (!strcmp(namec, "aux1prod")) {
+  aux1prod = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'aux1prod' from XML file.",valc); }
+}
+if (!strcmp(namec, "aux1prodmeso")) {
+  aux1prodmeso = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'aux1prodmeso' from XML file.",valc); }
+}
+if (!strcmp(namec, "aux1decay")) {
+  aux1decay = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'aux1decay' from XML file.",valc); }
+}
+if (!strcmp(namec, "aux1decaymeso")) {
+  aux1decaymeso = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'aux1decaymeso' from XML file.",valc); }
+}
+if (!strcmp(namec, "aux1transport")) {
+  aux1transport = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'aux1transport' from XML file.",valc); }
+}
+if (!strcmp(namec, "aux_cons")) {
+  aux_cons = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'aux_cons' from XML file.",valc); }
+}
+if (!strcmp(namec, "aux_breakdown")) {
+  aux_breakdown = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'aux_breakdown' from XML file.",valc); }
+}
+if (!strcmp(namec, "kaux1")) {
+  kaux1 = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'kaux1' from XML file.",valc); }
+}
+if (!strcmp(namec, "kap")) {
+  kap = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'kap' from XML file.",valc); }
+}
+if (!strcmp(namec, "leaf_tip_source")) {
+  leaf_tip_source = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'leaf_tip_source' from XML file.",valc); }
+}
+if (!strcmp(namec, "sam_efflux")) {
+  sam_efflux = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'sam_efflux' from XML file.",valc); }
+}
+if (!strcmp(namec, "sam_auxin")) {
+  sam_auxin = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'sam_auxin' from XML file.",valc); }
+}
+if (!strcmp(namec, "sam_auxin_breakdown")) {
+  sam_auxin_breakdown = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'sam_auxin_breakdown' from XML file.",valc); }
+}
+if (!strcmp(namec, "van3prod")) {
+  van3prod = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'van3prod' from XML file.",valc); }
+}
+if (!strcmp(namec, "van3autokat")) {
+  van3autokat = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'van3autokat' from XML file.",valc); }
+}
+if (!strcmp(namec, "van3sat")) {
+  van3sat = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'van3sat' from XML file.",valc); }
+}
+if (!strcmp(namec, "k2van3")) {
+  k2van3 = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'k2van3' from XML file.",valc); }
+}
+if (!strcmp(namec, "dt")) {
+  dt = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'dt' from XML file.",valc); }
+}
+if (!strcmp(namec, "rd_dt")) {
+  rd_dt = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'rd_dt' from XML file.",valc); }
+}
+if (!strcmp(namec, "datadir")) {
+  if (datadir) { free(datadir); }
+  datadir=strdup(valc);
+  datadir = AppendHomeDirIfPathRelative(datadir);
+}
+if (!strcmp(namec, "movie")) {
+movie = strtobool(valc);
+}
+if (!strcmp(namec, "nit")) {
+  nit = standardlocale.toInt(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'nit' from XML file.",valc); }
+}
+if (!strcmp(namec, "maxt")) {
+  maxt = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'maxt' from XML file.",valc); }
+}
+if (!strcmp(namec, "storage_stride")) {
+  storage_stride = standardlocale.toInt(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'storage_stride' from XML file.",valc); }
+}
+if (!strcmp(namec, "xml_storage_stride")) {
+  xml_storage_stride = standardlocale.toInt(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'xml_storage_stride' from XML file.",valc); }
+}
+if (!strcmp(namec, "rseed")) {
+  rseed = standardlocale.toInt(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'rseed' from XML file.",valc); }
+}
+if (!strcmp(namec, "constituous_expansion_limit")) {
+  constituous_expansion_limit = standardlocale.toInt(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'constituous_expansion_limit' from XML file.",valc); }
+}
+if (!strcmp(namec, "vessel_inh_level")) {
+  vessel_inh_level = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'vessel_inh_level' from XML file.",valc); }
+}
+if (!strcmp(namec, "vessel_expansion_rate")) {
+  vessel_expansion_rate = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'vessel_expansion_rate' from XML file.",valc); }
+}
+if (!strcmp(namec, "d")) {
+  d = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'd' from XML file.",valc); }
+}
+if (!strcmp(namec, "e")) {
+  e = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'e' from XML file.",valc); }
+}
+if (!strcmp(namec, "f")) {
+  f = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'f' from XML file.",valc); }
+}
+if (!strcmp(namec, "c")) {
+  c = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'c' from XML file.",valc); }
+}
+if (!strcmp(namec, "mu")) {
+  mu = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'mu' from XML file.",valc); }
+}
+if (!strcmp(namec, "nu")) {
+  nu = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'nu' from XML file.",valc); }
+}
+if (!strcmp(namec, "rho0")) {
+  rho0 = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'rho0' from XML file.",valc); }
+}
+if (!strcmp(namec, "rho1")) {
+  rho1 = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'rho1' from XML file.",valc); }
+}
+if (!strcmp(namec, "c0")) {
+  c0 = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'c0' from XML file.",valc); }
+}
+if (!strcmp(namec, "gamma")) {
+  gamma = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'gamma' from XML file.",valc); }
+}
+if (!strcmp(namec, "eps")) {
+  eps = standardlocale.toDouble(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to double while reading parameter 'eps' from XML file.",valc); }
+}
+if (!strcmp(namec, "i1")) {
+  i1 = standardlocale.toInt(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'i1' from XML file.",valc); }
+}
+if (!strcmp(namec, "i2")) {
+  i2 = standardlocale.toInt(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'i2' from XML file.",valc); }
+}
+if (!strcmp(namec, "i3")) {
+  i3 = standardlocale.toInt(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'i3' from XML file.",valc); }
+}
+if (!strcmp(namec, "i4")) {
+  i4 = standardlocale.toInt(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'i4' from XML file.",valc); }
+}
+if (!strcmp(namec, "i5")) {
+  i5 = standardlocale.toInt(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'i5' from XML file.",valc); }
+}
+if (!strcmp(namec, "s1")) {
+  if (s1) { free(s1); }
+  s1=strdup(valc);
+}
+if (!strcmp(namec, "s2")) {
+  if (s2) { free(s2); }
+  s2=strdup(valc);
+}
+if (!strcmp(namec, "s3")) {
+  if (s3) { free(s3); }
+  s3=strdup(valc);
+}
+if (!strcmp(namec, "b1")) {
+b1 = strtobool(valc);
+}
+if (!strcmp(namec, "b2")) {
+b2 = strtobool(valc);
+}
+if (!strcmp(namec, "b3")) {
+b3 = strtobool(valc);
+}
+if (!strcmp(namec, "b4")) {
+b4 = strtobool(valc);
+}
+if (!strcmp(namec, "dir1")) {
+  if (dir1) { free(dir1); }
+  dir1=strdup(valc);
+}
+if (!strcmp(namec, "dir2")) {
+  if (dir2) { free(dir2); }
+  dir2=strdup(valc);
+}
+if (!strcmp(namec, "export_interval")) {
+  export_interval = standardlocale.toInt(valc, &ok);
+  if (!ok) { MyWarning::error("Read error: cannot convert string \"%s\" to integer while reading parameter 'export_interval' from XML file.",valc); }
+}
+if (!strcmp(namec, "export_fn_prefix")) {
+  if (export_fn_prefix) { free(export_fn_prefix); }
+  export_fn_prefix=strdup(valc);
+}
+}
+void Parameter::AssignValArrayToPar(const char *namec, vector<double> valarray) {
+if (!strcmp(namec, "D")) {
+  int i=0;
+  vector<double>::const_iterator v=valarray.begin();
+  while (v!=valarray.end() && i <= 14 ) {
+     D[i++]=*(v++);
   }
 }
-void Parameter::AssignValArrayToPar(const char *namec, vector<double> valarray) {
-  if (!strcmp(namec, "D")) {
-    int i=0;
-    vector<double>::const_iterator v=valarray.begin();
-    while (v!=valarray.end() && i <= 14 ) {
-      D[i++]=*(v++);
-    }
+if (!strcmp(namec, "initval")) {
+  int i=0;
+  vector<double>::const_iterator v=valarray.begin();
+  while (v!=valarray.end() && i <= 14 ) {
+     initval[i++]=*(v++);
   }
-  if (!strcmp(namec, "initval")) {
-    int i=0;
-    vector<double>::const_iterator v=valarray.begin();
-    while (v!=valarray.end() && i <= 14 ) {
-      initval[i++]=*(v++);
-    }
+}
+if (!strcmp(namec, "k")) {
+  int i=0;
+  vector<double>::const_iterator v=valarray.begin();
+  while (v!=valarray.end() && i <= 14 ) {
+     k[i++]=*(v++);
   }
-  if (!strcmp(namec, "k")) {
-    int i=0;
-    vector<double>::const_iterator v=valarray.begin();
-    while (v!=valarray.end() && i <= 14 ) {
-      k[i++]=*(v++);
-    }
-  }
+}
 }
 
 ostream &operator<<(ostream &os, Parameter &p) {
-  p.Write(os);
-  return os;
+    p.Write(os);
+    return os;
 }
 
 /* finis */
diff --git a/src/parameter.h b/src/parameter.h
--- a/src/parameter.h
+++ b/src/parameter.h
@@ -21,6 +21,9 @@
  *
  */
 
+// WARNING: This file is automatically generated by make_parameter_source.pl. Do not edit.
+// All edits will be discarded.
+
 #ifndef _PARAMETER_H_
 #define _PARAMETER_H_
 #include "vector.h"
@@ -139,6 +142,8 @@
   bool b4;
   char * dir1;
   char * dir2;
+  int export_interval;
+  char * export_fn_prefix;
  private:
  };
 
@@ -147,5 +152,3 @@
 
 
 #endif
-
-/* finis */
diff --git a/src/pardialog.cpp b/src/pardialog.cpp
--- a/src/pardialog.cpp
+++ b/src/pardialog.cpp
@@ -19,6 +19,8 @@
  *
  */
 
+// Do not edit. All edits will be discarded.
+
 #include "pardialog.h"
 #include "parameter.h"
 #include <cstring>
@@ -30,7 +32,7 @@
 static const std::string _module_id("$Id$");
 
 ParameterDialog::ParameterDialog(QWidget *parent, const char *name, Qt::WindowFlags f) : QDialog(parent,name,false,f) {
-  extern Parameter par;
+    extern Parameter par;
   arrowcolor_edit = new QLineEdit( QString("%1").arg(par.arrowcolor), this, "arrowcolor_edit" );
   arrowsize_edit = new QLineEdit( QString("%1").arg(par.arrowsize), this, "arrowsize_edit" );
   textcolor_edit = new QLineEdit( QString("%1").arg(par.textcolor), this, "textcolor_edit" );
@@ -135,12 +137,14 @@ ParameterDialog::ParameterDialog(QWidget
   b4_edit = new QLineEdit( QString("%1").arg(sbool(par.b4)), this, "b4_edit" );
   dir1_edit = new QLineEdit( QString("%1").arg(par.dir1), this, "dir1_edit" );
   dir2_edit = new QLineEdit( QString("%1").arg(par.dir2), this, "dir2_edit" );
-  // make a 1x1 grid; it will auto-expand
-  QGridLayout *grid = new QGridLayout( this, 1, 1 );
+  export_interval_edit = new QLineEdit( QString("%1").arg(par.export_interval), this, "export_interval_edit" );
+  export_fn_prefix_edit = new QLineEdit( QString("%1").arg(par.export_fn_prefix), this, "export_fn_prefix_edit" );
+// make a 1x1 grid; it will auto-expand
+QGridLayout *grid = new QGridLayout( this, 1, 1 );
     
-  // add the first four widgets with (row, column) addressing
-  setWindowTitle( QString( " Parameter values for VirtualLeaf") );
-  grid->addWidget( new QLabel( "<h3> Parameter values for VirtualLeaf</h3>",this), 0, 0, 1, -1, Qt::AlignCenter);
+// add the first four widgets with (row, column) addressing
+  setWindowTitle( QString( " Parameter values for The Virtual Leaf") );
+  grid->addWidget( new QLabel( "<h3> Parameter values for The Virtual Leaf</h3>",this), 0, 0, 1, -1, Qt::AlignCenter);
   grid->addWidget( new QLabel( "", this), 0+1, 0, 1, -1);
   grid->addWidget( new QLabel( " <b>Visualization</b>", this), 3, 0, 1, 2 );
   grid->addWidget( new QLabel( "arrowcolor", this ),4, 0 );
@@ -349,123 +353,129 @@ ParameterDialog::ParameterDialog(QWidget
   grid->addWidget( dir1_edit, 29, 6+1  );
   grid->addWidget( new QLabel( "dir2", this ),3, 8 );
   grid->addWidget( dir2_edit, 3, 8+1  );
-  QPushButton *pb = new QPushButton( "&Write", this );
-  grid->addWidget(pb, 31, 6 );
-  connect( pb, SIGNAL( clicked() ), this, SLOT( write() ) );
-  QPushButton *pb2 = new QPushButton( "&Close", this );
-  grid->addWidget(pb2,31, 6+1 );
-  connect( pb2, SIGNAL( clicked() ), this, SLOT( close() ) );
-  QPushButton *pb3 = new QPushButton( "&Reset", this );
-  grid->addWidget(pb3, 31, 6+2 );
-  connect( pb3, SIGNAL( clicked() ), this, SLOT( Reset() ) );
-  show();
+  grid->addWidget( new QLabel( "export_interval", this ),4, 8 );
+  grid->addWidget( export_interval_edit, 4, 8+1  );
+  grid->addWidget( new QLabel( "export_fn_prefix", this ),5, 8 );
+  grid->addWidget( export_fn_prefix_edit, 5, 8+1  );
+QPushButton *pb = new QPushButton( "&Write", this );
+grid->addWidget(pb, 31, 6 );
+connect( pb, SIGNAL( clicked() ), this, SLOT( write() ) );
+QPushButton *pb2 = new QPushButton( "&Close", this );
+grid->addWidget(pb2,31, 6+1 );
+connect( pb2, SIGNAL( clicked() ), this, SLOT( close() ) );
+QPushButton *pb3 = new QPushButton( "&Reset", this );
+grid->addWidget(pb3, 31, 6+2 );
+connect( pb3, SIGNAL( clicked() ), this, SLOT( Reset() ) );
+show();
 };
 
 ParameterDialog::~ParameterDialog(void) {
-  delete arrowcolor_edit;
-  delete arrowsize_edit;
-  delete textcolor_edit;
-  delete cellnumsize_edit;
-  delete nodenumsize_edit;
-  delete node_mag_edit;
-  delete outlinewidth_edit;
-  delete cell_outline_color_edit;
-  delete resize_stride_edit;
-  delete T_edit;
-  delete lambda_length_edit;
-  delete lambda_celllength_edit;
-  delete target_length_edit;
-  delete cell_expansion_rate_edit;
-  delete cell_div_expansion_rate_edit;
-  delete auxin_dependent_growth_edit;
-  delete ode_accuracy_edit;
-  delete mc_stepsize_edit;
-  delete mc_cell_stepsize_edit;
-  delete energy_threshold_edit;
-  delete bend_lambda_edit;
-  delete alignment_lambda_edit;
-  delete rel_cell_div_threshold_edit;
-  delete rel_perimeter_stiffness_edit;
-  delete collapse_node_threshold_edit;
-  delete morphogen_div_threshold_edit;
-  delete morphogen_expansion_threshold_edit;
-  delete copy_wall_edit;
-  delete source_edit;
-  delete D_edit;
-  delete initval_edit;
-  delete k1_edit;
-  delete k2_edit;
-  delete r_edit;
-  delete kr_edit;
-  delete km_edit;
-  delete Pi_tot_edit;
-  delete transport_edit;
-  delete ka_edit;
-  delete pin_prod_edit;
-  delete pin_prod_in_epidermis_edit;
-  delete pin_breakdown_edit;
-  delete pin_breakdown_internal_edit;
-  delete aux1prod_edit;
-  delete aux1prodmeso_edit;
-  delete aux1decay_edit;
-  delete aux1decaymeso_edit;
-  delete aux1transport_edit;
-  delete aux_cons_edit;
-  delete aux_breakdown_edit;
-  delete kaux1_edit;
-  delete kap_edit;
-  delete leaf_tip_source_edit;
-  delete sam_efflux_edit;
-  delete sam_auxin_edit;
-  delete sam_auxin_breakdown_edit;
-  delete van3prod_edit;
-  delete van3autokat_edit;
-  delete van3sat_edit;
-  delete k2van3_edit;
-  delete dt_edit;
-  delete rd_dt_edit;
-  delete datadir_edit;
-  delete movie_edit;
-  delete nit_edit;
-  delete maxt_edit;
-  delete storage_stride_edit;
-  delete xml_storage_stride_edit;
-  delete rseed_edit;
-  delete constituous_expansion_limit_edit;
-  delete vessel_inh_level_edit;
-  delete vessel_expansion_rate_edit;
-  delete d_edit;
-  delete e_edit;
-  delete f_edit;
-  delete c_edit;
-  delete mu_edit;
-  delete nu_edit;
-  delete rho0_edit;
-  delete rho1_edit;
-  delete c0_edit;
-  delete gamma_edit;
-  delete eps_edit;
-  delete k_edit;
-  delete i1_edit;
-  delete i2_edit;
-  delete i3_edit;
-  delete i4_edit;
-  delete i5_edit;
-  delete s1_edit;
-  delete s2_edit;
-  delete s3_edit;
-  delete b1_edit;
-  delete b2_edit;
-  delete b3_edit;
-  delete b4_edit;
-  delete dir1_edit;
-  delete dir2_edit;
+delete arrowcolor_edit;
+delete arrowsize_edit;
+delete textcolor_edit;
+delete cellnumsize_edit;
+delete nodenumsize_edit;
+delete node_mag_edit;
+delete outlinewidth_edit;
+delete cell_outline_color_edit;
+delete resize_stride_edit;
+delete T_edit;
+delete lambda_length_edit;
+delete lambda_celllength_edit;
+delete target_length_edit;
+delete cell_expansion_rate_edit;
+delete cell_div_expansion_rate_edit;
+delete auxin_dependent_growth_edit;
+delete ode_accuracy_edit;
+delete mc_stepsize_edit;
+delete mc_cell_stepsize_edit;
+delete energy_threshold_edit;
+delete bend_lambda_edit;
+delete alignment_lambda_edit;
+delete rel_cell_div_threshold_edit;
+delete rel_perimeter_stiffness_edit;
+delete collapse_node_threshold_edit;
+delete morphogen_div_threshold_edit;
+delete morphogen_expansion_threshold_edit;
+delete copy_wall_edit;
+delete source_edit;
+delete D_edit;
+delete initval_edit;
+delete k1_edit;
+delete k2_edit;
+delete r_edit;
+delete kr_edit;
+delete km_edit;
+delete Pi_tot_edit;
+delete transport_edit;
+delete ka_edit;
+delete pin_prod_edit;
+delete pin_prod_in_epidermis_edit;
+delete pin_breakdown_edit;
+delete pin_breakdown_internal_edit;
+delete aux1prod_edit;
+delete aux1prodmeso_edit;
+delete aux1decay_edit;
+delete aux1decaymeso_edit;
+delete aux1transport_edit;
+delete aux_cons_edit;
+delete aux_breakdown_edit;
+delete kaux1_edit;
+delete kap_edit;
+delete leaf_tip_source_edit;
+delete sam_efflux_edit;
+delete sam_auxin_edit;
+delete sam_auxin_breakdown_edit;
+delete van3prod_edit;
+delete van3autokat_edit;
+delete van3sat_edit;
+delete k2van3_edit;
+delete dt_edit;
+delete rd_dt_edit;
+delete datadir_edit;
+delete movie_edit;
+delete nit_edit;
+delete maxt_edit;
+delete storage_stride_edit;
+delete xml_storage_stride_edit;
+delete rseed_edit;
+delete constituous_expansion_limit_edit;
+delete vessel_inh_level_edit;
+delete vessel_expansion_rate_edit;
+delete d_edit;
+delete e_edit;
+delete f_edit;
+delete c_edit;
+delete mu_edit;
+delete nu_edit;
+delete rho0_edit;
+delete rho1_edit;
+delete c0_edit;
+delete gamma_edit;
+delete eps_edit;
+delete k_edit;
+delete i1_edit;
+delete i2_edit;
+delete i3_edit;
+delete i4_edit;
+delete i5_edit;
+delete s1_edit;
+delete s2_edit;
+delete s3_edit;
+delete b1_edit;
+delete b2_edit;
+delete b3_edit;
+delete b4_edit;
+delete dir1_edit;
+delete dir2_edit;
+delete export_interval_edit;
+delete export_fn_prefix_edit;
 }
 
 void ParameterDialog::write(void) {
     
-  extern Parameter par;
-  QString tmpval;
+extern Parameter par;
+QString tmpval;
   par.arrowcolor = strdup((const char *)arrowcolor_edit->text());
   par.arrowsize = arrowsize_edit->text().toDouble();
   par.textcolor = strdup((const char *)textcolor_edit->text());
@@ -486,7 +496,7 @@ void ParameterDialog::write(void) {
   else if (tmpval == "false" || tmpval == "no") par.auxin_dependent_growth = false;
   else {
     if (QMessageBox::question(this, "Syntax error", tr("Value %1 of parameter %2 is not recognized as Boolean.\nDo you mean TRUE or FALSE?").arg(tmpval).arg("auxin_dependent_growth"),"True","False", QString::null, 0, 1)==0) par.auxin_dependent_growth=true;
-    else par.auxin_dependent_growth=false;
+      else par.auxin_dependent_growth=false;
   }
   par.ode_accuracy = ode_accuracy_edit->text().toDouble();
   par.mc_stepsize = mc_stepsize_edit->text().toDouble();
@@ -504,7 +514,7 @@ void ParameterDialog::write(void) {
   else if (tmpval == "false" || tmpval == "no") par.copy_wall = false;
   else {
     if (QMessageBox::question(this, "Syntax error", tr("Value %1 of parameter %2 is not recognized as Boolean.\nDo you mean TRUE or FALSE?").arg(tmpval).arg("copy_wall"),"True","False", QString::null, 0, 1)==0) par.copy_wall=true;
-    else par.copy_wall=false;
+      else par.copy_wall=false;
   }
   par.source = source_edit->text().toDouble();
   tmpval = D_edit->text().section(',', 0, 0);
@@ -604,7 +614,7 @@ void ParameterDialog::write(void) {
   else if (tmpval == "false" || tmpval == "no") par.movie = false;
   else {
     if (QMessageBox::question(this, "Syntax error", tr("Value %1 of parameter %2 is not recognized as Boolean.\nDo you mean TRUE or FALSE?").arg(tmpval).arg("movie"),"True","False", QString::null, 0, 1)==0) par.movie=true;
-    else par.movie=false;
+      else par.movie=false;
   }
   par.nit = nit_edit->text().toInt();
   par.maxt = maxt_edit->text().toDouble();
@@ -668,32 +678,34 @@ void ParameterDialog::write(void) {
   else if (tmpval == "false" || tmpval == "no") par.b1 = false;
   else {
     if (QMessageBox::question(this, "Syntax error", tr("Value %1 of parameter %2 is not recognized as Boolean.\nDo you mean TRUE or FALSE?").arg(tmpval).arg("b1"),"True","False", QString::null, 0, 1)==0) par.b1=true;
-    else par.b1=false;
+      else par.b1=false;
   }
   tmpval = b2_edit->text().stripWhiteSpace();
   if (tmpval == "true" || tmpval == "yes" ) par.b2 = true;
   else if (tmpval == "false" || tmpval == "no") par.b2 = false;
   else {
     if (QMessageBox::question(this, "Syntax error", tr("Value %1 of parameter %2 is not recognized as Boolean.\nDo you mean TRUE or FALSE?").arg(tmpval).arg("b2"),"True","False", QString::null, 0, 1)==0) par.b2=true;
-    else par.b2=false;
+      else par.b2=false;
   }
   tmpval = b3_edit->text().stripWhiteSpace();
   if (tmpval == "true" || tmpval == "yes" ) par.b3 = true;
   else if (tmpval == "false" || tmpval == "no") par.b3 = false;
   else {
     if (QMessageBox::question(this, "Syntax error", tr("Value %1 of parameter %2 is not recognized as Boolean.\nDo you mean TRUE or FALSE?").arg(tmpval).arg("b3"),"True","False", QString::null, 0, 1)==0) par.b3=true;
-    else par.b3=false;
+      else par.b3=false;
   }
   tmpval = b4_edit->text().stripWhiteSpace();
   if (tmpval == "true" || tmpval == "yes" ) par.b4 = true;
   else if (tmpval == "false" || tmpval == "no") par.b4 = false;
   else {
     if (QMessageBox::question(this, "Syntax error", tr("Value %1 of parameter %2 is not recognized as Boolean.\nDo you mean TRUE or FALSE?").arg(tmpval).arg("b4"),"True","False", QString::null, 0, 1)==0) par.b4=true;
-    else par.b4=false;
+      else par.b4=false;
   }
   par.dir1 = strdup((const char *)dir1_edit->text());
   par.dir2 = strdup((const char *)dir2_edit->text());
-  Reset();
+  par.export_interval = export_interval_edit->text().toInt();
+  par.export_fn_prefix = strdup((const char *)export_fn_prefix_edit->text());
+Reset();
 
 }
 void ParameterDialog::Reset(void) {
@@ -802,6 +814,8 @@ void ParameterDialog::Reset(void) {
   b4_edit->setText( QString("%1").arg(sbool(par.b4)));
   dir1_edit->setText( QString("%1").arg(par.dir1) );
   dir2_edit->setText( QString("%1").arg(par.dir2) );
+  export_interval_edit->setText( QString("%1").arg(par.export_interval) );
+  export_fn_prefix_edit->setText( QString("%1").arg(par.export_fn_prefix) );
 }
 
 /* finis */
diff --git a/src/pardialog.h b/src/pardialog.h
--- a/src/pardialog.h
+++ b/src/pardialog.h
@@ -1,3 +1,29 @@
+/*
+ *
+ *  $Id$
+ *
+ *  This file is part of the Virtual Leaf.
+ *
+ *  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.
+ *
+ *  VirtualLeaf is distributed in the hope that it will be useful,
+ *  but WITHOUT ANY WARRANTY; without even the implied warranty of
+ *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ *  GNU General Public License for more details.
+ *
+ *  You should have received a copy of the GNU General Public License
+ *  along with the Virtual Leaf.  If not, see <http://www.gnu.org/licenses/>.
+ *
+ *  Copyright 2010 Roeland Merks.
+ *
+ */
+
+// WARNING: This file is automatically generated by make_parameter_source.pl. Do not edit.
+// All edits will be discarded.
+
 #ifndef PARAMETER_DIALOG_H
 #define PARAMETER_DIALOG_H
 #include <qdialog.h>
@@ -118,5 +144,8 @@ class ParameterDialog : public QDialog {
   QLineEdit *b4_edit;
   QLineEdit *dir1_edit;
   QLineEdit *dir2_edit;
+  QLineEdit *export_interval_edit;
+  QLineEdit *export_fn_prefix_edit;
 };
 #endif
+
diff --git a/src/perl/ChangeLog b/src/perl/ChangeLog
--- a/src/perl/ChangeLog
+++ b/src/perl/ChangeLog
@@ -1,3 +1,8 @@
+2010-10-08    <guravage@caterpie.sen.cwi.nl>
+
+	* make_pardialog_source.pl:
+	* make_parameter_source.pl: Added GPL boilerplate and mercurial keywords.
+
 2010-06-25    <guravage@caterpie.sen.cwi.nl>
 
 	* make_parameter_source.pl: Added datadir changes.
diff --git a/src/perl/make_parameter_source.pl b/src/perl/make_parameter_source.pl
--- a/src/perl/make_parameter_source.pl
+++ b/src/perl/make_parameter_source.pl
@@ -88,6 +88,27 @@ while (<parfile>) {
 $lines=$i;
 
 print cppfile <<END_HEADER;
+/*
+ *
+ *  This file is part of the Virtual Leaf.
+ *
+ *  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.
+ *
+ *  VirtualLeaf is distributed in the hope that it will be useful,
+ *  but WITHOUT ANY WARRANTY; without even the implied warranty of
+ *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ *  GNU General Public License for more details.
+ *
+ *  You should have received a copy of the GNU General Public License
+ *  along with the Virtual Leaf.  If not, see <http://www.gnu.org/licenses/>.
+ *
+ *  Copyright 2010 Roeland Merks.
+ *
+ */
+
 // WARNING: This file is automatically generated by make_parameter_source.pl.
 // Do not edit. All edits will be discarded.
 
@@ -108,6 +129,8 @@ print cppfile <<END_HEADER;
 
 using namespace std;
 
+static const std::string _module_id("\$Id\$");
+
 Parameter::Parameter() {
 END_HEADER
 
@@ -373,8 +396,31 @@ END_TRAILER
 
 open hfile, ">parameter.h";
 print hfile <<END_HEADER2;
+/*
+ *
+ *  \$Id\$
+ *
+ *  This file is part of the Virtual Leaf.
+ *
+ *  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.
+ *
+ *  VirtualLeaf is distributed in the hope that it will be useful,
+ *  but WITHOUT ANY WARRANTY; without even the implied warranty of
+ *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ *  GNU General Public License for more details.
+ *
+ *  You should have received a copy of the GNU General Public License
+ *  along with the Virtual Leaf.  If not, see <http://www.gnu.org/licenses/>.
+ *
+ *  Copyright 2010 Roeland Merks.
+ *
+ */
+
 // WARNING: This file is automatically generated by make_parameter_source.pl. Do not edit.
-// All edits will be discarded.
+// Do not edit. All edits will be discarded.
 
 #ifndef _PARAMETER_H_
 #define _PARAMETER_H_
diff --git a/src/perl/make_pardialog_source.pl b/src/perl/make_pardialog_source.pl
--- a/src/perl/make_pardialog_source.pl
+++ b/src/perl/make_pardialog_source.pl
@@ -88,6 +88,30 @@ while (<parfile>) {
 
 
 print cppfile <<END_HEADER;
+/*
+ *
+ *  This file is part of the Virtual Leaf.
+ *
+ *  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.
+ *
+ *  VirtualLeaf is distributed in the hope that it will be useful,
+ *  but WITHOUT ANY WARRANTY; without even the implied warranty of
+ *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ *  GNU General Public License for more details.
+ *
+ *  You should have received a copy of the GNU General Public License
+ *  along with the Virtual Leaf.  If not, see <http://www.gnu.org/licenses/>.
+ *
+ *  Copyright 2010 Roeland Merks.
+ *
+ */
+
+// WARNING: This file is automatically generated by make_parameter_source.pl. Do not edit.
+// Do not edit. All edits will be discarded.
+
 #include "$basename.h"
 #include "parameter.h"
 #include <cstring>
@@ -96,6 +120,8 @@ print cppfile <<END_HEADER;
 #include <qlineedit.h>
 #include <qmessagebox.h>
 
+static const std::string _module_id("\$Id\$");
+
 ParameterDialog::ParameterDialog(QWidget *parent, const char *name, Qt::WindowFlags f) : QDialog(parent,name,false,f) {
     extern Parameter par;
 END_HEADER
@@ -265,6 +291,32 @@ print cppfile "}\n\n";
 
 open hfile, ">$basename.h";
 print hfile <<END_HEADER2;
+/*
+ *
+ *  \$Id\$
+ *
+ *  This file is part of the Virtual Leaf.
+ *
+ *  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.
+ *
+ *  VirtualLeaf is distributed in the hope that it will be useful,
+ *  but WITHOUT ANY WARRANTY; without even the implied warranty of
+ *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ *  GNU General Public License for more details.
+ *
+ *  You should have received a copy of the GNU General Public License
+ *  along with the Virtual Leaf.  If not, see <http://www.gnu.org/licenses/>.
+ *
+ *  Copyright 2010 Roeland Merks.
+ *
+ */
+
+// WARNING: This file is automatically generated by make_parameter_source.pl. Do not edit.
+// Do not edit. All edits will be discarded.
+
 #ifndef PARAMETER_DIALOG_H
 #define PARAMETER_DIALOG_H
 #include <qdialog.h>