Changeset 042f82 for src/config.cpp

Timestamp:

Jul 23, 2009, 2:21:57 PM (16 years ago)

Author:

Frederik Heber <heber@…>

Branches:

Action_Thermostats, Add_AtomRandomPerturbation, Add_FitFragmentPartialChargesAction, Add_RotateAroundBondAction, Add_SelectAtomByNameAction, Added_ParseSaveFragmentResults, AddingActions_SaveParseParticleParameters, Adding_Graph_to_ChangeBondActions, Adding_MD_integration_tests, Adding_ParticleName_to_Atom, Adding_StructOpt_integration_tests, AtomFragments, Automaking_mpqc_open, AutomationFragmentation_failures, Candidate_v1.5.4, Candidate_v1.6.0, Candidate_v1.6.1, ChangeBugEmailaddress, ChangingTestPorts, ChemicalSpaceEvaluator, CombiningParticlePotentialParsing, Combining_Subpackages, Debian_Package_split, Debian_package_split_molecuildergui_only, Disabling_MemDebug, Docu_Python_wait, EmpiricalPotential_contain_HomologyGraph, EmpiricalPotential_contain_HomologyGraph_documentation, Enable_parallel_make_install, Enhance_userguide, Enhanced_StructuralOptimization, Enhanced_StructuralOptimization_continued, Example_ManyWaysToTranslateAtom, Exclude_Hydrogens_annealWithBondGraph, FitPartialCharges_GlobalError, Fix_BoundInBox_CenterInBox_MoleculeActions, Fix_ChargeSampling_PBC, Fix_ChronosMutex, Fix_FitPartialCharges, Fix_FitPotential_needs_atomicnumbers, Fix_ForceAnnealing, Fix_IndependentFragmentGrids, Fix_ParseParticles, Fix_ParseParticles_split_forward_backward_Actions, Fix_PopActions, Fix_QtFragmentList_sorted_selection, Fix_Restrictedkeyset_FragmentMolecule, Fix_StatusMsg, Fix_StepWorldTime_single_argument, Fix_Verbose_Codepatterns, Fix_fitting_potentials, Fixes, ForceAnnealing_goodresults, ForceAnnealing_oldresults, ForceAnnealing_tocheck, ForceAnnealing_with_BondGraph, ForceAnnealing_with_BondGraph_continued, ForceAnnealing_with_BondGraph_continued_betteresults, ForceAnnealing_with_BondGraph_contraction-expansion, FragmentAction_writes_AtomFragments, FragmentMolecule_checks_bonddegrees, GeometryObjects, Gui_Fixes, Gui_displays_atomic_force_velocity, ImplicitCharges, IndependentFragmentGrids, IndependentFragmentGrids_IndividualZeroInstances, IndependentFragmentGrids_IntegrationTest, IndependentFragmentGrids_Sole_NN_Calculation, JobMarket_RobustOnKillsSegFaults, JobMarket_StableWorkerPool, JobMarket_unresolvable_hostname_fix, MoreRobust_FragmentAutomation, ODR_violation_mpqc_open, PartialCharges_OrthogonalSummation, PdbParser_setsAtomName, PythonUI_with_named_parameters, QtGui_reactivate_TimeChanged_changes, Recreated_GuiChecks, Rewrite_FitPartialCharges, RotateToPrincipalAxisSystem_UndoRedo, SaturateAtoms_findBestMatching, SaturateAtoms_singleDegree, StoppableMakroAction, Subpackage_CodePatterns, Subpackage_JobMarket, Subpackage_LinearAlgebra, Subpackage_levmar, Subpackage_mpqc_open, Subpackage_vmg, Switchable_LogView, ThirdParty_MPQC_rebuilt_buildsystem, TrajectoryDependenant_MaxOrder, TremoloParser_IncreasedPrecision, TremoloParser_MultipleTimesteps, TremoloParser_setsAtomName, Ubuntu_1604_changes, stable

Children:

36ec71

Parents:

205ccd

Message:

fixed indentation from tabs to two spaces.

File:

• src/config.cpp (modified) (14 diffs)

src/config.cpp

@@ -163 +163 @@
 config::config()
 {
-        mainname = (char *) MallocString(sizeof(char)*MAXSTRINGSIZE,"config constructor: mainname");
-        defaultpath = (char *) MallocString(sizeof(char)*MAXSTRINGSIZE,"config constructor: defaultpath");
-        pseudopotpath = (char *) MallocString(sizeof(char)*MAXSTRINGSIZE,"config constructor: pseudopotpath");
+  mainname = (char *) MallocString(sizeof(char)*MAXSTRINGSIZE,"config constructor: mainname");
+  defaultpath = (char *) MallocString(sizeof(char)*MAXSTRINGSIZE,"config constructor: defaultpath");
+  pseudopotpath = (char *) MallocString(sizeof(char)*MAXSTRINGSIZE,"config constructor: pseudopotpath");
   databasepath = (char *) MallocString(sizeof(char)*MAXSTRINGSIZE,"config constructor: databasepath");
-        configpath = (char *) MallocString(sizeof(char)*MAXSTRINGSIZE,"config constructor: configpath");
-        configname = (char *) MallocString(sizeof(char)*MAXSTRINGSIZE,"config constructor: configname");
-        strcpy(mainname,"pcp");
-        strcpy(defaultpath,"not specified");
-        strcpy(pseudopotpath,"not specified");
-        configpath[0]='\0';
-        configname[0]='\0';
-        basis = "3-21G";
-
-        FastParsing = false;
-        ProcPEGamma=8;
-        ProcPEPsi=1;
-        DoOutVis=0;
-        DoOutMes=1;
-        DoOutNICS=0;
-        DoOutOrbitals=0;
-        DoOutCurrent=0;
-        DoPerturbation=0;
-        DoFullCurrent=0;
-        DoWannier=0;
-        CommonWannier=0;
-        SawtoothStart=0.01;
-        VectorPlane=0;
-        VectorCut=0;
-        UseAddGramSch=1;
-        Seed=1;
-
-        MaxOuterStep=0;
-        Deltat=1;
-        OutVisStep=10;
-        OutSrcStep=5;
-        TargetTemp=0.00095004455;
-        ScaleTempStep=25;
-        MaxPsiStep=0;
-        EpsWannier=1e-7;
-
-        MaxMinStep=100;
-        RelEpsTotalEnergy=1e-7;
-        RelEpsKineticEnergy=1e-5;
-        MaxMinStopStep=1;
-        MaxMinGapStopStep=0;
-        MaxInitMinStep=100;
-        InitRelEpsTotalEnergy=1e-5;
-        InitRelEpsKineticEnergy=1e-4;
-        InitMaxMinStopStep=1;
-        InitMaxMinGapStopStep=0;
-
-        //BoxLength[NDIM*NDIM];
-
-        ECut=128.;
-        MaxLevel=5;
-        RiemannTensor=0;
-        LevRFactor=0;
-        RiemannLevel=0;
-        Lev0Factor=2;
-        RTActualUse=0;
-        PsiType=0;
-        MaxPsiDouble=0;
-        PsiMaxNoUp=0;
-        PsiMaxNoDown=0;
-        AddPsis=0;
-
-        RCut=20.;
-        StructOpt=0;
-        IsAngstroem=1;
-        RelativeCoord=0;
-        MaxTypes=0;
+  configpath = (char *) MallocString(sizeof(char)*MAXSTRINGSIZE,"config constructor: configpath");
+  configname = (char *) MallocString(sizeof(char)*MAXSTRINGSIZE,"config constructor: configname");
+  strcpy(mainname,"pcp");
+  strcpy(defaultpath,"not specified");
+  strcpy(pseudopotpath,"not specified");
+  configpath[0]='\0';
+  configname[0]='\0';
+  basis = "3-21G";
+
+  FastParsing = false;
+  ProcPEGamma=8;
+  ProcPEPsi=1;
+  DoOutVis=0;
+  DoOutMes=1;
+  DoOutNICS=0;
+  DoOutOrbitals=0;
+  DoOutCurrent=0;
+  DoPerturbation=0;
+  DoFullCurrent=0;
+  DoWannier=0;
+  CommonWannier=0;
+  SawtoothStart=0.01;
+  VectorPlane=0;
+  VectorCut=0;
+  UseAddGramSch=1;
+  Seed=1;
+
+  MaxOuterStep=0;
+  Deltat=1;
+  OutVisStep=10;
+  OutSrcStep=5;
+  TargetTemp=0.00095004455;
+  ScaleTempStep=25;
+  MaxPsiStep=0;
+  EpsWannier=1e-7;
+
+  MaxMinStep=100;
+  RelEpsTotalEnergy=1e-7;
+  RelEpsKineticEnergy=1e-5;
+  MaxMinStopStep=1;
+  MaxMinGapStopStep=0;
+  MaxInitMinStep=100;
+  InitRelEpsTotalEnergy=1e-5;
+  InitRelEpsKineticEnergy=1e-4;
+  InitMaxMinStopStep=1;
+  InitMaxMinGapStopStep=0;
+
+  //BoxLength[NDIM*NDIM];
+
+  ECut=128.;
+  MaxLevel=5;
+  RiemannTensor=0;
+  LevRFactor=0;
+  RiemannLevel=0;
+  Lev0Factor=2;
+  RTActualUse=0;
+  PsiType=0;
+  MaxPsiDouble=0;
+  PsiMaxNoUp=0;
+  PsiMaxNoDown=0;
+  AddPsis=0;
+
+  RCut=20.;
+  StructOpt=0;
+  IsAngstroem=1;
+  RelativeCoord=0;
+  MaxTypes=0;
 };
 
@@ -241 +241 @@
 config::~config()
 {
-        Free((void **)&mainname, "config::~config: *mainname");
-        Free((void **)&defaultpath, "config::~config: *defaultpath");
-        Free((void **)&pseudopotpath, "config::~config: *pseudopotpath");
+  Free((void **)&mainname, "config::~config: *mainname");
+  Free((void **)&defaultpath, "config::~config: *defaultpath");
+  Free((void **)&pseudopotpath, "config::~config: *pseudopotpath");
   Free((void **)&databasepath, "config::~config: *databasepath");
-        Free((void **)&configpath, "config::~config: *configpath");
-        Free((void **)&configname, "config::~config: *configname");
+  Free((void **)&configpath, "config::~config: *configpath");
+  Free((void **)&configname, "config::~config: *configname");
 };
 
