#!/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"