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

Reply via email to