Excellent.
If it helps anyone else, here's a new version incorporating Steven's
suggestions.
Best Regards,
Matt
(define-param unit 1e-3) ; base unit in meters (1e-3 = 1
mm)
(define-param fh 3e9) ; frequency in hertz
(define-param res_ppw 20) ; resolution in pixels per
wavelength
(define-param num_periods 10) ; number of periods to let
simulation run for
(define lam (/ 3e8 fh)) ; free space wavelength (meters)
(define lamu (/ lam unit)) ; wavelength in base units
(define fcen (/ unit lam)) ; source center frequency
(define-param csize (* lamu 8)) ; make the cell size 8*lambda
(set! geometry-lattice (make lattice (size csize csize no-size)))
(set! geometry (list
(make block (center 0 0) (size csize csize no-size) (material
air))))
(set! pml-layers (list (make pml (thickness 100))))
(set-param! resolution (/ res_ppw lamu))
; set temporal width equal to 2 wavelength
(set! sources (list
(make source
(src (make continuous-src (frequency fcen) (width 20)))
(component Ez) (center 0 0))))
(use-output-directory)
(run-until (* num_periods lamu) (at-every (/ lamu 10) (output-png Ez "-S 2 -Zc
dkbluered")))
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
