On Sun, Apr 08, 2012 at 12:42:33AM +0200, Filip Dominec wrote: > Dear Meep users, > I wrote a Python-MEEP simulation of a gaussian pulse passing through a simple > structure. My script can be found at http://fzu.cz/~dominecf/sphere.py; more > information is in this file and I hope it may even be useful for somebody. > > At the input and output of the simulation volume, there are flux planes > measuring the intensity transmission and reflection spectra, T(f) and R(f), > which are plotted succesfully. > > However, I would like to obtain the complex _amplitude_ spectra of > transmission and reflection, say t(f) and r(f). I believe that Meep knows > these values all the time and only at the end it simply calculates > R(f)=|r(f)|^2 and T(f)=|t(f)|^2. > > But I need to know the complex r(f), t(f) to calculate also the phase of the > passing wave. > > Any hints are appreciated. Thank you in advance!
The trick to this is to do the fourier transform yourself. Instead of your typical reflection/transmission monitors, output the electric field at these points at each timestep, then take the FFT of the resulting data. See here for an example: http://falsecolour.com/aw/meep_ptest.tar.gz In the file getsparams.m. But be warned! I was never quite able to work out all the kinks with this method. There's also post on meep-discuss where Prof. Johnson describes this kind of thing, but again the results are suspicious. -- Aaron Webster _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
