diff --git a/src/libnptm/tra_subs.cpp b/src/libnptm/tra_subs.cpp
index e9122e7dbd277bf18f3a876b1c833a648855a11e..1903f9d313246a982ef5d485d3abd903cab44e04 100644
--- a/src/libnptm/tra_subs.cpp
+++ b/src/libnptm/tra_subs.cpp
@@ -190,7 +190,7 @@ void ffrt(
 
   ctqce = new complex<double>[3]();
   ctqcs = new complex<double>[3]();
-  for (int l60 = 1; l60 < cil->le; l60++) {
+  for (int l60 = 1; l60 <= cil->le; l60++) {
     int lpo = l60 + 1;
     int il = l60 * lpo;
     int ltpo = l60 + lpo;