AJ Cressman wrote:
> AJ Cressman wrote:
>> Hello,
>> I am working on simulating a cylindrical waveguide, and thought that I
>> could change the k-point of the system so that I only have fields
>> traveling to the right, and then starting at the left again once they
>> hit the right edge. I got this error when using change-k-point! :
>>
>> meep: error on line 546 of ../../../src/h5file.cpp: file data is
>> inconsistent size for subsequent extend_data
>>
>> I cannot seem to find any information on this error. If you need more
>> information, or would like to see the source code, let me know.
>> Thanks,
>> AJ Cressman
>>
>> PS- Also, since I am working in cylindrical coordinates, when I
>> specify the direction of the PML would I use R instead of X, etc.?
Here is the complete source code; some values might not make sense,
since I've been tweaking every little piece I could. Any help is much
appreciated.
--Code Starts--
(set! dimensions CYLINDRICAL)
(define-param outputdirectory "./waveguide-out/")
(set! eps-averaging? #f)
(set! geometry-lattice (make lattice (size 3 0 12)))
(set! resolution 10)
(set! k-point #f)
(set! pml-layers (list (make pml (direction R) (side High) (thickness
1.0))))
(define-param lam 1.0) ; pulse wavelength
(define-param fcen (/ 1 lam)) ; pulse center frequency
(define-param df (* 0.85 fcen)) ; pulse width (in frequency)
(define-param s .75) ; beam sigma (gaussian beam width)
(define ((my-amp sig) x)
(exp (- (/ (vector3-dot x x) (* 2 sig sig)))))
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(amp-func (my-amp s))
(component Ez)
(center 0 0 -5) (size 1 0 0))))
(define right-side (volume (center 0 0 9) (size 3 0 6)))
(define left-side (volume (center 0 0 3) (size 3 0 6)))
(define mystep (lambda () (change-k-point! (vector3 0 0 0))))
(use-output-directory)
(run-until 10
(at-time 5 mystep)
(to-appended "ez" (at-every 0.1 output-efield-z)))
--Code Finished--
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
