Hi,
Another thing that I don't understand is that when you specify the frequency of
a gaussian source, it actually as a different frequency. In the attached ctl
file, frequency is set to 0.15, but the transmitted flux in air shows 0.1636 !
I don't understand why.
Cheers.
--
; test the frequency range of a certain pulse
(define-param sx 10)
(define-param sy 10)
(set! geometry-lattice (make lattice (size sx sy no-size)))
(set! geometry
(list
(make block
(center 0 0)
(size 0 0 0)
(material nothing))))
(define-param fcen 0.15) ; pulse center frequency
(define-param df 0.3) ; pulse width (in frequency)
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df) ))
(component Ex)
(center 0 (- (/ sy 2) 2))
(size sx 0))))
(set! pml-layers (list (make pml (direction Y) (thickness 1.0))))
(set! resolution 10)
(define-param nfreq 100) ; number of frequencies at which to compute flux
(define trans ; transmitted flux
(add-flux fcen df nfreq
(make flux-region
(center 0 (+ (/ sy -2) 2))
(size sx 0))))
(run-until 200
(at-beginning output-epsilon)
(to-appended "ex" (at-every 0.5 output-efield-x))
(at-end output-efield-x))
(display-fluxes trans)
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