@@ -255 +255 @@
 void config::Edit()
 {
-        char choice;
-
-        do {
-                cout << Verbose(0) << "===========EDIT CONFIGURATION============================" << endl;
-                cout << Verbose(0) << " A - mainname (prefix for all runtime files)" << endl;
-                cout << Verbose(0) << " B - Default path (for runtime files)" << endl;
-                cout << Verbose(0) << " C - Path of pseudopotential files" << endl;
-                cout << Verbose(0) << " D - Number of coefficient sharing processes" << endl;
-                cout << Verbose(0) << " E - Number of wave function sharing processes" << endl;
-                cout << Verbose(0) << " F - 0: Don't output density for OpenDX, 1: do" << endl;
-                cout << Verbose(0) << " G - 0: Don't output physical data, 1: do" << endl;
-                cout << Verbose(0) << " H - 0: Don't output densities of each unperturbed orbital for OpenDX, 1: do" << endl;
-                cout << Verbose(0) << " I - 0: Don't output current density for OpenDX, 1: do" << endl;
-                cout << Verbose(0) << " J - 0: Don't do the full current calculation, 1: do" << endl;
-                cout << Verbose(0) << " K - 0: Don't do perturbation calculation to obtain susceptibility and shielding, 1: do" << endl;
-                cout << Verbose(0) << " L - 0: Wannier centres as calculated, 1: common centre for all, 2: unite centres according to spread, 3: cell centre, 4: shifted to nearest grid point" << endl;
-                cout << Verbose(0) << " M - Absolute begin of unphysical sawtooth transfer for position operator within cell" << endl;
-                cout << Verbose(0) << " N - (0,1,2) x,y,z-plane to do two-dimensional current vector cut" << endl;
-                cout << Verbose(0) << " O - Absolute position along vector cut axis for cut plane" << endl;
-                cout << Verbose(0) << " P - Additional Gram-Schmidt-Orthonormalization to stabilize numerics" << endl;
-                cout << Verbose(0) << " Q - Initial integer value of random number generator" << endl;
-                cout << Verbose(0) << " R - for perturbation 0, for structure optimization defines upper limit of iterations" << endl;
-                cout << Verbose(0) << " T - Output visual after ...th step" << endl;
-                cout << Verbose(0) << " U - Output source densities of wave functions after ...th step" << endl;
-                cout << Verbose(0) << " X - minimization iterations per wave function, if unsure leave at default value 0" << endl;
-                cout << Verbose(0) << " Y - tolerance value for total spread in iterative Jacobi diagonalization" << endl;
-                cout << Verbose(0) << " Z - Maximum number of minimization iterations" << endl;
-                cout << Verbose(0) << " a - Relative change in total energy to stop min. iteration" << endl;
-                cout << Verbose(0) << " b - Relative change in kinetic energy to stop min. iteration" << endl;
-                cout << Verbose(0) << " c - Check stop conditions every ..th step during min. iteration" << endl;
-                cout << Verbose(0) << " e - Maximum number of minimization iterations during initial level" << endl;
-                cout << Verbose(0) << " f - Relative change in total energy to stop min. iteration during initial level" << endl;
-                cout << Verbose(0) << " g - Relative change in kinetic energy to stop min. iteration during initial level" << endl;
-                cout << Verbose(0) << " h - Check stop conditions every ..th step during min. iteration during initial level" << endl;
-//              cout << Verbose(0) << " j - six lower diagonal entries of matrix, defining the unit cell" << endl;
-                cout << Verbose(0) << " k - Energy cutoff of plane wave basis in Hartree" << endl;
-                cout << Verbose(0) << " l - Maximum number of levels in multi-level-ansatz" << endl;
-                cout << Verbose(0) << " m - Factor by which grid nodes increase between standard and upper level" << endl;
-                cout << Verbose(0) << " n - 0: Don't use RiemannTensor, 1: Do" << endl;
-                cout << Verbose(0) << " o - Factor by which grid nodes increase between Riemann and standard(?) level" << endl;
-                cout << Verbose(0) << " p - Number of Riemann levels" << endl;
-                cout << Verbose(0) << " r - 0: Don't Use RiemannTensor, 1: Do" << endl;
-                cout << Verbose(0) << " s - 0: Doubly occupied orbitals, 1: Up-/Down-Orbitals" << endl;
-                cout << Verbose(0) << " t - Number of orbitals (depends pn SpinType)" << endl;
-                cout << Verbose(0) << " u - Number of SpinUp orbitals (depends on SpinType)" << endl;
-                cout << Verbose(0) << " v - Number of SpinDown orbitals (depends on SpinType)" << endl;
-                cout << Verbose(0) << " w - Number of additional, unoccupied orbitals" << endl;
-                cout << Verbose(0) << " x - radial cutoff for ewald summation in Bohrradii" << endl;
-                cout << Verbose(0) << " y - 0: Don't do structure optimization beforehand, 1: Do" << endl;
-                cout << Verbose(0) << " z - 0: Units are in Bohr radii, 1: units are in Aengstrom" << endl;
-                cout << Verbose(0) << " i - 0: Coordinates given in file are absolute, 1: ... are relative to unit cell" << endl;
-                cout << Verbose(0) << "=========================================================" << endl;
-                cout << Verbose(0) << "INPUT: ";
-                cin >> choice;
-
-                switch (choice) {
-                                case 'A': // mainname
-                                        cout << Verbose(0) << "Old: " << config::mainname << "\t new: ";
-                                        cin >> config::mainname;
-                                        break;
-                                case 'B': // defaultpath
-                                        cout << Verbose(0) << "Old: " << config::defaultpath << "\t new: ";
-                                        cin >> config::defaultpath;
-                                        break;
-                                case 'C': // pseudopotpath
-                                        cout << Verbose(0) << "Old: " << config::pseudopotpath << "\t new: ";
-                                        cin >> config::pseudopotpath;
-                                        break;
-
-                                case 'D': // ProcPEGamma
-                                        cout << Verbose(0) << "Old: " << config::ProcPEGamma << "\t new: ";
-                                        cin >> config::ProcPEGamma;
-                                        break;
-                                case 'E': // ProcPEPsi
-                                        cout << Verbose(0) << "Old: " << config::ProcPEPsi << "\t new: ";
-                                        cin >> config::ProcPEPsi;
-                                        break;
-                                case 'F': // DoOutVis
-                                        cout << Verbose(0) << "Old: " << config::DoOutVis << "\t new: ";
-                                        cin >> config::DoOutVis;
-                                        break;
-                                case 'G': // DoOutMes
-                                        cout << Verbose(0) << "Old: " << config::DoOutMes << "\t new: ";
-                                        cin >> config::DoOutMes;
-                                        break;
-                                case 'H': // DoOutOrbitals
-                                        cout << Verbose(0) << "Old: " << config::DoOutOrbitals << "\t new: ";
-                                        cin >> config::DoOutOrbitals;
-                                        break;
-                                case 'I': // DoOutCurrent
-                                        cout << Verbose(0) << "Old: " << config::DoOutCurrent << "\t new: ";
-                                        cin >> config::DoOutCurrent;
-                                        break;
-                                case 'J': // DoFullCurrent
-                                        cout << Verbose(0) << "Old: " << config::DoFullCurrent << "\t new: ";
-                                        cin >> config::DoFullCurrent;
-                                        break;
-                                case 'K': // DoPerturbation
-                                        cout << Verbose(0) << "Old: " << config::DoPerturbation << "\t new: ";
-                                        cin >> config::DoPerturbation;
-                                        break;
-                                case 'L': // CommonWannier
-                                        cout << Verbose(0) << "Old: " << config::CommonWannier << "\t new: ";
-                                        cin >> config::CommonWannier;
-                                        break;
-                                case 'M': // SawtoothStart
-                                        cout << Verbose(0) << "Old: " << config::SawtoothStart << "\t new: ";
-                                        cin >> config::SawtoothStart;
-                                        break;
-                                case 'N': // VectorPlane
-                                        cout << Verbose(0) << "Old: " << config::VectorPlane << "\t new: ";
-                                        cin >> config::VectorPlane;
-                                        break;
-                                case 'O': // VectorCut
-                                        cout << Verbose(0) << "Old: " << config::VectorCut << "\t new: ";
-                                        cin >> config::VectorCut;
-                                        break;
-                                case 'P': // UseAddGramSch
-                                        cout << Verbose(0) << "Old: " << config::UseAddGramSch << "\t new: ";
-                                        cin >> config::UseAddGramSch;
-                                        break;
-                                case 'Q': // Seed
-                                        cout << Verbose(0) << "Old: " << config::Seed << "\t new: ";
-                                        cin >> config::Seed;
-                                        break;
-
-                                case 'R': // MaxOuterStep
-                                        cout << Verbose(0) << "Old: " << config::MaxOuterStep << "\t new: ";
-                                        cin >> config::MaxOuterStep;
-                                        break;
-                                case 'T': // OutVisStep
-                                        cout << Verbose(0) << "Old: " << config::OutVisStep << "\t new: ";
-                                        cin >> config::OutVisStep;
-                                        break;
-                                case 'U': // OutSrcStep
-                                        cout << Verbose(0) << "Old: " << config::OutSrcStep << "\t new: ";
-                                        cin >> config::OutSrcStep;
-                                        break;
-                                case 'X': // MaxPsiStep
-                                        cout << Verbose(0) << "Old: " << config::MaxPsiStep << "\t new: ";
-                                        cin >> config::MaxPsiStep;
-                                        break;
-                                case 'Y': // EpsWannier
-                                        cout << Verbose(0) << "Old: " << config::EpsWannier << "\t new: ";
-                                        cin >> config::EpsWannier;
-                                        break;
-
-                                case 'Z': // MaxMinStep
-                                        cout << Verbose(0) << "Old: " << config::MaxMinStep << "\t new: ";
-                                        cin >> config::MaxMinStep;
-                                        break;
-                                case 'a': // RelEpsTotalEnergy
-                                        cout << Verbose(0) << "Old: " << config::RelEpsTotalEnergy << "\t new: ";
-                                        cin >> config::RelEpsTotalEnergy;
-                                        break;
-                                case 'b': // RelEpsKineticEnergy
-                                        cout << Verbose(0) << "Old: " << config::RelEpsKineticEnergy << "\t new: ";
-                                        cin >> config::RelEpsKineticEnergy;
-                                        break;
-                                case 'c': // MaxMinStopStep
-                                        cout << Verbose(0) << "Old: " << config::MaxMinStopStep << "\t new: ";
-                                        cin >> config::MaxMinStopStep;
-                                        break;
-                                case 'e': // MaxInitMinStep
-                                        cout << Verbose(0) << "Old: " << config::MaxInitMinStep << "\t new: ";
-                                        cin >> config::MaxInitMinStep;
-                                        break;
-                                case 'f': // InitRelEpsTotalEnergy
-                                        cout << Verbose(0) << "Old: " << config::InitRelEpsTotalEnergy << "\t new: ";
-                                        cin >> config::InitRelEpsTotalEnergy;
-                                        break;
-                                case 'g': // InitRelEpsKineticEnergy
-                                        cout << Verbose(0) << "Old: " << config::InitRelEpsKineticEnergy << "\t new: ";
-                                        cin >> config::InitRelEpsKineticEnergy;
-                                        break;
-                                case 'h': // InitMaxMinStopStep
-                                        cout << Verbose(0) << "Old: " << config::InitMaxMinStopStep << "\t new: ";
-                                        cin >> config::InitMaxMinStopStep;
-                                        break;
-
-//                              case 'j': // BoxLength
-//                                      cout << Verbose(0) << "enter lower triadiagonalo form of basis matrix" << endl << endl;
-//                                      for (int i=0;i<6;i++) {
-//                                              cout << Verbose(0) << "Cell size" << i << ": ";
-//                                              cin >> mol->cell_size[i];
-//                                      }
-//                                      break;
-
-                                case 'k': // ECut
-                                        cout << Verbose(0) << "Old: " << config::ECut << "\t new: ";
-                                        cin >> config::ECut;
-                                        break;
-                                case 'l': // MaxLevel
-                                        cout << Verbose(0) << "Old: " << config::MaxLevel << "\t new: ";
-                                        cin >> config::MaxLevel;
-                                        break;
-                                case 'm': // RiemannTensor
-                                        cout << Verbose(0) << "Old: " << config::RiemannTensor << "\t new: ";
-                                        cin >> config::RiemannTensor;
-                                        break;
-                                case 'n': // LevRFactor
-                                        cout << Verbose(0) << "Old: " << config::LevRFactor << "\t new: ";
-                                        cin >> config::LevRFactor;
-                                        break;
-                                case 'o': // RiemannLevel
-                                        cout << Verbose(0) << "Old: " << config::RiemannLevel << "\t new: ";
-                                        cin >> config::RiemannLevel;
-                                        break;
-                                case 'p': // Lev0Factor
-                                        cout << Verbose(0) << "Old: " << config::Lev0Factor << "\t new: ";
-                                        cin >> config::Lev0Factor;
-                                        break;
-                                case 'r': // RTActualUse
-                                        cout << Verbose(0) << "Old: " << config::RTActualUse << "\t new: ";
-                                        cin >> config::RTActualUse;
-                                        break;
-                                case 's': // PsiType
-                                        cout << Verbose(0) << "Old: " << config::PsiType << "\t new: ";
-                                        cin >> config::PsiType;
-                                        break;
-                                case 't': // MaxPsiDouble
-                                        cout << Verbose(0) << "Old: " << config::MaxPsiDouble << "\t new: ";
-                                        cin >> config::MaxPsiDouble;
-                                        break;
-                                case 'u': // PsiMaxNoUp
-                                        cout << Verbose(0) << "Old: " << config::PsiMaxNoUp << "\t new: ";
-                                        cin >> config::PsiMaxNoUp;
-                                        break;
-                                case 'v': // PsiMaxNoDown
-                                        cout << Verbose(0) << "Old: " << config::PsiMaxNoDown << "\t new: ";
-                                        cin >> config::PsiMaxNoDown;
-                                        break;
-                                case 'w': // AddPsis
-                                        cout << Verbose(0) << "Old: " << config::AddPsis << "\t new: ";
-                                        cin >> config::AddPsis;
-                                        break;
-
-                                case 'x': // RCut
-                                        cout << Verbose(0) << "Old: " << config::RCut << "\t new: ";
-                                        cin >> config::RCut;
-                                        break;
-                                case 'y': // StructOpt
-                                        cout << Verbose(0) << "Old: " << config::StructOpt << "\t new: ";
-                                        cin >> config::StructOpt;
-                                        break;
-                                case 'z': // IsAngstroem
-                                        cout << Verbose(0) << "Old: " << config::IsAngstroem << "\t new: ";
-                                        cin >> config::IsAngstroem;
-                                        break;
-                                case 'i': // RelativeCoord
-                                        cout << Verbose(0) << "Old: " << config::RelativeCoord << "\t new: ";
-                                        cin >> config::RelativeCoord;
-                                        break;
-                };
-        } while (choice != 'q');
+  char choice;
+
+  do {
+    cout << Verbose(0) << "===========EDIT CONFIGURATION============================" << endl;
+    cout << Verbose(0) << " A - mainname (prefix for all runtime files)" << endl;
+    cout << Verbose(0) << " B - Default path (for runtime files)" << endl;
+    cout << Verbose(0) << " C - Path of pseudopotential files" << endl;
+    cout << Verbose(0) << " D - Number of coefficient sharing processes" << endl;
+    cout << Verbose(0) << " E - Number of wave function sharing processes" << endl;
+    cout << Verbose(0) << " F - 0: Don't output density for OpenDX, 1: do" << endl;
+    cout << Verbose(0) << " G - 0: Don't output physical data, 1: do" << endl;
+    cout << Verbose(0) << " H - 0: Don't output densities of each unperturbed orbital for OpenDX, 1: do" << endl;
+    cout << Verbose(0) << " I - 0: Don't output current density for OpenDX, 1: do" << endl;
+    cout << Verbose(0) << " J - 0: Don't do the full current calculation, 1: do" << endl;
+    cout << Verbose(0) << " K - 0: Don't do perturbation calculation to obtain susceptibility and shielding, 1: do" << endl;
+    cout << Verbose(0) << " L - 0: Wannier centres as calculated, 1: common centre for all, 2: unite centres according to spread, 3: cell centre, 4: shifted to nearest grid point" << endl;
+    cout << Verbose(0) << " M - Absolute begin of unphysical sawtooth transfer for position operator within cell" << endl;
+    cout << Verbose(0) << " N - (0,1,2) x,y,z-plane to do two-dimensional current vector cut" << endl;
+    cout << Verbose(0) << " O - Absolute position along vector cut axis for cut plane" << endl;
+    cout << Verbose(0) << " P - Additional Gram-Schmidt-Orthonormalization to stabilize numerics" << endl;
+    cout << Verbose(0) << " Q - Initial integer value of random number generator" << endl;
+    cout << Verbose(0) << " R - for perturbation 0, for structure optimization defines upper limit of iterations" << endl;
+    cout << Verbose(0) << " T - Output visual after ...th step" << endl;
+    cout << Verbose(0) << " U - Output source densities of wave functions after ...th step" << endl;
+    cout << Verbose(0) << " X - minimization iterations per wave function, if unsure leave at default value 0" << endl;
+    cout << Verbose(0) << " Y - tolerance value for total spread in iterative Jacobi diagonalization" << endl;
+    cout << Verbose(0) << " Z - Maximum number of minimization iterations" << endl;
+    cout << Verbose(0) << " a - Relative change in total energy to stop min. iteration" << endl;
+    cout << Verbose(0) << " b - Relative change in kinetic energy to stop min. iteration" << endl;
+    cout << Verbose(0) << " c - Check stop conditions every ..th step during min. iteration" << endl;
+    cout << Verbose(0) << " e - Maximum number of minimization iterations during initial level" << endl;
+    cout << Verbose(0) << " f - Relative change in total energy to stop min. iteration during initial level" << endl;
+    cout << Verbose(0) << " g - Relative change in kinetic energy to stop min. iteration during initial level" << endl;
+    cout << Verbose(0) << " h - Check stop conditions every ..th step during min. iteration during initial level" << endl;
+//    cout << Verbose(0) << " j - six lower diagonal entries of matrix, defining the unit cell" << endl;
+    cout << Verbose(0) << " k - Energy cutoff of plane wave basis in Hartree" << endl;
+    cout << Verbose(0) << " l - Maximum number of levels in multi-level-ansatz" << endl;
+    cout << Verbose(0) << " m - Factor by which grid nodes increase between standard and upper level" << endl;
+    cout << Verbose(0) << " n - 0: Don't use RiemannTensor, 1: Do" << endl;
+    cout << Verbose(0) << " o - Factor by which grid nodes increase between Riemann and standard(?) level" << endl;
+    cout << Verbose(0) << " p - Number of Riemann levels" << endl;
+    cout << Verbose(0) << " r - 0: Don't Use RiemannTensor, 1: Do" << endl;
+    cout << Verbose(0) << " s - 0: Doubly occupied orbitals, 1: Up-/Down-Orbitals" << endl;
+    cout << Verbose(0) << " t - Number of orbitals (depends pn SpinType)" << endl;
+    cout << Verbose(0) << " u - Number of SpinUp orbitals (depends on SpinType)" << endl;
+    cout << Verbose(0) << " v - Number of SpinDown orbitals (depends on SpinType)" << endl;
+    cout << Verbose(0) << " w - Number of additional, unoccupied orbitals" << endl;
+    cout << Verbose(0) << " x - radial cutoff for ewald summation in Bohrradii" << endl;
+    cout << Verbose(0) << " y - 0: Don't do structure optimization beforehand, 1: Do" << endl;
+    cout << Verbose(0) << " z - 0: Units are in Bohr radii, 1: units are in Aengstrom" << endl;
+    cout << Verbose(0) << " i - 0: Coordinates given in file are absolute, 1: ... are relative to unit cell" << endl;
+    cout << Verbose(0) << "=========================================================" << endl;
+    cout << Verbose(0) << "INPUT: ";
+    cin >> choice;
+
+    switch (choice) {
+        case 'A': // mainname
+          cout << Verbose(0) << "Old: " << config::mainname << "\t new: ";
+          cin >> config::mainname;
+          break;
+        case 'B': // defaultpath
+          cout << Verbose(0) << "Old: " << config::defaultpath << "\t new: ";
+          cin >> config::defaultpath;
+          break;
+        case 'C': // pseudopotpath
+          cout << Verbose(0) << "Old: " << config::pseudopotpath << "\t new: ";
+          cin >> config::pseudopotpath;
+          break;
+
+        case 'D': // ProcPEGamma
+          cout << Verbose(0) << "Old: " << config::ProcPEGamma << "\t new: ";
+          cin >> config::ProcPEGamma;
+          break;
+        case 'E': // ProcPEPsi
+          cout << Verbose(0) << "Old: " << config::ProcPEPsi << "\t new: ";
+          cin >> config::ProcPEPsi;
+          break;
+        case 'F': // DoOutVis
+          cout << Verbose(0) << "Old: " << config::DoOutVis << "\t new: ";
+          cin >> config::DoOutVis;
+          break;
+        case 'G': // DoOutMes
+          cout << Verbose(0) << "Old: " << config::DoOutMes << "\t new: ";
+          cin >> config::DoOutMes;
+          break;
+        case 'H': // DoOutOrbitals
+          cout << Verbose(0) << "Old: " << config::DoOutOrbitals << "\t new: ";
+          cin >> config::DoOutOrbitals;
+          break;
+        case 'I': // DoOutCurrent
+          cout << Verbose(0) << "Old: " << config::DoOutCurrent << "\t new: ";
+          cin >> config::DoOutCurrent;
+          break;
+        case 'J': // DoFullCurrent
+          cout << Verbose(0) << "Old: " << config::DoFullCurrent << "\t new: ";
+          cin >> config::DoFullCurrent;
+          break;
+        case 'K': // DoPerturbation
+          cout << Verbose(0) << "Old: " << config::DoPerturbation << "\t new: ";
+          cin >> config::DoPerturbation;
+          break;
+        case 'L': // CommonWannier
+          cout << Verbose(0) << "Old: " << config::CommonWannier << "\t new: ";
+          cin >> config::CommonWannier;
+          break;
+        case 'M': // SawtoothStart
+          cout << Verbose(0) << "Old: " << config::SawtoothStart << "\t new: ";
+          cin >> config::SawtoothStart;
+          break;
+        case 'N': // VectorPlane
+          cout << Verbose(0) << "Old: " << config::VectorPlane << "\t new: ";
+          cin >> config::VectorPlane;
+          break;
+        case 'O': // VectorCut
+          cout << Verbose(0) << "Old: " << config::VectorCut << "\t new: ";
+          cin >> config::VectorCut;
+          break;
+        case 'P': // UseAddGramSch
+          cout << Verbose(0) << "Old: " << config::UseAddGramSch << "\t new: ";
+          cin >> config::UseAddGramSch;
+          break;
+        case 'Q': // Seed
+          cout << Verbose(0) << "Old: " << config::Seed << "\t new: ";
+          cin >> config::Seed;
+          break;
+
+        case 'R': // MaxOuterStep
+          cout << Verbose(0) << "Old: " << config::MaxOuterStep << "\t new: ";
+          cin >> config::MaxOuterStep;
+          break;
+        case 'T': // OutVisStep
+          cout << Verbose(0) << "Old: " << config::OutVisStep << "\t new: ";
+          cin >> config::OutVisStep;
+          break;
+        case 'U': // OutSrcStep
+          cout << Verbose(0) << "Old: " << config::OutSrcStep << "\t new: ";
+          cin >> config::OutSrcStep;
+          break;
+        case 'X': // MaxPsiStep
+          cout << Verbose(0) << "Old: " << config::MaxPsiStep << "\t new: ";
+          cin >> config::MaxPsiStep;
+          break;
+        case 'Y': // EpsWannier
+          cout << Verbose(0) << "Old: " << config::EpsWannier << "\t new: ";
+          cin >> config::EpsWannier;
+          break;
+
+        case 'Z': // MaxMinStep
+          cout << Verbose(0) << "Old: " << config::MaxMinStep << "\t new: ";
+          cin >> config::MaxMinStep;
+          break;
+        case 'a': // RelEpsTotalEnergy
+          cout << Verbose(0) << "Old: " << config::RelEpsTotalEnergy << "\t new: ";
+          cin >> config::RelEpsTotalEnergy;
+          break;
+        case 'b': // RelEpsKineticEnergy
+          cout << Verbose(0) << "Old: " << config::RelEpsKineticEnergy << "\t new: ";
+          cin >> config::RelEpsKineticEnergy;
+          break;
+        case 'c': // MaxMinStopStep
+          cout << Verbose(0) << "Old: " << config::MaxMinStopStep << "\t new: ";
+          cin >> config::MaxMinStopStep;
+          break;
+        case 'e': // MaxInitMinStep
+          cout << Verbose(0) << "Old: " << config::MaxInitMinStep << "\t new: ";
+          cin >> config::MaxInitMinStep;
+          break;
+        case 'f': // InitRelEpsTotalEnergy
+          cout << Verbose(0) << "Old: " << config::InitRelEpsTotalEnergy << "\t new: ";
+          cin >> config::InitRelEpsTotalEnergy;
+          break;
+        case 'g': // InitRelEpsKineticEnergy
+          cout << Verbose(0) << "Old: " << config::InitRelEpsKineticEnergy << "\t new: ";
+          cin >> config::InitRelEpsKineticEnergy;
+          break;
+        case 'h': // InitMaxMinStopStep
+          cout << Verbose(0) << "Old: " << config::InitMaxMinStopStep << "\t new: ";
+          cin >> config::InitMaxMinStopStep;
+          break;
+
+//        case 'j': // BoxLength
+//          cout << Verbose(0) << "enter lower triadiagonalo form of basis matrix" << endl << endl;
+//          for (int i=0;i<6;i++) {
+//            cout << Verbose(0) << "Cell size" << i << ": ";
+//            cin >> mol->cell_size[i];
+//          }
+//          break;
+
+        case 'k': // ECut
+          cout << Verbose(0) << "Old: " << config::ECut << "\t new: ";
+          cin >> config::ECut;
+          break;
+        case 'l': // MaxLevel
+          cout << Verbose(0) << "Old: " << config::MaxLevel << "\t new: ";
+          cin >> config::MaxLevel;
+          break;
+        case 'm': // RiemannTensor
+          cout << Verbose(0) << "Old: " << config::RiemannTensor << "\t new: ";
+          cin >> config::RiemannTensor;
+          break;
+        case 'n': // LevRFactor
+          cout << Verbose(0) << "Old: " << config::LevRFactor << "\t new: ";
+          cin >> config::LevRFactor;
+          break;
+        case 'o': // RiemannLevel
+          cout << Verbose(0) << "Old: " << config::RiemannLevel << "\t new: ";
+          cin >> config::RiemannLevel;
+          break;
+        case 'p': // Lev0Factor
+          cout << Verbose(0) << "Old: " << config::Lev0Factor << "\t new: ";
+          cin >> config::Lev0Factor;
+          break;
+        case 'r': // RTActualUse
+          cout << Verbose(0) << "Old: " << config::RTActualUse << "\t new: ";
+          cin >> config::RTActualUse;
+          break;
+        case 's': // PsiType
+          cout << Verbose(0) << "Old: " << config::PsiType << "\t new: ";
+          cin >> config::PsiType;
+          break;
+        case 't': // MaxPsiDouble
+          cout << Verbose(0) << "Old: " << config::MaxPsiDouble << "\t new: ";
+          cin >> config::MaxPsiDouble;
+          break;
+        case 'u': // PsiMaxNoUp
+          cout << Verbose(0) << "Old: " << config::PsiMaxNoUp << "\t new: ";
+          cin >> config::PsiMaxNoUp;
+          break;
+        case 'v': // PsiMaxNoDown
+          cout << Verbose(0) << "Old: " << config::PsiMaxNoDown << "\t new: ";
+          cin >> config::PsiMaxNoDown;
+          break;
+        case 'w': // AddPsis
+          cout << Verbose(0) << "Old: " << config::AddPsis << "\t new: ";
+          cin >> config::AddPsis;
+          break;
+
+        case 'x': // RCut
+          cout << Verbose(0) << "Old: " << config::RCut << "\t new: ";
+          cin >> config::RCut;
+          break;
+        case 'y': // StructOpt
+          cout << Verbose(0) << "Old: " << config::StructOpt << "\t new: ";
+          cin >> config::StructOpt;
+          break;
+        case 'z': // IsAngstroem
+          cout << Verbose(0) << "Old: " << config::IsAngstroem << "\t new: ";
+          cin >> config::IsAngstroem;
+          break;
+        case 'i': // RelativeCoord
+          cout << Verbose(0) << "Old: " << config::RelativeCoord << "\t new: ";
+          cin >> config::RelativeCoord;
+          break;
+    };
+  } while (choice != 'q');
 };
 
