lda0 1[um] "Center wavelength"
f0 c_const/lda0 "Center frequency"
omega0 2*pi*f0 "Center angular frequency"
n0 1 "Refractive index, superstrate"
nH 2.5 "Refractive index, high-index layers"
nL 1.5 "Refractive index, low-index layers"
ns 1.5 "Refractive index, substrate"
E0 1[V/m] "Electric field amplitude"
T0 1/f0 "Period for center frequency"
df0 f0*4/pi*asin((nH-nL)/(nH+nL)) "Reflector stop band"
d0 lda0/(4*n0) "Thickness, superstrate"
dH lda0/(4*nH) "Thickness, high-index layers"
dL lda0/(4*nL) "Thickness, low-index layers"
dPeriod dH+dL "Thickness, pair of low- and high-index layers"
ds lda0/(4*ns) "Thickness, substrate"
NPeriod 20 "Number of periods"
L_tot NPeriod*dPeriod+dH "Total thickness"
deltaf 1[THz] "Frequency resolution"
Tend 1/deltaf "End time"
Tc 3*Td "Delay time"
Td 1/(2*f0) "Pulse duration"
Tmax max(200*Td,Tend) "Simulation end time"