Help about CPML 2D code...

Hi boys... I programmed this code for CPML 2D simulation based on Taflove 's fortran code for 3D. But when I plot relative error for UPML (alpha max=0) and CPML (both simulations with the same caracteristics), the relative error is smaller than UPML than for CPML. I leave my code. I think it's wrong. Can somebody check it and correct it? , please.

I need it to finish my thesis...
Many thanks

Mutasili