I will try to do what you have suggested.
I am also posting the entire ctl file.
Kindly let me know if there is anything in the problem that is causing such a
behaviour.
Thanks a lot!
Sanchit
-------------------------------ctl file-------------------------------------
(define-param l 8)
(define-param w 0.1)
(define-param lsub 4)
(define-param wsub 4)
(define-param ewvg 12)
(define-param esub 4)
(define-param pad 2)
(define-param dpml 2)
(define sxy 32)
(set! geometry-lattice (make lattice (size sxy sxy no-size)))
(define wvg-ycen 0)
(define wvg-xcen (+ (* -0.5 (- sxy (* 2 pad) (* 2 dpml))) 4) )
(define sub-ycen -2.05)
(define sub-xcen (- wvg-xcen 2))
(set! geometry
(list
(make block
(center wvg-xcen wvg-ycen)
(size l w infinity)
(material (make dielectric (epsilon ewvg))))
(make block
(center sub-xcen sub-ycen)
(size lsub wsub infinity)
(material (make dielectric (epsilon esub))))))
(set! sources (list
(make source
(src (make continuous-src (wavelength 1.55)))
(component Ez)
(center (* -0.5 (- sxy (* 2 pad) (* 2 dpml))) wvg-ycen)
(size 0 w))))
(set! pml-layers (list (make pml (thickness dpml))))
(set-param! resolution 10)
(define (myoutput)
(print (meep-time) ","" Incident, " (get-field-point Ez (vector3
-12 0))
","" Transmitted," (get-field-point Ez (vector3 -3.5 0))"\n"))
(run-until 200 (at-every 0.5 myoutput))
-----------------------------------------------------------------------------
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
