Dear Prof. Steven and meep users! I would like to ask you about the Amplitude in (run-k-points ...) simulation. What I want to say... If I do the simulation of light propagation in 2D structure (planar dielectric waveguide) with the cw source, I set the amplitude A, for example, and I can obtain the Power (there were a lot of Prof. Steven's explanation about how to do this in this mailing list). So, If I do the simulation of gaussian pulse propagation in such structure, I also can obtain the something like (flux-in-box X (volume (center 0 0 0) (size 0 sy 0))) I have I wave, propagating in positive X direction. So I can obtain some value of power over again.
If I do the simulation in order to obtain the dispersion characteristic (omega(k) dependency, omega - frequency, k - longitudal wavenumber, for example ) I use the (run-k-point ...) in meep. I write smt like (run-until 500 (at-beginning output-epsilon) (at-every dt flux-plane) (at-every dt energy-box) (to-appended "ez" (at-every 0.5 output-efield-z)) (after-sources (harminv Ez (vector3 0.1234 0) fcen df)) ) where (define (flux-plane) (print "flux-poynting-plane:, " (meep-time) ", " (flux-in-box X (volume (center 0 0 0) (size 0 sy 0))) "\n")) So I obtain the power characteristics as well - and righ away I make a comparison and see that there is a differences in power calculation between the gaussian calculation and (run-k-points ...) calculation. The difference is big. What is the reason for... you will ask me) This is the attempt to obtain a dispersion for nonlinear structure. So I've obtained the dispersion for nonlinear periodic structure and now I would like to estimate the specific power for the specific value of nonlinear frequency shift. Or if it is possible, how to find the detailed information about how does (run-k-point ...) work ? f (maybe only the meep's source code) Cordially, Alexandr. PhD Student, Russia. _______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss