#!/usr/local/bin/gnuplot -persist
# set terminal svg size 600,400 dynamic enhanced font 'arial,10' mousing name "map_projection_1" butt dashlength 1.0 
# set output 'map_projection.1.svg'
unset border
set angles degrees
set key at screen 0.98, 0.95 right top vertical Right noreverse enhanced autotitle nobox
set key noinvert samplen 0.3 spacing 1 width 0 height 0 
set size ratio 0.5 1,1
set style data lines
unset xtics
unset ytics
set title "'{/:Bold Winkel tripel}' map projection" 
set xrange [ -180.000 : 180.000 ] noreverse nowriteback
set x2range [ * : * ] noreverse writeback
set yrange [ -90.0000 : 90.0000 ] noreverse nowriteback
set y2range [ * : * ] noreverse writeback
set zrange [ * : * ] noreverse writeback
set cbrange [ * : * ] noreverse writeback
set rrange [ * : * ] noreverse writeback
set lmargin  1
set rmargin  1
set colorbox vertical origin screen 0.9, 0.2 size screen 0.05, 0.6 front  noinvert bdefault
α(λ,φ) = acos(cos(φ) * cos(λ/2.))
x_W3(λ, φ) = 0.5 * (λ*cos(φ1) + (2*cos(φ)*sin(λ/2.))/sinc(α(λ,φ)))
y_W3(λ, φ) = 0.5 * (φ + sin(φ)/sinc(α(λ,φ)))
NO_ANIMATION = 1
π = 3.14159265358979
φ1 = 0.880689235420357
## Last datafile plotted: "world.dat"
plot for [λ=-180:180:10] [φ=-90:90] '+' using (x_W3(λ,φ)):(y_W3(λ,φ)) with lines lc "cyan" lw .5 notitle,      for [φ = -90:90:30] [λ=-180:180:10] '+' using (x_W3(λ,φ)):(y_W3(λ,φ)) with lines lc "cyan" lw .5 notitle,      'world.dat' using (x_W3($1,$2)):(y_W3($1,$2)) with filledcurve fc "brown" title "fill",      'world.dat' using (x_W3($1,$2)):(y_W3($1,$2)) with lines lc "black" title "outline"