# This file calculates eye safety for Ronja Inferno.
# Worst case scenario is assumed, where the lens is accidentally misfocused
# so the infrared light is concentrated into the viewer's eye or
# binoculars

set terminal postscript portrait enhanced color solid
set output
set style line 1 lt 1 lw 3
set style line 2 lt 2 lw 3
pi=acos(-1)

# Wavelength
lambda=875.0 

# R_lambda is constant from page 7606 row "700-1049" for wavelength
# lambda.
R_lambda=10**((700-lambda)/500)

# K is from page 7605.
K=1.0e4

# Maximum radiant intensities of the LEDs from datasheet
rpmax_4220=0.06 # W/sr
rpmax_4230=0.131 # W/sr

# Length of slide path to calculate misfocused operation
pathlen=0.04 # 4cm

# Lens gains. Measured on real lens samples using
# Ronja benchpress, which is a measurement bench for
# lens optical gain measurements.
# See http://ronja.twibright.com/benchpress
# In dB, HSDL 4220 LED
gain_60_4220=25.56
gain_75_4220=22.51
gain_80_4220=27.3
gain_90_4220=27.2
gain_100_4220=28.81
gain_130_4220=29.41

# In dB, HSDL 4230 LED
gain_60_4230=22.82
gain_75_4230=19.54
gain_80_4230=24.02
gain_90_4230=27.24
gain_100_4230=25.91
gain_130_4230=27.21

# Focal lengths of lens samples, measured.
# In meters
flen_60=0.211
flen_75=0.146
flen_80=0.231
flen_90=0.244
flen_100=0.241
flen_130=0.312

# Diameters of lens samples, measured.
# In meters
dia_60=0.06
dia_75=0.075
dia_80=0.08
dia_90=0.09
dia_100=0.1
dia_130=0.13

h=3600 # 1 hour
day=24*h
week=7*day
month=30*day
year=day*365.25

# Magnification of binoculars. 1 means naked eye.
mag=1

# Minimum distance into whis it is possible to misfocus
# the transmitter
mindist_true(f)=1/(1/f-1/(f+pathlen))

mindist(f)=floor(mindist_true(f)*100+0.5)/100.0

# Area of the lens
area(dummy)=pi*aperture**2/4

# Radiant intensity after Ronja output lens
intensity(dummy)=rpmax*10**(gain/10)

# We use radiance here and not spectral radiance, because the source
# has narrow spectrum.
# For definition of radiance see http://en.wikipedia.org/wiki/Radiance
radiance(dummy)=intensity(0)/area(0)

# Viewing angle of the source on binocular output
angle(l)=aperture/l*mag

# Maximum permitted time according to 3 (b), page 7605
tmax(l)=K**2/((angle(l)*radiance(0)*R_lambda)**2)

# For definition of irradiance see http://en.wikipedia.org/wiki/Irradiance
irradiance(l)=intensity(0)/l**2*mag**2

set xlabel "Distance of observer from transmitter [m]"
set ylabel "Time [s]"
set xtics nomirror
set grid xtics mxtics ytics mytics lt 7, lt 0
set yrange[1e3:1e11]
set y2range [3:11]
set y2tics ("hour" log10(h), \
"2 hours" log10(2*h), "1/2 hour" log10(h/2), "day" log10(day), \
"6 hours" log10(6*h), "12 hours" log10(12*h), "week" log10(week), \
"month" log10(month), "year" log10(year), "10 years" log10(10*year), \
"100 years" log10(100*year), "minute" log10(60), "1/4 hour" log10(h/4), \
"10 minutes" log10(600), "5 minutes" log10(300), "2 minutes" log10(120), \
"1000 years" log10(1000*year))
set logscale xy
set mytics 10
set tics out
set ytics nomirror

rpmax=rpmax_4220
gain=gain_60_4220
flen=flen_60 # In meters
aperture=dia_60
set title "Par. 3b): Max. allowed exposure time for naked eye\n60mm loupe, HSDL4220"
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)

plot tmax(x) title "Max time" ls 1, \
h title "" ls 2, \
h/2 title "" ls 2, \
h/4 title "" ls 2, \
day title "" ls 2, \
2*h title "" ls 2, \
6*h title "" ls 2, \
12*h title "" ls 2, \
week title "" ls 2, \
month title "" ls 2, \
year title "" ls 2, \
10*year title "" ls 2, \
100*year title "" ls 2, \
1000*year title "" ls 2, \
60 title "" ls 2, \
120 title "" ls 2, \
600 title "" ls 2, \
300 title "" ls 2

