Correct normalization factor for Electric field

- Use 1 instead of -IQ(NP) to convert E field from V/cm to 1/cm. - Use RM instead of hard-coded 0.511 for the electron rest mass.
nrc-cnrc · Jan 20, 2016 · 9e18712 · 9e18712 · ftessier · Jan 28, 2016
1 parent aadebb9
commit 9e18712
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/HEN_HOUSE/src/emf_macros.mortran b/HEN_HOUSE/src/emf_macros.mortran
@@ -420,10 +420,10 @@ REPLACE {$GET-POTENTIAL(#,#,#,#);} WITH {;
 "***********************************************"
 REPLACE {$GET-EM-FIELD(#,#,#,#,#);} WITH {;
 "fixed e-field normalization below--factor of 10^6 must be in denominator!--BW"
-fbtemp=-IQ(NP)/(RM*10**6);
+fbtemp=  1.0/(RM*10**6);
 {P1}X0=ExIN*fbtemp;{P1}Y0=EyIN*fbtemp;{P1}Z0=EzIN*fbtemp;
-fbtemp=3.0/0.511;
-{P2}X0=BxIN*fbtemp;{P2}Y0=ByIN*fbtemp;{P2}Z0=BzIN*fbtemp;
+fbtemp=3.0/RM;
+{P2}X0=Bx*fbtemp;{P2}Y0=By*fbtemp;{P2}Z0=Bz*fbtemp;
 }
 ;
 "REPLACE {$SET-TUSTEP-EM-FIELD;} WITH {;}"