@@ -520 +520 @@
 int config::TestSyntax(char *filename, periodentafel *periode, molecule *mol)
 {
-        int test;
-        ifstream file(filename);
-
-        // search file for keyword: ProcPEGamma (new syntax)
-        if (ParseForParameter(1,&file,"ProcPEGamma", 0, 1, 1, int_type, &test, 1, optional)) {
-                file.close();
-                return 1;
-        }
-        // search file for keyword: ProcsGammaPsi (old syntax)
-        if (ParseForParameter(1,&file,"ProcsGammaPsi", 0, 1, 1, int_type, &test, 1, optional)) {
-                file.close();
-                return 0;
-        }
-        file.close();
-        return -1;
+  int test;
+  ifstream file(filename);
+
+  // search file for keyword: ProcPEGamma (new syntax)
+  if (ParseForParameter(1,&file,"ProcPEGamma", 0, 1, 1, int_type, &test, 1, optional)) {
+    file.close();
+    return 1;
+  }
+  // search file for keyword: ProcsGammaPsi (old syntax)
+  if (ParseForParameter(1,&file,"ProcsGammaPsi", 0, 1, 1, int_type, &test, 1, optional)) {
+    file.close();
+    return 0;
+  }
+  file.close();
+  return -1;
 }
 
@@ -542 +542 @@
 bool config::GetIsAngstroem() const
 {
-        return (IsAngstroem == 1);
+  return (IsAngstroem == 1);
 };
 
@@ -550 +550 @@
 char * config::GetDefaultPath() const
 {
-        return defaultpath;
+  return defaultpath;
 };
 
@@ -559 +559 @@
 void config::SetDefaultPath(const char *path)
 {
-        strcpy(defaultpath, path);
+  strcpy(defaultpath, path);
 };
 
@@ -567 +567 @@
 void config::RetrieveConfigPathAndName(string filename)
 {
-        char *ptr = NULL;
-        char *buffer = new char[MAXSTRINGSIZE];
-        strncpy(buffer, filename.c_str(), MAXSTRINGSIZE);
-        int last = -1;
-        for(last=MAXSTRINGSIZE;last--;) {
-                if (buffer[last] == '/')
-                        break;
-        }
-        if (last == -1) { // no path in front, set to local directory.
-                strcpy(configpath, "./");
-                ptr = buffer;
-        } else {
-                strncpy(configpath, buffer, last+1);
-                ptr = &buffer[last+1];
-                if (last < 254)
-                        configpath[last+1]='\0';
-        }
-        strcpy(configname, ptr);
-        cout << "Found configpath: " << configpath << ", dir slash was found at " << last << ", config name is " << configname << "." << endl;
-        delete[](buffer);
+  char *ptr = NULL;
+  char *buffer = new char[MAXSTRINGSIZE];
+  strncpy(buffer, filename.c_str(), MAXSTRINGSIZE);
+  int last = -1;
+  for(last=MAXSTRINGSIZE;last--;) {
+    if (buffer[last] == '/')
+      break;
+  }
+  if (last == -1) { // no path in front, set to local directory.
+    strcpy(configpath, "./");
+    ptr = buffer;
+  } else {
+    strncpy(configpath, buffer, last+1);
+    ptr = &buffer[last+1];
+    if (last < 254)
+      configpath[last+1]='\0';
+  }
+  strcpy(configname, ptr);
+  cout << "Found configpath: " << configpath << ", dir slash was found at " << last << ", config name is " << configname << "." << endl;
+  delete[](buffer);
 };
 
