Advertisement
el_matt

Working-Meep-Code-Gratings

Jan 23rd, 2013
252
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Scheme 2.78 KB | None | 0 0
  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. )
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement