          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 @@
+@@ -1152,6 +1160,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,6 +1168,8 @@
         y_Hz = self.y
         x_Hz_FDTD = EMpy.utils.centered1d(x)
         y_Hz_FDTD = EMpy.utils.centered1d(y)

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]
