mponent pt decay-by
В· 50 - dT - через які проміжки часу потрібно проводити вимірювання величини ...
В· Ez - component
В· (if no-bend?
(vector3 (- dmin 2) l_place)
(vector3 (+ dmax 2) l_place)) - pt - ... щоб її значення в одній з розрахункових точок (за наявності або відсутності решітки) ...
В· 1e-3)) - decay-by - ... було б менше заданої величини (у нашому випадку Ez <0.001)
У другому способі, можна вказати скільки разів треба провести розрахунки поля після вимкнення джерела випромінювання. Приклад:
(run-sources + 500) - розрахунок припиниться через 500 тимчасових інтервалів після вимкнення джерела.
6) Всі розрахунки проведені і можна виконати висновок результатів.
В· (display-fluxes trans)
Результат ::, 0.35, 5.58516844105163e-5:, 0.351010101010101, 8.45812001141722e-5:, 0.352020202020202, 1.27883318656805e-4:, 0.353030303030303, 1.90004998820637e-4:, 0.354040404040404, 2.75747528407135e-4 :, 0.355050505050505, 3.92351676416799e-4:, 0.356060606060606, 5.52562540438976e-4:, 0.357070707070707, 7.76885293373528e-4:, 0.358080808080808, 0.001092440293566721:, 0.359090909090909, 0.00152809295946778
...
У першій колонці задається частота випромінювання, у другій - кількість минулої енергії.
Розраховувати можна не тільки минуле крізь структуру кількість енергії. Також, наприклад, можна розрахувати величину відбитої від структури енергії або величину втрат випромінювання. Все залежить від поставлених цілей і вимагає деяких додаткових математичних операцій
7) Зробимо аналіз отриманих результатів:
а) Запуск програми: Запуск програми проводиться два рази. У перший раз без кристалічної структури:
meep no-bend? = true flux1.ctl | tee bend0.out
Вдруге - зі структурою:
meep flux1.ctl | tee bend.out
В обох випадках результати перенаправляються у файл. Bend0.out - зберігає значення енергій падаючого випромінювання; bend.out - зберігає значення енергій минулого випромінювання. p align="justify"> Програма запускається два рази, щоб не відбулося накладення хвиль, що залишилися після першого запуску, на хвилі при другому запуску поля.
б) Далі переведемо результати програми у формат dat, який сприймається такими додатками як Excel і Matlab:
grep flux1: bend0.out> bend0.datflux1: bend.out> b...
