From 9741b654e05e1f11d89c9bb3839d7a472d7fa863 Mon Sep 17 00:00:00 2001
From: Thorsten Liebig <Thorsten.Liebig@gmx.de>
Date: Sun, 26 Feb 2017 18:37:58 +0100
Subject: [PATCH] nf2ff: fix P_rad calculation with radius != 1m

Signed-off-by: Thorsten Liebig <Thorsten.Liebig@gmx.de>
---
 nf2ff/nf2ff_calc.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/nf2ff/nf2ff_calc.cpp b/nf2ff/nf2ff_calc.cpp
index 335dd9a..2bfd3b5 100644
--- a/nf2ff/nf2ff_calc.cpp
+++ b/nf2ff/nf2ff_calc.cpp
@@ -496,7 +496,7 @@ bool nf2ff_calc::AddSinglePlane(float **lines, unsigned int* numLines, complex<f
 			m_H_theta[tn][pn] += factor*(Np[tn][pn] - Lt[tn][pn]/Z0);
 			m_H_phi[tn][pn] -= factor*(Nt[tn][pn] + Lp[tn][pn]/Z0);
 
-			m_P_rad[tn][pn] = m_radius*m_radius/(2*fZ0) * abs((m_E_theta[tn][pn]*conj(m_E_theta[tn][pn])+m_E_phi[tn][pn]*conj(m_E_phi[tn][pn])));
+			m_P_rad[tn][pn] = abs((m_E_theta[tn][pn]*conj(m_E_theta[tn][pn])+m_E_phi[tn][pn]*conj(m_E_phi[tn][pn])))/(2*fZ0);
 			if (m_P_rad[tn][pn]>P_max)
 				P_max = m_P_rad[tn][pn];
 		}