try to setup an extended drude

to be tested

tests/seteps/setepsdriver-f90-new-1.f90

@@ -17,7 +17,7 @@ program setepsdriver
   double precision, allocatable :: tmp_osc(:,:)
   integer, allocatable :: tmp_oscType(:)
   
-  integer :: i, j, k, l, noPoints
+  integer :: i, j, k, l, m, noPoints
 
   write(*,*) 'setepsdriver-f90-new-1 start'
   
@@ -72,6 +72,14 @@ program setepsdriver
           osc(2, jos) = osc(3, jos)
           osc(3, jos) = 0
         endif
+        if (osc(1, jos) < 1e-7) then
+          m = nos(l) / 2
+          oscType(jos) = extendedDrude
+          osc(1, jos) = osc(1, jos + m)
+          osc(2, jos) = osc(3, jos)
+          osc(3, jos) = osc(3, jos + m)
+          nos(l) = nos(l) - 1
+        endif
       enddo
     endif
   enddo