Hi Young, Can you please explain more on what result do you expect for your structure?
I am working on the same problem as you. The only difference is that I perform a get-field-point action on a particular point and on the origin location of the source. Next I perform a fourier transform on these two points and normalize. This is what Steven advised previously and I hope it is useful for you. Chan-Hoe -----Original Message----- From: [EMAIL PROTECTED] on behalf of Young Chul Jun Sent: Thu 4/19/2007 2:08 AM To: [email protected] Subject: [Meep-discuss] decay rate calculation Dear Steven and other Meep users, I'm trying to calculate J*E at the source position to invetigate the decay rate variation (or LDOS) in my cavity structure. First, I located a point source having a fixed frequency and polarization Ex and waited for some time until transients died away. And I tried to get the magnitude of the electric field (Ex) at the same source position. (I used complex field to get the amplitude of the electric field at the chosen frequency - not time-varying electric fields) Then, I ran the simulation again without my structure and obtained the field value again at the source position to get the normalized decay rate. I repeated simulations with different point source positions (pointsrc_y), but my result is quite different from what I expected for my structure.. Do you have any ideas or comments for my code? Thanks very much. Any comment will be really appreciated.. Best Regards, Young (set-param! force-complex-fields? true) (define-param lambda 600) (define wcen (/ 1 lambda)) (define-param pointsrc_y 100) ... (set! sources (list (make source (src (make continuous-src (frequency wcen))) (size 0 0) (component Ex) (center 0 pointsrc_y)))) (define (print-field-at-point) (print "At y=" pointsrc_y " : " (magnitude (get-field-point Ex (vector3 0 pointsrc_y))) "\n")) (run-until (* 10 lambda) (at-end print-field-at-point)) _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
