correct the index of the oscillators in the log

bug showed up with more than one layer

tests/SetEpsTestCasesFromScratch/f77/eels.for

@@ -194,11 +194,12 @@ C                                                                       ABTI0209
 
 C **** log modification start
       write (99, '(i5, i5)')         neps, layers
-      write (99, '(a)')              mode
+      j1 = 0
       do i = 1, neps
         write (99, '(a, g15.7, i5)') name(i), epsinf(i), nos(i)
         do j = 1, nos(i)
-          write (99, '(3g15.7)')     osc(1,j), osc(2,j), osc(3,j)
+          j1 = j1 + 1
+          write (99, '(3g15.7)')     osc(1,j1), osc(2,j1), osc(3,j1)
         enddo
       enddo
       write (99, *)

tests/SetEpsTestCasesFromScratch/f90/doeels.f90

@@ -63,7 +63,7 @@ subroutine doeels (e0, theta, phia, phib, wmin, wmax, dw, comment, comment_size,
   double complex, allocatable :: eps(:)
 
 ! **** log modification start
-  integer :: j
+  integer :: j, j1
 ! **** log modification end
 
   external fun, qrat
@@ -196,10 +196,12 @@ subroutine doeels (e0, theta, phia, phib, wmin, wmax, dw, comment, comment_size,
 ! **** log modification start
   write (99, '(i5, i5)')         neps, layers
   write (99, '(a)')              mode
+  j1 = 0
   do i = 1, neps
     write (99, '(a, g15.7, i5)') name(i), epsinf(i), nos(i)
     do j = 1, nos(i)
-      write (99, '(3g15.7)')     osc(1,j), osc(2,j), osc(3,j)
+      j1 = j1 + 1
+      write (99, '(3g15.7)')     osc(1,j1), osc(2,j1), osc(3,j1)
     enddo
   enddo
   write (99, *)

tests/SetEpsTestCasesFromScratch/f90/epsLog/seteps006.log

@@ -3,7 +3,7 @@ No-layers
 MnO          4.950000        1
    269.0000       16.00000      0.5000000E-01
 Platinum     8.900000        1
-   269.0000       16.00000      0.5000000E-01
+   185541.0      -1.000000      0.1200000    
 
    50.00000       21.52103      0.1595170      -60.54409       30923.34    
    52.00000       21.56943      0.1668695      -60.54407       29733.97    

tests/SetEpsTestCasesFromScratch/original/eels.for

@@ -187,11 +187,12 @@ C                                                                       ABTI0209
 
 C **** log modification start
       write (99, '(i5, i5)')         neps, layers
-      write (99, '(a)')              mode
+      j1 = 0
       do i = 1, neps
         write (99, '(a, g15.7, i5)') name, epsinf(i), nos(i)
         do j = 1, nos(i)
-          write (99, '(3g15.7)')     osc(1,j), osc(2,j), osc(3,j)
+          j1 = j1 + 1
+          write (99, '(3g15.7)')     osc(1,j1), osc(2,j1), osc(3,j1)
         enddo
       enddo
       write (99, *)