gain=gain_60_4230
rpmax=rpmax_4230
set title "Par. 3b): Max. allowed exposure time for naked eye\n\
60mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_75_4220
flen=flen_75 # In meters
aperture=dia_75
set title "Par. 3b): Max. allowed exposure time for naked eye\n\
75mm loupe, HSDL4220"
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)
replot

gain=gain_75_4230
rpmax=rpmax_4230
set title "Par. 3b): Max. allowed exposure time for naked eye\n\
75mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_80_4220
flen=flen_80 # In meters
aperture=dia_80
set title "Par. 3b): Max. allowed exposure time for naked eye\n\
80mm loupe, HSDL4220"
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)
replot

gain=gain_80_4230
rpmax=rpmax_4230
set title "Par. 3b): Max. allowed exposure time for naked eye\n\
80mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_90_4220
flen=flen_90 # In meters
aperture=dia_90
set title "Par. 3b): Max. allowed exposure time for naked eye\n\
90mm loupe, HSDL4220"
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)
replot

gain=gain_90_4230
rpmax=rpmax_4230
set title "Par. 3b): Max. allowed exposure time for naked eye\n\
90mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_100_4220
flen=flen_100 # In meters
aperture=dia_100
set title "Par. 3b): Max. allowed exposure time for naked eye\n\
100mm loupe, HSDL4220"
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)
replot

gain=gain_100_4230
rpmax=rpmax_4230
set title "Par. 3b): Max. allowed exposure time for naked eye\n\
100mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_130_4220
flen=flen_130 # In meters
aperture=dia_130
set title "Par. 3b): Max. allowed exposure time for naked eye\n\
130mm loupe, HSDL4220"
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)
replot

gain=gain_130_4230
rpmax=rpmax_4230
set title "Par. 3b): Max. allowed exposure time for naked eye\n\
130mm loupe, HSDL4230"
replot

# 20x magnifying telescope
mag=20

set xrange [10:2000]
set x2range [10:2000]
set x2tics (2000)
set yrange [1:1e9]
set y2range [0:9]

rpmax=rpmax_4220
gain=gain_60_4220
flen=flen_60 # In meters
aperture=dia_60
set title "Par. 3b): Max. allowed exposure time for 20x telescope\n\
60mm loupe, HSDL4220"
replot

gain=gain_60_4230
rpmax=rpmax_4230
set title "Par. 3b): Max. allowed exposure time for 20x telescope\n\
60mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_75_4220
flen=flen_75 # In meters
aperture=dia_75
set title "Par. 3b): Max. allowed exposure time for 20x telescope\n\
75mm loupe, HSDL4220"
replot

gain=gain_75_4230
rpmax=rpmax_4230
set title "Par. 3b): Max. allowed exposure time for 20x telescope\n\
75mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_80_4220
flen=flen_80 # In meters
aperture=dia_80
set title "Par. 3b): Max. allowed exposure time for 20x telescope\n\
80mm loupe, HSDL4220"
replot

gain=gain_80_4230
rpmax=rpmax_4230
set title "Par. 3b): Max. allowed exposure time for 20x telescope\n\
80mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_90_4220
flen=flen_90 # In meters
aperture=dia_90
set title "Par. 3b): Max. allowed exposure time for 20x telescope\n\
90mm loupe, HSDL4220"
replot

gain=gain_90_4230
rpmax=rpmax_4230
set title "Par. 3b): Max. allowed exposure time for 20x telescope\n\
90mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_100_4220
flen=flen_100 # In meters
aperture=dia_100
set title "Par. 3b): Max. allowed exposure time for 20x telescope\n\
100mm loupe, HSDL4220"
replot

gain=gain_100_4230
rpmax=rpmax_4230
set title "Par. 3b): Max. allowed exposure time for 20x telescope\n\
100mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_130_4220
flen=flen_130 # In meters
aperture=dia_130
set title "Par. 3b): Max. allowed exposure time for 20x telescope\n\
130mm loupe, HSDL4220"
replot

gain=gain_130_4230
rpmax=rpmax_4230
set title "Par. 3b): Max. allowed exposure time for 20x telescope\n\
130mm loupe, HSDL4230"
replot

mag=1
set yrange [1e-6:1e3]
set xrange [mindist(flen):]
set format y "%g"
set noy2tics
set tics in
set ytics mirror
set ylabel "Radiant flux density [W*m^{-2}]

