diff --git a/Common/engine_interface_base.h b/Common/engine_interface_base.h
index 15be488..8b101a6 100644
--- a/Common/engine_interface_base.h
+++ b/Common/engine_interface_base.h
@@ -22,9 +22,9 @@
 
 class Operator_Base;
 
-//! This is the abstact base for all Engine Interface classes.
+//! This is the abstract base for all Engine Interface classes.
 /*!
-	 This is the abstact base for all Engine Interface classes. It will provide unified access to the field information of the corresponding engine.
+	 This is the abstract base for all Engine Interface classes. It will provide unified access to the field information of the corresponding engine.
 	 All processing methods should only access this base class.
 */
 class Engine_Interface_Base
diff --git a/Common/operator_base.h b/Common/operator_base.h
index cfbd40b..5c98305 100644
--- a/Common/operator_base.h
+++ b/Common/operator_base.h
@@ -71,7 +71,7 @@ public:
 	//! Get the length of an FDTD edge (unit is meter).
 	virtual double GetEdgeLength(int ny, const unsigned int pos[3], bool dualMesh = false) const =0;
 
-	//! Get the area around an edge for a given direction \a n and a given mesh posisition \a pos
+	//! Get the area around an edge for a given direction \a n and a given mesh position \a pos
 	/*!
 		This will return the area around an edge with a given direction, measured at the middle of the edge.
 		In a cartesian mesh this is equal to the NodeArea, may be different in other coordinate systems.
@@ -81,7 +81,7 @@ public:
 	//! Get the volume of an FDTD cell
 	virtual double GetCellVolume(const unsigned int pos[3], bool dualMesh = false) const =0;
 
-	//! Snap the given coodinates to mesh indices, return box dimension
+	//! Snap the given coordinates to mesh indices, return box dimension
 	virtual bool SnapToMesh(const double* coord, unsigned int* uicoord, bool dualMesh=false, bool fullMesh=false, bool* inside=NULL) const =0;
 
 	//! Snap a given box to the operator mesh, uiStart will be always <= uiStop
diff --git a/Common/processfields.h b/Common/processfields.h
index d954484..1f9e7f5 100644
--- a/Common/processfields.h
+++ b/Common/processfields.h
@@ -60,9 +60,9 @@ public:
 
 	//! Define the Dump-Mode
 	void SetDumpMode(Engine_Interface_Base::InterpolationType mode);
-	//! This methode will dump all fields on a main cell node using 2 E-field and 4 H-fields per direction.
+	//! This method will dump all fields on a main cell node using 2 E-field and 4 H-fields per direction.
 	void SetDumpMode2Node() {SetDumpMode(Engine_Interface_Base::NODE_INTERPOLATE);}
-	//! This methode will dump all fields in the center of a main cell (dual-node) using 4 E-field and 2 H-fields per direction.
+	//! This method will dump all fields in the center of a main cell (dual-node) using 4 E-field and 2 H-fields per direction.
 	void SetDumpMode2Cell() {SetDumpMode(Engine_Interface_Base::CELL_INTERPOLATE);}
 
 	//! Set dump type: 0 for E-fields, 1 for H-fields, 2 for D-fields, 3 for B-fields, 4 for J-fields, etc...
diff --git a/Common/processing.h b/Common/processing.h
index ab20f78..38d507b 100644
--- a/Common/processing.h
+++ b/Common/processing.h
@@ -102,7 +102,7 @@ public:
 	//! Set the dump precision
 	void SetPrecision(unsigned int val) {m_precision = val;}
 
-	//! Dump probe geometry to file (will obay main or dual mesh property)
+	//! Dump probe geometry to file (will obey main or dual mesh property)
 	virtual void DumpBox2File(std::string vtkfilenameprefix) const {DumpBox2File(vtkfilenameprefix,m_dualMesh);}
 
 	//! Dump probe geometry to file
@@ -184,7 +184,7 @@ public:
 
 	void Reset();
 
-	//! Deletes all given processing's, can be helpful, but use carefull!!!
+	//! Deletes all given processing's, can be helpful, but use carefully!!!
 	void DeleteAll();
 
 	//! Invoke PreProcess() on all Processings.
diff --git a/Common/processmodematch.cpp b/Common/processmodematch.cpp
index 85d2d58..e7e4b5d 100644
--- a/Common/processmodematch.cpp
+++ b/Common/processmodematch.cpp
@@ -74,7 +74,7 @@ void ProcessModeMatch::InitProcess()
 
 	if (m_Eng_Interface==NULL)
 	{
-		cerr << "ProcessModeMatch::InitProcess: Error, Engine_Interface is NULL, abort mode mathcing..." << endl;
+		cerr << "ProcessModeMatch::InitProcess: Error, Engine_Interface is NULL, abort mode matching..." << endl;
 		Enabled=false;
 		return;
 	}
@@ -123,7 +123,7 @@ void ProcessModeMatch::InitProcess()
 		int res = m_ModeParser[n]->Parse(m_ModeFunction[ny], "x,y,z,rho,a,r,t");
 		if (res >= 0)
 		{
-			cerr << "ProcessModeMatch::InitProcess(): Warning, an error occured parsing the mode matching function (see below) ..." << endl;
+			cerr << "ProcessModeMatch::InitProcess(): Warning, an error occurred parsing the mode matching function (see below) ..." << endl;
 			cerr << m_ModeFunction[ny] << "\n" << string(res, ' ') << "^\n" << m_ModeParser[n]->ErrorMsg() << "\n";
 			SetEnable(false);
 			Reset();
@@ -138,7 +138,7 @@ void ProcessModeMatch::InitProcess()
 	bool dualMesh = m_ModeFieldType==1;
 	unsigned int pos[3] = {0,0,0};
 	double discLine[3] = {0,0,0};
-	double gridDelta = 1; // 1 -> mode-matching function is definied in drawing units...
+	double gridDelta = 1; // 1 -> mode-matching function is defined in drawing units...
 	double var[7];
 	pos[m_ny] = start[m_ny];
 	discLine[m_ny] = Op->GetDiscLine(m_ny,pos[m_ny],dualMesh);
diff --git a/Common/processmodematch.h b/Common/processmodematch.h
index 0bb03b7..6105886 100644
--- a/Common/processmodematch.h
+++ b/Common/processmodematch.h
@@ -24,7 +24,7 @@ class CSFunctionParser;
 
 //! Processing class to match a mode to a given analytic function and return the integral value.
 /*!
-  The analytric function has to be definied in drawing units.
+  The analytric function has to be defined in drawing units.
   It will return the integral value and the mode purity as a secondary value.
   */
 class ProcessModeMatch : public ProcessIntegral
diff --git a/FDTD/excitation.h b/FDTD/excitation.h
index b43ad28..6f09c9a 100644
--- a/FDTD/excitation.h
+++ b/FDTD/excitation.h
@@ -64,7 +64,7 @@ public:
 	//! Get the length of the excitation signal
 	unsigned int GetLength() const {return Length;}
 
-	//! Get the max frequeny excited by this signal
+	//! Get the max frequency excited by this signal
 	double GetMaxFrequency() const {return m_f_max;}
 
 	//! Get the frequency of interest
diff --git a/FDTD/extensions/engine_ext_cylindermultigrid.cpp b/FDTD/extensions/engine_ext_cylindermultigrid.cpp
index eca2c9e..a9b23cb 100644
--- a/FDTD/extensions/engine_ext_cylindermultigrid.cpp
+++ b/FDTD/extensions/engine_ext_cylindermultigrid.cpp
@@ -55,7 +55,7 @@ void Engine_Ext_CylinderMultiGrid::DoPreVoltageUpdates()
 	if (!m_IsBase)
 	{
 		//cerr << "child: volt wait on base " << endl;
-		m_WaitOnBase->wait(); //wait on base to finisch current sync and/or to finisch voltage updates, than start child voltage updates
+		m_WaitOnBase->wait(); //wait on base to finish current sync and/or to finish voltage updates, than start child voltage updates
 	}
 }
 
@@ -69,14 +69,14 @@ void Engine_Ext_CylinderMultiGrid::Apply2Voltages()
 	if (m_IsBase)
 	{
 		m_WaitOnBase->wait();	//base voltage updates are done, tell child to start its voltage updates
-		m_WaitOnChild->wait();	//wait for child to finisch its updates
-		SyncVoltages();			//child is finisch, run sync and go to current updates next
+		m_WaitOnChild->wait();	//wait for child to finish its updates
+		SyncVoltages();			//child is finish, run sync and go to current updates next
 		m_WaitOnSync->wait();	//sync is done... move on and tell child to move on...
 	}
 	else
 	{
 		m_WaitOnChild->wait();	//child is finished voltage updates, will tell base to run sync
-		m_WaitOnSync->wait();	//wait for base to finisch sync before going to wait for current updates
+		m_WaitOnSync->wait();	//wait for base to finish sync before going to wait for current updates
 	}
 }
 
@@ -125,7 +125,7 @@ void Engine_Ext_CylinderMultiGrid::DoPreCurrentUpdates()
 	if (!m_IsBase)
 	{
 		//cerr << "child: curr wait on base " << endl;
-		m_WaitOnBase->wait(); //wait on base to finisch voltage sync and current updates, than start child current updates
+		m_WaitOnBase->wait(); //wait on base to finish voltage sync and current updates, than start child current updates
 	}
 }
 
@@ -140,14 +140,14 @@ void Engine_Ext_CylinderMultiGrid::Apply2Current()
 	{
 		//cerr << "Base: curr wait on base done, wait on sync" << endl;
 		m_WaitOnBase->wait();	//base current updates are done, tell child to start its current updates
-		m_WaitOnChild->wait();	//wait for child to finisch its updates
-		SyncCurrents();			//child is finisch, run sync and go to voltage updates next
+		m_WaitOnChild->wait();	//wait for child to finish its updates
+		SyncCurrents();			//child is finish, run sync and go to voltage updates next
 		m_WaitOnSync->wait();	//sync is done... move on and tell child to move on...
 	}
 	else
 	{
 		m_WaitOnChild->wait();	//child is finished current updates, will tell base to run sync...
-		m_WaitOnSync->wait();	//wait for base to finisch sync before going to wait for next voltage updates
+		m_WaitOnSync->wait();	//wait for base to finish sync before going to wait for next voltage updates
 		//cerr << "Child: curr done, wait on sync" << endl;
 	}
 }
diff --git a/FDTD/extensions/engine_extension.h b/FDTD/extensions/engine_extension.h
index 12ef712..7dfdb00 100644
--- a/FDTD/extensions/engine_extension.h
+++ b/FDTD/extensions/engine_extension.h
@@ -41,27 +41,27 @@ public:
 
 	virtual void SetNumberOfThreads(int nrThread);
 
-	//! This methode will be called __before__ the main engine does the usual voltage updates. This methode may __not__ change the engine voltages!!!
+	//! This method will be called __before__ the main engine does the usual voltage updates. This method may __not__ change the engine voltages!!!
 	virtual void DoPreVoltageUpdates() {}
 	virtual void DoPreVoltageUpdates(int threadID);
-	//! This methode will be called __after__ the main engine does the usual voltage updates. This methode may __not__ change the engine voltages!!!
+	//! This method will be called __after__ the main engine does the usual voltage updates. This method may __not__ change the engine voltages!!!
 	virtual void DoPostVoltageUpdates() {}
 	virtual void DoPostVoltageUpdates(int threadID);
-	//! This methode will be called __after__ all updates to the voltages and extensions and may add/set its results to the engine voltages, but may __not__ rely on the current value of the engine voltages!!!
+	//! This method will be called __after__ all updates to the voltages and extensions and may add/set its results to the engine voltages, but may __not__ rely on the current value of the engine voltages!!!
 	virtual void Apply2Voltages() {}
 	virtual void Apply2Voltages(int threadID);
 
-	//! This methode will be called __before__ the main engine does the usual current updates. This methode may __not__ change the engine current!!!
+	//! This method will be called __before__ the main engine does the usual current updates. This method may __not__ change the engine current!!!
 	virtual void DoPreCurrentUpdates() {}
 	virtual void DoPreCurrentUpdates(int threadID);
-	//! This methode will be called __after__ the main engine does the usual current updates. This methode may __not__ change the engine current!!!
+	//! This method will be called __after__ the main engine does the usual current updates. This method may __not__ change the engine current!!!
 	virtual void DoPostCurrentUpdates() {}
 	virtual void DoPostCurrentUpdates(int threadID);
-	//! This methode will be called __after__ all updates to the current and extensions and may add/set its results to the engine current, but may __not__ rely on the current value of the engine current!!!
+	//! This method will be called __after__ all updates to the current and extensions and may add/set its results to the engine current, but may __not__ rely on the current value of the engine current!!!
 	virtual void Apply2Current() {}
 	virtual void Apply2Current(int threadID);
 
-	//! Set the Engine to this extention. This will usually done automatically by Engine::AddExtension
+	//! Set the Engine to this extension. This will usually done automatically by Engine::AddExtension
 	virtual void SetEngine(Engine* eng) {m_Eng=eng;}
 
 	//! Get the priority for this extension
diff --git a/FDTD/extensions/operator_ext_lorentzmaterial.cpp b/FDTD/extensions/operator_ext_lorentzmaterial.cpp
index 21ad184..6c261c8 100644
--- a/FDTD/extensions/operator_ext_lorentzmaterial.cpp
+++ b/FDTD/extensions/operator_ext_lorentzmaterial.cpp
@@ -314,7 +314,7 @@ bool Operator_Ext_LorentzMaterial::BuildExtension()
 							if (L_D[n]>0)
 							{
 								v_int[n].push_back((2.0*L_D[n]-dT*R_D[n])/(2.0*L_D[n]+dT*R_D[n]));
-								// check for r==0 in clyindrical coords and get special VI cooefficient
+								// check for r==0 in clyindrical coords and get special VI coefficient
 								if (m_CC_R0_included && n==2 && pos[0]==0)
 									v_ext[n].push_back(dT/(L_D[n]+dT*R_D[n]/2.0)*m_Op_Cyl->m_Cyl_Ext->vi_R0[pos[2]]);
 								else
diff --git a/FDTD/extensions/operator_ext_upml.cpp b/FDTD/extensions/operator_ext_upml.cpp
index e61dbdf..6683459 100644
--- a/FDTD/extensions/operator_ext_upml.cpp
+++ b/FDTD/extensions/operator_ext_upml.cpp
@@ -281,7 +281,7 @@ bool Operator_Ext_UPML::SetGradingFunction(string func)
 	int res = m_GradingFunction->Parse(m_GradFunc.c_str(), "D,dl,W,Z,N");
 	if (res < 0) return true;
 
-	cerr << "Operator_Ext_UPML::SetGradingFunction: Warning, an error occured parsing the pml grading function (see below) ..." << endl;
+	cerr << "Operator_Ext_UPML::SetGradingFunction: Warning, an error occurred parsing the pml grading function (see below) ..." << endl;
 	cerr << func << "\n" << string(res, ' ') << "^\n" << m_GradingFunction->ErrorMsg() << "\n";
 	return false;
 }
diff --git a/FDTD/extensions/operator_ext_upml.h b/FDTD/extensions/operator_ext_upml.h
index 5256655..0649d70 100644
--- a/FDTD/extensions/operator_ext_upml.h
+++ b/FDTD/extensions/operator_ext_upml.h
@@ -23,11 +23,11 @@
 
 class FunctionParser;
 
