our %vxname2pos; # $vxname2pos{VXNAME} = Math::GSL::Vector
our %bb; # $r{Bounding...}[0..3] = points
our $maxpaper = 'a3';
+my $a1sfactor = 1/sqrt(8);
sub TAU { M_PI * 2.0; }
sub MM2PT { 72.0 / 25.4; }
sub prs2r ($) {
# bodgery for the sea and land
- return 'NZ' if m/ \| / && !grep { !m{^(?:NZ|L2?)$} } split / \| /, $_[0];
+ return $c{Sea} if m/ \| / && !grep { !m{^$c{OuterfaceRe}$} } split / \| /, $_[0];
$prs2region{$_[0]} // confess "@_ ?"
}
sub prs2rr ($) { $region{prs2r($_[0])} }
my @edges;
my $process_rr = sub {
my $last_ai;
- if ($rr->{Name} eq 'NZ') {
- # We combined L and NZ; now we must split them again
- # in fact, we just throw away L and L2 entirely.
- # Here, when processing NZ, we keep only vertices that
- # are part of NZ.
+ if ($rr->{Name} eq $c{Sea}) {
+ # We combined land and sea; now we must split them again
+ # in fact, we just throw away all L entirely.
+ # Here, when processing sea, we keep only vertices that
+ # are part of the sea.
@edges = grep {
my $evxname = $_->{VxName};
my @eregions = split / \| ?/, $evxname;
- grep { $_ eq 'NZ' } @eregions
+ grep { $_ eq $c{Sea} } @eregions
} @edges;
}
$rr->{Polygon} = [ @edges ];
- if ($rr->{Name} ne 'NZ') {
+ if ($rr->{Name} ne $c{Sea}) {
for my $ei (0..$#edges) {
if (!defined $last_ai) {
my $ai;
}
$process_rr->();
- $region{$_}{Special} = 1 foreach qw(NZ L L2);
+ $region{$_}{Special} = 1 foreach @{ $c{SpecialRegions} };
}
#----- geometry mangling -----
$mid = ($lims[1] - $lims[0]) * 0.5;
}
-sub adjust_nz() {
- my $poly = $region{NZ}{Polygon};
+sub adjust_sea() {
+ # Adjusts the top and bottom edges of the sea
+ my $poly = $region{$c{Sea}}{Polygon};
my %occurs;
foreach my $v (@$poly) {
$occurs{$v->{VxName}}++
}
foreach my $v (@$poly) {
- $v->{NZ_Occurs} = $occurs{$v->{VxName}};
+ $v->{SeaOccurs} = $occurs{$v->{VxName}};
}
- #print STDERR " nz ".(scalar @$poly)."\n";
+ #print STDERR " sea ".(scalar @$poly)."\n";
my $occurs_m = sub {
my ($ei) = @_;
- $poly->[($ei+@$poly) % @$poly]{NZ_Occurs} > 1
+ $poly->[($ei+@$poly) % @$poly]{SeaOccurs} > 1
};
for my $ei (0..$#$poly) {
next if $occurs_m->($ei);
next unless $occurs_m->($ei-1);
next unless $occurs_m->($ei+1);
- #print STDERR " nz would adjust $poly->[$ei]{VxName}\n";
+ #print STDERR " sea would adjust $poly->[$ei]{VxName}\n";
# adjust coord outwards
}
for my $ei (0..$#$poly) {
next unless $occurs_m->($ei);
- #print STDERR " nz occurs_m $ei $poly->[$ei]{VxName}\n";
+ #print STDERR " sea occurs_m $ei $poly->[$ei]{VxName}\n";
next unless $occurs_m->($ei-1);
my $dleft = mkvec(-1000,0);
foreach my $which (0,-1) {
my $posj = \ $poly->[$ej]{Pos};
my $posk = \ $poly->[$ek]{Pos};
my $dout = $$posk - $$posj;
- #print STDERR " nz adj j=$ej ",prvec($$posj)," $poly->[$ej]{VxName}\n";
- #print STDERR " nz adj k=$ek ",prvec($$posk)," $poly->[$ek]{VxName}\n";
+ #print STDERR " sea adj j=$ej ",prvec($$posj)," $poly->[$ej]{VxName}\n";
+ #print STDERR " sea adj k=$ek ",prvec($$posk)," $poly->[$ek]{VxName}\n";
if ($poly->[$ej]{VxName} =~ /Noorderzijlvest/) {
- #print STDERR " nz adj extra\n";
+ #print STDERR " sea adj extra\n";
$$posk = $$posj + 10 * mkvec(100,-50);
}
$$posj = $$posj + $dleft + $dout*10;
@{ $bb{PrintableT} }[2..3],
];
} elsif ($maxpaper =~ m/^a1/) {
- my $offx = 0.5*($a1_w - $bb{Entire}[2]);
+ my $offx = ($a1_w - $bb{Entire}[2] * (1 + $a1sfactor)) / 3;
my $offy = 0.5*($a1_h - $bb{Entire}[3]);
$somepage->('P', -$offx,-$offy, $a1_w-$offx, $a1_h-$offy);
my $hairs = 30 * MM2PT;
/pgsave save def
&&EndPageSetup
all
+END
+ if ($maxpaper =~ m/^a1/) {
+ o(<<END);
+ $a1_w 0 translate
+ $a1sfactor dup scale
+ $a1_w neg 0 translate
+ 1 1 3 { pop all 0 $a1_h 0.8 mul translate } for
+END
+ }
+ o_amble(<<END);
pgsave restore
showpage
&&EOF
}
}
-sub afsluitdijk() {
- my ($adj) = grep { $_->{Name} eq 'NZ' } @{ $region{'ZZ'}{Adj} };
+sub afsluitdijk ($$) {
+ my ($inner, $comment) = @_;
+ my ($adj) = grep { $_->{Name} eq $c{Sea} } @{ $region{$inner}{Adj} };
my $dash = 20;
my (@pos) = adj_ends($adj);
my $mid = 0.5 * ($pos[0] + $pos[1]);
my $r = 0.5 * ($pos[1] - $pos[0])->norm();
- o(sprintf <<'END', facepscolour('sea'), psvec($mid), $r);
- %% Afsluitdijk
+ o(sprintf <<'END', $comment, facepscolour('sea'), psvec($mid), $r);
+ %% %s
%s %s %s 0 360 arc closepath fill
END
some_edge_faces(()); # just sets the colour
one_hs('purple', $mid);
}
+sub afsluitdijken() {
+ afsluitdijk($_->{Region}, $_->{Comment}) foreach @{ $c{DashedHsInners} };
+}
+
sub edge_faces() {
- some_edge_faces($region{'NZ'});
- afsluitdijk();
+ some_edge_faces($region{$c{Sea}});
+ afsluitdijken();
some_edge_faces(grep {
$_->{Name} !~ m{^(?:L2?|NZ|ZZ)$}
} values %region);
next if $rr->{Name} =~ m/^L\d?$/;
my $words = $rr->{DisplayName};
my $pos;
- if ($rr->{Name} eq 'NZ') {
+ if ($rr->{Name} eq $c{Sea}) {
$pos = mkvec(90, 1280);
} else {
$pos = $rr->{Centre} + mkvec( 0, 0.5 * $sz * @$words );
$pos += $sz * &mkvec(@$adjust);
}
my $w = $rr->{Water};
- $w = 2 if $rr->{Special} || $rr->{Name} eq 'ZZ';
+ $w = 2 if $rr->{Name} eq $c{Sea};
if ($w) {
#print STDERR " water $w in $rr->{Name}\n";
my $down = water_cubes_sz(2) - water_cubes_sz(1) + $cube_gap;
$pos -= mkvec(0, -0.5 * $down);
my $cube_y = -($sz * @$words);
- if ($rr->{Name} eq 'Wieringermeer') {
+ if (grep { $_ eq $rr->{Name} } @{ $c{CubeAbove} }) {
$pos -= mkvec(0,+$sz);
$cube_y = 0;
$down = 0;
read_faces();
calculate_centres();
transform_coordinates();
-adjust_nz();
+adjust_sea();
computeboundings();
preamble();