Una fuente puntual


;; From the Meep tutorial: plotting permittivity and fields of a straight waveguide

(set! geometry-lattice (make lattice (size 100 100 no-size)))

(set! geometry 
 (list
  (make block 
   (center -9 -9) 
   (size 0.1 0.1)
   (material (make medium (epsilon 1.1)))
  )
 )
)

(set! sources 
 (list
  (make source
   (src (make continuous-src (frequency 0.1)))
   (component Ez)
   (center 0 0)
  )
 )
)

(set! pml-layers (list (make pml (thickness 1.0))))

(set! resolution 10)

(run-until 50
 (at-beginning output-epsilon)
 (at-end output-efield-z)
)
meep meep01.ctl >  meep01.out
h5topng -S4 -Zc dkbluered -a yarg -A meep01-eps-000050.00.h5 meep01-ez-000200.00.h5

________________________________________

(set! geometry-lattice (make lattice (size 100 100 no-size)))

(set! geometry 
 (list
  (make block 
   (center -9 -9) 
   (size 0.1 0.1)
   (material (make medium (epsilon 1.1)))
  )
 )
)

(set! sources 
 (list
  (make source
   (src (make continuous-src (frequency 0.5)))
   (component Ez)
   (center 0 0)
  )
 )
)

(set! pml-layers (list (make pml (thickness 1.0))))

(set! resolution 10)

(run-until 50
 (at-beginning output-epsilon)
 (at-end output-efield-z)
)
jmanza@peje:~/2019/08_2019/fdtd$ meep meep02.ctl > meep02.out
jmanza@peje:~/2019/08_2019/fdtd$ h5topng -S4 -Zc dkbluered -a yarg -A meep02-eps-000000.00.h5 meep02-ez-000050.00.h5

———————————————–

(set! geometry-lattice (make lattice (size 100 100 no-size)))

(set! geometry 
 (list
  (make block 
   (center -9 -9) 
   (size 0.1 0.1)
   (material (make medium (epsilon 1.1)))
  )
 )
)

(set! sources 
 (list
  (make source
   (src (make continuous-src (frequency 0.1)))
   (component Ez)
   (center 0 0)
  )
 )
)

(set! pml-layers (list (make pml (thickness 1.0))))

(set! resolution 10)




(run-until 50
 (at-beginning output-epsilon)
 (to-appended "ez" (at-every 1.0 output-efield-z))
)
meep meep11.ctl > meep11.out
h5topng -t 0:49 -R -Zc dkbluered -a yarg -A meep11-eps-000000.00.h5 meep11-ez.h5
convert meep11*.png meep11.gif

———————————————

(set! geometry-lattice (make lattice (size 100 100 no-size)))

(set! geometry 
 (list
  (make block 
   (center -9 -9) 
   (size 0.1 0.1)
   (material (make medium (epsilon 1.1)))
  )
 )
)

(set! sources 
 (list
  (make source
   (src (make continuous-src (frequency 0.1)))
   (component Ez)
   (center 0 0)
  )
 )
)

(set! pml-layers (list (make pml (thickness 1.0))))

(set! resolution 10)




(run-until 50
 (at-beginning output-epsilon)
 (to-appended "ez" (at-every 1.0 output-efield-z))
)
meep meep12.ctl > meep12.out
jmanza@peje:~/2019/08_2019/fdtd$ h5topng -t 0:49 -R -Zc dkbluered -a yarg -A meep12-eps-000000.00.h5 meep12-ez.h5
jmanza@peje:~/2019/08_2019/fdtd$ convert meep12*.png meep12.gif