diff --git a/Common/processfields_fd.cpp b/Common/processfields_fd.cpp
index 4057535..dd362d0 100644
--- a/Common/processfields_fd.cpp
+++ b/Common/processfields_fd.cpp
@@ -92,7 +92,9 @@ int ProcessFieldsFD::Process()
 	unsigned int pos[3];
 	for (size_t n = 0; n<m_FD_Samples.size(); ++n)
 	{
-		std::complex<float> exp_jwt = std::exp( (std::complex<float>)(-2.0 * _I * M_PI * m_FD_Samples.at(n) * T) );
+		std::complex<float> exp_jwt_2_dt = std::exp( (std::complex<float>)(-2.0 * _I * M_PI * m_FD_Samples.at(n) * T) );
+		exp_jwt_2_dt *= 2; // *2 for single-sided spectrum
+		exp_jwt_2_dt *= Op->GetTimestep(); // multiply with timestep
 		field_fd = m_FD_Fields.at(n);
 		for (pos[0]=0; pos[0]<numLines[0]; ++pos[0])
 		{
@@ -100,9 +102,9 @@ int ProcessFieldsFD::Process()
 			{
 				for (pos[2]=0; pos[2]<numLines[2]; ++pos[2])
 				{
-					field_fd[0][pos[0]][pos[1]][pos[2]] += field_td[0][pos[0]][pos[1]][pos[2]] * exp_jwt;
-					field_fd[1][pos[0]][pos[1]][pos[2]] += field_td[1][pos[0]][pos[1]][pos[2]] * exp_jwt;
-					field_fd[2][pos[0]][pos[1]][pos[2]] += field_td[2][pos[0]][pos[1]][pos[2]] * exp_jwt;
+					field_fd[0][pos[0]][pos[1]][pos[2]] += field_td[0][pos[0]][pos[1]][pos[2]] * exp_jwt_2_dt;
+					field_fd[1][pos[0]][pos[1]][pos[2]] += field_td[1][pos[0]][pos[1]][pos[2]] * exp_jwt_2_dt;
+					field_fd[2][pos[0]][pos[1]][pos[2]] += field_td[2][pos[0]][pos[1]][pos[2]] * exp_jwt_2_dt;
 				}
 			}
 		}
@@ -147,9 +149,9 @@ void ProcessFieldsFD::DumpFDData()
 					{
 						for (pos[2]=0; pos[2]<numLines[2]; ++pos[2])
 						{
-							field[0][pos[0]][pos[1]][pos[2]] = real(field_fd[0][pos[0]][pos[1]][pos[2]] * exp_jwt)/(float)m_FD_SampleCount;
-							field[1][pos[0]][pos[1]][pos[2]] = real(field_fd[1][pos[0]][pos[1]][pos[2]] * exp_jwt)/(float)m_FD_SampleCount;
-							field[2][pos[0]][pos[1]][pos[2]] = real(field_fd[2][pos[0]][pos[1]][pos[2]] * exp_jwt)/(float)m_FD_SampleCount;
+							field[0][pos[0]][pos[1]][pos[2]] = real(field_fd[0][pos[0]][pos[1]][pos[2]] * exp_jwt);
+							field[1][pos[0]][pos[1]][pos[2]] = real(field_fd[1][pos[0]][pos[1]][pos[2]] * exp_jwt);
+							field[2][pos[0]][pos[1]][pos[2]] = real(field_fd[2][pos[0]][pos[1]][pos[2]] * exp_jwt);
 						}
 					}
 				}
@@ -170,9 +172,9 @@ void ProcessFieldsFD::DumpFDData()
 					{
 						for (pos[2]=0; pos[2]<numLines[2]; ++pos[2])
 						{
-							field[0][pos[0]][pos[1]][pos[2]] = abs(field_fd[0][pos[0]][pos[1]][pos[2]])/(float)m_FD_SampleCount;
-							field[1][pos[0]][pos[1]][pos[2]] = abs(field_fd[1][pos[0]][pos[1]][pos[2]])/(float)m_FD_SampleCount;
-							field[2][pos[0]][pos[1]][pos[2]] = abs(field_fd[2][pos[0]][pos[1]][pos[2]])/(float)m_FD_SampleCount;
+							field[0][pos[0]][pos[1]][pos[2]] = abs(field_fd[0][pos[0]][pos[1]][pos[2]]);
+							field[1][pos[0]][pos[1]][pos[2]] = abs(field_fd[1][pos[0]][pos[1]][pos[2]]);
+							field[2][pos[0]][pos[1]][pos[2]] = abs(field_fd[2][pos[0]][pos[1]][pos[2]]);
 						}
 					}
 				}
@@ -218,7 +220,7 @@ void ProcessFieldsFD::DumpFDData()
 		{
 			stringstream ss;
 			ss << "f" << n;
-			DumpVectorArray2HDF5(m_filename.c_str(), "/FieldData/FD", ss.str(), m_FD_Fields.at(n),numLines,1.0/(float)m_FD_SampleCount,m_FD_Samples.at(n));
+			DumpVectorArray2HDF5(m_filename.c_str(), "/FieldData/FD", ss.str(), m_FD_Fields.at(n),numLines,1.0,m_FD_Samples.at(n));
 		}
 		return;
 	}
diff --git a/Common/processintegral.cpp b/Common/processintegral.cpp
index 3da7107..9aa37b4 100644
--- a/Common/processintegral.cpp
+++ b/Common/processintegral.cpp
@@ -76,7 +76,7 @@ void ProcessIntegral::InitProcess()
 void ProcessIntegral::FlushData()
 {
 	if (m_FD_Samples.size())
-		Dump_FD_Data(1.0/(double)m_FD_SampleCount,m_filename + "_FD");
+		Dump_FD_Data(1.0,m_filename + "_FD");
 }
 
 
@@ -149,7 +149,7 @@ int ProcessIntegral::Process()
 			for (size_t n=0; n<m_FD_Samples.size(); ++n)
 			{
 				for (int i=0; i<NrInt; ++i)
-					m_FD_Results[i].at(n) += (double)m_Results[i] * m_weight * std::exp( -2.0 * _I * M_PI * m_FD_Samples.at(n) * time );
+					m_FD_Results[i].at(n) += (double)m_Results[i] * m_weight * std::exp( -2.0 * _I * M_PI * m_FD_Samples.at(n) * time ) * 2.0 * Op->GetTimestep();
 			}
 			++m_FD_SampleCount;
 			if (m_Flush)