@@ -597 +597 @@
 void config::Load(char *filename, periodentafel *periode, molecule *mol)
 {
-        RetrieveConfigPathAndName(filename);
-
-        // ParseParameterFile
+  RetrieveConfigPathAndName(filename);
+
+  // ParseParameterFile
   struct ConfigFileBuffer *FileBuffer = new ConfigFileBuffer(filename);
-        FileBuffer->InitMapping();
-
-        /* Oeffne Hauptparameterdatei */
-        int di;
-        double BoxLength[9];
-        string zeile;
-        string dummy;
-        element *elementhash[MAX_ELEMENTS];
-        char name[MAX_ELEMENTS];
-        char keyword[MAX_ELEMENTS];
-        int Z, No[MAX_ELEMENTS];
-        int verbose = 0;
-        double value[3];
-
-        /* Namen einlesen */
-
-        ParseForParameter(verbose,FileBuffer, "mainname", 0, 1, 1, string_type, (config::mainname), 1, critical);
-        ParseForParameter(verbose,FileBuffer, "defaultpath", 0, 1, 1, string_type, (config::defaultpath), 1, critical);
-        ParseForParameter(verbose,FileBuffer, "pseudopotpath", 0, 1, 1, string_type, (config::pseudopotpath), 1, critical);
-        ParseForParameter(verbose,FileBuffer,"ProcPEGamma", 0, 1, 1, int_type, &(config::ProcPEGamma), 1, critical);
-        ParseForParameter(verbose,FileBuffer,"ProcPEPsi", 0, 1, 1, int_type, &(config::ProcPEPsi), 1, critical);
-
-        if (!ParseForParameter(verbose,FileBuffer,"Seed", 0, 1, 1, int_type, &(config::Seed), 1, optional))
-                config::Seed = 1;
-
-        if(!ParseForParameter(verbose,FileBuffer,"DoOutOrbitals", 0, 1, 1, int_type, &(config::DoOutOrbitals), 1, optional)) {
-                config::DoOutOrbitals = 0;
-        } else {
-                if (config::DoOutOrbitals < 0) config::DoOutOrbitals = 0;
-                if (config::DoOutOrbitals > 1) config::DoOutOrbitals = 1;
-        }
-        ParseForParameter(verbose,FileBuffer,"DoOutVis", 0, 1, 1, int_type, &(config::DoOutVis), 1, critical);
-        if (config::DoOutVis < 0) config::DoOutVis = 0;
-        if (config::DoOutVis > 1) config::DoOutVis = 1;
-        if (!ParseForParameter(verbose,FileBuffer,"VectorPlane", 0, 1, 1, int_type, &(config::VectorPlane), 1, optional))
-                config::VectorPlane = -1;
-        if (!ParseForParameter(verbose,FileBuffer,"VectorCut", 0, 1, 1, double_type, &(config::VectorCut), 1, optional))
-                config::VectorCut = 0.;
-        ParseForParameter(verbose,FileBuffer,"DoOutMes", 0, 1, 1, int_type, &(config::DoOutMes), 1, critical);
-        if (config::DoOutMes < 0) config::DoOutMes = 0;
-        if (config::DoOutMes > 1) config::DoOutMes = 1;
-        if (!ParseForParameter(verbose,FileBuffer,"DoOutCurr", 0, 1, 1, int_type, &(config::DoOutCurrent), 1, optional))
-                config::DoOutCurrent = 0;
-        if (config::DoOutCurrent < 0) config::DoOutCurrent = 0;
-        if (config::DoOutCurrent > 1) config::DoOutCurrent = 1;
-        ParseForParameter(verbose,FileBuffer,"AddGramSch", 0, 1, 1, int_type, &(config::UseAddGramSch), 1, critical);
-        if (config::UseAddGramSch < 0) config::UseAddGramSch = 0;
-        if (config::UseAddGramSch > 2) config::UseAddGramSch = 2;
-        if(!ParseForParameter(verbose,FileBuffer,"DoWannier", 0, 1, 1, int_type, &(config::DoWannier), 1, optional)) {
-                config::DoWannier = 0;
-        } else {
-                if (config::DoWannier < 0) config::DoWannier = 0;
-                if (config::DoWannier > 1) config::DoWannier = 1;
-        }
-        if(!ParseForParameter(verbose,FileBuffer,"CommonWannier", 0, 1, 1, int_type, &(config::CommonWannier), 1, optional)) {
-                config::CommonWannier = 0;
-        } else {
-                if (config::CommonWannier < 0) config::CommonWannier = 0;
-                if (config::CommonWannier > 4) config::CommonWannier = 4;
-        }
-        if(!ParseForParameter(verbose,FileBuffer,"SawtoothStart", 0, 1, 1, double_type, &(config::SawtoothStart), 1, optional)) {
-                config::SawtoothStart = 0.01;
-        } else {
-                if (config::SawtoothStart < 0.) config::SawtoothStart = 0.;
-                if (config::SawtoothStart > 1.) config::SawtoothStart = 1.;
-        }
-
-        ParseForParameter(verbose,FileBuffer,"MaxOuterStep", 0, 1, 1, int_type, &(config::MaxOuterStep), 1, critical);
-        if (!ParseForParameter(verbose,FileBuffer,"Deltat", 0, 1, 1, double_type, &(config::Deltat), 1, optional))
-                config::Deltat = 1;
-        ParseForParameter(verbose,FileBuffer,"OutVisStep", 0, 1, 1, int_type, &(config::OutVisStep), 1, optional);
-        ParseForParameter(verbose,FileBuffer,"OutSrcStep", 0, 1, 1, int_type, &(config::OutSrcStep), 1, optional);
-        ParseForParameter(verbose,FileBuffer,"TargetTemp", 0, 1, 1, double_type, &(config::TargetTemp), 1, optional);
-        //ParseForParameter(verbose,FileBuffer,"Thermostat", 0, 1, 1, int_type, &(config::ScaleTempStep), 1, optional);
-        if (!ParseForParameter(verbose,FileBuffer,"EpsWannier", 0, 1, 1, double_type, &(config::EpsWannier), 1, optional))
-                config::EpsWannier = 1e-8;
-
-        // stop conditions
-        //if (config::MaxOuterStep <= 0) config::MaxOuterStep = 1;
-        ParseForParameter(verbose,FileBuffer,"MaxPsiStep", 0, 1, 1, int_type, &(config::MaxPsiStep), 1, critical);
-        if (config::MaxPsiStep <= 0) config::MaxPsiStep = 3;
-
-        ParseForParameter(verbose,FileBuffer,"MaxMinStep", 0, 1, 1, int_type, &(config::MaxMinStep), 1, critical);
-        ParseForParameter(verbose,FileBuffer,"RelEpsTotalE", 0, 1, 1, double_type, &(config::RelEpsTotalEnergy), 1, critical);
-        ParseForParameter(verbose,FileBuffer,"RelEpsKineticE", 0, 1, 1, double_type, &(config::RelEpsKineticEnergy), 1, critical);
-        ParseForParameter(verbose,FileBuffer,"MaxMinStopStep", 0, 1, 1, int_type, &(config::MaxMinStopStep), 1, critical);
-        ParseForParameter(verbose,FileBuffer,"MaxMinGapStopStep", 0, 1, 1, int_type, &(config::MaxMinGapStopStep), 1, critical);
-        if (config::MaxMinStep <= 0) config::MaxMinStep = config::MaxPsiStep;
-        if (config::MaxMinStopStep < 1) config::MaxMinStopStep = 1;
-        if (config::MaxMinGapStopStep < 1) config::MaxMinGapStopStep = 1;
-
-        ParseForParameter(verbose,FileBuffer,"MaxInitMinStep", 0, 1, 1, int_type, &(config::MaxInitMinStep), 1, critical);
-        ParseForParameter(verbose,FileBuffer,"InitRelEpsTotalE", 0, 1, 1, double_type, &(config::InitRelEpsTotalEnergy), 1, critical);
-        ParseForParameter(verbose,FileBuffer,"InitRelEpsKineticE", 0, 1, 1, double_type, &(config::InitRelEpsKineticEnergy), 1, critical);
-        ParseForParameter(verbose,FileBuffer,"InitMaxMinStopStep", 0, 1, 1, int_type, &(config::InitMaxMinStopStep), 1, critical);
-        ParseForParameter(verbose,FileBuffer,"InitMaxMinGapStopStep", 0, 1, 1, int_type, &(config::InitMaxMinGapStopStep), 1, critical);
-        if (config::MaxInitMinStep <= 0) config::MaxInitMinStep = config::MaxPsiStep;
-        if (config::InitMaxMinStopStep < 1) config::InitMaxMinStopStep = 1;
-        if (config::InitMaxMinGapStopStep < 1) config::InitMaxMinGapStopStep = 1;
-
-        // Unit cell and magnetic field
-        ParseForParameter(verbose,FileBuffer, "BoxLength", 0, 3, 3, lower_trigrid, BoxLength, 1, critical); /* Lattice->RealBasis */
-        mol->cell_size[0] = BoxLength[0];
-        mol->cell_size[1] = BoxLength[3];
-        mol->cell_size[2] = BoxLength[4];
-        mol->cell_size[3] = BoxLength[6];
-        mol->cell_size[4] = BoxLength[7];
-        mol->cell_size[5] = BoxLength[8];
-        //if (1) fprintf(stderr,"\n");
-
-        ParseForParameter(verbose,FileBuffer,"DoPerturbation", 0, 1, 1, int_type, &(config::DoPerturbation), 1, optional);
-        ParseForParameter(verbose,FileBuffer,"DoOutNICS", 0, 1, 1, int_type, &(config::DoOutNICS), 1, optional);
-        if (!ParseForParameter(verbose,FileBuffer,"DoFullCurrent", 0, 1, 1, int_type, &(config::DoFullCurrent), 1, optional))
-                config::DoFullCurrent = 0;
-        if (config::DoFullCurrent < 0) config::DoFullCurrent = 0;
-        if (config::DoFullCurrent > 2) config::DoFullCurrent = 2;
-        if (config::DoOutNICS < 0) config::DoOutNICS = 0;
-        if (config::DoOutNICS > 2) config::DoOutNICS = 2;
-        if (config::DoPerturbation == 0) {
-                config::DoFullCurrent = 0;
-                config::DoOutNICS = 0;
-        }
-
-        ParseForParameter(verbose,FileBuffer,"ECut", 0, 1, 1, double_type, &(config::ECut), 1, critical);
-        ParseForParameter(verbose,FileBuffer,"MaxLevel", 0, 1, 1, int_type, &(config::MaxLevel), 1, critical);
-        ParseForParameter(verbose,FileBuffer,"Level0Factor", 0, 1, 1, int_type, &(config::Lev0Factor), 1, critical);
-        if (config::Lev0Factor < 2) {
-                config::Lev0Factor = 2;
-        }
-        ParseForParameter(verbose,FileBuffer,"RiemannTensor", 0, 1, 1, int_type, &di, 1, critical);
-        if (di >= 0 && di < 2) {
-                config::RiemannTensor = di;
-        } else {
-                fprintf(stderr, "0 <= RiemanTensor < 2: 0 UseNotRT, 1 UseRT");
-                exit(1);
-        }
-        switch (config::RiemannTensor) {
-                case 0: //UseNoRT
-                        if (config::MaxLevel < 2) {
-                                config::MaxLevel = 2;
-                        }
-                        config::LevRFactor = 2;
-                        config::RTActualUse = 0;
-                        break;
-                case 1: // UseRT
-                        if (config::MaxLevel < 3) {
-                                config::MaxLevel = 3;
-                        }
-                        ParseForParameter(verbose,FileBuffer,"RiemannLevel", 0, 1, 1, int_type, &(config::RiemannLevel), 1, critical);
-                        if (config::RiemannLevel < 2) {
-                                config::RiemannLevel = 2;
-                        }
-                        if (config::RiemannLevel > config::MaxLevel-1) {
-                                config::RiemannLevel = config::MaxLevel-1;
-                        }
-                        ParseForParameter(verbose,FileBuffer,"LevRFactor", 0, 1, 1, int_type, &(config::LevRFactor), 1, critical);
-                        if (config::LevRFactor < 2) {
-                                config::LevRFactor = 2;
-                        }
-                        config::Lev0Factor = 2;
-                        config::RTActualUse = 2;
-                        break;
-        }
-        ParseForParameter(verbose,FileBuffer,"PsiType", 0, 1, 1, int_type, &di, 1, critical);
-        if (di >= 0 && di < 2) {
-                config::PsiType = di;
-        } else {
-                fprintf(stderr, "0 <= PsiType < 2: 0 UseSpinDouble, 1 UseSpinUpDown");
-                exit(1);
-        }
-        switch (config::PsiType) {
-        case 0: // SpinDouble
-                ParseForParameter(verbose,FileBuffer,"MaxPsiDouble", 0, 1, 1, int_type, &(config::MaxPsiDouble), 1, critical);
-                ParseForParameter(verbose,FileBuffer,"AddPsis", 0, 1, 1, int_type, &(config::AddPsis), 1, optional);
-                break;
-        case 1: // SpinUpDown
-                if (config::ProcPEGamma % 2) config::ProcPEGamma*=2;
-                ParseForParameter(verbose,FileBuffer,"PsiMaxNoUp", 0, 1, 1, int_type, &(config::PsiMaxNoUp), 1, critical);
-                ParseForParameter(verbose,FileBuffer,"PsiMaxNoDown", 0, 1, 1, int_type, &(config::PsiMaxNoDown), 1, critical);
-                ParseForParameter(verbose,FileBuffer,"AddPsis", 0, 1, 1, int_type, &(config::AddPsis), 1, optional);
-                break;
-        }
-
-        // IonsInitRead
-
-        ParseForParameter(verbose,FileBuffer,"RCut", 0, 1, 1, double_type, &(config::RCut), 1, critical);
-        ParseForParameter(verbose,FileBuffer,"IsAngstroem", 0, 1, 1, int_type, &(config::IsAngstroem), 1, critical);
-        ParseForParameter(verbose,FileBuffer,"MaxTypes", 0, 1, 1, int_type, &(config::MaxTypes), 1, critical);
-        if (!ParseForParameter(verbose,FileBuffer,"RelativeCoord", 0, 1, 1, int_type, &(config::RelativeCoord) , 1, optional))
-                config::RelativeCoord = 0;
-        if (!ParseForParameter(verbose,FileBuffer,"StructOpt", 0, 1, 1, int_type, &(config::StructOpt), 1, optional))
-                config::StructOpt = 0;
-        if (MaxTypes == 0) {
-                cerr << "There are no atoms according to MaxTypes in this config file." << endl;
-        } else {
-                // prescan number of ions per type
-                cout << Verbose(0) << "Prescanning ions per type: " << endl;
-                int NoAtoms = 0;
-                for (int i=0; i < config::MaxTypes; i++) {
-                        sprintf(name,"Ion_Type%i",i+1);
-                        ParseForParameter(verbose,FileBuffer, (const char*)name, 0, 1, 1, int_type, &No[i], 1, critical);
-                        ParseForParameter(verbose,FileBuffer, name, 0, 2, 1, int_type, &Z, 1, critical);
-                        elementhash[i] = periode->FindElement(Z);
-                        cout << Verbose(1) << i << ". Z = " << elementhash[i]->Z << " with " << No[i] << " ions." << endl;
-                        NoAtoms += No[i];
-                }
-                int repetition = 0; // which repeated keyword shall be read
-
-                // sort the lines via the LineMapping
+  FileBuffer->InitMapping();
+
+  /* Oeffne Hauptparameterdatei */
+  int di;
+  double BoxLength[9];
+  string zeile;
+  string dummy;
+  element *elementhash[MAX_ELEMENTS];
+  char name[MAX_ELEMENTS];
+  char keyword[MAX_ELEMENTS];
+  int Z, No[MAX_ELEMENTS];
+  int verbose = 0;
+  double value[3];
+
+  /* Namen einlesen */
+
+  ParseForParameter(verbose,FileBuffer, "mainname", 0, 1, 1, string_type, (config::mainname), 1, critical);
+  ParseForParameter(verbose,FileBuffer, "defaultpath", 0, 1, 1, string_type, (config::defaultpath), 1, critical);
+  ParseForParameter(verbose,FileBuffer, "pseudopotpath", 0, 1, 1, string_type, (config::pseudopotpath), 1, critical);
+  ParseForParameter(verbose,FileBuffer,"ProcPEGamma", 0, 1, 1, int_type, &(config::ProcPEGamma), 1, critical);
+  ParseForParameter(verbose,FileBuffer,"ProcPEPsi", 0, 1, 1, int_type, &(config::ProcPEPsi), 1, critical);
+
+  if (!ParseForParameter(verbose,FileBuffer,"Seed", 0, 1, 1, int_type, &(config::Seed), 1, optional))
+    config::Seed = 1;
+
+  if(!ParseForParameter(verbose,FileBuffer,"DoOutOrbitals", 0, 1, 1, int_type, &(config::DoOutOrbitals), 1, optional)) {
+    config::DoOutOrbitals = 0;
+  } else {
+    if (config::DoOutOrbitals < 0) config::DoOutOrbitals = 0;
+    if (config::DoOutOrbitals > 1) config::DoOutOrbitals = 1;
+  }
+  ParseForParameter(verbose,FileBuffer,"DoOutVis", 0, 1, 1, int_type, &(config::DoOutVis), 1, critical);
+  if (config::DoOutVis < 0) config::DoOutVis = 0;
+  if (config::DoOutVis > 1) config::DoOutVis = 1;
+  if (!ParseForParameter(verbose,FileBuffer,"VectorPlane", 0, 1, 1, int_type, &(config::VectorPlane), 1, optional))
+    config::VectorPlane = -1;
+  if (!ParseForParameter(verbose,FileBuffer,"VectorCut", 0, 1, 1, double_type, &(config::VectorCut), 1, optional))
+    config::VectorCut = 0.;
+  ParseForParameter(verbose,FileBuffer,"DoOutMes", 0, 1, 1, int_type, &(config::DoOutMes), 1, critical);
+  if (config::DoOutMes < 0) config::DoOutMes = 0;
+  if (config::DoOutMes > 1) config::DoOutMes = 1;
+  if (!ParseForParameter(verbose,FileBuffer,"DoOutCurr", 0, 1, 1, int_type, &(config::DoOutCurrent), 1, optional))
+    config::DoOutCurrent = 0;
+  if (config::DoOutCurrent < 0) config::DoOutCurrent = 0;
+  if (config::DoOutCurrent > 1) config::DoOutCurrent = 1;
+  ParseForParameter(verbose,FileBuffer,"AddGramSch", 0, 1, 1, int_type, &(config::UseAddGramSch), 1, critical);
+  if (config::UseAddGramSch < 0) config::UseAddGramSch = 0;
+  if (config::UseAddGramSch > 2) config::UseAddGramSch = 2;
+  if(!ParseForParameter(verbose,FileBuffer,"DoWannier", 0, 1, 1, int_type, &(config::DoWannier), 1, optional)) {
+    config::DoWannier = 0;
+  } else {
+    if (config::DoWannier < 0) config::DoWannier = 0;
+    if (config::DoWannier > 1) config::DoWannier = 1;
+  }
+  if(!ParseForParameter(verbose,FileBuffer,"CommonWannier", 0, 1, 1, int_type, &(config::CommonWannier), 1, optional)) {
+    config::CommonWannier = 0;
+  } else {
+    if (config::CommonWannier < 0) config::CommonWannier = 0;
+    if (config::CommonWannier > 4) config::CommonWannier = 4;
+  }
+  if(!ParseForParameter(verbose,FileBuffer,"SawtoothStart", 0, 1, 1, double_type, &(config::SawtoothStart), 1, optional)) {
+    config::SawtoothStart = 0.01;
+  } else {
+    if (config::SawtoothStart < 0.) config::SawtoothStart = 0.;
+    if (config::SawtoothStart > 1.) config::SawtoothStart = 1.;
+  }
+
+  ParseForParameter(verbose,FileBuffer,"MaxOuterStep", 0, 1, 1, int_type, &(config::MaxOuterStep), 1, critical);
+  if (!ParseForParameter(verbose,FileBuffer,"Deltat", 0, 1, 1, double_type, &(config::Deltat), 1, optional))
+    config::Deltat = 1;
+  ParseForParameter(verbose,FileBuffer,"OutVisStep", 0, 1, 1, int_type, &(config::OutVisStep), 1, optional);
+  ParseForParameter(verbose,FileBuffer,"OutSrcStep", 0, 1, 1, int_type, &(config::OutSrcStep), 1, optional);
+  ParseForParameter(verbose,FileBuffer,"TargetTemp", 0, 1, 1, double_type, &(config::TargetTemp), 1, optional);
+  //ParseForParameter(verbose,FileBuffer,"Thermostat", 0, 1, 1, int_type, &(config::ScaleTempStep), 1, optional);
+  if (!ParseForParameter(verbose,FileBuffer,"EpsWannier", 0, 1, 1, double_type, &(config::EpsWannier), 1, optional))
+    config::EpsWannier = 1e-8;
+
+  // stop conditions
+  //if (config::MaxOuterStep <= 0) config::MaxOuterStep = 1;
+  ParseForParameter(verbose,FileBuffer,"MaxPsiStep", 0, 1, 1, int_type, &(config::MaxPsiStep), 1, critical);
+  if (config::MaxPsiStep <= 0) config::MaxPsiStep = 3;
+
+  ParseForParameter(verbose,FileBuffer,"MaxMinStep", 0, 1, 1, int_type, &(config::MaxMinStep), 1, critical);
+  ParseForParameter(verbose,FileBuffer,"RelEpsTotalE", 0, 1, 1, double_type, &(config::RelEpsTotalEnergy), 1, critical);
+  ParseForParameter(verbose,FileBuffer,"RelEpsKineticE", 0, 1, 1, double_type, &(config::RelEpsKineticEnergy), 1, critical);
+  ParseForParameter(verbose,FileBuffer,"MaxMinStopStep", 0, 1, 1, int_type, &(config::MaxMinStopStep), 1, critical);
+  ParseForParameter(verbose,FileBuffer,"MaxMinGapStopStep", 0, 1, 1, int_type, &(config::MaxMinGapStopStep), 1, critical);
+  if (config::MaxMinStep <= 0) config::MaxMinStep = config::MaxPsiStep;
+  if (config::MaxMinStopStep < 1) config::MaxMinStopStep = 1;
+  if (config::MaxMinGapStopStep < 1) config::MaxMinGapStopStep = 1;
+
+  ParseForParameter(verbose,FileBuffer,"MaxInitMinStep", 0, 1, 1, int_type, &(config::MaxInitMinStep), 1, critical);
+  ParseForParameter(verbose,FileBuffer,"InitRelEpsTotalE", 0, 1, 1, double_type, &(config::InitRelEpsTotalEnergy), 1, critical);
+  ParseForParameter(verbose,FileBuffer,"InitRelEpsKineticE", 0, 1, 1, double_type, &(config::InitRelEpsKineticEnergy), 1, critical);
+  ParseForParameter(verbose,FileBuffer,"InitMaxMinStopStep", 0, 1, 1, int_type, &(config::InitMaxMinStopStep), 1, critical);
+  ParseForParameter(verbose,FileBuffer,"InitMaxMinGapStopStep", 0, 1, 1, int_type, &(config::InitMaxMinGapStopStep), 1, critical);
+  if (config::MaxInitMinStep <= 0) config::MaxInitMinStep = config::MaxPsiStep;
+  if (config::InitMaxMinStopStep < 1) config::InitMaxMinStopStep = 1;
+  if (config::InitMaxMinGapStopStep < 1) config::InitMaxMinGapStopStep = 1;
+
+  // Unit cell and magnetic field
+  ParseForParameter(verbose,FileBuffer, "BoxLength", 0, 3, 3, lower_trigrid, BoxLength, 1, critical); /* Lattice->RealBasis */
+  mol->cell_size[0] = BoxLength[0];
+  mol->cell_size[1] = BoxLength[3];
+  mol->cell_size[2] = BoxLength[4];
+  mol->cell_size[3] = BoxLength[6];
+  mol->cell_size[4] = BoxLength[7];
+  mol->cell_size[5] = BoxLength[8];
+  //if (1) fprintf(stderr,"\n");
+
+  ParseForParameter(verbose,FileBuffer,"DoPerturbation", 0, 1, 1, int_type, &(config::DoPerturbation), 1, optional);
+  ParseForParameter(verbose,FileBuffer,"DoOutNICS", 0, 1, 1, int_type, &(config::DoOutNICS), 1, optional);
+  if (!ParseForParameter(verbose,FileBuffer,"DoFullCurrent", 0, 1, 1, int_type, &(config::DoFullCurrent), 1, optional))
+    config::DoFullCurrent = 0;
+  if (config::DoFullCurrent < 0) config::DoFullCurrent = 0;
+  if (config::DoFullCurrent > 2) config::DoFullCurrent = 2;
+  if (config::DoOutNICS < 0) config::DoOutNICS = 0;
+  if (config::DoOutNICS > 2) config::DoOutNICS = 2;
+  if (config::DoPerturbation == 0) {
+    config::DoFullCurrent = 0;
+    config::DoOutNICS = 0;
+  }
+
+  ParseForParameter(verbose,FileBuffer,"ECut", 0, 1, 1, double_type, &(config::ECut), 1, critical);
+  ParseForParameter(verbose,FileBuffer,"MaxLevel", 0, 1, 1, int_type, &(config::MaxLevel), 1, critical);
+  ParseForParameter(verbose,FileBuffer,"Level0Factor", 0, 1, 1, int_type, &(config::Lev0Factor), 1, critical);
+  if (config::Lev0Factor < 2) {
+    config::Lev0Factor = 2;
+  }
+  ParseForParameter(verbose,FileBuffer,"RiemannTensor", 0, 1, 1, int_type, &di, 1, critical);
+  if (di >= 0 && di < 2) {
+    config::RiemannTensor = di;
+  } else {
+    fprintf(stderr, "0 <= RiemanTensor < 2: 0 UseNotRT, 1 UseRT");
+    exit(1);
+  }
+  switch (config::RiemannTensor) {
+    case 0: //UseNoRT
+      if (config::MaxLevel < 2) {
+        config::MaxLevel = 2;
+      }
+      config::LevRFactor = 2;
+      config::RTActualUse = 0;
+      break;
+    case 1: // UseRT
+      if (config::MaxLevel < 3) {
+        config::MaxLevel = 3;
+      }
+      ParseForParameter(verbose,FileBuffer,"RiemannLevel", 0, 1, 1, int_type, &(config::RiemannLevel), 1, critical);
+      if (config::RiemannLevel < 2) {
+        config::RiemannLevel = 2;
+      }
+      if (config::RiemannLevel > config::MaxLevel-1) {
+        config::RiemannLevel = config::MaxLevel-1;
+      }
+      ParseForParameter(verbose,FileBuffer,"LevRFactor", 0, 1, 1, int_type, &(config::LevRFactor), 1, critical);
+      if (config::LevRFactor < 2) {
+        config::LevRFactor = 2;
+      }
+      config::Lev0Factor = 2;
+      config::RTActualUse = 2;
+      break;
+  }
+  ParseForParameter(verbose,FileBuffer,"PsiType", 0, 1, 1, int_type, &di, 1, critical);
+  if (di >= 0 && di < 2) {
+    config::PsiType = di;
+  } else {
+    fprintf(stderr, "0 <= PsiType < 2: 0 UseSpinDouble, 1 UseSpinUpDown");
+    exit(1);
+  }
+  switch (config::PsiType) {
+  case 0: // SpinDouble
+    ParseForParameter(verbose,FileBuffer,"MaxPsiDouble", 0, 1, 1, int_type, &(config::MaxPsiDouble), 1, critical);
+    ParseForParameter(verbose,FileBuffer,"AddPsis", 0, 1, 1, int_type, &(config::AddPsis), 1, optional);
+    break;
+  case 1: // SpinUpDown
+    if (config::ProcPEGamma % 2) config::ProcPEGamma*=2;
+    ParseForParameter(verbose,FileBuffer,"PsiMaxNoUp", 0, 1, 1, int_type, &(config::PsiMaxNoUp), 1, critical);
+    ParseForParameter(verbose,FileBuffer,"PsiMaxNoDown", 0, 1, 1, int_type, &(config::PsiMaxNoDown), 1, critical);
+    ParseForParameter(verbose,FileBuffer,"AddPsis", 0, 1, 1, int_type, &(config::AddPsis), 1, optional);
+    break;
+  }
+
+  // IonsInitRead
+
+  ParseForParameter(verbose,FileBuffer,"RCut", 0, 1, 1, double_type, &(config::RCut), 1, critical);
+  ParseForParameter(verbose,FileBuffer,"IsAngstroem", 0, 1, 1, int_type, &(config::IsAngstroem), 1, critical);
+  ParseForParameter(verbose,FileBuffer,"MaxTypes", 0, 1, 1, int_type, &(config::MaxTypes), 1, critical);
+  if (!ParseForParameter(verbose,FileBuffer,"RelativeCoord", 0, 1, 1, int_type, &(config::RelativeCoord) , 1, optional))
+    config::RelativeCoord = 0;
+  if (!ParseForParameter(verbose,FileBuffer,"StructOpt", 0, 1, 1, int_type, &(config::StructOpt), 1, optional))
+    config::StructOpt = 0;
+  if (MaxTypes == 0) {
+    cerr << "There are no atoms according to MaxTypes in this config file." << endl;
+  } else {
+    // prescan number of ions per type
+    cout << Verbose(0) << "Prescanning ions per type: " << endl;
+    int NoAtoms = 0;
+    for (int i=0; i < config::MaxTypes; i++) {
+      sprintf(name,"Ion_Type%i",i+1);
+      ParseForParameter(verbose,FileBuffer, (const char*)name, 0, 1, 1, int_type, &No[i], 1, critical);
+      ParseForParameter(verbose,FileBuffer, name, 0, 2, 1, int_type, &Z, 1, critical);
+      elementhash[i] = periode->FindElement(Z);
+      cout << Verbose(1) << i << ". Z = " << elementhash[i]->Z << " with " << No[i] << " ions." << endl;
+      NoAtoms += No[i];
+    }
+    int repetition = 0; // which repeated keyword shall be read
+
+    // sort the lines via the LineMapping
     sprintf(name,"Ion_Type%i",config::MaxTypes);
-                if (!ParseForParameter(verbose,FileBuffer, (const char*)name, 1, 1, 1, int_type, &value[0], 1, critical)) {
-                  cerr << "There are no atoms in the config file!" << endl;
+    if (!ParseForParameter(verbose,FileBuffer, (const char*)name, 1, 1, 1, int_type, &value[0], 1, critical)) {
+      cerr << "There are no atoms in the config file!" << endl;
       return;
     }
-                FileBuffer->CurrentLine++;
-                //cout << FileBuffer->buffer[ FileBuffer->LineMapping[FileBuffer->CurrentLine]];
-                FileBuffer->MapIonTypesInBuffer(NoAtoms);
-                //for (int i=0; i<(NoAtoms < 100 ? NoAtoms : 100 < 100 ? NoAtoms : 100);++i) {
-                //      cout << FileBuffer->buffer[ FileBuffer->LineMapping[FileBuffer->CurrentLine+i]];
-                //}
-
-                map<int, atom *> AtomList[config::MaxTypes];
-                map<int, atom *> LinearList;
-                if (!FastParsing) {
-                        // parse in trajectories
-                        bool status = true;
-                        atom *neues = NULL;
-                        while (status) {
-                                cout << "Currently parsing MD step " << repetition << "." << endl;
-                                for (int i=0; i < config::MaxTypes; i++) {
-                                        sprintf(name,"Ion_Type%i",i+1);
-                                        for(int j=0;j<No[i];j++) {
-                                                sprintf(keyword,"%s_%i",name, j+1);
-                                                if (repetition == 0) {
-                                                        neues = new atom();
-                                                        AtomList[i][j] = neues;
-                                                        LinearList[FileBuffer->CurrentLine] = neues;
-                                                        neues->type = elementhash[i]; // find element type
-                                                } else
-                                                        neues = AtomList[i][j];
+    FileBuffer->CurrentLine++;
+    //cout << FileBuffer->buffer[ FileBuffer->LineMapping[FileBuffer->CurrentLine]];
+    FileBuffer->MapIonTypesInBuffer(NoAtoms);
+    //for (int i=0; i<(NoAtoms < 100 ? NoAtoms : 100 < 100 ? NoAtoms : 100);++i) {
+    //  cout << FileBuffer->buffer[ FileBuffer->LineMapping[FileBuffer->CurrentLine+i]];
+    //}
+
+    map<int, atom *> AtomList[config::MaxTypes];
+    map<int, atom *> LinearList;
+    if (!FastParsing) {
+      // parse in trajectories
+      bool status = true;
+      atom *neues = NULL;
+      while (status) {
+        cout << "Currently parsing MD step " << repetition << "." << endl;
+        for (int i=0; i < config::MaxTypes; i++) {
+          sprintf(name,"Ion_Type%i",i+1);
+          for(int j=0;j<No[i];j++) {
+            sprintf(keyword,"%s_%i",name, j+1);
+            if (repetition == 0) {
+              neues = new atom();
+              AtomList[i][j] = neues;
+              LinearList[FileBuffer->CurrentLine] = neues;
+              neues->type = elementhash[i]; // find element type
+            } else
+              neues = AtomList[i][j];
             status = (status &&
-                                                                                ParseForParameter(verbose,FileBuffer, keyword, 0, 1, 1, double_type, &neues->x.x[0], 1, (repetition == 0) ? critical : optional) &&
-                                                                                ParseForParameter(verbose,FileBuffer, keyword, 0, 2, 1, double_type, &neues->x.x[1], 1, (repetition == 0) ? critical : optional) &&
-                                                                                ParseForParameter(verbose,FileBuffer, keyword, 0, 3, 1, double_type, &neues->x.x[2], 1, (repetition == 0) ? critical : optional) &&
-                                                                                ParseForParameter(verbose,FileBuffer, keyword, 0, 4, 1, int_type, &neues->FixedIon, 1, (repetition == 0) ? critical : optional));
-                                                if (!status) break;
-
-                                                // check size of vectors
-                                                if (mol->Trajectories[neues].R.size() <= (unsigned int)(repetition)) {
-                                                        //cout << "Increasing size for trajectory array of " << keyword << " to " << (repetition+10) << "." << endl;
-                                                        mol->Trajectories[neues].R.resize(repetition+10);
-                                                        mol->Trajectories[neues].U.resize(repetition+10);
-                                                        mol->Trajectories[neues].F.resize(repetition+10);
-                                                }
-
-                                                // put into trajectories list
-                                                for (int d=0;d<NDIM;d++)
-                                                        mol->Trajectories[neues].R.at(repetition).x[d] = neues->x.x[d];
-
-                                                // parse velocities if present
-                                                if(!ParseForParameter(verbose,FileBuffer, keyword, 0, 5, 1, double_type, &neues->v.x[0], 1,optional))
-                                                        neues->v.x[0] = 0.;
-                                                if(!ParseForParameter(verbose,FileBuffer, keyword, 0, 6, 1, double_type, &neues->v.x[1], 1,optional))
-                                                        neues->v.x[1] = 0.;
-                                                if(!ParseForParameter(verbose,FileBuffer, keyword, 0, 7, 1, double_type, &neues->v.x[2], 1,optional))
-                                                        neues->v.x[2] = 0.;
-                                                for (int d=0;d<NDIM;d++)
-                                                        mol->Trajectories[neues].U.at(repetition).x[d] = neues->v.x[d];
-
-                                                // parse forces if present
-                                                if(!ParseForParameter(verbose,FileBuffer, keyword, 0, 8, 1, double_type, &value[0], 1,optional))
-                                                        value[0] = 0.;
-                                                if(!ParseForParameter(verbose,FileBuffer, keyword, 0, 9, 1, double_type, &value[1], 1,optional))
-                                                        value[1] = 0.;
-                                                if(!ParseForParameter(verbose,FileBuffer, keyword, 1, 10, 1, double_type, &value[2], 1,optional))
-                                                        value[2] = 0.;
-                                                for (int d=0;d<NDIM;d++)
-                                                        mol->Trajectories[neues].F.at(repetition).x[d] = value[d];
-
-        //                                              cout << "Parsed position of step " << (repetition) << ": (";
-        //                                              for (int d=0;d<NDIM;d++)
-        //                                                      cout << mol->Trajectories[neues].R.at(repetition).x[d] << " ";                                  // next step
-        //                                              cout << ")\t(";
-        //                                              for (int d=0;d<NDIM;d++)
-        //                                                      cout << mol->Trajectories[neues].U.at(repetition).x[d] << " ";                                  // next step
-        //                                              cout << ")\t(";
-        //                                              for (int d=0;d<NDIM;d++)
-        //                                                      cout << mol->Trajectories[neues].F.at(repetition).x[d] << " ";                                  // next step
-        //                                              cout << ")" << endl;
-                                        }
-                                }
-                                repetition++;
-                        }
-                        // put atoms into the molecule in their original order
-                        for(map<int, atom*>::iterator runner = LinearList.begin(); runner != LinearList.end(); ++runner) {
-                                mol->AddAtom(runner->second);
-                        }
-                        repetition--;
-                        cout << "Found " << repetition << " trajectory steps." << endl;
-                        mol->MDSteps = repetition;
-                } else {
-                        // find the maximum number of MD steps so that we may parse last one (Ion_Type1_1 must always be present, because is the first atom)
-                        repetition = 0;
-                        while ( ParseForParameter(verbose,FileBuffer, "Ion_Type1_1", 0, 1, 1, double_type, &value[0], repetition, (repetition == 0) ? critical : optional) &&
-                                                        ParseForParameter(verbose,FileBuffer, "Ion_Type1_1", 0, 2, 1, double_type, &value[1], repetition, (repetition == 0) ? critical : optional) &&
-                                                        ParseForParameter(verbose,FileBuffer, "Ion_Type1_1", 0, 3, 1, double_type, &value[2], repetition, (repetition == 0) ? critical : optional))
-                                repetition++;
-                        cout << "I found " << repetition << " times the keyword Ion_Type1_1." << endl;
-                        // parse in molecule coordinates
-                        for (int i=0; i < config::MaxTypes; i++) {
-                                sprintf(name,"Ion_Type%i",i+1);
-                                for(int j=0;j<No[i];j++) {
-                                        sprintf(keyword,"%s_%i",name, j+1);
-                                        atom *neues = new atom();
-                                        // then parse for each atom the coordinates as often as present
-                                        ParseForParameter(verbose,FileBuffer, keyword, 0, 1, 1, double_type, &neues->x.x[0], repetition,critical);
-                                        ParseForParameter(verbose,FileBuffer, keyword, 0, 2, 1, double_type, &neues->x.x[1], repetition,critical);
-                                        ParseForParameter(verbose,FileBuffer, keyword, 0, 3, 1, double_type, &neues->x.x[2], repetition,critical);
-                                        ParseForParameter(verbose,FileBuffer, keyword, 0, 4, 1, int_type, &neues->FixedIon, repetition,critical);
-                                        if(!ParseForParameter(verbose,FileBuffer, keyword, 0, 5, 1, double_type, &neues->v.x[0], repetition,optional))
-                                                neues->v.x[0] = 0.;
-                                        if(!ParseForParameter(verbose,FileBuffer, keyword, 0, 6, 1, double_type, &neues->v.x[1], repetition,optional))
-                                                neues->v.x[1] = 0.;
-                                        if(!ParseForParameter(verbose,FileBuffer, keyword, 0, 7, 1, double_type, &neues->v.x[2], repetition,optional))
-                                                neues->v.x[2] = 0.;
-                                        // here we don't care if forces are present (last in trajectories is always equal to current position)
-                                        neues->type = elementhash[i]; // find element type
-                                        mol->AddAtom(neues);
-                                }
-                        }
-                }
-        }
-        delete(FileBuffer);
+                    ParseForParameter(verbose,FileBuffer, keyword, 0, 1, 1, double_type, &neues->x.x[0], 1, (repetition == 0) ? critical : optional) &&
+                    ParseForParameter(verbose,FileBuffer, keyword, 0, 2, 1, double_type, &neues->x.x[1], 1, (repetition == 0) ? critical : optional) &&
+                    ParseForParameter(verbose,FileBuffer, keyword, 0, 3, 1, double_type, &neues->x.x[2], 1, (repetition == 0) ? critical : optional) &&
+                    ParseForParameter(verbose,FileBuffer, keyword, 0, 4, 1, int_type, &neues->FixedIon, 1, (repetition == 0) ? critical : optional));
+            if (!status) break;
+
+            // check size of vectors
+            if (mol->Trajectories[neues].R.size() <= (unsigned int)(repetition)) {
+              //cout << "Increasing size for trajectory array of " << keyword << " to " << (repetition+10) << "." << endl;
+              mol->Trajectories[neues].R.resize(repetition+10);
+              mol->Trajectories[neues].U.resize(repetition+10);
+              mol->Trajectories[neues].F.resize(repetition+10);
+            }
+
+            // put into trajectories list
+            for (int d=0;d<NDIM;d++)
+              mol->Trajectories[neues].R.at(repetition).x[d] = neues->x.x[d];
+
+            // parse velocities if present
+            if(!ParseForParameter(verbose,FileBuffer, keyword, 0, 5, 1, double_type, &neues->v.x[0], 1,optional))
+              neues->v.x[0] = 0.;
+            if(!ParseForParameter(verbose,FileBuffer, keyword, 0, 6, 1, double_type, &neues->v.x[1], 1,optional))
+              neues->v.x[1] = 0.;
+            if(!ParseForParameter(verbose,FileBuffer, keyword, 0, 7, 1, double_type, &neues->v.x[2], 1,optional))
+              neues->v.x[2] = 0.;
+            for (int d=0;d<NDIM;d++)
+              mol->Trajectories[neues].U.at(repetition).x[d] = neues->v.x[d];
+
+            // parse forces if present
+            if(!ParseForParameter(verbose,FileBuffer, keyword, 0, 8, 1, double_type, &value[0], 1,optional))
+              value[0] = 0.;
+            if(!ParseForParameter(verbose,FileBuffer, keyword, 0, 9, 1, double_type, &value[1], 1,optional))
+              value[1] = 0.;
+            if(!ParseForParameter(verbose,FileBuffer, keyword, 1, 10, 1, double_type, &value[2], 1,optional))
+              value[2] = 0.;
+            for (int d=0;d<NDIM;d++)
+              mol->Trajectories[neues].F.at(repetition).x[d] = value[d];
+
+  //            cout << "Parsed position of step " << (repetition) << ": (";
+  //            for (int d=0;d<NDIM;d++)
+  //              cout << mol->Trajectories[neues].R.at(repetition).x[d] << " ";          // next step
+  //            cout << ")\t(";
+  //            for (int d=0;d<NDIM;d++)
+  //              cout << mol->Trajectories[neues].U.at(repetition).x[d] << " ";          // next step
+  //            cout << ")\t(";
+  //            for (int d=0;d<NDIM;d++)
+  //              cout << mol->Trajectories[neues].F.at(repetition).x[d] << " ";          // next step
+  //            cout << ")" << endl;
+          }
+        }
+        repetition++;
+      }
+      // put atoms into the molecule in their original order
+      for(map<int, atom*>::iterator runner = LinearList.begin(); runner != LinearList.end(); ++runner) {
+        mol->AddAtom(runner->second);
+      }
+      repetition--;
+      cout << "Found " << repetition << " trajectory steps." << endl;
+      mol->MDSteps = repetition;
+    } else {
+      // find the maximum number of MD steps so that we may parse last one (Ion_Type1_1 must always be present, because is the first atom)
+      repetition = 0;
+      while ( ParseForParameter(verbose,FileBuffer, "Ion_Type1_1", 0, 1, 1, double_type, &value[0], repetition, (repetition == 0) ? critical : optional) &&
+              ParseForParameter(verbose,FileBuffer, "Ion_Type1_1", 0, 2, 1, double_type, &value[1], repetition, (repetition == 0) ? critical : optional) &&
+              ParseForParameter(verbose,FileBuffer, "Ion_Type1_1", 0, 3, 1, double_type, &value[2], repetition, (repetition == 0) ? critical : optional))
+        repetition++;
+      cout << "I found " << repetition << " times the keyword Ion_Type1_1." << endl;
+      // parse in molecule coordinates
+      for (int i=0; i < config::MaxTypes; i++) {
+        sprintf(name,"Ion_Type%i",i+1);
+        for(int j=0;j<No[i];j++) {
+          sprintf(keyword,"%s_%i",name, j+1);
+          atom *neues = new atom();
+          // then parse for each atom the coordinates as often as present
+          ParseForParameter(verbose,FileBuffer, keyword, 0, 1, 1, double_type, &neues->x.x[0], repetition,critical);
+          ParseForParameter(verbose,FileBuffer, keyword, 0, 2, 1, double_type, &neues->x.x[1], repetition,critical);
+          ParseForParameter(verbose,FileBuffer, keyword, 0, 3, 1, double_type, &neues->x.x[2], repetition,critical);
+          ParseForParameter(verbose,FileBuffer, keyword, 0, 4, 1, int_type, &neues->FixedIon, repetition,critical);
+          if(!ParseForParameter(verbose,FileBuffer, keyword, 0, 5, 1, double_type, &neues->v.x[0], repetition,optional))
+            neues->v.x[0] = 0.;
+          if(!ParseForParameter(verbose,FileBuffer, keyword, 0, 6, 1, double_type, &neues->v.x[1], repetition,optional))
+            neues->v.x[1] = 0.;
+          if(!ParseForParameter(verbose,FileBuffer, keyword, 0, 7, 1, double_type, &neues->v.x[2], repetition,optional))
+            neues->v.x[2] = 0.;
+          // here we don't care if forces are present (last in trajectories is always equal to current position)
+          neues->type = elementhash[i]; // find element type
+          mol->AddAtom(neues);
+        }
+      }
+    }
+  }
+  delete(FileBuffer);
 };
 
