Dear Meep users,
I have problem in reflection flux from a slab of
dielectric. I changed my file into very simple example however, it is
not working yet. Could you please let me to know where the problem is?
I have always this comment after run.
creating output file "./1-refl-flux.h5"...
HDF5-DIAG: Error detected in HDF5 (1.8.4) thread 0:
#000: ../../../src/H5S.c line 1331 in H5Screate_simple(): zero sized
dimension for non-unlimited dimension
major: Invalid arguments to routine
minor: Bad value
HDF5-DIAG: Error detected in HDF5 (1.8.4) thread 0:
#000: ../../../src/H5Ddeprec.c line 160 in H5Dcreate1(): not a dataspace ID
major: Invalid arguments to routine
minor: Inappropriate type
meep: Error creating dataset
my ctl file is:
(define-param nSi 4.67 )
(define-param nSiO2 1.432)
(define-param nNb2O5 2.222)
(define-param epsSi (* nSi nSi))
(define-param epsSiO2 (* nSiO2 nSiO2))
(define-param epsGlass (* 1.5 1.5))
(define-param epsNb2O5 (* nNb2O5 nNb2O5))
(define-param w11 .129775)
(define-param w12 .045909) ;Si
(define-param w21 0.038672) ; Nb2O5
(define-param w22 0.060468)
;(define-param w3 2.830)
(define-param N1 19)
(define-param N2 79)
(define-param sy 10000)
(define-param pad 5)
(define-param dpml 2.5) ; 1/0.4(micrometer)
(define-param sx (+ 14 (* 2 ( + pad dpml)))) ; 29
; the computational cell is:
(set! geometry-lattice (make lattice (size sx sy no-size )))
; the block centers:
(define-param c11 -7)
(define-param c2 -2.247)
(define-param c3 3.084)
(define-param c1 -5.331)
(define-param c22 7)
; (define-param c22 -.832)
; the PC geometry:
(define-param PC? false) ; if true, have PC mirror
(set! geometry
( if PC?
(make block (center (+ c11 (/ w12 2)) 0) (size w12 infinity infinity)
(material (make dielectric (epsilon epsSi))))
;(set! geometry '())
(list (make block (center 0 0) (size 14 infinity infinity)
(material
air)))
)
)
(set! pml-layers (list (make pml (thickness dpml))))
(set-param! resolution 10)
; Define source:
(define-param fcen 1.5)
(define-param df 1)
(define-param nfreq 100)
(define trans 0)
(define refl 0)
;(map (lambda (fcen)
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component
Ex)
(center (- dpml (* 0.5 sx) 2.5) 0)
(size 0 sy))))
(set! symmetries (list (make mirror-sym (direction Y) (phase -1))))
(set! trans
(add-flux fcen df nfreq
(make flux-region
(center (+ (* -0.5 sx) dpml 3.5) 0) (size 0 sy))))
(set! refl ; reflected
flux
(add-flux fcen df nfreq
(make flux-region
(center (- dpml (* 0.5 sx) 3) 0) (size 0 sy))))
(if PC? (load-minus-flux "refl-flux" refl))
(run-sources+ 1 (at-beginning output-epsilon))
(if (not PC?) (save-flux "refl-flux" refl))
(print "fcen," fcen " \n")
(display-fluxes trans
refl)
(reset-meep)
Best regards
Mahdieh
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