rpmax=rpmax_4220
gain=gain_60_4220
flen=flen_60 # In meters
aperture=dia_60
set title "Par. 3c1): Radiant flux density for naked eye\n\
60mm loupe, HSDL4220"
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)
plot irradiance(x) title "Radiant flux density" ls 1, \
 100 title "Par. 3c1) limit" ls 2

rpmax=rpmax_4230
gain=gain_60_4230
set title "Par. 3c1): Radiant flux density for naked eye\n\
60mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_75_4220
flen=flen_75 # In meters
aperture=dia_75
set title "Par. 3c1): Radiant flux density for naked eye\n\
75mm loupe, HSDL4220"
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)
replot

rpmax=rpmax_4230
gain=gain_75_4230
set title "Par. 3c1): Radiant flux density for naked eye\n\
75mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_80_4220
flen=flen_80 # In meters
aperture=dia_80
set title "Par. 3c1): Radiant flux density for naked eye\n\
80mm loupe, HSDL4220"
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)
replot

rpmax=rpmax_4230
gain=gain_80_4230
set title "Par. 3c1): Radiant flux density for naked eye\n\
80mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_90_4220
flen=flen_90 # In meters
aperture=dia_90
set title "Par. 3c1): Radiant flux density for naked eye\n\
90mm loupe, HSDL4220"
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)
replot

rpmax=rpmax_4230
gain=gain_90_4230
set title "Par. 3c1): Radiant flux density for naked eye\n\
90mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_100_4220
flen=flen_100 # In meters
aperture=dia_100
set title "Par. 3c1): Radiant flux density for naked eye\n\
100mm loupe, HSDL4220"
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)
replot

rpmax=rpmax_4230
gain=gain_100_4230
set title "Par. 3c1): Radiant flux density for naked eye\n\
100mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_130_4220
flen=flen_130 # In meters
aperture=dia_130
set title "Par. 3c1): Radiant flux density for naked eye\n\
130mm loupe, HSDL4220"
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)
replot

rpmax=rpmax_4230
gain=gain_130_4230
set title "Par. 3c1): Radiant flux density for naked eye\n\
130mm loupe, HSDL4230"
replot

mag=20
set xrange [10:2000]
set x2range [10:2000]
set x2tics (2000)
set yrange [1e-3:1e3]

rpmax=rpmax_4220
gain=gain_60_4220
flen=flen_60 # In meters
aperture=dia_60
set title "Par. 3c1): Radiant flux density for 20x telescope\n\
60mm loupe, HSDL4220"
plot irradiance(x) title "Radiant flux density" ls 1, \
 100 title "Par. 3c1) limit" ls 2

rpmax=rpmax_4230
gain=gain_60_4230
set title "Par. 3c1): Radiant flux density for 20x telescope\n\
60mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_75_4220
flen=flen_75 # In meters
aperture=dia_75
set title "Par. 3c1): Radiant flux density for 20x telescope\n\
75mm loupe, HSDL4220"
replot

rpmax=rpmax_4230
gain=gain_75_4230
set title "Par. 3c1): Radiant flux density for 20x telescope\n\
75mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_80_4220
flen=flen_80 # In meters
aperture=dia_80
set title "Par. 3c1): Radiant flux density for 20x telescope\n\
80mm loupe, HSDL4220"
replot

rpmax=rpmax_4230
gain=gain_80_4230
set title "Par. 3c1): Radiant flux density for 20x telescope\n\
80mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_90_4220
flen=flen_90 # In meters
aperture=dia_90
set title "Par. 3c1): Radiant flux density for 20x telescope\n\
90mm loupe, HSDL4220"
replot

rpmax=rpmax_4230
gain=gain_90_4230
set title "Par. 3c1): Radiant flux density for 20x telescope\n\
90mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_100_4220
flen=flen_100 # In meters
aperture=dia_100
set title "Par. 3c1): Radiant flux density for 20x telescope\n\
100mm loupe, HSDL4220"
replot

rpmax=rpmax_4230
gain=gain_100_4230
set title "Par. 3c1): Radiant flux density for 20x telescope\n\
100mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_130_4220
flen=flen_130 # In meters
aperture=dia_130
set title "Par. 3c1): Radiant flux density for 20x telescope\n\
130mm loupe, HSDL4220"
replot

rpmax=rpmax_4230
gain=gain_130_4230
set title "Par. 3c1): Radiant flux density for 20x telescope\n\
130mm loupe, HSDL4230"
replot

set ylabel "Radiance [W*m^{-2}*sr^{-1}]"
set yrange [1e3:1e9]
limit_3c2(l)=0.6e4/angle(l)
mag=1

