Dear Steven and others, My ctl is as follows, **********************************************************************************************************
(define-param sy 100) ; size of cell in y direction
(define-param sx 50) ; padding between last hole and PML edge
(define-param dpml 2) ; PML thickness
(set! geometry-lattice (make lattice (size sx sy no-size)))
(set! pml-layers (list (make pml (thickness dpml))))
(set-param! resolution 10) ;
(use-output-directory)
(define-param fcen 0.3) ;
(define-param df 0.1) ;
(define ( (my-amp sigma k) x)
(exp (- (* 0+2i pi (vector3-dot k x))
(/ (vector3-dot x x) (* 2 sigma sigma)))))
(set! sources (list
(make source
(src (make continuous-src (frequency fcen) (fwidth df)))
(component Hz)
(center 0 -48) ;
(size 10 0)
(amp-func (my-amp 1 (vector3 0))))))
(run-until 400
(at-beginning output-epsilon)
(at-every 1.2 (output-png Hz "-Zc gray")))
**********************************************************************************************************
I expect that a gaussian beam propagating along the y direction with the beam
waist of 1 unit will be generated. But my anticipation is not consistent with
the result which is attached. So how to generate a guassian beam?
Best Regards,
Li
PS: If I want to have a gaussian beam propagating along the bisector between x
and y, how to specify k? (vector3 1 1)?
hz-000154.90.png
Description: PNG image
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