@@ -943 +943 @@
 void config::LoadOld(char *filename, periodentafel *periode, molecule *mol)
 {
-        ifstream *file = new ifstream(filename);
-        if (file == NULL) {
-                cerr << "ERROR: config file " << filename << " missing!" << endl;
-                return;
-        }
-        RetrieveConfigPathAndName(filename);
-        // ParseParameters
-
-        /* Oeffne Hauptparameterdatei */
-        int l, i, di;
-        double a,b;
-        double BoxLength[9];
-        string zeile;
-        string dummy;
-        element *elementhash[128];
-        int Z, No, AtomNo, found;
-        int verbose = 0;
-
-        /* Namen einlesen */
-
-        ParseForParameter(verbose,file, "mainname", 0, 1, 1, string_type, (config::mainname), 1, critical);
-        ParseForParameter(verbose,file, "defaultpath", 0, 1, 1, string_type, (config::defaultpath), 1, critical);
-        ParseForParameter(verbose,file, "pseudopotpath", 0, 1, 1, string_type, (config::pseudopotpath), 1, critical);
-        ParseForParameter(verbose,file, "ProcsGammaPsi", 0, 1, 1, int_type, &(config::ProcPEGamma), 1, critical);
-        ParseForParameter(verbose,file, "ProcsGammaPsi", 0, 2, 1, int_type, &(config::ProcPEPsi), 1, critical);
-        config::Seed = 1;
-        config::DoOutOrbitals = 0;
-        ParseForParameter(verbose,file,"DoOutVis", 0, 1, 1, int_type, &(config::DoOutVis), 1, critical);
-        if (config::DoOutVis < 0) config::DoOutVis = 0;
-        if (config::DoOutVis > 1) config::DoOutVis = 1;
-        config::VectorPlane = -1;
-        config::VectorCut = 0.;
-        ParseForParameter(verbose,file,"DoOutMes", 0, 1, 1, int_type, &(config::DoOutMes), 1, critical);
-        if (config::DoOutMes < 0) config::DoOutMes = 0;
-        if (config::DoOutMes > 1) config::DoOutMes = 1;
-        config::DoOutCurrent = 0;
-        ParseForParameter(verbose,file,"AddGramSch", 0, 1, 1, int_type, &(config::UseAddGramSch), 1, critical);
-        if (config::UseAddGramSch < 0) config::UseAddGramSch = 0;
-        if (config::UseAddGramSch > 2) config::UseAddGramSch = 2;
-        config::CommonWannier = 0;
-        config::SawtoothStart = 0.01;
-
-        ParseForParameter(verbose,file,"MaxOuterStep", 0, 1, 1, double_type, &(config::MaxOuterStep), 1, critical);
-        ParseForParameter(verbose,file,"Deltat", 0, 1, 1, double_type, &(config::Deltat), 1, optional);
-        ParseForParameter(verbose,file,"VisOuterStep", 0, 1, 1, int_type, &(config::OutVisStep), 1, optional);
-        ParseForParameter(verbose,file,"VisSrcOuterStep", 0, 1, 1, int_type, &(config::OutSrcStep), 1, optional);
-        ParseForParameter(verbose,file,"TargetTemp", 0, 1, 1, double_type, &(config::TargetTemp), 1, optional);
-        ParseForParameter(verbose,file,"ScaleTempStep", 0, 1, 1, int_type, &(config::ScaleTempStep), 1, optional);
-        config::EpsWannier = 1e-8;
-
-        // stop conditions
-        //if (config::MaxOuterStep <= 0) config::MaxOuterStep = 1;
-        ParseForParameter(verbose,file,"MaxPsiStep", 0, 1, 1, int_type, &(config::MaxPsiStep), 1, critical);
-        if (config::MaxPsiStep <= 0) config::MaxPsiStep = 3;
-
-        ParseForParameter(verbose,file,"MaxMinStep", 0, 1, 1, int_type, &(config::MaxMinStep), 1, critical);
-        ParseForParameter(verbose,file,"MaxMinStep", 0, 2, 1, double_type, &(config::RelEpsTotalEnergy), 1, critical);
-        ParseForParameter(verbose,file,"MaxMinStep", 0, 3, 1, double_type, &(config::RelEpsKineticEnergy), 1, critical);
-        ParseForParameter(verbose,file,"MaxMinStep", 0, 4, 1, int_type, &(config::MaxMinStopStep), 1, critical);
-        if (config::MaxMinStep <= 0) config::MaxMinStep = config::MaxPsiStep;
-        if (config::MaxMinStopStep < 1) config::MaxMinStopStep = 1;
-        config::MaxMinGapStopStep = 1;
-
-        ParseForParameter(verbose,file,"MaxInitMinStep", 0, 1, 1, int_type, &(config::MaxInitMinStep), 1, critical);
-        ParseForParameter(verbose,file,"MaxInitMinStep", 0, 2, 1, double_type, &(config::InitRelEpsTotalEnergy), 1, critical);
-        ParseForParameter(verbose,file,"MaxInitMinStep", 0, 3, 1, double_type, &(config::InitRelEpsKineticEnergy), 1, critical);
-        ParseForParameter(verbose,file,"MaxInitMinStep", 0, 4, 1, int_type, &(config::InitMaxMinStopStep), 1, critical);
-        if (config::MaxInitMinStep <= 0) config::MaxInitMinStep = config::MaxPsiStep;
-        if (config::InitMaxMinStopStep < 1) config::InitMaxMinStopStep = 1;
-        config::InitMaxMinGapStopStep = 1;
-
-        ParseForParameter(verbose,file, "BoxLength", 0, 3, 3, lower_trigrid, BoxLength, 1, critical); /* Lattice->RealBasis */
-        mol->cell_size[0] = BoxLength[0];
-        mol->cell_size[1] = BoxLength[3];
-        mol->cell_size[2] = BoxLength[4];
-        mol->cell_size[3] = BoxLength[6];
-        mol->cell_size[4] = BoxLength[7];
-        mol->cell_size[5] = BoxLength[8];
-        if (1) fprintf(stderr,"\n");
-        config::DoPerturbation = 0;
-        config::DoFullCurrent = 0;
-
-        ParseForParameter(verbose,file,"ECut", 0, 1, 1, double_type, &(config::ECut), 1, critical);
-        ParseForParameter(verbose,file,"MaxLevel", 0, 1, 1, int_type, &(config::MaxLevel), 1, critical);
-        ParseForParameter(verbose,file,"Level0Factor", 0, 1, 1, int_type, &(config::Lev0Factor), 1, critical);
-        if (config::Lev0Factor < 2) {
-                config::Lev0Factor = 2;
-        }
-        ParseForParameter(verbose,file,"RiemannTensor", 0, 1, 1, int_type, &di, 1, critical);
-        if (di >= 0 && di < 2) {
-                config::RiemannTensor = di;
-        } else {
-                fprintf(stderr, "0 <= RiemanTensor < 2: 0 UseNotRT, 1 UseRT");
-                exit(1);
-        }
-        switch (config::RiemannTensor) {
-                case 0: //UseNoRT
-                        if (config::MaxLevel < 2) {
-                                config::MaxLevel = 2;
-                        }
-                        config::LevRFactor = 2;
-                        config::RTActualUse = 0;
-                        break;
-                case 1: // UseRT
-                        if (config::MaxLevel < 3) {
-                                config::MaxLevel = 3;
-                        }
-                        ParseForParameter(verbose,file,"RiemannLevel", 0, 1, 1, int_type, &(config::RiemannLevel), 1, critical);
-                        if (config::RiemannLevel < 2) {
-                                config::RiemannLevel = 2;
-                        }
-                        if (config::RiemannLevel > config::MaxLevel-1) {
-                                config::RiemannLevel = config::MaxLevel-1;
-                        }
-                        ParseForParameter(verbose,file,"LevRFactor", 0, 1, 1, int_type, &(config::LevRFactor), 1, critical);
-                        if (config::LevRFactor < 2) {
-                                config::LevRFactor = 2;
-                        }
-                        config::Lev0Factor = 2;
-                        config::RTActualUse = 2;
-                        break;
-        }
-        ParseForParameter(verbose,file,"PsiType", 0, 1, 1, int_type, &di, 1, critical);
-        if (di >= 0 && di < 2) {
-                config::PsiType = di;
-        } else {
-                fprintf(stderr, "0 <= PsiType < 2: 0 UseSpinDouble, 1 UseSpinUpDown");
-                exit(1);
-        }
-        switch (config::PsiType) {
-        case 0: // SpinDouble
-                ParseForParameter(verbose,file,"MaxPsiDouble", 0, 1, 1, int_type, &(config::MaxPsiDouble), 1, critical);
-                config::AddPsis = 0;
-                break;
-        case 1: // SpinUpDown
-                if (config::ProcPEGamma % 2) config::ProcPEGamma*=2;
-                ParseForParameter(verbose,file,"MaxPsiUp", 0, 1, 1, int_type, &(config::PsiMaxNoUp), 1, critical);
-                ParseForParameter(verbose,file,"MaxPsiDown", 0, 1, 1, int_type, &(config::PsiMaxNoDown), 1, critical);
-                config::AddPsis = 0;
-                break;
-        }
-
-        // IonsInitRead
-
-        ParseForParameter(verbose,file,"RCut", 0, 1, 1, double_type, &(config::RCut), 1, critical);
-        ParseForParameter(verbose,file,"IsAngstroem", 0, 1, 1, int_type, &(config::IsAngstroem), 1, critical);
-        config::RelativeCoord = 0;
-        config::StructOpt = 0;
-
-        // Routine from builder.cpp
-
-
-        for (i=MAX_ELEMENTS;i--;)
-                elementhash[i] = NULL;
-        cout << Verbose(0) << "Parsing Ions ..." << endl;
-        No=0;
-        found = 0;
-        while (getline(*file,zeile,'\n')) {
-                if (zeile.find("Ions_Data") == 0) {
-                        cout << Verbose(1) << "found Ions_Data...begin parsing" << endl;
-                        found ++;
-                }
-                if (found > 0) {
-                        if (zeile.find("Ions_Data") == 0)
-                                getline(*file,zeile,'\n'); // read next line and parse this one
-                        istringstream input(zeile);
-                        input >> AtomNo;        // number of atoms
-                        input >> Z;                      // atomic number
-                        input >> a;
-                        input >> l;
-                        input >> l;
-                        input >> b;              // element mass
-                        elementhash[No] = periode->FindElement(Z);
-                        cout << Verbose(1) << "AtomNo: " << AtomNo << "\tZ: " << Z << "\ta:" << a << "\tl:"     << l << "\b:" << b << "\tElement:" << elementhash[No] << "\t:" << endl;
-                        for(i=0;i<AtomNo;i++) {
-                                if (!getline(*file,zeile,'\n')) {// parse on and on
-                                        cout << Verbose(2) << "Error: Too few items in ionic list of element" << elementhash[No] << "." << endl << "Exiting." << endl;
-                                        // return 1;
-                                } else {
-                                        //cout << Verbose(2) << "Reading line: " << zeile << endl;
-                                }
-                                istringstream input2(zeile);
-                                atom *neues = new atom();
-                                input2 >> neues->x.x[0]; // x
-                                input2 >> neues->x.x[1]; // y
-                                input2 >> neues->x.x[2]; // z
-                                input2 >> l;
-                                neues->type = elementhash[No]; // find element type
-                                mol->AddAtom(neues);
-                        }
-                        No++;
-                }
-        }
-        file->close();
-        delete(file);
+  ifstream *file = new ifstream(filename);
+  if (file == NULL) {
+    cerr << "ERROR: config file " << filename << " missing!" << endl;
+    return;
+  }
+  RetrieveConfigPathAndName(filename);
+  // ParseParameters
+
+  /* Oeffne Hauptparameterdatei */
+  int l, i, di;
+  double a,b;
+  double BoxLength[9];
+  string zeile;
+  string dummy;
+  element *elementhash[128];
+  int Z, No, AtomNo, found;
+  int verbose = 0;
+
+  /* Namen einlesen */
+
+  ParseForParameter(verbose,file, "mainname", 0, 1, 1, string_type, (config::mainname), 1, critical);
+  ParseForParameter(verbose,file, "defaultpath", 0, 1, 1, string_type, (config::defaultpath), 1, critical);
+  ParseForParameter(verbose,file, "pseudopotpath", 0, 1, 1, string_type, (config::pseudopotpath), 1, critical);
+  ParseForParameter(verbose,file, "ProcsGammaPsi", 0, 1, 1, int_type, &(config::ProcPEGamma), 1, critical);
+  ParseForParameter(verbose,file, "ProcsGammaPsi", 0, 2, 1, int_type, &(config::ProcPEPsi), 1, critical);
+  config::Seed = 1;
+  config::DoOutOrbitals = 0;
+  ParseForParameter(verbose,file,"DoOutVis", 0, 1, 1, int_type, &(config::DoOutVis), 1, critical);
+  if (config::DoOutVis < 0) config::DoOutVis = 0;
+  if (config::DoOutVis > 1) config::DoOutVis = 1;
+  config::VectorPlane = -1;
+  config::VectorCut = 0.;
+  ParseForParameter(verbose,file,"DoOutMes", 0, 1, 1, int_type, &(config::DoOutMes), 1, critical);
+  if (config::DoOutMes < 0) config::DoOutMes = 0;
+  if (config::DoOutMes > 1) config::DoOutMes = 1;
+  config::DoOutCurrent = 0;
+  ParseForParameter(verbose,file,"AddGramSch", 0, 1, 1, int_type, &(config::UseAddGramSch), 1, critical);
+  if (config::UseAddGramSch < 0) config::UseAddGramSch = 0;
+  if (config::UseAddGramSch > 2) config::UseAddGramSch = 2;
+  config::CommonWannier = 0;
+  config::SawtoothStart = 0.01;
+
+  ParseForParameter(verbose,file,"MaxOuterStep", 0, 1, 1, double_type, &(config::MaxOuterStep), 1, critical);
+  ParseForParameter(verbose,file,"Deltat", 0, 1, 1, double_type, &(config::Deltat), 1, optional);
+  ParseForParameter(verbose,file,"VisOuterStep", 0, 1, 1, int_type, &(config::OutVisStep), 1, optional);
+  ParseForParameter(verbose,file,"VisSrcOuterStep", 0, 1, 1, int_type, &(config::OutSrcStep), 1, optional);
+  ParseForParameter(verbose,file,"TargetTemp", 0, 1, 1, double_type, &(config::TargetTemp), 1, optional);
+  ParseForParameter(verbose,file,"ScaleTempStep", 0, 1, 1, int_type, &(config::ScaleTempStep), 1, optional);
+  config::EpsWannier = 1e-8;
+
+  // stop conditions
+  //if (config::MaxOuterStep <= 0) config::MaxOuterStep = 1;
+  ParseForParameter(verbose,file,"MaxPsiStep", 0, 1, 1, int_type, &(config::MaxPsiStep), 1, critical);
+  if (config::MaxPsiStep <= 0) config::MaxPsiStep = 3;
+
+  ParseForParameter(verbose,file,"MaxMinStep", 0, 1, 1, int_type, &(config::MaxMinStep), 1, critical);
+  ParseForParameter(verbose,file,"MaxMinStep", 0, 2, 1, double_type, &(config::RelEpsTotalEnergy), 1, critical);
+  ParseForParameter(verbose,file,"MaxMinStep", 0, 3, 1, double_type, &(config::RelEpsKineticEnergy), 1, critical);
+  ParseForParameter(verbose,file,"MaxMinStep", 0, 4, 1, int_type, &(config::MaxMinStopStep), 1, critical);
+  if (config::MaxMinStep <= 0) config::MaxMinStep = config::MaxPsiStep;
+  if (config::MaxMinStopStep < 1) config::MaxMinStopStep = 1;
+  config::MaxMinGapStopStep = 1;
+
+  ParseForParameter(verbose,file,"MaxInitMinStep", 0, 1, 1, int_type, &(config::MaxInitMinStep), 1, critical);
+  ParseForParameter(verbose,file,"MaxInitMinStep", 0, 2, 1, double_type, &(config::InitRelEpsTotalEnergy), 1, critical);
+  ParseForParameter(verbose,file,"MaxInitMinStep", 0, 3, 1, double_type, &(config::InitRelEpsKineticEnergy), 1, critical);
+  ParseForParameter(verbose,file,"MaxInitMinStep", 0, 4, 1, int_type, &(config::InitMaxMinStopStep), 1, critical);
+  if (config::MaxInitMinStep <= 0) config::MaxInitMinStep = config::MaxPsiStep;
+  if (config::InitMaxMinStopStep < 1) config::InitMaxMinStopStep = 1;
+  config::InitMaxMinGapStopStep = 1;
+
+  ParseForParameter(verbose,file, "BoxLength", 0, 3, 3, lower_trigrid, BoxLength, 1, critical); /* Lattice->RealBasis */
+  mol->cell_size[0] = BoxLength[0];
+  mol->cell_size[1] = BoxLength[3];
+  mol->cell_size[2] = BoxLength[4];
+  mol->cell_size[3] = BoxLength[6];
+  mol->cell_size[4] = BoxLength[7];
+  mol->cell_size[5] = BoxLength[8];
+  if (1) fprintf(stderr,"\n");
+  config::DoPerturbation = 0;
+  config::DoFullCurrent = 0;
+
+  ParseForParameter(verbose,file,"ECut", 0, 1, 1, double_type, &(config::ECut), 1, critical);
+  ParseForParameter(verbose,file,"MaxLevel", 0, 1, 1, int_type, &(config::MaxLevel), 1, critical);
+  ParseForParameter(verbose,file,"Level0Factor", 0, 1, 1, int_type, &(config::Lev0Factor), 1, critical);
+  if (config::Lev0Factor < 2) {
+    config::Lev0Factor = 2;
+  }
+  ParseForParameter(verbose,file,"RiemannTensor", 0, 1, 1, int_type, &di, 1, critical);
+  if (di >= 0 && di < 2) {
+    config::RiemannTensor = di;
+  } else {
+    fprintf(stderr, "0 <= RiemanTensor < 2: 0 UseNotRT, 1 UseRT");
+    exit(1);
+  }
+  switch (config::RiemannTensor) {
+    case 0: //UseNoRT
+      if (config::MaxLevel < 2) {
+        config::MaxLevel = 2;
+      }
+      config::LevRFactor = 2;
+      config::RTActualUse = 0;
+      break;
+    case 1: // UseRT
+      if (config::MaxLevel < 3) {
+        config::MaxLevel = 3;
+      }
+      ParseForParameter(verbose,file,"RiemannLevel", 0, 1, 1, int_type, &(config::RiemannLevel), 1, critical);
+      if (config::RiemannLevel < 2) {
+        config::RiemannLevel = 2;
+      }
+      if (config::RiemannLevel > config::MaxLevel-1) {
+        config::RiemannLevel = config::MaxLevel-1;
+      }
+      ParseForParameter(verbose,file,"LevRFactor", 0, 1, 1, int_type, &(config::LevRFactor), 1, critical);
+      if (config::LevRFactor < 2) {
+        config::LevRFactor = 2;
+      }
+      config::Lev0Factor = 2;
+      config::RTActualUse = 2;
+      break;
+  }
+  ParseForParameter(verbose,file,"PsiType", 0, 1, 1, int_type, &di, 1, critical);
+  if (di >= 0 && di < 2) {
+    config::PsiType = di;
+  } else {
+    fprintf(stderr, "0 <= PsiType < 2: 0 UseSpinDouble, 1 UseSpinUpDown");
+    exit(1);
+  }
+  switch (config::PsiType) {
+  case 0: // SpinDouble
+    ParseForParameter(verbose,file,"MaxPsiDouble", 0, 1, 1, int_type, &(config::MaxPsiDouble), 1, critical);
+    config::AddPsis = 0;
+    break;
+  case 1: // SpinUpDown
+    if (config::ProcPEGamma % 2) config::ProcPEGamma*=2;
+    ParseForParameter(verbose,file,"MaxPsiUp", 0, 1, 1, int_type, &(config::PsiMaxNoUp), 1, critical);
+    ParseForParameter(verbose,file,"MaxPsiDown", 0, 1, 1, int_type, &(config::PsiMaxNoDown), 1, critical);
+    config::AddPsis = 0;
+    break;
+  }
+
+  // IonsInitRead
+
+  ParseForParameter(verbose,file,"RCut", 0, 1, 1, double_type, &(config::RCut), 1, critical);
+  ParseForParameter(verbose,file,"IsAngstroem", 0, 1, 1, int_type, &(config::IsAngstroem), 1, critical);
+  config::RelativeCoord = 0;
+  config::StructOpt = 0;
+
+  // Routine from builder.cpp
+
+
+  for (i=MAX_ELEMENTS;i--;)
+    elementhash[i] = NULL;
+  cout << Verbose(0) << "Parsing Ions ..." << endl;
+  No=0;
+  found = 0;
+  while (getline(*file,zeile,'\n')) {
+    if (zeile.find("Ions_Data") == 0) {
+      cout << Verbose(1) << "found Ions_Data...begin parsing" << endl;
+      found ++;
+    }
+    if (found > 0) {
+      if (zeile.find("Ions_Data") == 0)
+        getline(*file,zeile,'\n'); // read next line and parse this one
+      istringstream input(zeile);
+      input >> AtomNo;  // number of atoms
+      input >> Z;       // atomic number
+      input >> a;
+      input >> l;
+      input >> l;
+      input >> b;     // element mass
+      elementhash[No] = periode->FindElement(Z);
+      cout << Verbose(1) << "AtomNo: " << AtomNo << "\tZ: " << Z << "\ta:" << a << "\tl:"  << l << "\b:" << b << "\tElement:" << elementhash[No] << "\t:" << endl;
+      for(i=0;i<AtomNo;i++) {
+        if (!getline(*file,zeile,'\n')) {// parse on and on
+          cout << Verbose(2) << "Error: Too few items in ionic list of element" << elementhash[No] << "." << endl << "Exiting." << endl;
+          // return 1;
+        } else {
+          //cout << Verbose(2) << "Reading line: " << zeile << endl;
+        }
+        istringstream input2(zeile);
+        atom *neues = new atom();
+        input2 >> neues->x.x[0]; // x
+        input2 >> neues->x.x[1]; // y
+        input2 >> neues->x.x[2]; // z
+        input2 >> l;
+        neues->type = elementhash[No]; // find element type
+        mol->AddAtom(neues);
+      }
+      No++;
+    }
+  }
+  file->close();
+  delete(file);
 };
 