-//! Operator extension implemention an uniaxial perfectly matched layer (upml)
+//! Operator extension implementation an uniaxial perfectly matched layer (upml)
 /*
   The priority for this extension should be the highest of all extensions since this operator will use the main engine to perform vital parts in the upml implementation.
   Therefore the voltages and currents as well as the operator are replaced during these update process.
-  This extension is propably incompatible with the most other extensions operating in the same regions.
+  This extension is probably incompatible with the most other extensions operating in the same regions.
   */
 class Operator_Ext_UPML : public Operator_Extension
 {
diff --git a/FDTD/extensions/operator_extension.h b/FDTD/extensions/operator_extension.h
index 21e42af..503f1fe 100644
--- a/FDTD/extensions/operator_extension.h
+++ b/FDTD/extensions/operator_extension.h
@@ -40,7 +40,7 @@ public:
 	//! Create a clone of this extension, will return NULL if this is impossible
 	/*!
 		Create a clone of this extension, will return NULL if this is impossible (e.g. derived extension has no clone method and copy-constructor)...
-		BuildExtension has to be called separatly!
+		BuildExtension has to be called separately!
 	*/
 	virtual Operator_Extension* Clone(Operator* op) {UNUSED(op); return NULL;}
 
diff --git a/FDTD/openems_fdtd_mpi.cpp b/FDTD/openems_fdtd_mpi.cpp
index 964bdd4..3fa41d9 100644
--- a/FDTD/openems_fdtd_mpi.cpp
+++ b/FDTD/openems_fdtd_mpi.cpp
@@ -334,7 +334,7 @@ unsigned int openEMS_FDTD_MPI::GetNextStep()
 
 	int local_step=step;
 
-	//find the smallest next step requestes by all processings
+	//find the smallest next step requests by all processings
 	MPI_Reduce(&local_step, &step, 1, MPI_INT, MPI_MIN, 0, MPI_COMM_WORLD);
 	//send the smallest next step to all
 	MPI_Bcast(&step, 1, MPI_INT, 0, MPI_COMM_WORLD);
@@ -397,7 +397,7 @@ bool openEMS_FDTD_MPI::SetupProcessing()
 			if (dynamic_cast<ProcessIntegral*>(proc)!=NULL)
 			{
 				//type is integral processing --> disable! Needs to be fixed!
-				cerr << "openEMS_FDTD_MPI::SetupProcessing(): Warning: Processing: " << proc->GetName() << " occures multiple times and is being deactivated..." << endl;
+				cerr << "openEMS_FDTD_MPI::SetupProcessing(): Warning: Processing: " << proc->GetName() << " occurs multiple times and is being deactivated..." << endl;
 				cerr << "openEMS_FDTD_MPI::SetupProcessing(): Note: Processing: Make sure that there are no splits inside probes or sources." << endl;
 				deactivate = true;
 				rename = false;
@@ -405,12 +405,12 @@ bool openEMS_FDTD_MPI::SetupProcessing()
 			if (dynamic_cast<ProcessFields*>(proc)!=NULL)
 			{
 				//type is field processing --> renameing! Needs to be fixed!
-				cerr << "openEMS_FDTD_MPI::SetupProcessing(): Warning: Processing: " << proc->GetName() << " occures multiple times and is being renamed..." << endl;
+				cerr << "openEMS_FDTD_MPI::SetupProcessing(): Warning: Processing: " << proc->GetName() << " occurs multiple times and is being renamed..." << endl;
 				deactivate = false;
 				rename = true;
 			}
 		}
-		//broadcast informations to all
+		//broadcast information to all
 		MPI_Bcast(&deactivate, 1, MPI::BOOL, 0, MPI_COMM_WORLD);
 		MPI_Bcast(&rename, 1, MPI::BOOL, 0, MPI_COMM_WORLD);
 		if (deactivate)
diff --git a/FDTD/operator.cpp b/FDTD/operator.cpp
index 259b1fb..d435c45 100644
--- a/FDTD/operator.cpp
+++ b/FDTD/operator.cpp
@@ -373,7 +373,7 @@ int Operator::SnapLine2Mesh(const double* start, const double* stop, unsigned in
 		return ret;
 
 	//fixme, do we need to do something about start or stop being outside the field domain?
-	//maybe caclulate the intersection point and snap to that?
+	//maybe calculate the intersection point and snap to that?
 	//it seems to work like this as well...
 
 	return ret;
@@ -450,7 +450,7 @@ Grid_Path Operator::FindPath(double start[], double stop[])
 			currPos[minDir]+=-1;
 			minPos[minDir]-=1;
 		}
-		//check validity of current postion
+		//check validity of current position
 		for (int n=0;n<3;++n)
 			if (currPos[n]>=numLines[n])
 			{
@@ -1397,7 +1397,7 @@ bool Operator::AverageMatCellCenter(int ny, const unsigned int* pos, double* Eff
 	for (int n=0; n<4; ++n)
 		if (std::isnan(EffMat[n]) || std::isinf(EffMat[n]))
 		{
-			cerr << "Operator::" << __func__ << ": Error, an effective material parameter is not a valid result, this should NOT have happend... exit..." << endl;
+			cerr << "Operator::" << __func__ << ": Error, an effective material parameter is not a valid result, this should NOT have happened... exit..." << endl;
 			cerr << ny << "@" << n << " : " << pos[0] << "," << pos[1] << ","  << pos[2] << endl;
 			exit(0);
 		}
@@ -1510,7 +1510,7 @@ bool Operator::AverageMatQuarterCell(int ny, const unsigned int* pos, double* Ef
 	for (int n=0; n<4; ++n)
 		if (std::isnan(EffMat[n]) || std::isinf(EffMat[n]))
 		{
-			cerr << "Operator::" << __func__ << ": Error, An effective material parameter is not a valid result, this should NOT have happend... exit..." << endl;
+			cerr << "Operator::" << __func__ << ": Error, An effective material parameter is not a valid result, this should NOT have happened... exit..." << endl;
 			cerr << ny << "@" << n << " : " << pos[0] << "," << pos[1] << ","  << pos[2] << endl;
 			exit(0);
 		}
@@ -1696,7 +1696,7 @@ bool Operator::Calc_LumpedElements()
 
 			}
 			else
-				cerr << "Operator::Calc_LumpedElements(): Warning: Primitves other than boxes are not supported for lumped elements! skipping "
+				cerr << "Operator::Calc_LumpedElements(): Warning: Primitives other than boxes are not supported for lumped elements! skipping "
 						<< prims.at(bn)->GetTypeName() << " ID: " << prims.at(bn)->GetID() << " @ Property: " << PLE->GetName() << endl;
 		}
 	}
diff --git a/FDTD/operator.h b/FDTD/operator.h
index 8fc0409..233ecf6 100644
--- a/FDTD/operator.h
+++ b/FDTD/operator.h
@@ -134,7 +134,7 @@ public:
 	//! Get the volume of an FDTD cell
 	virtual double GetCellVolume(const unsigned int pos[3], bool dualMesh = false) const;
 
-	//! Get the area around an edge for a given direction \a n and a given mesh posisition \a pos
+	//! Get the area around an edge for a given direction \a n and a given mesh position \a pos
 	/*!
 		This will return the area around an edge with a given direction, measured at the middle of the edge.
 		In a cartesian mesh this is equal to the NodeArea, may be different in other coordinate systems.
@@ -143,7 +143,7 @@ public:
 
 	virtual unsigned int SnapToMeshLine(int ny, double coord, bool &inside, bool dualMesh=false, bool fullMesh=false) const;
 
-	//! Snap the given coodinates to mesh indices
+	//! Snap the given coordinates to mesh indices
 	virtual bool SnapToMesh(const double* coord, unsigned int* uicoord, bool dualMesh=false, bool fullMesh=false, bool* inside=NULL) const;
 
 	//! Snap a given box to the FDTD mesh
@@ -223,7 +223,7 @@ protected:
 	/*!
 	  Get the raw disc delta for a given position and direction.
 	  The result will be positive if a disc delta inside the simulation domain is requested.
-	  The result will be the negative value of the first or last disc delta respectivly if the position is outside the field domain.
+	  The result will be the negative value of the first or last disc delta respectively if the position is outside the field domain.
 	  */
 	virtual double GetRawDiscDelta(int ny, const int pos) const;
 
diff --git a/FDTD/operator_cylinder.h b/FDTD/operator_cylinder.h
index 152f618..b924b72 100644
--- a/FDTD/operator_cylinder.h
+++ b/FDTD/operator_cylinder.h
@@ -48,14 +48,14 @@ public:
 	//! Get the coordinates for a given node index and component, according to the cylindrical yee-algorithm. Returns true if inside the FDTD domain.
 	virtual bool GetYeeCoords(int ny, unsigned int pos[3], double* coords, bool dualMesh) const;
 
-	//! Get the node width for a given direction \a n and a given mesh posisition \a pos
+	//! Get the node width for a given direction \a n and a given mesh position \a pos
 	virtual double GetNodeWidth(int ny, const unsigned int pos[3], bool dualMesh = false) const;
-	//! Get the node width for a given direction \a n and a given mesh posisition \a pos
+	//! Get the node width for a given direction \a n and a given mesh position \a pos
 	virtual double GetNodeWidth(int ny, const int pos[3], bool dualMesh = false) const;
 
-	//! Get the node area for a given direction \a n and a given mesh posisition \a pos
+	//! Get the node area for a given direction \a n and a given mesh position \a pos
 	virtual double GetNodeArea(int n, const unsigned int* pos, bool dualMesh=false) const;
-	//! Get the node area for a given direction \a n and a given mesh posisition \a pos
+	//! Get the node area for a given direction \a n and a given mesh position \a pos
 	virtual double GetNodeArea(int ny, const int pos[3], bool dualMesh = false) const;
 
 	//! Get the length of an FDTD edge, including radius corrected alpha-mesh width.
@@ -64,7 +64,7 @@ public:
 	//! Get the volume of an FDTD cell
 	virtual double GetCellVolume(const unsigned int pos[3], bool dualMesh = false) const;
 
-	//! Get the area around an edge for a given direction \a n and a given mesh posisition \a pos
+	//! Get the area around an edge for a given direction \a n and a given mesh position \a pos
 	/*!
 		This will return the area around an edge with a given direction, measured at the middle of the edge.
 		In a cartesian mesh this is equal to the NodeArea, may be different in other coordinate systems.
diff --git a/FDTD/operator_mpi.cpp b/FDTD/operator_mpi.cpp
index 432f5d7..6fa8d1c 100644
--- a/FDTD/operator_mpi.cpp
+++ b/FDTD/operator_mpi.cpp
@@ -53,7 +53,7 @@ double Operator_MPI::CalcTimestep()
 		return ret;
 
 	double local_dT = dT;
-	//find the smallest time-step requestes by all processings
+	//find the smallest time-step requests by all processings
 	MPI_Reduce(&local_dT, &dT, 1, MPI_DOUBLE, MPI_MIN, 0, MPI_COMM_WORLD);
 	//send the smallest time-step to all
 	MPI_Bcast(&dT, 1, MPI_DOUBLE, 0, MPI_COMM_WORLD);
diff --git a/INSTALL b/INSTALL
index 6678190..4035c28 100644
--- a/INSTALL
+++ b/INSTALL
@@ -22,6 +22,6 @@ Install instructions for openEMS:
     make install (may require root)
 
     Note:
-      - all path informations may be stored in a localConfig.cmake
+      - all path information may be stored in a localConfig.cmake
       - the default "prefix" is /usr/local
 
diff --git a/TODO b/TODO
index a6e6d2c..df3f09f 100644
--- a/TODO
+++ b/TODO
@@ -3,10 +3,10 @@
 
 todo / mandatory for v0.1.0:
 - more examples and lots of testing...
-- improvments and testing for MPI engine
+- improvements and testing for MPI engine
 
 wishes:
-- location dependend excitation given by a file (e.g. mode-profile simulated with comsol)
+- location dependent excitation given by a file (e.g. mode-profile simulated with comsol)
 - more import filter (e.g. gerber import)
 
 outlook:
diff --git a/linux/debian.changelog b/linux/debian.changelog
index 755d0f4..506840b 100644
--- a/linux/debian.changelog
+++ b/linux/debian.changelog
@@ -17,7 +17,7 @@ openems (0.0.25-4) stable; urgency=low
   * Upstream fixes
  -- Sebastian Held <sebastian.held@gmx.de>  Sun, 25 Dec 2011 21:23:38 +0100
 openems (0.0.25-3) stable; urgency=low
-  * Changed dependancy on boost to enable build on oneiric
+  * Changed dependency on boost to enable build on oneiric
  -- Sebastian Held <sebastian.held@gmx.de>  Sun, 18 Dec 2011 21:39:38 +0100
 openems (0.0.25-2) stable; urgency=low
   * New upstream release
diff --git a/matlab/AR_estimate.m b/matlab/AR_estimate.m
index ab6e293..b5eb39b 100644
--- a/matlab/AR_estimate.m
+++ b/matlab/AR_estimate.m
@@ -20,7 +20,7 @@ function [val_ar t_ar f_val_ar EC] = AR_estimate( t, val, freq, nu, mu, expand_f
 %           0 --> no error
 %           1 --> input error: t and val mismatch
 %           2 --> input error: mu has to be larger than 2*nu
-%           3 --> inout error: expand_factor has to be larger than 1
+%           3 --> input error: expand_factor has to be larger than 1
 %           10 --> AR error: signal is to short for AR estimate --> decrease AR order
 %           11 --> AR error: estimated signal appears to be unstable --> use a different mu
 %
diff --git a/matlab/Add2Queue.m b/matlab/Add2Queue.m
index a80bdc9..646566c 100644
--- a/matlab/Add2Queue.m
+++ b/matlab/Add2Queue.m
@@ -1,7 +1,7 @@
 function [queue] = Add2Queue(queue,func_name, func_args, varargin)
 % function [queue] = Add2Queue(queue,func_name, func_args, varargin)
 %
-% Use this function to add a funtion to the queue.
+% Use this function to add a function to the queue.
 %
 % For more details see: InitQueue
 %
diff --git a/matlab/AddCPWPort.m b/matlab/AddCPWPort.m
index fefebf8..9308eab 100644
--- a/matlab/AddCPWPort.m
+++ b/matlab/AddCPWPort.m
@@ -68,7 +68,7 @@ evec0 = evec ./ sum(evec); % evec0 is a unit vector
 
 %set defaults
 feed_shift = 0;
-feed_R = inf; %(default is open, no resitance)
+feed_R = inf; %(default is open, no resistance)
 excite = false;
 measplanepos = nan;
 PortNamePrefix = '';
@@ -277,7 +277,7 @@ if excite
     CSX = AddBox( CSX, [PortNamePrefix 'port_excite_2_' num2str(portnr)], prio, ex_start+width_add_start, ex_stop+width_add_stop );
 end
 
-%% CPW resitance at start of CPW line
+%% CPW resistance at start of CPW line
 ex_start(idx_prop) = start(idx_prop);
 ex_stop(idx_prop) = ex_start(idx_prop);
 
@@ -292,6 +292,6 @@ elseif feed_R == 0
     CSX = AddBox( CSX, materialname, prio, ex_start-width_add_start, ex_stop-width_add_stop );
     CSX = AddBox( CSX, materialname, prio, ex_start+width_add_start, ex_stop+width_add_stop );
 else
-    error('openEMS:AddCPWPort','CPW port with resitance <= 0 it not possible');
+    error('openEMS:AddCPWPort','CPW port with resistance <= 0 it not possible');
 end
 end
diff --git a/matlab/AddCoaxialPort.m b/matlab/AddCoaxialPort.m
index 7897a6b..aae5a92 100644
--- a/matlab/AddCoaxialPort.m
+++ b/matlab/AddCoaxialPort.m
@@ -50,7 +50,7 @@ dir = DirChar2Int(dir);
 
 %set defaults
 feed_shift = 0;
-feed_R = inf; %(default is open, no resitance)
+feed_R = inf; %(default is open, no resistance)
 excite_amp = 0;
 measplanepos = nan;
 PortNamePrefix = '';
diff --git a/matlab/AddMSLPort.m b/matlab/AddMSLPort.m
index 7fa179a..8170bde 100644
--- a/matlab/AddMSLPort.m
+++ b/matlab/AddMSLPort.m
@@ -69,7 +69,7 @@ evec0 = evec ./ sum(evec); % evec0 is a unit vector
 
 %set defaults
 feed_shift = 0;
-feed_R = inf; %(default is open, no resitance)
+feed_R = inf; %(default is open, no resistance)
 excite = false;
 measplanepos = nan;
 PortNamePrefix = '';
diff --git a/matlab/AddStripLinePort.m b/matlab/AddStripLinePort.m
index 030c4ed..83ec79b 100644
--- a/matlab/AddStripLinePort.m
+++ b/matlab/AddStripLinePort.m
@@ -17,7 +17,7 @@ function [CSX,port] = AddStripLinePort( CSX, prio, portnr, materialname, start,
 %               is false)
 % 'FeedShift'   shift to port from start by a given distance in drawing
 %               units. Default is 0. Only active if 'ExcitePort' is set!
-% 'Feed_R'      Specifiy a lumped port resistance. Default is no lumped
+% 'Feed_R'      Specify a lumped port resistance. Default is no lumped
 %               port resistance --> port has to end in an ABC.
 % 'MeasPlaneShift'  Shift the measurement plane from start t a given distance
 %               in drawing units. Default is the middle of start/stop.
diff --git a/matlab/AddWaveGuidePort.m b/matlab/AddWaveGuidePort.m
index 1a17db2..336bc1c 100644
--- a/matlab/AddWaveGuidePort.m
+++ b/matlab/AddWaveGuidePort.m
@@ -74,7 +74,7 @@ for n=1:2:numel(varargin_tmp)
     end
 end
 
-% matlab adressing
+% matlab addressing
 dir = dir + 1;
 dir_sign = sign(stop(dir) - start(dir));
 if (dir_sign==0)
diff --git a/matlab/CalcNF2FF.m b/matlab/CalcNF2FF.m
index c80b52e..b6e0227 100644
--- a/matlab/CalcNF2FF.m
+++ b/matlab/CalcNF2FF.m
@@ -14,7 +14,7 @@ function nf2ff = CalcNF2FF(nf2ff, Sim_Path, freq, theta, phi, varargin)
 % freq:     array of frequencies to analyse
 % theta,phi: spherical coordinates to evaluate the far-field on (in radians)
 %
-% optional paramater:
+% optional parameter:
 % 'Center': nf2ff phase center, default is [0 0 0]
 %           !! Make sure the center is never outside of your nf2ff box!!
 %           Definition is the correct coordinate system necessary
diff --git a/matlab/ConvertHDF5_VTK.m b/matlab/ConvertHDF5_VTK.m
index ff9af8f..1cf6cd2 100644
--- a/matlab/ConvertHDF5_VTK.m
+++ b/matlab/ConvertHDF5_VTK.m
@@ -56,7 +56,7 @@ end
 
 if (do_FD_dump)
     if (~isfield(field,'FD'))
-        warning('openEMS:ConvertHDF5_VTK','no FD data found skipping frequency domian vtk dump...');
+        warning('openEMS:ConvertHDF5_VTK','no FD data found skipping frequency domain vtk dump...');
     else
         %set weighting
         if (numel(weight)~=numel(field.FD.frequency))
@@ -88,7 +88,7 @@ end
 
 if (do_TD_dump)
     if (~isfield(field,'TD'))
-        warning('openEMS:ConvertHDF5_VTK','no TD data found skipping time domian vtk dump...');
+        warning('openEMS:ConvertHDF5_VTK','no TD data found skipping time domain vtk dump...');
     else
         disp('dumping time domain data...')
         acc = ['%0' int2str(ceil(log10(numel(field.TD.time)+1))) 'd'];
diff --git a/matlab/Dump2VTK.m b/matlab/Dump2VTK.m
index e57b153..84f086c 100644
--- a/matlab/Dump2VTK.m
+++ b/matlab/Dump2VTK.m
@@ -1,7 +1,7 @@
 function Dump2VTK(filename, fields, mesh, fieldname, varargin)
 % Dump2VTK(filename, fields, mesh, fieldname, varargin)
 %
-%   Dump fields extraced from an hdf5 file to a vtk file format
+%   Dump fields extracted from an hdf5 file to a vtk file format
 %
 %   possible arguments:
 %       'NativeDump': 0 (default) / 1, dump in native coordinate system
diff --git a/matlab/FindFreeSSH.m b/matlab/FindFreeSSH.m
index 8fd6c22..1879e87 100644
--- a/matlab/FindFreeSSH.m
+++ b/matlab/FindFreeSSH.m
@@ -44,7 +44,7 @@ end
 if ischar(host_list)
     fid=fopen(host_list);
     if (fid==-1)
-        error('FindFreeSSH: cant open host file');
+        error('FindFreeSSH: cannot open host file');
     end
     clear host_list;
     host_list = {};
diff --git a/matlab/InitFDTD.m b/matlab/InitFDTD.m
index dba5120..2c94e9c 100644
--- a/matlab/InitFDTD.m
+++ b/matlab/InitFDTD.m
@@ -1,7 +1,7 @@
 function FDTD = InitFDTD(varargin)
 % function FDTD = InitFDTD(varargin)
 %
-% Inititalize the FDTD data-structure.
+% Initialize the FDTD data-structure.
 %
 % optional field arguments for usage with openEMS:
 %   NrTS:           max. number of timesteps to simulate (e.g. default=1e9)
diff --git a/matlab/SetBoundaryCond.m b/matlab/SetBoundaryCond.m
index fbfbd4c..eed4eab 100644
--- a/matlab/SetBoundaryCond.m
+++ b/matlab/SetBoundaryCond.m
@@ -11,7 +11,7 @@ function FDTD = SetBoundaryCond(FDTD, BC, varargin)
 % 
 % example:
 % BC = [ 1     1     0     0     2     3     ]  %using numbers or
-% BC = {'PMC' 'PMC' 'PEC' 'PEC' 'MUR' 'PML_8'}  %usign equivalent strings
+% BC = {'PMC' 'PMC' 'PEC' 'PEC' 'MUR' 'PML_8'}  %using equivalent strings
 %
 % mur-abc definitions
 % define a phase-velocity to be used by the mur-abc
diff --git a/matlab/Tutorials/Bent_Patch_Antenna.m b/matlab/Tutorials/Bent_Patch_Antenna.m
index 33d5cac..b7f10b0 100644
--- a/matlab/Tutorials/Bent_Patch_Antenna.m
+++ b/matlab/Tutorials/Bent_Patch_Antenna.m
@@ -1,7 +1,7 @@
 %
 % Tutorials / bent patch antenna
 %
-% Describtion at:
+% Description at:
 % http://openems.de/index.php/Tutorial:_Bent_Patch_Antenna
 %
 % Tested with
@@ -76,7 +76,7 @@ start = [patch.radius+substrate.thickness -substr_ang_width/2 -substrate.length/
 stop  = [patch.radius+substrate.thickness +substr_ang_width/2  substrate.length/2];
 CSX = AddBox( CSX, 'Jt_patch', 0, start, stop );
 
-%% create ground (not really necessary, only for esthetic reasons)
+%% create ground (not really necessary, only for aesthetic reasons)
 CSX = AddMetal( CSX, 'gnd' ); % create a perfect electric conductor (PEC)
 start = [patch.radius -substr_ang_width/2 -substrate.length/2];
 stop  = [patch.radius +substr_ang_width/2 +substrate.length/2];
@@ -158,7 +158,7 @@ ylabel( 'reflection coefficient |S_{11}|' );
 
 drawnow
 
-%find resonance frequncy from s11
+%find resonance frequency from s11
 f_res_ind = find(s11==min(s11));
 f_res = freq(f_res_ind);
 
diff --git a/matlab/Tutorials/CRLH_Extraction.m b/matlab/Tutorials/CRLH_Extraction.m
index e0bc361..921271a 100644
--- a/matlab/Tutorials/CRLH_Extraction.m
+++ b/matlab/Tutorials/CRLH_Extraction.m
@@ -1,7 +1,7 @@
 %
 % Tutorials / CRLH_Extraction
 %
-% Describtion at:
+% Description at:
 % http://openems.de/index.php/Tutorial:_CRLH_Extraction
 %
 % Tested with
diff --git a/matlab/Tutorials/CRLH_LeakyWaveAnt.m b/matlab/Tutorials/CRLH_LeakyWaveAnt.m
index a8359b7..e5f8aa0 100644
--- a/matlab/Tutorials/CRLH_LeakyWaveAnt.m
+++ b/matlab/Tutorials/CRLH_LeakyWaveAnt.m
@@ -1,7 +1,7 @@
 %
 % Tutorials / CRLH_LeakyWaveAnt
 %
-% Describtion at:
+% Description at:
 % http://openems.de/index.php/Tutorial:_CRLH_Leaky_Wave_Antenna
 %
 % Tested with
diff --git a/matlab/Tutorials/Circ_Waveguide.m b/matlab/Tutorials/Circ_Waveguide.m
index d61a713..6126521 100644
--- a/matlab/Tutorials/Circ_Waveguide.m
+++ b/matlab/Tutorials/Circ_Waveguide.m
@@ -1,7 +1,7 @@
 %
 % Tutorials / Circ_Waveguide
 %
-% Describtion at:
+% Description at:
 % http://openems.de/index.php/Tutorial:_Circular_Waveguide
 %
 % Tested with
diff --git a/matlab/Tutorials/Conical_Horn_Antenna.m b/matlab/Tutorials/Conical_Horn_Antenna.m
index f266993..dc640ef 100644
--- a/matlab/Tutorials/Conical_Horn_Antenna.m
+++ b/matlab/Tutorials/Conical_Horn_Antenna.m
@@ -1,7 +1,7 @@
 %
 % Tutorials / conical horn antenna
 %
-% Describtion at:
+% Description at:
 % http://openems.de/index.php/Tutorial:_Conical_Horn_Antenna
 %
 % Tested with
diff --git a/matlab/Tutorials/CylindricalWave_CC.m b/matlab/Tutorials/CylindricalWave_CC.m
index b81d2d4..33b63ea 100644
--- a/matlab/Tutorials/CylindricalWave_CC.m
+++ b/matlab/Tutorials/CylindricalWave_CC.m
@@ -1,7 +1,7 @@
 %
 % Tutorials / CylindricalWave_CC
 %
-% Describtion at:
+% Description at:
 % http://openems.de/index.php/Tutorial:_2D_Cylindrical_Wave
 %
 % Tested with
diff --git a/matlab/Tutorials/Dipole_SAR.m b/matlab/Tutorials/Dipole_SAR.m
index a172117..12ca996 100644
--- a/matlab/Tutorials/Dipole_SAR.m
+++ b/matlab/Tutorials/Dipole_SAR.m
@@ -1,7 +1,7 @@
 %
 % Tutorials / Dipole SAR + Power budget
 %
-% Describtion at:
+% Description at:
 % http://openems.de/index.php/Tutorial:_Dipole_SAR
 %
 % Tested with
diff --git a/matlab/Tutorials/Helical_Antenna.m b/matlab/Tutorials/Helical_Antenna.m
index b94ebd6..394dae2 100644
--- a/matlab/Tutorials/Helical_Antenna.m
+++ b/matlab/Tutorials/Helical_Antenna.m
@@ -1,7 +1,7 @@
 %
 % Tutorials / helical antenna
 %
-% Describtion at:
+% Description at:
 % http://openems.de/index.php/Tutorial:_Helical_Antenna
 %
 % Tested with
diff --git a/matlab/Tutorials/Horn_Antenna.m b/matlab/Tutorials/Horn_Antenna.m
index d463734..6b4e911 100644
--- a/matlab/Tutorials/Horn_Antenna.m
+++ b/matlab/Tutorials/Horn_Antenna.m
@@ -1,7 +1,7 @@
 %
 % Tutorials / horn antenna
 %
-% Describtion at:
+% Description at:
 % http://openems.de/index.php/Tutorial:_Horn_Antenna
 %
 % Tested with
diff --git a/matlab/Tutorials/MRI_LP_Birdcage.m b/matlab/Tutorials/MRI_LP_Birdcage.m
index 2434f03..9d98048 100644
--- a/matlab/Tutorials/MRI_LP_Birdcage.m
+++ b/matlab/Tutorials/MRI_LP_Birdcage.m
@@ -1,7 +1,7 @@
 %
 % Tutorials / 3T MRI Low Pass Birdcage coil
 %
-% Describtion at:
+% Description at:
 % http://openems.de/index.php/Tutorial:_MRI_LP_Birdcage
 %
 % Estimated time to run:    ~7h @ ~65MC/s
@@ -69,7 +69,7 @@ body_model_transform = {'Rotate_X',pi,'Rotate_Z',pi, ...
                         'Translate',[0,5,-720]};
 
 %% some internal parameter
-physical_constants % load important physical constans
+physical_constants % load important physical constants
 end_crit = 1e-5;    %abort simulation at -50dB energy drop
 unit = 1e-3;        %drawing unit used
 
@@ -164,7 +164,7 @@ for n=1:BC.N_rungs
     stop  = [BC.rad a0+da_Segs/2+da_Strip/2 -0.5*BC.portlength];
     CSX = AddBox(CSX,'metal',1, start, stop);
     
-    % some additonal mesh lines
+    % some additional mesh lines
     mesh.a = [mesh.a a0+da_Segs/2];
     
     a0 = a0 + da_Segs;
diff --git a/matlab/Tutorials/MRI_Loop_Coil.m b/matlab/Tutorials/MRI_Loop_Coil.m
index 088f8f3..30cdfe7 100644
--- a/matlab/Tutorials/MRI_Loop_Coil.m
+++ b/matlab/Tutorials/MRI_Loop_Coil.m
@@ -1,7 +1,7 @@
 %
 % Tutorials / 7T MRI Loop Coil
 %
-% Describtion at:
+% Description at:
 % http://openems.de/index.php/Tutorial:_MRI_Loop_Coil
 %
 % Tested with
@@ -145,7 +145,7 @@ CSX = AddBox(CSX, 'body_model', 0, body_box.start, body_box.stop);
 % create loop mesh
 mesh = DetectEdges(CSX);
 
-% add a dense homegeneous mesh inside the human body model
+% add a dense homogeneous mesh inside the human body model
 mesh.x = [mesh.x mesh_box.start(1) mesh_box.stop(1)];
 mesh.y = [mesh.y mesh_box.start(2) mesh_box.stop(2)];
 mesh.z = [mesh.z mesh_box.start(3) mesh_box.stop(3)];
@@ -178,7 +178,7 @@ CSX = AddBox(CSX,'SAR_xz',0, body_box.start.*[1 0 1], body_box.stop.*[1 0 1]);
 %% boundary conditions
 mesh = AddPML(mesh, 10);
 
-%% finaly define the FDTD mesh grid
+%% finally define the FDTD mesh grid
 disp(['number of cells: ' num2str(1e-6*numel(mesh.x)*numel(mesh.y)*numel(mesh.z)) ' Mcells'])
 CSX = DefineRectGrid( CSX, unit, mesh );
 
diff --git a/matlab/Tutorials/MSL_NotchFilter.m b/matlab/Tutorials/MSL_NotchFilter.m
index 612456c..ecca91d 100644
--- a/matlab/Tutorials/MSL_NotchFilter.m
+++ b/matlab/Tutorials/MSL_NotchFilter.m
@@ -1,7 +1,7 @@
 %
 % Tutorials / MSL_NotchFilter
 %
-% Describtion at:
+% Description at:
 % http://openems.de/index.php/Tutorial:_Microstrip_Notch_Filter
 %
 % Tested with
diff --git a/matlab/Tutorials/Parallel_Plate_Waveguide.m b/matlab/Tutorials/Parallel_Plate_Waveguide.m
index cbf76c0..36a1a49 100644
--- a/matlab/Tutorials/Parallel_Plate_Waveguide.m
+++ b/matlab/Tutorials/Parallel_Plate_Waveguide.m
@@ -1,7 +1,7 @@
 %
 % Tutorials / Parallel_Plate_Waveguide
 %
-% Describtion at:
+% Description at:
 % http://openems.de/index.php/Tutorial:_Parallel_Plate_Waveguide
 %
 % Tested with
diff --git a/matlab/Tutorials/Patch_Antenna_Array.m b/matlab/Tutorials/Patch_Antenna_Array.m
index ed12c67..87e98cb 100644
--- a/matlab/Tutorials/Patch_Antenna_Array.m
+++ b/matlab/Tutorials/Patch_Antenna_Array.m
@@ -6,10 +6,10 @@ function [port nf2ff] = Patch_Antenna_Array(Sim_Path, postproc_only, show_struct
 %
 % Sim_Path: Simulation path
 % postproc_only: set to post process only 0/1
-% show_structure: show the strucuture in AppCSXCAD 0/1
+% show_structure: show the structure in AppCSXCAD 0/1
 % xpos: the x-position for each antenna is defined
 % caps: the port capacity (will override active port)
-% resist: port resitance
+% resist: port resistance
 % active: switch port active
 %
 % References:
diff --git a/matlab/Tutorials/Patch_Antenna_Phased_Array.m b/matlab/Tutorials/Patch_Antenna_Phased_Array.m
index 873f345..2f4d0a9 100644
--- a/matlab/Tutorials/Patch_Antenna_Phased_Array.m
+++ b/matlab/Tutorials/Patch_Antenna_Phased_Array.m
@@ -1,7 +1,7 @@
 %
 % Tutorials / Patch Antenna Phased Array
 %
-% Describtion at:
+% Description at:
 %
 % Tested with
 %  - Matlab 2011a
@@ -122,7 +122,7 @@ disp(['I2/I1: Matlab: ' num2str(I_out(2)/I_out(1))])
 disp(['I3/I1: Matlab: ' num2str(I_out(3)/I_out(1))])
 
 
-%% do a referenc simulation for the given C2/C3 values
+%% do a reference simulation for the given C2/C3 values
 if (do_reference_simulation)
     active = [1 0 0];
     caps = [0 C2 C3];
@@ -132,7 +132,7 @@ if (do_reference_simulation)
     port_ref = calcPort( port_ref, Sim_Path, f, 'RefImpedance', 50);
     nf2ff_ref = CalcNF2FF(nf2ff_ref, Sim_Path, f0, [-180:2:180]*pi/180, 0);
 
-    % extract currents from referenc simulation
+    % extract currents from reference simulation
     for p=1:3
         I_ref(p,1) = interp1(f, port_ref{p}.if.tot,f0);
     end
@@ -141,7 +141,7 @@ if (do_reference_simulation)
     disp(['I3/I1: openEMS: ' num2str(I_ref(3)/I_ref(1))])
 end
 
-%% calculate and apply weighting cooefficients [3]
+%% calculate and apply weighting coefficients [3]
 % calculate
 coeff = I\I_out;
 
diff --git a/matlab/Tutorials/RCS_Sphere.m b/matlab/Tutorials/RCS_Sphere.m
index b8bb5bc..26c2081 100644
--- a/matlab/Tutorials/RCS_Sphere.m
+++ b/matlab/Tutorials/RCS_Sphere.m
@@ -1,7 +1,7 @@
 %
 % Tutorials / radar cross section of a metal sphere
 %
-% Describtion at:
+% Description at:
 % http://openems.de/index.php/Tutorial:_RCS_Sphere
 %
 % Tested with
diff --git a/matlab/Tutorials/RadarUWBTutorial.m b/matlab/Tutorials/RadarUWBTutorial.m
index 6fd2f5b..0f07bfc 100644
--- a/matlab/Tutorials/RadarUWBTutorial.m
+++ b/matlab/Tutorials/RadarUWBTutorial.m
@@ -27,7 +27,7 @@ physical_constants;
 % in different channels.
 % 
 % The Gaussian excitation is set to the same 3dB bandwidth as the
-% channels of the IEEE 802.15.4 UWB PHY. One exeption is channel4twice
+% channels of the IEEE 802.15.4 UWB PHY. One exception is channel4twice
 % which has the double bandwidth of channel 4. It can be seen that the
 % delay is larger and the fidelity is smaller in the vicinity of the
 % (undesired) resonances of the antenna. Note that for a real UWB system
@@ -117,7 +117,7 @@ CSX = AddBox(CSX, 'Patch', 2, [-patchsize/2, gap, 0], [patchsize/2, gap + patchs
 mesh.x = [];
 mesh.y = [];
 
-% two mesh lines for the metal coatings of teh substrate
+% two mesh lines for the metal coatings of the substrate
 mesh.z = linspace(-substrate.height, 0, substrate.cells +1);
 
 % find optimal mesh lines for the patch and ground, not yes the microstrip line
@@ -139,7 +139,7 @@ mesh.x = [mesh.x meshline.x];
 % and only the top of the substrate, the other edges are covered by the ground plane
 mesh.y = [mesh.y, meshsubstrate.y(end)]; % top of substrate
 
-% for now we have only the edges, now calculate mesh lines inbetween
+% for now we have only the edges, now calculate mesh lines in between
 mesh = SmoothMesh(mesh, fineResolution);
 
 % add the outer boundary
@@ -228,7 +228,7 @@ figure %("visible", "off");
 polarFF(nf2ff, 'xaxis', 'theta', 'logscale', [minplot, maxplot]);
 title(["fidelity ", suffix, " / %"]);
 
-% save the plots in order to compare them afer simulating the different channels
+% save the plots in order to compare them after simulating the different channels
 print(1, ["s11_", suffix, ".png"]);
 print(2, ["farfield_", suffix, ".png"]);
 print(3, ["delay_mm_", suffix, ".png"]);
diff --git a/matlab/Tutorials/Rect_Waveguide.m b/matlab/Tutorials/Rect_Waveguide.m
index 3b1a79a..505be6e 100644
--- a/matlab/Tutorials/Rect_Waveguide.m
+++ b/matlab/Tutorials/Rect_Waveguide.m
@@ -1,7 +1,7 @@
 %
 % Tutorials / Rect_Waveguide
 %
-% Describtion at:
+% Description at:
 % http://openems.de/index.php/Tutorial:_Rectangular_Waveguide
 %
 % Tested with
@@ -21,7 +21,7 @@ unit = 1e-6; %drawing unit in um
 % waveguide dimensions
 % WR42
 a = 10700;   %waveguide width
-b = 4300;    %waveguide heigth
+b = 4300;    %waveguide height
 length = 50000;
 
 % frequency range of interest
diff --git a/matlab/Tutorials/Simple_Patch_Antenna.m b/matlab/Tutorials/Simple_Patch_Antenna.m
index 1e24c4c..1f34b19 100644
--- a/matlab/Tutorials/Simple_Patch_Antenna.m
+++ b/matlab/Tutorials/Simple_Patch_Antenna.m
@@ -1,6 +1,6 @@
 %% Simple Patch Antenna Tutorial
 %
-% Describtion at:
+% Description at:
 % <http://openems.de/index.php/Tutorial:_Simple_Patch_Antenna>
 %
 % Tested with
@@ -147,7 +147,7 @@ ylabel( 'reflection coefficient |S_{11}|' );
 drawnow
 
 %% NFFF Plots
-%find resonance frequncy from s11
+%find resonance frequency from s11
 f_res_ind = find(s11==min(s11));
 f_res = freq(f_res_ind);
 
diff --git a/matlab/Tutorials/StripLine2MSL.m b/matlab/Tutorials/StripLine2MSL.m
index 64eeaa1..7fb0024 100644
--- a/matlab/Tutorials/StripLine2MSL.m
+++ b/matlab/Tutorials/StripLine2MSL.m
@@ -1,7 +1,7 @@
 %
 % Stripline to Microstrip Line Transition
 %
-% Describtion at:
+% Description at:
 % <http://openems.de/index.php/Tutorial:_Stripline_to_MSL_Transition>
 %
 % Tested with
@@ -68,7 +68,7 @@ CSX = AddBox( CSX, 'RO4350B', 0, start, stop );
 CSX = AddMetal( CSX, 'gnd' ); 
 CSX = AddMetal( CSX, 'metal' ); 
 
-% Create strip line port (incl. metal stip line)
+% Create strip line port (incl. metal strip line)
 start = [-line_length -strip_width/2  0];
 stop  = [0            +strip_width/2  0];
 [CSX,port{1}] = AddStripLinePort( CSX, 100, 1, 'metal', start, stop, strip_substrate_thickness, 'x', [0 0 -1], 'ExcitePort', true, 'FeedShift', feed_shift, 'MeasPlaneShift', meas_shift );
diff --git a/matlab/Tutorials/readme b/matlab/Tutorials/readme
index a80dfa6..a8c6811 100644
--- a/matlab/Tutorials/readme
+++ b/matlab/Tutorials/readme
@@ -1 +1 @@
-* Find the tutorial describtions at http://openems.de/index.php/Tutorials
+* Find the tutorial descriptions at http://openems.de/index.php/Tutorials
diff --git a/matlab/calcPort.m b/matlab/calcPort.m
index 9b56955..867f0f2 100644
--- a/matlab/calcPort.m
+++ b/matlab/calcPort.m
@@ -32,7 +32,7 @@ function [port] = calcPort( port, SimDir, f, varargin)
 %   port.f                  the given frequency fector
 %   port.uf.tot/inc/ref     total, incoming and reflected voltage
 %   port.if.tot/inc/ref     total, incoming and reflected current
-%   port.ZL_ref             used refernce impedance
+%   port.ZL_ref             used reference impedance
 %
 %   port.P_inc              incoming power
 %   port.P_ref              reflected power
diff --git a/matlab/calcTLPort.m b/matlab/calcTLPort.m
index 68912c6..e294489 100644
--- a/matlab/calcTLPort.m
+++ b/matlab/calcTLPort.m
@@ -34,7 +34,7 @@ function [port] = calcTLPort( port, SimDir, f, varargin)
 %   port.if.tot/inc/ref     total, incoming and reflected current
 %   port.beta:              propagation constant
 %   port.ZL:                characteristic line impedance
-%   port.ZL_ref             used refernce impedance
+%   port.ZL_ref             used reference impedance
 %
 % example:
 %   port{1} = calcTLPort( port{1}, Sim_Path, f, 'RefImpedance', 50);
diff --git a/matlab/examples/antennas/inverted_f.m b/matlab/examples/antennas/inverted_f.m
index 175f94c..d2dc473 100644
--- a/matlab/examples/antennas/inverted_f.m
+++ b/matlab/examples/antennas/inverted_f.m
@@ -184,7 +184,7 @@ ylabel( 'reflection coefficient |S_{11}|' );
 drawnow
 
 %% NFFF contour plots %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-%find resonance frequncy from s11
+%find resonance frequency from s11
 f_res_ind = find(s11==min(s11));
 f_res = freq(f_res_ind);
 
diff --git a/matlab/examples/optimizer/optimizer_simfun.m b/matlab/examples/optimizer/optimizer_simfun.m
index 14152de..31bb0d1 100644
--- a/matlab/examples/optimizer/optimizer_simfun.m
+++ b/matlab/examples/optimizer/optimizer_simfun.m
@@ -37,7 +37,7 @@ disp( ['DONE. L = ' num2str(L(1)/1e-6) ' uH'] );
 goal   = 2e-6; % specify the goal: 2 uH
 result = abs(goal - L(1)); % costs must not be negative
 
-% restore curent folder
+% restore current folder
 cd( oldpwd );
 
 
diff --git a/matlab/examples/other/Metamaterial_PlaneWave_Drude.m b/matlab/examples/other/Metamaterial_PlaneWave_Drude.m
index db72b3e..feea063 100644
--- a/matlab/examples/other/Metamaterial_PlaneWave_Drude.m
+++ b/matlab/examples/other/Metamaterial_PlaneWave_Drude.m
@@ -22,7 +22,7 @@ pic_size = round([1400 1400/4]); %define the animation picture size
 length = 500;
 width = 10;
 mesh_res = 0.5;        % mesh resolution
-height = 3*mesh_res; % hight is ony 3 lines with PEC (top/bottom) --> quasi 2D
+height = 3*mesh_res; % height is only 3 lines with PEC (top/bottom) --> quasi 2D
 
 %FDTD setup
 f0 = 5e9;         %center frequency
diff --git a/matlab/examples/transmission_lines/CPW_Line.m b/matlab/examples/transmission_lines/CPW_Line.m
index 6dac636..7d2fdaf 100644
--- a/matlab/examples/transmission_lines/CPW_Line.m
+++ b/matlab/examples/transmission_lines/CPW_Line.m
@@ -1,7 +1,7 @@
 %
 % Tutorials / CPW_Line
 %
-% Describtion at:
+% Description at:
 %
 % Tested with
 %  - Octave 3.8.1
@@ -33,7 +33,7 @@ if 1
   feed_shift_cells = 0;
   x_spacing = air_spacing;
 else % or use a waveguide with start/end in a pml
-  feed_R = inf; % CPW ends in a pml --> disable termination resitance 
+  feed_R = inf; % CPW ends in a pml --> disable termination resistance 
   feed_shift_cells = 10; % CPW ends in an 8 cells thick pml --> shift feed 10 cells
   pml_add_cells = [0 0 8 8 8 8]; % do not add air-space in x-direction
   x_spacing = 0; % do not add air-space in x-direction
diff --git a/matlab/examples/transmission_lines/MSL.m b/matlab/examples/transmission_lines/MSL.m
index b6ec0b3..3bc8447 100644
--- a/matlab/examples/transmission_lines/MSL.m
+++ b/matlab/examples/transmission_lines/MSL.m
@@ -88,7 +88,7 @@ CSX = SetMaterialWeight( CSX, 'fakepml', 'Sigma', ['pow(z-' num2str(length-abs_l
 start = [mesh.x(1)  mesh.y(1)  length-abs_length];
 stop  = [mesh.x(end) mesh.y(end) length];
 % the geometric priority is set to 0, which is lower than the priority
-% of the MSL, thus the MSL (copper) has precendence
+% of the MSL, thus the MSL (copper) has precedence
 priority = 0;
 CSX = AddBox( CSX, 'fakepml', priority, start, stop );
 
diff --git a/matlab/examples/waveguide/Rect_Waveguide.m b/matlab/examples/waveguide/Rect_Waveguide.m
index a9601ad..8d57196 100644
--- a/matlab/examples/waveguide/Rect_Waveguide.m
+++ b/matlab/examples/waveguide/Rect_Waveguide.m
@@ -28,7 +28,7 @@ numTS = 50000; %max. number of timesteps
 % waveguide dimensions
 length = 1000;
 a = 1000;   %waveguide width
-b = 600;    %waveguide heigth
+b = 600;    %waveguide height
 
 %waveguide TE-mode definition
 m = 1;
diff --git a/openems.cpp b/openems.cpp
index a3ee256..1a40200 100644
--- a/openems.cpp
+++ b/openems.cpp
@@ -1218,7 +1218,7 @@ bool openEMS::DumpStatistics(const string& filename, double time)
 	stat_file << FDTD_Op->GetNumberCells() << "\t% number of cells" << endl;
 	stat_file << FDTD_Op->GetTimestep() << "\t% timestep (s)" << endl;
 	stat_file << FDTD_Eng->GetNumberOfTimesteps() << "\t% number of iterations" << endl;
-	stat_file << FDTD_Eng->GetNumberOfTimesteps()*FDTD_Op->GetTimestep() << "\t% total numercial time (s)" << endl;
+	stat_file << FDTD_Eng->GetNumberOfTimesteps()*FDTD_Op->GetTimestep() << "\t% total numerical time (s)" << endl;
 	stat_file << time << "\t% simulation time (s)" << endl;
 	stat_file << (double)FDTD_Op->GetNumberCells()*(double)FDTD_Eng->GetNumberOfTimesteps()/time << "\t% speed (cells/s)" << endl;
 
diff --git a/openems.h b/openems.h
index b308b13..b76b113 100644
--- a/openems.h
+++ b/openems.h
@@ -90,7 +90,7 @@ public:
 	int Get_PML_Size(int idx);
 	void Set_Mur_PhaseVel(int idx, double val);
 
-	//! Get informations about external libs used by openEMS
+	//! Get information about external libs used by openEMS
 	static std::string GetExtLibsInfo(std::string prefix="\t");
 
 	//! Get welcome screen for openEMS
diff --git a/python/Tutorials/CRLH_Extraction.py b/python/Tutorials/CRLH_Extraction.py
index 4e1850e..00984d3 100644
--- a/python/Tutorials/CRLH_Extraction.py
+++ b/python/Tutorials/CRLH_Extraction.py
@@ -2,7 +2,7 @@
 """
  Tutorials / CRLH_Extraction
 
- Describtion at:
+ Description at:
  http://openems.de/index.php/Tutorial:_CRLH_Extraction
 
  Tested with
diff --git a/python/Tutorials/MSL_NotchFilter.py b/python/Tutorials/MSL_NotchFilter.py
index f036e0f..0448aa0 100644
--- a/python/Tutorials/MSL_NotchFilter.py
+++ b/python/Tutorials/MSL_NotchFilter.py
@@ -2,7 +2,7 @@
 """
  Microstrip Notch Filter Tutorial
 
- Describtion at:
+ Description at:
  http://openems.de/doc/openEMS/Tutorials.html#microstrip-notch-filter
 
  Tested with
diff --git a/python/Tutorials/RCS_Sphere.py b/python/Tutorials/RCS_Sphere.py
index 6daaebd..db3eae2 100644
--- a/python/Tutorials/RCS_Sphere.py
+++ b/python/Tutorials/RCS_Sphere.py
@@ -86,7 +86,7 @@ if not post_proc_only:
     FDTD.Run(Sim_Path, verbose=3, cleanup=True)
 
 ### Postprocessing & plotting
-# get Gaussian pulse stength at frequency f0
+# get Gaussian pulse strength at frequency f0
 ef = UI_data('et', Sim_Path, freq=f0)
 
 Pin = 0.5*norm(E_dir)**2/Z0 * abs(ef.ui_f_val[0])**2
diff --git a/python/Tutorials/Rect_Waveguide.py b/python/Tutorials/Rect_Waveguide.py
index 5d38115..226c2e2 100644
--- a/python/Tutorials/Rect_Waveguide.py
+++ b/python/Tutorials/Rect_Waveguide.py
@@ -2,7 +2,7 @@
 """
  Rectangular Waveguide Tutorial
 
- Describtion at:
+ Description at:
  http://openems.de/doc/openEMS/Tutorials.html#rectangular-waveguide
 
  Tested with
@@ -30,7 +30,7 @@ unit = 1e-6; #drawing unit in um
 # waveguide dimensions
 # WR42
 a = 10700;   #waveguide width
-b = 4300;    #waveguide heigth
+b = 4300;    #waveguide height
 length = 50000;
 
 # frequency range of interest
diff --git a/python/openEMS/nf2ff.py b/python/openEMS/nf2ff.py
index d0ed198..6f852f1 100644
--- a/python/openEMS/nf2ff.py
+++ b/python/openEMS/nf2ff.py
@@ -61,7 +61,7 @@ class nf2ff:
         self.dump_type = 0  # default Et/Ht
         self.dump_mode = 1  # default cell interpolated
 
-        self.freq = None   # broadband recording by defualt
+        self.freq = None   # broadband recording by default
         if 'frequency' in kw:
             self.freq = kw['frequency']
             del kw['frequency']
@@ -147,13 +147,13 @@ class nf2ff:
         result = nf2ff_results(fn)
         if result.phi is not None:
             if not np.abs((result.r-radius)/radius)<1e-6:
-                raise Exception('Radius does not match. Did you read an invalid chached result? Try "read_cached=False"')
+                raise Exception('Radius does not match. Did you read an invalid cached result? Try "read_cached=False"')
             if not utilities.Check_Array_Equal(np.rad2deg(result.theta), self.theta, 1e-4):
-                raise Exception('Theta array does not match. Did you read an invalid chached result? Try "read_cached=False"')
+                raise Exception('Theta array does not match. Did you read an invalid cached result? Try "read_cached=False"')
             if not utilities.Check_Array_Equal(np.rad2deg(result.phi), self.phi, 1e-4):
-                raise Exception('Phi array does not match. Did you read an invalid chached result? Try "read_cached=False"')
+                raise Exception('Phi array does not match. Did you read an invalid cached result? Try "read_cached=False"')
             if not utilities.Check_Array_Equal(result.freq, self.freq, 1e-6, relative=True):
-                raise Exception('Frequency array does not match. Did you read an invalid chached result? Try "read_cached=False"')
+                raise Exception('Frequency array does not match. Did you read an invalid cached result? Try "read_cached=False"')
         return result
 
 class nf2ff_results:
diff --git a/python/openEMS/openEMS.pyx b/python/openEMS/openEMS.pyx
index d5d129f..e68fc07 100644
--- a/python/openEMS/openEMS.pyx
+++ b/python/openEMS/openEMS.pyx
@@ -122,7 +122,7 @@ cdef class openEMS:
     def SetEndCriteria(self, val):
         """ SetEndCriteria(val)
 
-        Set the end critera value. E.g. 1e-6 for -60dB
+        Set the end criteria value. E.g. 1e-6 for -60dB
         """
         self.thisptr.SetEndCriteria(val)
 
@@ -265,7 +265,7 @@ cdef class openEMS:
     def AddLumpedPort(self, port_nr, R, start, stop, p_dir, excite=0, **kw):
         """ AddLumpedPort(port_nr, R, start, stop, p_dir, excite=0, **kw)
 
-        Add a lumped port wit the given values and location.
+        Add a lumped port with the given values and location.
 
         See Also
         --------
@@ -426,7 +426,7 @@ cdef class openEMS:
         Run the openEMS FDTD simulation.
 
         :param sim_path: str -- path to run in and create result data
-        :param cleanup: bool -- remove exisiting sim_path to cleanup old results
+        :param cleanup: bool -- remove existing sim_path to cleanup old results
         :param setup_only: bool -- only perform FDTD setup, do not run simulation
         :param verbose: int -- set the openEMS verbosity level 0..3
 
diff --git a/python/openEMS/ports.py b/python/openEMS/ports.py
index 6185ac4..5a934f4 100644
--- a/python/openEMS/ports.py
+++ b/python/openEMS/ports.py
@@ -58,7 +58,7 @@ class Port(object):
     :param excite: float -- port excitation amplitude
     :param priority: int -- priority of all contained primtives
     :param PortNamePrefix: str -- a prefix for all ports-names
-    :param delay: float -- a positiv delay value to e.g. emulate a phase shift
+    :param delay: float -- a positive delay value to e.g. emulate a phase shift
     """
     def __init__(self, CSX, port_nr, start, stop, excite, **kw):
         self.CSX      = CSX
diff --git a/tools/AdrOp.h b/tools/AdrOp.h
index e4a690b..c9ad239 100644
--- a/tools/AdrOp.h
+++ b/tools/AdrOp.h
@@ -42,7 +42,7 @@ public:
 	///Deconstructor
 	virtual ~AdrOp();
 	///Set the current n-dim position, get 1-dim array position as return value
-	/*!A position has to be set or all other methodes will case error! \n The methode will exit with error message if invalid position is set! \sa ErrorMsg */
+	/*!A position has to be set or all other methods will case error! \n The method will exit with error message if invalid position is set! \sa ErrorMsg */
 	unsigned int SetPos(unsigned int muiIpos, unsigned int muiJpos, unsigned int muiKpos=0, unsigned int muiLpos=0);
 
 	bool SetPosChecked(unsigned int muiIpos, unsigned int muiJpos, unsigned int muiKpos=0, unsigned int muiLpos=0);
@@ -60,25 +60,25 @@ public:
 
 	double GetIndexWidth(int ny, int index);
 	double GetIndexCoord(int ny, int index);
-	///Get the gird delta at the given  index of direction ny. (if index<0 return negative value as index=0 would give, if index>=max-1 returns negative value as index=max-2 would give)
+	///Get the grid delta at the given index of direction ny. (if index<0 return negative value as index=0 would give, if index>=max-1 returns negative value as index=max-2 would give)
 	double GetIndexDelta(int ny, int index);
 
 //	double GetCellVolume(unsigned int uiCell);
 
 	unsigned int GetPosFromNode(int ny, unsigned int uiNode);
 	///Set a shift in ny direction (e.g. 0 for i-direction)
-	/*!Shift set by this methode will be ignored by methode GetPos*/
+	/*!Shift set by this method will be ignored by method GetPos*/
 	unsigned int Shift(int ny, int step);
 	///Set a checked shift in ny direction (e.g. 0 for i-direction)
-	/*!Shift set by this methode will be ignored by methode GetPos*/
+	/*!Shift set by this method will be ignored by method GetPos*/
 	bool CheckShift(int ny, int step);
-	///Returns the current 1-dim position including shift by methode "Shift" + additional (transitory) shift
+	///Returns the current 1-dim position including shift by method "Shift" + additional (transitory) shift
 	unsigned int GetShiftedPos(int ny=-1, int step=0);
-	///Reset shift set by "Shift"-methode
+	///Reset shift set by "Shift"-method
 	void ResetShift();
 	///Iterates through AdrOp; --- obsolete ---
 	unsigned int Iterate(int jump=1);
-	///Retruns size of array
+	///Returns size of array
 	unsigned int GetSize();
 	///Set mode to reflect by node
 	/*!1D-example (6 nodes): \image html node_reflect.PNG order: 0,1,2,3,4,5,4,3,...*/
diff --git a/tools/ErrorMsg.cpp b/tools/ErrorMsg.cpp
index 9a1974c..ff01d49 100644
--- a/tools/ErrorMsg.cpp
+++ b/tools/ErrorMsg.cpp
@@ -60,14 +60,14 @@ void ErrorMsg::Error(unsigned int nr,char *chAddMsg)
 	if ((nr>0) && (nr<=NoMsg))
 	{
 		if (Msg[nr-1]!=NULL) fprintf(stderr,"%s",Msg[nr-1]);
-		else fprintf(stderr,"unkown error occured!! Error code: %d exiting...",nr);
+		else fprintf(stderr,"unknown error occurred!! Error code: %d exiting...",nr);
 		if (chAddMsg!=NULL) fprintf(stderr,"%s",chAddMsg);
 		getchar();
 		exit(nr);
 	}
 	else
 	{
-		fprintf(stderr,"unkown error occured!! Error code: %d exiting...",nr);
+		fprintf(stderr,"unknown error occurred!! Error code: %d exiting...",nr);
 		getchar();
 		exit(nr);
 	}
@@ -78,14 +78,14 @@ void ErrorMsg::Error(unsigned int nr,int addNr)
 	if ((nr>0) && (nr<=NoMsg))
 	{
 		if (Msg[nr-1]!=NULL) fprintf(stderr,"%s",Msg[nr-1]);
-		else fprintf(stderr,"unkown error occured!! Error code: %d exiting...",nr);
+		else fprintf(stderr,"unknown error occurred!! Error code: %d exiting...",nr);
 		fprintf(stderr,"%d",addNr);
 		getchar();
 		exit(nr);
 	}
 	else
 	{
-		fprintf(stderr,"unkown error occured!! Error code: %d exiting...",nr);
+		fprintf(stderr,"unknown error occurred!! Error code: %d exiting...",nr);
 		getchar();
 		exit(nr);
 	}
@@ -93,6 +93,6 @@ void ErrorMsg::Error(unsigned int nr,int addNr)
 
 void ErrorMsg::ownError(void)
 {
-	fprintf(stdout," Error occured by using Error Message class!! ... exiting...");
+	fprintf(stdout," Error occurred by using Error Message class!! ... exiting...");
 	exit(-1);
 }
diff --git a/tools/ErrorMsg.h b/tools/ErrorMsg.h
index fa63c93..b98aaa5 100644
--- a/tools/ErrorMsg.h
+++ b/tools/ErrorMsg.h
@@ -32,7 +32,7 @@ public:
 	ErrorMsg(unsigned int NoMessage=0);
 	///Deconstructor
 	virtual ~ErrorMsg();
-	///Methode for defining error messages
+	///Method for defining error messages
 	/*! \param nr Number of defining error message \param *Message Set error message string  \sa Error */
 	void SetMsg(unsigned int nr, const char *Message);
 	///Call an error message. Will exit the program!
diff --git a/tools/hdf5_file_writer.cpp b/tools/hdf5_file_writer.cpp
index 56d3214..f989cc0 100644
--- a/tools/hdf5_file_writer.cpp
+++ b/tools/hdf5_file_writer.cpp
@@ -455,7 +455,7 @@ bool HDF5_File_Writer::WriteAtrribute(std::string locName, std::string attr_name
 	hid_t attribute_id = H5Acreate(loc, attr_name.c_str(), mem_type, dataspace_id, H5P_DEFAULT, H5P_DEFAULT);
 	if (attribute_id<0)
 	{
-		cerr << "HDF5_File_Writer::WriteAtrribute: Error, failed to create the attrbute" << endl;
+		cerr << "HDF5_File_Writer::WriteAtrribute: Error, failed to create the attribute" << endl;
 		H5Sclose(dataspace_id);
 		H5Oclose(loc);
 		H5Fclose(hdf5_file);
@@ -465,7 +465,7 @@ bool HDF5_File_Writer::WriteAtrribute(std::string locName, std::string attr_name
 	/* Write the attribute data. */
 	if (H5Awrite(attribute_id, mem_type, value)<0)
 	{
-		cerr << "HDF5_File_Writer::WriteAtrribute: Error, failed to write the attrbute" << endl;
+		cerr << "HDF5_File_Writer::WriteAtrribute: Error, failed to write the attribute" << endl;
 		H5Aclose(attribute_id);
 		H5Sclose(dataspace_id);
 		H5Oclose(loc);
diff --git a/tools/sar_calculation.cpp b/tools/sar_calculation.cpp
index 791e842..099ea00 100644
--- a/tools/sar_calculation.cpp
+++ b/tools/sar_calculation.cpp
@@ -267,16 +267,16 @@ int SAR_Calculation::FindFittingCubicalMass(unsigned int pos[3], float box_size,
 			return 0;
 		}
 
-		// if no valid or finally invalid cube is found, calculate an alternaive cube size
+		// if no valid or finally invalid cube is found, calculate an alternative cube size
 		if (mass_iterations==0)
 		{
-			// on first interation, try a relative resize
+			// on first iteration, try a relative resize
 			old_box_size=box_size;
 			box_size*=pow(m_avg_mass/mass,1.0/3.0);
 		}
 		else
 		{
-			// on later interations, try a newton approach
+			// on later iterations, try a newton approach
 			float new_box_size = box_size - (mass-m_avg_mass)/(mass-old_mass)*(box_size-old_box_size);
 			old_box_size = box_size;
 			box_size = new_box_size;
diff --git a/tools/useful.cpp b/tools/useful.cpp
index 209781c..b15b27f 100644
--- a/tools/useful.cpp
+++ b/tools/useful.cpp
@@ -63,7 +63,7 @@ std::vector<unsigned int> AssignJobs2Threads(unsigned int jobs, unsigned int nrT
 	}
 
 	if (jobs>0)
-		std::cerr << "AssignJobs2Threads: Error, " << jobs << " remain to be assigned, this should not have happend..." << std::endl;
+		std::cerr << "AssignJobs2Threads: Error, " << jobs << " remain to be assigned, this should not have happened..." << std::endl;
 
 	if (RemoveEmpty)
 	{
diff --git a/tools/vtk_file_writer.cpp b/tools/vtk_file_writer.cpp
index 79c40f3..be6586f 100644
--- a/tools/vtk_file_writer.cpp
+++ b/tools/vtk_file_writer.cpp
@@ -70,7 +70,7 @@ void VTK_File_Writer::SetMeshLines(double const* const* lines, unsigned int cons
 		vtkRectilinearGrid* RectGrid = dynamic_cast<vtkRectilinearGrid*>(m_GridData);
 		if (RectGrid==NULL)
 		{
-			cerr << "VTK_File_Writer::SetMeshLines: Error, grid invalid, this should not have happend! " << endl;
+			cerr << "VTK_File_Writer::SetMeshLines: Error, grid invalid, this should not have happened! " << endl;
 			exit(1);
 		}
 		RectGrid->SetDimensions(count[0],count[1],count[2]);
@@ -97,7 +97,7 @@ void VTK_File_Writer::SetMeshLines(double const* const* lines, unsigned int cons
 		vtkStructuredGrid* StructGrid = dynamic_cast<vtkStructuredGrid*>(m_GridData);
 		if (StructGrid==NULL)
 		{
-			cerr << "VTK_File_Writer::SetMeshLines: Error, grid invalid, this should not have happend! " << endl;
+			cerr << "VTK_File_Writer::SetMeshLines: Error, grid invalid, this should not have happened! " << endl;
 			exit(1);
 		}