Dear meepers,
I just started to use Meep and I'm a bit confused with the code.
I would to change the position of the waist in a Gaussian beam. For that, I'm
trying to introduce a nonlinear phase shift in y² (in accordance to the
wavefront curvature). But I can't figure out where and how this phase term can
be introduced!
Any idea?
Thank you very much!
My code so far looks like:
-----------------------------------------------------------------------------------------------------
(set! geometry-lattice (make lattice (size 400 80 no-size)))
(define ((my-amp sigma) x)
(exp (- (/ (vector3-dot x x) (* 2 sigma sigma))))
)
(set! sources (list
(make source
(src (make continuous-src (wavelength 2) (width 20) ))
(component Ez) (center -180 20 0) (size 0 78 99)
(amp-func (my-amp 4))
)))
(set! pml-layers (list (make pml (thickness 1.0))))
(set! resolution 8)
(set! force-complex-fields? true)
(run-until 4000
(at-beginning output-epsilon)
(at-beginning (output-png Ez "-a /usr/share/h5utils/colormaps/yarg
-A $EPS -Zc /usr/share/h5utils/colormaps/dkbluered"))
(at-end output-efield-z)
(at-every 50 (output-png Ez "-Zc
/usr/share/h5utils/colormaps/dkbluered"))
(at-every 50 output-efield-z)
(at-end output-efield)
(at-end output-hfield))
-----------------------------------------------------------------------------------------------------------------
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