rpmax=rpmax_4220
gain=gain_60_4220
flen=flen_60 # In meters
aperture=dia_60
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)
set title "Par. 3c2) Radiance for naked eye\n\
60mm loupe, HSDL4220"
plot radiance(0) title "Radiance" ls 1, \
limit_3c2(x) title "Limit by par. 3c2)" ls 2

rpmax=rpmax_4230
gain=gain_60_4230
set title "Par. 3c2) Radiance for naked eye\n\
60mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_75_4220
flen=flen_75 # In meters
aperture=dia_75
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)
set title "Par. 3c2) Radiance for naked eye\n\
75mm loupe, HSDL4220"
replot

rpmax=rpmax_4230
gain=gain_75_4230
set title "Par. 3c2) Radiance for naked eye\n\
75mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_80_4220
flen=flen_80 # In meters
aperture=dia_80
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)
set title "Par. 3c2) Radiance for naked eye\n\
80mm loupe, HSDL4220"
replot

rpmax=rpmax_4230
gain=gain_80_4230
set title "Par. 3c2) Radiance for naked eye\n\
80mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_90_4220
flen=flen_90 # In meters
aperture=dia_90
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)
set title "Par. 3c2) Radiance for naked eye\n\
90mm loupe, HSDL4220"
replot

rpmax=rpmax_4230
gain=gain_90_4230
set title "Par. 3c2) Radiance for naked eye\n\
90mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_100_4220
flen=flen_100 # In meters
aperture=dia_100
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)
set title "Par. 3c2) Radiance for naked eye\n\
100mm loupe, HSDL4220"
replot

rpmax=rpmax_4230
gain=gain_100_4230
set title "Par. 3c2) Radiance for naked eye\n\
100mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_130_4220
flen=flen_130 # In meters
aperture=dia_130
set xrange [mindist(flen):2000]
set x2range [mindist(flen):2000]
set x2tics (mindist(flen), 2000)
set title "Par. 3c2) Radiance for naked eye\n\
130mm loupe, HSDL4220"
replot

rpmax=rpmax_4230
gain=gain_130_4230
set title "Par. 3c2) Radiance for naked eye\n\
130mm loupe, HSDL4230"
replot

mag=20
set xrange [10:2000]
set x2range [10:2000]
set x2tics (2000)
set yrange [:1e8]

rpmax=rpmax_4220
gain=gain_60_4220
flen=flen_60 # In meters
aperture=dia_60
set title "Par. 3c2) Radiance for 20x telescope\n\
60mm loupe, HSDL4220"
plot radiance(0) title "Radiance" ls 1, \
limit_3c2(x) title "Limit by par. 3c2)" ls 2

rpmax=rpmax_4230
gain=gain_60_4230
set title "Par. 3c2) Radiance for 20x telescope\n\
60mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_75_4220
flen=flen_75 # In meters
aperture=dia_75
set title "Par. 3c2) Radiance for 20x telescope\n\
75mm loupe, HSDL4220"
replot

rpmax=rpmax_4230
gain=gain_75_4230
set title "Par. 3c2) Radiance for 20x telescope\n\
75mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_80_4220
flen=flen_80 # In meters
aperture=dia_80
set title "Par. 3c2) Radiance for 20x telescope\n\
80mm loupe, HSDL4220"
replot

rpmax=rpmax_4230
gain=gain_80_4230
set title "Par. 3c2) Radiance for 20x telescope\n\
80mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_90_4220
flen=flen_90 # In meters
aperture=dia_90
set title "Par. 3c2) Radiance for 20x telescope\n\
90mm loupe, HSDL4220"
replot

rpmax=rpmax_4230
gain=gain_90_4230
set title "Par. 3c2) Radiance for 20x telescope\n\
90mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_100_4220
flen=flen_100 # In meters
aperture=dia_100
set title "Par. 3c2) Radiance for 20x telescope\n\
100mm loupe, HSDL4220"
replot

rpmax=rpmax_4230
gain=gain_100_4230
set title "Par. 3c2) Radiance for 20x telescope\n\
100mm loupe, HSDL4230"
replot

rpmax=rpmax_4220
gain=gain_130_4220
flen=flen_130 # In meters
aperture=dia_130
set title "Par. 3c2) Radiance for 20x telescope\n\
130mm loupe, HSDL4220"
replot

rpmax=rpmax_4230
gain=gain_130_4230
set title "Par. 3c2) Radiance for 20x telescope\n\
130mm loupe, HSDL4230"
replot