@@ -1147 +1147 @@
 bool config::Save(const char *filename, periodentafel *periode, molecule *mol) const
 {
-        bool result = true;
-        // bring MaxTypes up to date
-        mol->CountElements();
-        ofstream *output = NULL;
-        output = new ofstream(filename, ios::out);
-        if (output != NULL) {
-                *output << "# ParallelCarParinello - main configuration file - created with molecuilder" << endl;
-                *output << endl;
-                *output << "mainname\t" << config::mainname << "\t# programm name (for runtime files)" << endl;
-                *output << "defaultpath\t" << config::defaultpath << "\t# where to put files during runtime" << endl;
-                *output << "pseudopotpath\t" << config::pseudopotpath << "\t# where to find pseudopotentials" << endl;
-                *output << endl;
-                *output << "ProcPEGamma\t" << config::ProcPEGamma << "\t# for parallel computing: share constants" << endl;
-                *output << "ProcPEPsi\t" << config::ProcPEPsi << "\t# for parallel computing: share wave functions" << endl;
-                *output << "DoOutVis\t" << config::DoOutVis << "\t# Output data for OpenDX" << endl;
-                *output << "DoOutMes\t" << config::DoOutMes << "\t# Output data for measurements" << endl;
-                *output << "DoOutOrbitals\t" << config::DoOutOrbitals << "\t# Output all Orbitals" << endl;
-                *output << "DoOutCurr\t" << config::DoOutCurrent << "\t# Ouput current density for OpenDx" << endl;
-                *output << "DoOutNICS\t" << config::DoOutNICS << "\t# Output Nucleus independent current shieldings" << endl;
-                *output << "DoPerturbation\t" << config::DoPerturbation << "\t# Do perturbation calculate and determine susceptibility and shielding" << endl;
-                *output << "DoFullCurrent\t" << config::DoFullCurrent << "\t# Do full perturbation" << endl;
-                *output << "CommonWannier\t" << config::CommonWannier << "\t# Put virtual centers at indivual orbits, all common, merged by variance, to grid point, to cell center" << endl;
-                *output << "SawtoothStart\t" << config::SawtoothStart << "\t# Absolute value for smooth transition at cell border " << endl;
-                *output << "VectorPlane\t" << config::VectorPlane << "\t# Cut plane axis (x, y or z: 0,1,2) for two-dim current vector plot" << endl;
-                *output << "VectorCut\t" << config::VectorCut << "\t# Cut plane axis value" << endl;
-                *output << "AddGramSch\t" << config::UseAddGramSch << "\t# Additional GramSchmidtOrtogonalization to be safe" << endl;
-                *output << "Seed\t\t" << config::Seed << "\t# initial value for random seed for Psi coefficients" << endl;
-                *output << endl;
-                *output << "MaxOuterStep\t" << config::MaxOuterStep << "\t# number of MolecularDynamics/Structure optimization steps" << endl;
-                *output << "Deltat\t" << config::Deltat << "\t# time per MD step" << endl;
-                *output << "OutVisStep\t" << config::OutVisStep << "\t# Output visual data every ...th step" << endl;
-                *output << "OutSrcStep\t" << config::OutSrcStep << "\t# Output \"restart\" data every ..th step" << endl;
-                *output << "TargetTemp\t" << config::TargetTemp << "\t# Target temperature" << endl;
-                *output << "MaxPsiStep\t" << config::MaxPsiStep << "\t# number of Minimisation steps per state (0 - default)" << endl;
-                *output << "EpsWannier\t" << config::EpsWannier << "\t# tolerance value for spread minimisation of orbitals" << endl;
-                *output << endl;
-                *output << "# Values specifying when to stop" << endl;
-                *output << "MaxMinStep\t" << config::MaxMinStep << "\t# Maximum number of steps" << endl;
-                *output << "RelEpsTotalE\t" << config::RelEpsTotalEnergy << "\t# relative change in total energy" << endl;
-                *output << "RelEpsKineticE\t" << config::RelEpsKineticEnergy << "\t# relative change in kinetic energy" << endl;
-                *output << "MaxMinStopStep\t" << config::MaxMinStopStep << "\t# check every ..th steps" << endl;
-                *output << "MaxMinGapStopStep\t" << config::MaxMinGapStopStep << "\t# check every ..th steps" << endl;
-                *output << endl;
-                *output << "# Values specifying when to stop for INIT, otherwise same as above" << endl;
-                *output << "MaxInitMinStep\t" << config::MaxInitMinStep << "\t# Maximum number of steps" << endl;
-                *output << "InitRelEpsTotalE\t" << config::InitRelEpsTotalEnergy << "\t# relative change in total energy" << endl;
-                *output << "InitRelEpsKineticE\t" << config::InitRelEpsKineticEnergy << "\t# relative change in kinetic energy" << endl;
-                *output << "InitMaxMinStopStep\t" << config::InitMaxMinStopStep << "\t# check every ..th steps" << endl;
-                *output << "InitMaxMinGapStopStep\t" << config::InitMaxMinGapStopStep << "\t# check every ..th steps" << endl;
-                *output << endl;
-                *output << "BoxLength\t\t\t# (Length of a unit cell)" << endl;
-                *output << mol->cell_size[0] << "\t" << endl;
-                *output << mol->cell_size[1] << "\t" << mol->cell_size[2] << "\t" << endl;
-                *output << mol->cell_size[3] << "\t" << mol->cell_size[4] << "\t" << mol->cell_size[5] << "\t" << endl;
-                // FIXME
-                *output << endl;
-                *output << "ECut\t\t" << config::ECut << "\t# energy cutoff for discretization in Hartrees" << endl;
-                *output << "MaxLevel\t" << config::MaxLevel << "\t# number of different levels in the code, >=2" << endl;
-                *output << "Level0Factor\t" << config::Lev0Factor << "\t# factor by which node number increases from S to 0 level" << endl;
-                *output << "RiemannTensor\t" << config::RiemannTensor << "\t# (Use metric)" << endl;
-                switch (config::RiemannTensor) {
-                        case 0: //UseNoRT
-                                break;
-                        case 1: // UseRT
-                                *output << "RiemannLevel\t" << config::RiemannLevel << "\t# Number of Riemann Levels" << endl;
-                                *output << "LevRFactor\t" << config::LevRFactor << "\t# factor by which node number increases from 0 to R level from" << endl;
-                                break;
-                }
-                *output << "PsiType\t\t" << config::PsiType << "\t# 0 - doubly occupied, 1 - SpinUp,SpinDown" << endl;
-                // write out both types for easier changing afterwards
-        //      switch (PsiType) {
-        //              case 0:
-                                *output << "MaxPsiDouble\t" << config::MaxPsiDouble << "\t# here: specifying both maximum number of SpinUp- and -Down-states" << endl;
-        //                      break;
-        //              case 1:
-                                *output << "PsiMaxNoUp\t" << config::PsiMaxNoUp << "\t# here: specifying maximum number of SpinUp-states" << endl;
-                                *output << "PsiMaxNoDown\t" << config::PsiMaxNoDown << "\t# here: specifying maximum number of SpinDown-states" << endl;
-        //                      break;
-        //      }
-                *output << "AddPsis\t\t" << config::AddPsis << "\t# Additional unoccupied Psis for bandgap determination" << endl;
-                *output << endl;
-                *output << "RCut\t\t" << config::RCut << "\t# R-cut for the ewald summation" << endl;
-                *output << "StructOpt\t" << config::StructOpt << "\t# Do structure optimization beforehand" << endl;
-                *output << "IsAngstroem\t" << config::IsAngstroem << "\t# 0 - Bohr, 1 - Angstroem" << endl;
-                *output << "RelativeCoord\t" << config::RelativeCoord << "\t# whether ion coordinates are relative (1) or absolute (0)" << endl;
-                *output << "MaxTypes\t" << mol->ElementCount << "\t# maximum number of different ion types" << endl;
-                *output << endl;
-                result = result && mol->Checkout(output);
-                if (mol->MDSteps <=1 )
-                        result = result && mol->Output(output);
-                else
-                        result = result && mol->OutputTrajectories(output);
-                output->close();
-                output->clear();
-                delete(output);
-                return result;
-        } else
-                return false;
+  bool result = true;
+  // bring MaxTypes up to date
+  mol->CountElements();
+  ofstream *output = NULL;
+  output = new ofstream(filename, ios::out);
+  if (output != NULL) {
+    *output << "# ParallelCarParinello - main configuration file - created with molecuilder" << endl;
+    *output << endl;
+    *output << "mainname\t" << config::mainname << "\t# programm name (for runtime files)" << endl;
+    *output << "defaultpath\t" << config::defaultpath << "\t# where to put files during runtime" << endl;
+    *output << "pseudopotpath\t" << config::pseudopotpath << "\t# where to find pseudopotentials" << endl;
+    *output << endl;
+    *output << "ProcPEGamma\t" << config::ProcPEGamma << "\t# for parallel computing: share constants" << endl;
+    *output << "ProcPEPsi\t" << config::ProcPEPsi << "\t# for parallel computing: share wave functions" << endl;
+    *output << "DoOutVis\t" << config::DoOutVis << "\t# Output data for OpenDX" << endl;
+    *output << "DoOutMes\t" << config::DoOutMes << "\t# Output data for measurements" << endl;
+    *output << "DoOutOrbitals\t" << config::DoOutOrbitals << "\t# Output all Orbitals" << endl;
+    *output << "DoOutCurr\t" << config::DoOutCurrent << "\t# Ouput current density for OpenDx" << endl;
+    *output << "DoOutNICS\t" << config::DoOutNICS << "\t# Output Nucleus independent current shieldings" << endl;
+    *output << "DoPerturbation\t" << config::DoPerturbation << "\t# Do perturbation calculate and determine susceptibility and shielding" << endl;
+    *output << "DoFullCurrent\t" << config::DoFullCurrent << "\t# Do full perturbation" << endl;
+    *output << "CommonWannier\t" << config::CommonWannier << "\t# Put virtual centers at indivual orbits, all common, merged by variance, to grid point, to cell center" << endl;
+    *output << "SawtoothStart\t" << config::SawtoothStart << "\t# Absolute value for smooth transition at cell border " << endl;
+    *output << "VectorPlane\t" << config::VectorPlane << "\t# Cut plane axis (x, y or z: 0,1,2) for two-dim current vector plot" << endl;
+    *output << "VectorCut\t" << config::VectorCut << "\t# Cut plane axis value" << endl;
+    *output << "AddGramSch\t" << config::UseAddGramSch << "\t# Additional GramSchmidtOrtogonalization to be safe" << endl;
+    *output << "Seed\t\t" << config::Seed << "\t# initial value for random seed for Psi coefficients" << endl;
+    *output << endl;
+    *output << "MaxOuterStep\t" << config::MaxOuterStep << "\t# number of MolecularDynamics/Structure optimization steps" << endl;
+    *output << "Deltat\t" << config::Deltat << "\t# time per MD step" << endl;
+    *output << "OutVisStep\t" << config::OutVisStep << "\t# Output visual data every ...th step" << endl;
+    *output << "OutSrcStep\t" << config::OutSrcStep << "\t# Output \"restart\" data every ..th step" << endl;
+    *output << "TargetTemp\t" << config::TargetTemp << "\t# Target temperature" << endl;
+    *output << "MaxPsiStep\t" << config::MaxPsiStep << "\t# number of Minimisation steps per state (0 - default)" << endl;
+    *output << "EpsWannier\t" << config::EpsWannier << "\t# tolerance value for spread minimisation of orbitals" << endl;
+    *output << endl;
+    *output << "# Values specifying when to stop" << endl;
+    *output << "MaxMinStep\t" << config::MaxMinStep << "\t# Maximum number of steps" << endl;
+    *output << "RelEpsTotalE\t" << config::RelEpsTotalEnergy << "\t# relative change in total energy" << endl;
+    *output << "RelEpsKineticE\t" << config::RelEpsKineticEnergy << "\t# relative change in kinetic energy" << endl;
+    *output << "MaxMinStopStep\t" << config::MaxMinStopStep << "\t# check every ..th steps" << endl;
+    *output << "MaxMinGapStopStep\t" << config::MaxMinGapStopStep << "\t# check every ..th steps" << endl;
+    *output << endl;
+    *output << "# Values specifying when to stop for INIT, otherwise same as above" << endl;
+    *output << "MaxInitMinStep\t" << config::MaxInitMinStep << "\t# Maximum number of steps" << endl;
+    *output << "InitRelEpsTotalE\t" << config::InitRelEpsTotalEnergy << "\t# relative change in total energy" << endl;
+    *output << "InitRelEpsKineticE\t" << config::InitRelEpsKineticEnergy << "\t# relative change in kinetic energy" << endl;
+    *output << "InitMaxMinStopStep\t" << config::InitMaxMinStopStep << "\t# check every ..th steps" << endl;
+    *output << "InitMaxMinGapStopStep\t" << config::InitMaxMinGapStopStep << "\t# check every ..th steps" << endl;
+    *output << endl;
+    *output << "BoxLength\t\t\t# (Length of a unit cell)" << endl;
+    *output << mol->cell_size[0] << "\t" << endl;
+    *output << mol->cell_size[1] << "\t" << mol->cell_size[2] << "\t" << endl;
+    *output << mol->cell_size[3] << "\t" << mol->cell_size[4] << "\t" << mol->cell_size[5] << "\t" << endl;
+    // FIXME
+    *output << endl;
+    *output << "ECut\t\t" << config::ECut << "\t# energy cutoff for discretization in Hartrees" << endl;
+    *output << "MaxLevel\t" << config::MaxLevel << "\t# number of different levels in the code, >=2" << endl;
+    *output << "Level0Factor\t" << config::Lev0Factor << "\t# factor by which node number increases from S to 0 level" << endl;
+    *output << "RiemannTensor\t" << config::RiemannTensor << "\t# (Use metric)" << endl;
+    switch (config::RiemannTensor) {
+      case 0: //UseNoRT
+        break;
+      case 1: // UseRT
+        *output << "RiemannLevel\t" << config::RiemannLevel << "\t# Number of Riemann Levels" << endl;
+        *output << "LevRFactor\t" << config::LevRFactor << "\t# factor by which node number increases from 0 to R level from" << endl;
+        break;
+    }
+    *output << "PsiType\t\t" << config::PsiType << "\t# 0 - doubly occupied, 1 - SpinUp,SpinDown" << endl;
+    // write out both types for easier changing afterwards
+  //  switch (PsiType) {
+  //    case 0:
+        *output << "MaxPsiDouble\t" << config::MaxPsiDouble << "\t# here: specifying both maximum number of SpinUp- and -Down-states" << endl;
+  //      break;
+  //    case 1:
+        *output << "PsiMaxNoUp\t" << config::PsiMaxNoUp << "\t# here: specifying maximum number of SpinUp-states" << endl;
+        *output << "PsiMaxNoDown\t" << config::PsiMaxNoDown << "\t# here: specifying maximum number of SpinDown-states" << endl;
+  //      break;
+  //  }
+    *output << "AddPsis\t\t" << config::AddPsis << "\t# Additional unoccupied Psis for bandgap determination" << endl;
+    *output << endl;
+    *output << "RCut\t\t" << config::RCut << "\t# R-cut for the ewald summation" << endl;
+    *output << "StructOpt\t" << config::StructOpt << "\t# Do structure optimization beforehand" << endl;
+    *output << "IsAngstroem\t" << config::IsAngstroem << "\t# 0 - Bohr, 1 - Angstroem" << endl;
+    *output << "RelativeCoord\t" << config::RelativeCoord << "\t# whether ion coordinates are relative (1) or absolute (0)" << endl;
+    *output << "MaxTypes\t" << mol->ElementCount <<  "\t# maximum number of different ion types" << endl;
+    *output << endl;
+    result = result && mol->Checkout(output);
+    if (mol->MDSteps <=1 )
+      result = result && mol->Output(output);
+    else
+      result = result && mol->OutputTrajectories(output);
+    output->close();
+    output->clear();
+    delete(output);
+    return result;
+  } else
+    return false;
 };
 
