Sorry, I misunderstood "output".
It's fine, now.
----- Original Message ----
From: Son-Geun Kang <[EMAIL PROTECTED]>
To: meep-discuss@ab-initio.mit.edu; meep-discuss@ab-initio.mit.edu
Sent: Saturday, February 10, 2007 2:05:34 PM
Subject: what's wrong with at-every?
Hi,
It seems MEEP should output Ez field every 0.3 unit. Why I get only one? And I
find the file changes when MEEP is on.
Thank you in
advance.
Songeun_Kang
******************************************************************************************************
; Some parameters to describe the geometry:
(define-param eps 7.673) ; dielectric constant of waveguide
(define-param w 13.636) ; width of waveguide
(define-param r 0.437) ; radius of holes
(define-param N 15) ; number of holes on either side of defect
; The cell dimensions
(define-param sy 40) ; size of cell in y direction
(define-param sx 40) ; padding between last hole and PML edge
(define-param dpml 1) ; PML thickness
(set! geometry-lattice (make lattice (size sx sy no-size)))
(set! geometry
(append ; combine lists of objects:
(list (make block (center 0 0) (size w infinity infinity)
(material (make dielectric (epsilon eps)))))
(list (make block (center 10 0) (size 20 w infinity)
(material (make dielectric (epsilon eps)))))
(geometric-object-duplicates (vector3 0.707 -0.707) 0 (- N 1)
(make cylinder (center 0 0) (radius r) (height infinity)
(material air)))
(geometric-object-duplicates (vector3 -0.707 0.707) 0 (- N 1)
(make cylinder (center 0 0) (radius r) (height infinity)
(material air)))))
(set! pml-layers (list (make pml (thickness dpml))))
(set-param! resolution 20)
(define-param fcen 0.145) ; pulse center frequency
(define-param df 0.25) ; pulse width (in frequency)
(define-param nfreq 500) ; number of frequencies at which to compute flux
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ez)
(center 0 18.5)
(size 0 w))))
(define transX19 ;
(add-flux fcen df nfreq
(make flux-region
(center 19 0) (size 0 (* w 2)))))
(define transY-19 ;
(add-flux fcen df nfreq
(make flux-region
(center 0 -19) (size (* w 2) 0))))
(define refl18.25 ;
(add-flux fcen df nfreq
(make flux-region
(center 0 18.25) (size 0 (* w 2)))))
(define refl18.75 ;
(add-flux fcen df nfreq
(make flux-region
(center 0 18.75) (size 0 (* w 2)))))
(run-sources+ (stop-when-fields-decayed
50 Ez
(vector3 19 0)
1e-3)
(at-beginning output-epsilon)
(to-appended "ez" (at-every 0.3 output-efield-z)))
(display-fluxes transX19 transY-19 refl18.25 refl18.75) ; print out the flux
spectrum
____________________________________________________________________________________
The fish are biting.
Get more visitors on your site using Yahoo! Search Marketing.
http://searchmarketing.yahoo.com/arp/sponsoredsearch_v2.php_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
