Respected Sir I am trying to find the dispersion curve for Air medium.But for the given values of k it is unable to generate corresponding frequencies.
(set! default-material (make dielectric (epsilon 1))) (define-param sx 12) (define-param sy 12) (define-param PMLThickness 2) (define-param latticex (+ sx (* 2 PMLThickness))) (define-param latticey (+ sy (* 2 PMLThickness))) (define-param res 20) (define-param fcen 0.5) ; pulse center frequency (define-param df 1) ; pulse width (in frequency) (set! geometry-lattice (make lattice (size latticex latticey no-size))) (set! sources (list (make source (src (make gaussian-src (frequency fcen) (fwidth df))) (component Hz) (center 0 0)))) (set! pml-layers (list (make pml (thickness PMLThickness)))) (set! resolution res) (run-until (/ 3 fcen) (at-beginning output-epsilon) (at-every (/ 1 fcen 20) output-hfield-z)) (define-param k-points (list (vector3 0.0 0 0) (vector3 3 0 0))) ; define a series of k-points (define-param k-interp 19) (set! k-points (interpolate k-interp k-points)) ; calculates frequencies for these k-points (run-k-points 300 k-points) I donot understand what is the problem with my code>kindly help.I will be looking forward to your reply. Regards
_______________________________________________ meep-discuss mailing list meep-discuss@ab-initio.mit.edu http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss