Pastebin launched a little side project called VERYVIRAL.com, check it out ;-) Want more features on Pastebin? Sign Up, it's FREE!
Guest

Working-Meep-Code-Gratings

By: el_matt on Jan 23rd, 2013  |  syntax: Scheme  |  size: 2.78 KB  |  views: 107  |  expires: Never
download  |  raw  |  embed  |  report abuse  |  print
Text below is selected. Please press Ctrl+C to copy to your clipboard. (⌘+C on Mac)
  1. (use-output-directory "/home/mrea1g12/Dropbox/AR-Coatings/Meep-Files/outputs")
  2.  
  3. (define-param f 2) ; light frequency
  4. (define-param lx 50) ; lattice x limit
  5. (define-param ly 50) ; lattice y limit
  6. (define-param a 4) ; grating periodicity
  7. (define-param w 1) ; grating aperture width
  8. (define-param d 2); grating aperture depth
  9. (define-param T 100); runtime
  10.  
  11. (set! geometry-lattice (make lattice (size lx ly no-size)))
  12.  
  13. (set! geometry (list
  14.                 (make block
  15.                                  (center 25 25)
  16.                                  (size 50 50 no-size)
  17.                                  (material
  18.                                  (make dielectric
  19.                                   (epsilon 2.25)
  20.                                  )
  21.                                  )
  22.                             )
  23.                             (make block
  24.                                  (center 0 0)
  25.                                  (size 35.355 35.355 no-size)
  26.                                  (e1 1 1 0)
  27.                                  (e2 1 -1 0)
  28.                                  (material
  29.                                  (make dielectric
  30.                                   (epsilon 1)
  31.                                  )
  32.                                  )
  33.                             )
  34.                             (make block
  35.                                  (center 12.5 12.5)
  36.                                  (size d w no-size)
  37.                                  (e1 1 1 0)
  38.                                  (e2 1 -1 0)
  39.                                  (material
  40.                                   (make dielectric
  41.                                    (epsilon 1)
  42.                                   )
  43.                                  )
  44.                             )
  45.                             (make block
  46.                                  (center
  47.                                   (+ 12.5 a)
  48.                                   (- 12.5 a)
  49.                                  )
  50.                                  (size d w no-size)
  51.                                  (e1 1 1 0)
  52.                                  (e2 1 -1 0)
  53.                                  (material
  54.                                   (make dielectric
  55.                                    (epsilon 1)
  56.                                   )
  57.                                  )
  58.                             )
  59.                             (make block
  60.                                  (center
  61.                                   (- 12.5 a)
  62.                                   (+ 12.5 a)
  63.                                  )
  64.                                  (size d w no-size)
  65.                                  (e1 1 1 0)
  66.                                  (e2 1 -1 0)
  67.                                  (material
  68.                                   (make dielectric
  69.                                    (epsilon 1)
  70.                                   )
  71.                                  )
  72.                             )
  73.                             (make block
  74.                                  (center
  75.                                   (+ 12.5 (* 2 a))
  76.                                   (- 12.5 (* 2 a))
  77.                                  )
  78.                                  (size d w no-size)
  79.                                  (e1 1 1 0)
  80.                                  (e2 1 -1 0)
  81.                                  (material
  82.                                   (make dielectric
  83.                                    (epsilon 1)
  84.                                   )
  85.                                  )
  86.                             )
  87.                             (make block
  88.                                  (center
  89.                                   (- 12.5 (* 2 a))
  90.                                   (+ 12.5 (* 2 a))
  91.                                  )
  92.                                  (size d w no-size)
  93.                                  (e1 1 1 0)
  94.                                  (e2 1 -1 0)
  95.                                  (material
  96.                                   (make dielectric
  97.                                    (epsilon 1)
  98.                               )
  99.                              )
  100.                             )
  101.                                                     (make block
  102.                                  (center
  103.                                   (+ 12.5 (* 3 a))
  104.                                   (- 12.5 (* 3 a))
  105.                                  )
  106.                                  (size d w no-size)
  107.                                  (e1 1 1 0)
  108.                                  (e2 1 -1 0)
  109.                                  (material
  110.                                   (make dielectric
  111.                                    (epsilon 1)
  112.                                   )
  113.                                  )
  114.                             )
  115.                             (make block
  116.                                  (center
  117.                                   (- 12.5 (* 3 a))
  118.                                   (+ 12.5 (* 3 a))
  119.                                  )
  120.                                  (size d w no-size)
  121.                                  (e1 1 1 0)
  122.                                  (e2 1 -1 0)
  123.                                  (material
  124.                                   (make dielectric
  125.                                    (epsilon 1)
  126.                               )
  127.                              )
  128.                             )
  129.                            )
  130. )
  131.  
  132. (set! sources (list
  133.                                 (make source
  134.                                         (src
  135.                                          (make continuous-src
  136.                                           (frequency f)
  137.                                          )
  138.                                         )
  139.                                 (component Ez)
  140.                                 (center -24 10)
  141.                                 (size 0 24)
  142.                                 )
  143.                            )
  144. )
  145.                                        
  146. (set! pml-layers
  147.                         (list
  148.                          (make pml
  149.                          (thickness 1.0)
  150.                          )
  151.                         )
  152. )
  153.  
  154. (set! resolution 16)
  155.  
  156.  
  157. (run-until T
  158.                         (at-beginning output-epsilon)
  159.                         (at-end
  160.                          (output-png Ez "-Zc dkbluered")
  161.                         )
  162. )