From c0596dfb5f33a13dfdc4e8ebf5af9dc53e954a05 Mon Sep 17 00:00:00 2001
From: Karl-Michael Schindler <karl-michael.schindler@web.de>
Date: Wed, 4 Sep 2024 18:13:54 +0200
Subject: [PATCH] try to setup an extended drude

to be tested
---
 tests/seteps/setepsdriver-f90-new-1.f90 | 10 +++++++++-
 1 file changed, 9 insertions(+), 1 deletion(-)

diff --git a/tests/seteps/setepsdriver-f90-new-1.f90 b/tests/seteps/setepsdriver-f90-new-1.f90
index 60df3aa..1cf80f4 100644
--- a/tests/seteps/setepsdriver-f90-new-1.f90
+++ b/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
-- 
GitLab