Dear Qian, as for the problems in your code,
(1) the position of your source is in the centre and should actually be at (center 0 (/ sy -2)) (2) the correct kx should probably be in the medium of incidence, i.e. (define kx (* refractiveindex frec (sin theta))) (3) You may want to check whether the angle of incidence you specify is actually above the critical angle of total internal reflection. Yours is actually not. (4) the time you run your simulation is rather long, maybe you should check the output in between in order to see the evanescent wave forming. Have fun, Andreas. ----------------------------------------- Andreas Erbe Interface Spectroscopy Group Department of Interface Chemistry and Surface Engineering Max Planck Institute for Iron Research GmbH PO Box 140444 40074 Düsseldorf Germany --- Max-Planck-Institut für Eisenforschung GmbH Max-Planck-Straße 1 D-40237 Düsseldorf Handelsregister B 2533 Amtsgericht Düsseldorf Geschäftsführung Prof. Dr. Jörg Neugebauer Prof. Dr. Dierk Raabe Prof. Dr. Martin Stratmann Dipl.-Kfm. Herbert Wilk Ust.-Id.-Nr.: DE 11 93 58 514 Steuernummer: 105 5891 1000 ------------------------------------------------- _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
