Help FDTD+CPML urgent

Hello,

I'm working in FDTD+CPML and I need urgent help.

I′m simulating the situation of the paper that I send in attach.

In the case of PEC buildings, I can′t get the results of power delay profiles shown in the paper. I use a different mode than they, I use TEz mode, with the same characteristics dt and dx that they use.

When I calculate the power using the expression:

10log(((Ex^2+Ey^2)*lambda^2*1)/(120*pi*4*pi))

and using the excitation like:

Jmax =1000;
alpha = (4/(32*dt))^2;
source=zeros(1,nSteps);
for n=1:88%(9 e 30)
source(n)= Jmax*(exp(-alpha*(((dt*n)-(32*dt))^2)));
end

then I apply in the time loop the source like:

Hz(x,y)=source(x,y);

I′m using everything like they have in this document and I can′t reach the same results.

The plot aspect is similar but the power is bigger (50 dBw bigger like you can see in the file attached)

Please Help me.

Thanks a lot

Joao Roldao