-diff -ruN EMpy-1.1.45/EMpy/modesolvers/FD.py EMpy-1.1.45new/EMpy/modesolvers/FD.py
---- EMpy-1.1.45/EMpy/modesolvers/FD.py	2017-01-13 07:39:49.000000000 -0700
-+++ EMpy-1.1.45new/EMpy/modesolvers/FD.py	2017-01-27 18:25:37.000000000 -0700
-@@ -1121,30 +1120,38 @@
-             x = self.x
-         if y is None:
-             y = self.y
--
-+            
-         # Ex: ignores y = 0, max
-         x_Ex = EMpy.utils.centered1d(self.x)
+diff -ruN EMpy-1.1.45-orig/EMpy/modesolvers/FD.py EMpy-1.1.45/EMpy/modesolvers/FD.py
+--- EMpy-1.1.45-orig/EMpy/modesolvers/FD.py	2017-02-06 18:18:22.000000000 -0700
++++ EMpy-1.1.45/EMpy/modesolvers/FD.py	2017-02-06 18:34:23.000000000 -0700
+@@ -1127,24 +1127,32 @@
          y_Ex = EMpy.utils.centered1d(self.y)
          x_Ex_FDTD = EMpy.utils.centered1d(x)
          y_Ex_FDTD = y[1:-1]

         Hx_FDTD = EMpy.utils.interp2(
             x_Hx_FDTD, y_Hx_FDTD, x_Hx, y_Hx, self.Hx) / (-120. * numpy.pi)
         # Hy: ignores y = 0, max, /120pi, reverse direction
@@ -1152,6 +1159,8 @@
         y_Hy = self.y
         x_Hy_FDTD = EMpy.utils.centered1d(x)
         y_Hy_FDTD = y[1:-1]

         Hy_FDTD = EMpy.utils.interp2(
             x_Hy_FDTD, y_Hy_FDTD, x_Hy, y_Hy, self.Hy) / (-120. * numpy.pi)
         # Hz: /120pi, reverse direction
@@ -1159,102 +1169,104 @@
         y_Hz = self.y
         x_Hz_FDTD = EMpy.utils.centered1d(x)
         y_Hz_FDTD = EMpy.utils.centered1d(y)