@@ -1254 +1254 @@
 bool config::SaveMPQC(const char *filename, molecule *mol) const
 {
-        int ElementNo = 0;
-        int AtomNo;
-        atom *Walker = NULL;
-        element *runner = NULL;
-        Vector *center = NULL;
-        ofstream *output = NULL;
-        stringstream *fname = NULL;
-
-        // first without hessian
-        fname = new stringstream;
-        *fname << filename << ".in";
-        output = new ofstream(fname->str().c_str(), ios::out);
-        *output << "% Created by MoleCuilder" << endl;
-        *output << "mpqc: (" << endl;
-        *output << "\tsavestate = no" << endl;
-        *output << "\tdo_gradient = yes" << endl;
-        *output << "\tmole<MBPT2>: (" << endl;
-        *output << "\t\tmaxiter = 200" << endl;
-        *output << "\t\tbasis = $:basis" << endl;
-        *output << "\t\tmolecule = $:molecule" << endl;
-        *output << "\t\treference<CLHF>: (" << endl;
-        *output << "\t\t\tbasis = $:basis" << endl;
-        *output << "\t\t\tmolecule = $:molecule" << endl;
-        *output << "\t\t)" << endl;
-        *output << "\t)" << endl;
-        *output << ")" << endl;
-        *output << "molecule<Molecule>: (" << endl;
-        *output << "\tunit = " << (IsAngstroem ? "angstrom" : "bohr" ) << endl;
-        *output << "\t{ atoms geometry } = {" << endl;
-        center = mol->DetermineCenterOfAll(output);
-        // output of atoms
-        runner = mol->elemente->start;
-        while (runner->next != mol->elemente->end) { // go through every element
-                runner = runner->next;
-                if (mol->ElementsInMolecule[runner->Z]) { // if this element got atoms
-                        ElementNo++;
-                        AtomNo = 0;
-                        Walker = mol->start;
-                        while (Walker->next != mol->end) { // go through every atom of this element
-                                Walker = Walker->next;
-                                if (Walker->type == runner) { // if this atom fits to element
-                                        AtomNo++;
-                                        *output << "\t\t" << Walker->type->symbol << " [ " << Walker->x.x[0]-center->x[0] << "\t" << Walker->x.x[1]-center->x[1] << "\t" << Walker->x.x[2]-center->x[2] << " ]" << endl;
-                                }
-                        }
-                }
-        }
-        delete(center);
-        *output << "\t}" << endl;
-        *output << ")" << endl;
-        *output << "basis<GaussianBasisSet>: (" << endl;
-        *output << "\tname = \""<< basis << "\"" << endl;
-        *output << "\tmolecule = $:molecule" << endl;
-        *output << ")" << endl;
-        output->close();
-        delete(output);
-        delete(fname);
-
-        // second with hessian
-        fname = new stringstream;
-        *fname << filename << ".hess.in";
-        output = new ofstream(fname->str().c_str(), ios::out);
-        *output << "% Created by MoleCuilder" << endl;
-        *output << "mpqc: (" << endl;
-        *output << "\tsavestate = no" << endl;
-        *output << "\tdo_gradient = yes" << endl;
-        *output << "\tmole<CLHF>: (" << endl;
-        *output << "\t\tmaxiter = 200" << endl;
-        *output << "\t\tbasis = $:basis" << endl;
-        *output << "\t\tmolecule = $:molecule" << endl;
-        *output << "\t)" << endl;
-        *output << "\tfreq<MolecularFrequencies>: (" << endl;
-        *output << "\t\tmolecule=$:molecule" << endl;
-        *output << "\t)" << endl;
-        *output << ")" << endl;
-        *output << "molecule<Molecule>: (" << endl;
-        *output << "\tunit = " << (IsAngstroem ? "angstrom" : "bohr" ) << endl;
-        *output << "\t{ atoms geometry } = {" << endl;
-        center = mol->DetermineCenterOfAll(output);
-        // output of atoms
-        runner = mol->elemente->start;
-        while (runner->next != mol->elemente->end) { // go through every element
-                runner = runner->next;
-                if (mol->ElementsInMolecule[runner->Z]) { // if this element got atoms
-                        ElementNo++;
-                        AtomNo = 0;
-                        Walker = mol->start;
-                        while (Walker->next != mol->end) { // go through every atom of this element
-                                Walker = Walker->next;
-                                if (Walker->type == runner) { // if this atom fits to element
-                                        AtomNo++;
-                                        *output << "\t\t" << Walker->type->symbol << " [ " << Walker->x.x[0]-center->x[0] << "\t" << Walker->x.x[1]-center->x[1] << "\t" << Walker->x.x[2]-center->x[2] << " ]" << endl;
-                                }
-                        }
-                }
-        }
-        delete(center);
-        *output << "\t}" << endl;
-        *output << ")" << endl;
-        *output << "basis<GaussianBasisSet>: (" << endl;
-        *output << "\tname = \"3-21G\"" << endl;
-        *output << "\tmolecule = $:molecule" << endl;
-        *output << ")" << endl;
-        output->close();
-        delete(output);
-        delete(fname);
-
-        return true;
+  int ElementNo = 0;
+  int AtomNo;
+  atom *Walker = NULL;
+  element *runner = NULL;
+  Vector *center = NULL;
+  ofstream *output = NULL;
+  stringstream *fname = NULL;
+
+  // first without hessian
+  fname = new stringstream;
+  *fname << filename << ".in";
+  output = new ofstream(fname->str().c_str(), ios::out);
+  *output << "% Created by MoleCuilder" << endl;
+  *output << "mpqc: (" << endl;
+  *output << "\tsavestate = no" << endl;
+  *output << "\tdo_gradient = yes" << endl;
+  *output << "\tmole<MBPT2>: (" << endl;
+  *output << "\t\tmaxiter = 200" << endl;
+  *output << "\t\tbasis = $:basis" << endl;
+  *output << "\t\tmolecule = $:molecule" << endl;
+  *output << "\t\treference<CLHF>: (" << endl;
+  *output << "\t\t\tbasis = $:basis" << endl;
+  *output << "\t\t\tmolecule = $:molecule" << endl;
+  *output << "\t\t)" << endl;
+  *output << "\t)" << endl;
+  *output << ")" << endl;
+  *output << "molecule<Molecule>: (" << endl;
+  *output << "\tunit = " << (IsAngstroem ? "angstrom" : "bohr" ) << endl;
+  *output << "\t{ atoms geometry } = {" << endl;
+  center = mol->DetermineCenterOfAll(output);
+  // output of atoms
+  runner = mol->elemente->start;
+  while (runner->next != mol->elemente->end) { // go through every element
+    runner = runner->next;
+    if (mol->ElementsInMolecule[runner->Z]) { // if this element got atoms
+      ElementNo++;
+      AtomNo = 0;
+      Walker = mol->start;
+      while (Walker->next != mol->end) { // go through every atom of this element
+        Walker = Walker->next;
+        if (Walker->type == runner) { // if this atom fits to element
+          AtomNo++;
+          *output << "\t\t" << Walker->type->symbol << " [ " << Walker->x.x[0]-center->x[0] << "\t" << Walker->x.x[1]-center->x[1] << "\t" << Walker->x.x[2]-center->x[2] << " ]" << endl;
+        }
+      }
+    }
+  }
+  delete(center);
+  *output << "\t}" << endl;
+  *output << ")" << endl;
+  *output << "basis<GaussianBasisSet>: (" << endl;
+  *output << "\tname = \""<< basis << "\"" << endl;
+  *output << "\tmolecule = $:molecule" << endl;
+  *output << ")" << endl;
+  output->close();
+  delete(output);
+  delete(fname);
+
+  // second with hessian
+  fname = new stringstream;
+  *fname << filename << ".hess.in";
+  output = new ofstream(fname->str().c_str(), ios::out);
+  *output << "% Created by MoleCuilder" << endl;
+  *output << "mpqc: (" << endl;
+  *output << "\tsavestate = no" << endl;
+  *output << "\tdo_gradient = yes" << endl;
+  *output << "\tmole<CLHF>: (" << endl;
+  *output << "\t\tmaxiter = 200" << endl;
+  *output << "\t\tbasis = $:basis" << endl;
+  *output << "\t\tmolecule = $:molecule" << endl;
+  *output << "\t)" << endl;
+  *output << "\tfreq<MolecularFrequencies>: (" << endl;
+  *output << "\t\tmolecule=$:molecule" << endl;
+  *output << "\t)" << endl;
+  *output << ")" << endl;
+  *output << "molecule<Molecule>: (" << endl;
+  *output << "\tunit = " << (IsAngstroem ? "angstrom" : "bohr" ) << endl;
+  *output << "\t{ atoms geometry } = {" << endl;
+  center = mol->DetermineCenterOfAll(output);
+  // output of atoms
+  runner = mol->elemente->start;
+  while (runner->next != mol->elemente->end) { // go through every element
+    runner = runner->next;
+    if (mol->ElementsInMolecule[runner->Z]) { // if this element got atoms
+      ElementNo++;
+      AtomNo = 0;
+      Walker = mol->start;
+      while (Walker->next != mol->end) { // go through every atom of this element
+        Walker = Walker->next;
+        if (Walker->type == runner) { // if this atom fits to element
+          AtomNo++;
+          *output << "\t\t" << Walker->type->symbol << " [ " << Walker->x.x[0]-center->x[0] << "\t" << Walker->x.x[1]-center->x[1] << "\t" << Walker->x.x[2]-center->x[2] << " ]" << endl;
+        }
+      }
+    }
+  }
+  delete(center);
+  *output << "\t}" << endl;
+  *output << ")" << endl;
+  *output << "basis<GaussianBasisSet>: (" << endl;
+  *output << "\tname = \"3-21G\"" << endl;
+  *output << "\tmolecule = $:molecule" << endl;
+  *output << ")" << endl;
+  output->close();
+  delete(output);
+  delete(fname);
+
+  return true;
 };
 
@@ -1373 +1373 @@
  * \param name Name of value in file (at least 3 chars!)
  * \param sequential 1 - do not reset file pointer to begin of file, 0 - set to beginning
