The code is atached. It is 2D TM implementation. No Boundary conditions are implemented. The extension '.java' is not allowed in attachments. So I changed the extension to '.for'. A C version of the code I may be able to give later. (I am a bit busy for next 3 weeks).