Drude model in Transfer Matrix Method?

[Dungnv]

Base on described theory about TMM for photonic crystal J.B. Pendry, J. Mod. Opt.41, 209 (1994), Prof. Pendry also uploaded a program using Fortran language [attachment] to calculate transmission/reflection coefficients of complex structure (2D dimensions).
The problem is, i want to perform this code to calculate T&R spectra View attachment adag_v1_0.gzfrom cut-wire pair structure.
This is my subroutine to define CWP structure, but it provide wrong results when i compare with CST program:

---

Code Fortran - [expand]
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
subroutine DEFSTR(EPSPOT,IXMAX,IYMAX,IZMAX,EPS,EPSREF,EPSSUB)
c
      integer IXMAX,IYMAX,IZMAX,IX,IY,IZ
      complex*16 EPSPOT(IXMAX,IYMAX,IZMAX),EPS,EPSREF,EPSSUB
 
      do 5 IX = 1,IXMAX
         do 5 IY = 1,IYMAX
            do 5 IZ = 1,IZMAX
               if ((IZ .EQ. 2) .OR. (IZ .EQ. 4)) then
                   if ((IY .GE. 3).AND.(IY .LE. 16) .AND.
     &                (IX .GE. 4) .AND. (IX .LE. 6)) then
                        EPSPOT(IX,IY,IZ) = EPS
                  else
                        EPSPOT(IX,IY,IZ) = EPSREF
                  end if
               elseif ((IZ .EQ. 3)) then
                   EPSPOT(IX,IY,IZ) = EPSSUB
               else
                        EPSPOT(IX,IY,IZ) = EPSREF
                end if
 5    continue
      return
      end

---

and the second question is, when i try to use Drude's model for epsilon of metal part, then T and R become infinitive.
Anyone know this method, help me plz.