- *                              (if file is sequentially parsed this can be way faster! However, beware of multiple values per line, as whole line is read -
- *                               best approach: 0 0 0 1 (not resetted only on last value of line) - and of yth, which is now
- *                               counted from this unresetted position!)
+ *        (if file is sequentially parsed this can be way faster! However, beware of multiple values per line, as whole line is read -
+ *         best approach: 0 0 0 1 (not resetted only on last value of line) - and of yth, which is now
+ *         counted from this unresetted position!)
  * \param xth Which among a number of parameters it is (in grid case it's row number as grid is read as a whole!)
  * \param yth In grid case specifying column number, otherwise the yth \a name matching line
@@ -1386 +1386 @@
  */
 int config::ParseForParameter(int verbose, ifstream *file, const char *name, int sequential, int const xth, int const yth, int type, void *value, int repetition, int critical) {
-        int i,j;        // loop variables
-        int length = 0, maxlength = -1;
-        long file_position = file->tellg(); // mark current position
-        char *dummy1, *dummy, *free_dummy;              // pointers in the line that is read in per step
-        dummy1 = free_dummy = (char *) Malloc(256 * sizeof(char), "config::ParseForParameter: *free_dummy");
-
-        //fprintf(stderr,"Parsing for %s\n",name);
-        if (repetition == 0)
-                //Error(SomeError, "ParseForParameter(): argument repetition must not be 0!");
-                return 0;
-
-        int line = 0; // marks line where parameter was found
-        int found = (type >= grid) ? 0 : (-yth + 1);    // marks if yth parameter name was found
-        while((found != repetition)) {
-                dummy1 = dummy = free_dummy;
-                do {
-                        file->getline(dummy1, 256); // Read the whole line
-                        if (file->eof()) {
-                                if ((critical) && (found == 0)) {
-                                        Free((void **)&free_dummy, "config::ParseForParameter: *free_dummy");
-                                        //Error(InitReading, name);
-                                        fprintf(stderr,"Error:InitReading, critical %s not found\n", name);
-                                        exit(255);
-                                } else {
-                                        //if (!sequential)
-                                        file->clear();
-                                        file->seekg(file_position, ios::beg);   // rewind to start position
-                                        Free((void **)&free_dummy, "config::ParseForParameter: *free_dummy");
-                                        return 0;
-                                }
-                        }
-                        line++;
-                } while (dummy != NULL && dummy1 != NULL && ((dummy1[0] == '#') || (dummy1[0] == '\0'))); // skip commentary and empty lines
-
-                // C++ getline removes newline at end, thus re-add
-                if ((dummy1 != NULL) && (strchr(dummy1,'\n') == NULL)) {
-                        i = strlen(dummy1);
-                        dummy1[i] = '\n';
-                        dummy1[i+1] = '\0';
-                }
-                //fprintf(stderr,"line %i ends at %i, newline at %i\n", line, strlen(dummy1), strchr(dummy1,'\n')-free_dummy);
-
-                if (dummy1 == NULL) {
-                        if (verbose) fprintf(stderr,"Error reading line %i\n",line);
-                } else {
-                        //fprintf(stderr,"Now parsing the line %i: %s\n", line, dummy1);
-                }
-                // Seek for possible end of keyword on line if given ...
-                if (name != NULL) {
-                        dummy = strchr(dummy1,'\t');    // set dummy on first tab or space which ever's nearer
-                        if (dummy == NULL) {
-                                dummy = strchr(dummy1, ' ');    // if not found seek for space
-                                while ((dummy != NULL) && ((*dummy == '\t') || (*dummy == ' ')))                // skip some more tabs and spaces if necessary
-                                        dummy++;
-                        }
-                        if (dummy == NULL) {
-                                dummy = strchr(dummy1, '\n'); // set on line end then (whole line = keyword)
-                                //fprintf(stderr,"Error: Cannot find tabs or spaces on line %i in search for %s\n", line, name);
-                                //Free((void **)&free_dummy);
-                                //Error(FileOpenParams, NULL);
-                        } else {
-                                //fprintf(stderr,"found tab at %i\n",(char *)dummy-(char *)dummy1);
-                        }
-                } else dummy = dummy1;
-                // ... and check if it is the keyword!
-                //fprintf(stderr,"name %p, dummy %i/%c, dummy1 %i/%c, strlen(name) %i\n", &name, dummy, *dummy, dummy1, *dummy1, strlen(name));
-                if ((name == NULL) || (((dummy-dummy1 >= 3) && (strncmp(dummy1, name, strlen(name)) == 0)) && ((unsigned int)(dummy-dummy1) == strlen(name)))) {
-                        found++; // found the parameter!
-                        //fprintf(stderr,"found %s at line %i between %i and %i\n", name, line, dummy1, dummy);
-
-                        if (found == repetition) {
-                                for (i=0;i<xth;i++) { // i = rows
-                                        if (type >= grid) {
-                                                // grid structure means that grid starts on the next line, not right after keyword
-                                                dummy1 = dummy = free_dummy;
-                                                do {
-                                                        file->getline(dummy1, 256); // Read the whole line, skip commentary and empty ones
-                                                        if (file->eof()) {
-                                                                if ((critical) && (found == 0)) {
-                                                                        Free((void **)&free_dummy, "config::ParseForParameter: *free_dummy");
-                                                                        //Error(InitReading, name);
-                                                                        fprintf(stderr,"Error:InitReading, critical %s not found\n", name);
-                                                                        exit(255);
-                                                                } else {
-                                                                        //if (!sequential)
-                                                                        file->clear();
-                                                                        file->seekg(file_position, ios::beg);   // rewind to start position
-                                                                        Free((void **)&free_dummy, "config::ParseForParameter: *free_dummy");
-                                                                        return 0;
-                                                                }
-                                                        }
-                                                        line++;
-                                                } while ((dummy1[0] == '#') || (dummy1[0] == '\n'));
-                                                if (dummy1 == NULL){
-                                                        if (verbose) fprintf(stderr,"Error reading line %i\n", line);
-                                                } else {
-                                                        //fprintf(stderr,"Reading next line %i: %s\n", line, dummy1);
-                                                }
-                                        } else { // simple int, strings or doubles start in the same line
-                                                while ((*dummy == '\t') || (*dummy == ' '))      // skip interjacent tabs and spaces
-                                                        dummy++;
-                                        }
-                                        // C++ getline removes newline at end, thus re-add
-                                        if ((dummy1 != NULL) && (strchr(dummy1,'\n') == NULL)) {
-                                                j = strlen(dummy1);
-                                                dummy1[j] = '\n';
-                                                dummy1[j+1] = '\0';
-                                        }
-
-                                        int start = (type >= grid) ? 0 : yth-1 ;
-                                        for (j=start;j<yth;j++) { // j = columns
-                                                // check for lower triangular area and upper triangular area
-                                                if ( ((i > j) && (type == upper_trigrid)) || ((j > i) && (type == lower_trigrid))) {
-                                                        *((double *)value) = 0.0;
-                                                        fprintf(stderr,"%f\t",*((double *)value));
-                                                        value = (void *)((long)value + sizeof(double));
-                                                        //value += sizeof(double);
-                                                } else {
-                                                        // otherwise we must skip all interjacent tabs and spaces and find next value
-                                                        dummy1 = dummy;
-                                                        dummy = strchr(dummy1, '\t'); // seek for tab or space
-                                                        if (dummy == NULL)
-                                                                dummy = strchr(dummy1, ' ');    // if not found seek for space
-                                                        if (dummy == NULL) { // if still zero returned ...
-                                                                dummy = strchr(dummy1, '\n');   // ... at line end then
-                                                                if ((j < yth-1) && (type < 4)) {        // check if xth value or not yet
-                                                                        if (critical) {
-                                                                                if (verbose) fprintf(stderr,"Error: EoL at %i and still missing %i value(s) for parameter %s\n", line, yth-j, name);
-                                                                                Free((void **)&free_dummy, "config::ParseForParameter: *free_dummy");
-                                                                                //return 0;
-                                                                                exit(255);
-                                                                                //Error(FileOpenParams, NULL);
-                                                                        } else {
-                                                                                //if (!sequential)
-                                                                                file->clear();
-                                                                                file->seekg(file_position, ios::beg);   // rewind to start position
-                                                                                Free((void **)&free_dummy, "config::ParseForParameter: *free_dummy");
-                                                                                return 0;
-                                                                        }
-                                                                }
-                                                        } else {
-                                                                //fprintf(stderr,"found tab at %i\n",(char *)dummy-(char *)free_dummy);
-                                                        }
-                                                        if (*dummy1 == '#') {
-                                                                // found comment, skipping rest of line
-                                                                //if (verbose) fprintf(stderr,"Error: '#' at %i and still missing %i value(s) for parameter %s\n", line, yth-j, name);
-                                                                if (!sequential) { // here we need it!
-                                                                        file->seekg(file_position, ios::beg);   // rewind to start position
-                                                                }
-                                                                Free((void **)&free_dummy, "config::ParseForParameter: *free_dummy");
-                                                                return 0;
-                                                        }
-                                                        //fprintf(stderr,"value from %i to %i\n",(char *)dummy1-(char *)free_dummy,(char *)dummy-(char *)free_dummy);
-                                                        switch(type) {
-                                                                case (row_int):
-                                                                        *((int *)value) = atoi(dummy1);
-                                                                        if ((verbose) && (i==0) && (j==0)) fprintf(stderr,"%s = ", name);
-                                                                        if (verbose) fprintf(stderr,"%i\t",*((int *)value));
-                                                                                value = (void *)((long)value + sizeof(int));
-                                                                                //value += sizeof(int);
-                                                                        break;
-                                                                case(row_double):
-                                                                case(grid):
-                                                                case(lower_trigrid):
-                                                                case(upper_trigrid):
-                                                                        *((double *)value) = atof(dummy1);
-                                                                        if ((verbose) && (i==0) && (j==0)) fprintf(stderr,"%s = ", name);
-                                                                        if (verbose) fprintf(stderr,"%lg\t",*((double *)value));
-                                                                        value = (void *)((long)value + sizeof(double));
-                                                                        //value += sizeof(double);
-                                                                        break;
-                                                                case(double_type):
-                                                                        *((double *)value) = atof(dummy1);
-                                                                        if ((verbose) && (i == xth-1)) fprintf(stderr,"%s = %lg\n", name, *((double *) value));
-                                                                        //value += sizeof(double);
-                                                                        break;
-                                                                case(int_type):
-                                                                        *((int *)value) = atoi(dummy1);
-                                                                        if ((verbose) && (i == xth-1)) fprintf(stderr,"%s = %i\n", name, *((int *) value));
-                                                                        //value += sizeof(int);
-                                                                        break;
-                                                                default:
-                                                                case(string_type):
-                                                                        if (value != NULL) {
-                                                                                //if (maxlength == -1) maxlength = strlen((char *)value); // get maximum size of string array
-                                                                                maxlength = MAXSTRINGSIZE;
-                                                                                length = maxlength > (dummy-dummy1) ? (dummy-dummy1) : maxlength; // cap at maximum
-                                                                                strncpy((char *)value, dummy1, length); // copy as much
-                                                                                ((char *)value)[length] = '\0'; // and set end marker
-                                                                                if ((verbose) && (i == xth-1)) fprintf(stderr,"%s is '%s' (%i chars)\n",name,((char *) value), length);
-                                                                                //value += sizeof(char);
-                                                                        } else {
-                                                                        }
-                                                                break;
-                                                        }
-                                                }
-                                                while (*dummy == '\t')
-                                                        dummy++;
-                                        }
-                                }
-                        }
-                }
-        }
-        if ((type >= row_int) && (verbose)) fprintf(stderr,"\n");
-        Free((void **)&free_dummy, "config::ParseForParameter: *free_dummy");
-        if (!sequential) {
-                file->clear();
-                file->seekg(file_position, ios::beg);   // rewind to start position
-        }
-        //fprintf(stderr, "End of Parsing\n\n");
-
-        return (found); // true if found, false if not
+  int i,j;  // loop variables
+  int length = 0, maxlength = -1;
+  long file_position = file->tellg(); // mark current position
+  char *dummy1, *dummy, *free_dummy;    // pointers in the line that is read in per step
+  dummy1 = free_dummy = (char *) Malloc(256 * sizeof(char), "config::ParseForParameter: *free_dummy");
+
+  //fprintf(stderr,"Parsing for %s\n",name);
+  if (repetition == 0)
+    //Error(SomeError, "ParseForParameter(): argument repetition must not be 0!");
+    return 0;
+
+  int line = 0; // marks line where parameter was found
+  int found = (type >= grid) ? 0 : (-yth + 1);  // marks if yth parameter name was found
+  while((found != repetition)) {
+    dummy1 = dummy = free_dummy;
+    do {
+      file->getline(dummy1, 256); // Read the whole line
+      if (file->eof()) {
+        if ((critical) && (found == 0)) {
+          Free((void **)&free_dummy, "config::ParseForParameter: *free_dummy");
+          //Error(InitReading, name);
+          fprintf(stderr,"Error:InitReading, critical %s not found\n", name);
+          exit(255);
+        } else {
+          //if (!sequential)
+          file->clear();
+          file->seekg(file_position, ios::beg);  // rewind to start position
+          Free((void **)&free_dummy, "config::ParseForParameter: *free_dummy");
+          return 0;
+        }
+      }
+      line++;
+    } while (dummy != NULL && dummy1 != NULL && ((dummy1[0] == '#') || (dummy1[0] == '\0'))); // skip commentary and empty lines
+
+    // C++ getline removes newline at end, thus re-add
+    if ((dummy1 != NULL) && (strchr(dummy1,'\n') == NULL)) {
+      i = strlen(dummy1);
+      dummy1[i] = '\n';
+      dummy1[i+1] = '\0';
+    }
+    //fprintf(stderr,"line %i ends at %i, newline at %i\n", line, strlen(dummy1), strchr(dummy1,'\n')-free_dummy);
+
+    if (dummy1 == NULL) {
+      if (verbose) fprintf(stderr,"Error reading line %i\n",line);
+    } else {
+      //fprintf(stderr,"Now parsing the line %i: %s\n", line, dummy1);
+    }
+    // Seek for possible end of keyword on line if given ...
+    if (name != NULL) {
+      dummy = strchr(dummy1,'\t');  // set dummy on first tab or space which ever's nearer
+      if (dummy == NULL) {
+        dummy = strchr(dummy1, ' ');  // if not found seek for space
+        while ((dummy != NULL) && ((*dummy == '\t') || (*dummy == ' ')))    // skip some more tabs and spaces if necessary
+          dummy++;
+      }
+      if (dummy == NULL) {
+        dummy = strchr(dummy1, '\n'); // set on line end then (whole line = keyword)
+        //fprintf(stderr,"Error: Cannot find tabs or spaces on line %i in search for %s\n", line, name);
+        //Free((void **)&free_dummy);
+        //Error(FileOpenParams, NULL);
+      } else {
+        //fprintf(stderr,"found tab at %i\n",(char *)dummy-(char *)dummy1);
+      }
+    } else dummy = dummy1;
+    // ... and check if it is the keyword!
+    //fprintf(stderr,"name %p, dummy %i/%c, dummy1 %i/%c, strlen(name) %i\n", &name, dummy, *dummy, dummy1, *dummy1, strlen(name));
+    if ((name == NULL) || (((dummy-dummy1 >= 3) && (strncmp(dummy1, name, strlen(name)) == 0)) && ((unsigned int)(dummy-dummy1) == strlen(name)))) {
+      found++; // found the parameter!
+      //fprintf(stderr,"found %s at line %i between %i and %i\n", name, line, dummy1, dummy);
+
+      if (found == repetition) {
+        for (i=0;i<xth;i++) { // i = rows
+          if (type >= grid) {
+            // grid structure means that grid starts on the next line, not right after keyword
+            dummy1 = dummy = free_dummy;
+            do {
+              file->getline(dummy1, 256); // Read the whole line, skip commentary and empty ones
+              if (file->eof()) {
+                if ((critical) && (found == 0)) {
+                  Free((void **)&free_dummy, "config::ParseForParameter: *free_dummy");
+                  //Error(InitReading, name);
+                  fprintf(stderr,"Error:InitReading, critical %s not found\n", name);
+                  exit(255);
+                } else {
+                  //if (!sequential)
+                  file->clear();
+                  file->seekg(file_position, ios::beg);  // rewind to start position
+                  Free((void **)&free_dummy, "config::ParseForParameter: *free_dummy");
+                  return 0;
+                }
+              }
+              line++;
+            } while ((dummy1[0] == '#') || (dummy1[0] == '\n'));
+            if (dummy1 == NULL){
+              if (verbose) fprintf(stderr,"Error reading line %i\n", line);
+            } else {
+              //fprintf(stderr,"Reading next line %i: %s\n", line, dummy1);
+            }
+          } else { // simple int, strings or doubles start in the same line
+            while ((*dummy == '\t') || (*dummy == ' '))   // skip interjacent tabs and spaces
+              dummy++;
+          }
+          // C++ getline removes newline at end, thus re-add
+          if ((dummy1 != NULL) && (strchr(dummy1,'\n') == NULL)) {
+            j = strlen(dummy1);
+            dummy1[j] = '\n';
+            dummy1[j+1] = '\0';
+          }
+
+          int start = (type >= grid) ? 0 : yth-1 ;
+          for (j=start;j<yth;j++) { // j = columns
+            // check for lower triangular area and upper triangular area
+            if ( ((i > j) && (type == upper_trigrid)) || ((j > i) && (type == lower_trigrid))) {
+              *((double *)value) = 0.0;
+              fprintf(stderr,"%f\t",*((double *)value));
+              value = (void *)((long)value + sizeof(double));
+              //value += sizeof(double);
+            } else {
+              // otherwise we must skip all interjacent tabs and spaces and find next value
+              dummy1 = dummy;
+              dummy = strchr(dummy1, '\t'); // seek for tab or space
+              if (dummy == NULL)
+                dummy = strchr(dummy1, ' ');  // if not found seek for space
+              if (dummy == NULL) { // if still zero returned ...
+                dummy = strchr(dummy1, '\n');  // ... at line end then
+                if ((j < yth-1) && (type < 4)) {  // check if xth value or not yet
+                  if (critical) {
+                    if (verbose) fprintf(stderr,"Error: EoL at %i and still missing %i value(s) for parameter %s\n", line, yth-j, name);
+                    Free((void **)&free_dummy, "config::ParseForParameter: *free_dummy");
+                    //return 0;
+                    exit(255);
+                    //Error(FileOpenParams, NULL);
+                  } else {
+                    //if (!sequential)
+                    file->clear();
+                    file->seekg(file_position, ios::beg);  // rewind to start position
+                    Free((void **)&free_dummy, "config::ParseForParameter: *free_dummy");
+                    return 0;
+                  }
+                }
+              } else {
+                //fprintf(stderr,"found tab at %i\n",(char *)dummy-(char *)free_dummy);
+              }
+              if (*dummy1 == '#') {
+                // found comment, skipping rest of line
+                //if (verbose) fprintf(stderr,"Error: '#' at %i and still missing %i value(s) for parameter %s\n", line, yth-j, name);
+                if (!sequential) { // here we need it!
+                  file->seekg(file_position, ios::beg);  // rewind to start position
+                }
+                Free((void **)&free_dummy, "config::ParseForParameter: *free_dummy");
+                return 0;
+              }
+              //fprintf(stderr,"value from %i to %i\n",(char *)dummy1-(char *)free_dummy,(char *)dummy-(char *)free_dummy);
+              switch(type) {
+                case (row_int):
+                  *((int *)value) = atoi(dummy1);
+                  if ((verbose) && (i==0) && (j==0)) fprintf(stderr,"%s = ", name);
+                  if (verbose) fprintf(stderr,"%i\t",*((int *)value));
+                    value = (void *)((long)value + sizeof(int));
+                    //value += sizeof(int);
+                  break;
+                case(row_double):
+                case(grid):
+                case(lower_trigrid):
+                case(upper_trigrid):
+                  *((double *)value) = atof(dummy1);
+                  if ((verbose) && (i==0) && (j==0)) fprintf(stderr,"%s = ", name);
+                  if (verbose) fprintf(stderr,"%lg\t",*((double *)value));
+                  value = (void *)((long)value + sizeof(double));
+                  //value += sizeof(double);
+                  break;
+                case(double_type):
+                  *((double *)value) = atof(dummy1);
+                  if ((verbose) && (i == xth-1)) fprintf(stderr,"%s = %lg\n", name, *((double *) value));
+                  //value += sizeof(double);
+                  break;
+                case(int_type):
+                  *((int *)value) = atoi(dummy1);
+                  if ((verbose) && (i == xth-1)) fprintf(stderr,"%s = %i\n", name, *((int *) value));
+                  //value += sizeof(int);
+                  break;
+                default:
+                case(string_type):
+                  if (value != NULL) {
+                    //if (maxlength == -1) maxlength = strlen((char *)value); // get maximum size of string array
+                    maxlength = MAXSTRINGSIZE;
+                    length = maxlength > (dummy-dummy1) ? (dummy-dummy1) : maxlength; // cap at maximum
+                    strncpy((char *)value, dummy1, length);  // copy as much
+                    ((char *)value)[length] = '\0';  // and set end marker
+                    if ((verbose) && (i == xth-1)) fprintf(stderr,"%s is '%s' (%i chars)\n",name,((char *) value), length);
+                    //value += sizeof(char);
+                  } else {
+                  }
+                break;
+              }
+            }
+            while (*dummy == '\t')
+              dummy++;
+          }
+        }
+      }
+    }
+  }
+  if ((type >= row_int) && (verbose)) fprintf(stderr,"\n");
+  Free((void **)&free_dummy, "config::ParseForParameter: *free_dummy");
+  if (!sequential) {
+    file->clear();
+    file->seekg(file_position, ios::beg);  // rewind to start position
+  }
+  //fprintf(stderr, "End of Parsing\n\n");
+
+  return (found); // true if found, false if not
 }