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

Reply via email to