diff --git a/doc/figures/grids.eps b/doc/figures/grids.eps
new file mode 100644
index 0000000000000000000000000000000000000000..702a60dc9c3f0626d9a668e3fb634fa2dedfbb8d
--- /dev/null
+++ b/doc/figures/grids.eps
@@ -0,0 +1,435 @@
+%!PS-Adobe-3.0 EPSF-3.0
+%%Title: grids2.fig
+%%Creator: fig2dev Version 3.2 Patchlevel 5d
+%%CreationDate: Wed Apr 10 13:51:09 2013
+%%For: greta@mario.sen.cwi.nl (Margreet Nool)
+%%BoundingBox: 0 0 367 579
+%Magnification: 1.0000
+%%EndComments
+%%BeginProlog
+/$F2psDict 200 dict def
+$F2psDict begin
+$F2psDict /mtrx matrix put
+/col-1 {0 setgray} bind def
+/col0 {0.000 0.000 0.000 srgb} bind def
+/col1 {0.000 0.000 1.000 srgb} bind def
+/col2 {0.000 1.000 0.000 srgb} bind def
+/col3 {0.000 1.000 1.000 srgb} bind def
+/col4 {1.000 0.000 0.000 srgb} bind def
+/col5 {1.000 0.000 1.000 srgb} bind def
+/col6 {1.000 1.000 0.000 srgb} bind def
+/col7 {1.000 1.000 1.000 srgb} bind def
+/col8 {0.000 0.000 0.560 srgb} bind def
+/col9 {0.000 0.000 0.690 srgb} bind def
+/col10 {0.000 0.000 0.820 srgb} bind def
+/col11 {0.530 0.810 1.000 srgb} bind def
+/col12 {0.000 0.560 0.000 srgb} bind def
+/col13 {0.000 0.690 0.000 srgb} bind def
+/col14 {0.000 0.820 0.000 srgb} bind def
+/col15 {0.000 0.560 0.560 srgb} bind def
+/col16 {0.000 0.690 0.690 srgb} bind def
+/col17 {0.000 0.820 0.820 srgb} bind def
+/col18 {0.560 0.000 0.000 srgb} bind def
+/col19 {0.690 0.000 0.000 srgb} bind def
+/col20 {0.820 0.000 0.000 srgb} bind def
+/col21 {0.560 0.000 0.560 srgb} bind def
+/col22 {0.690 0.000 0.690 srgb} bind def
+/col23 {0.820 0.000 0.820 srgb} bind def
+/col24 {0.500 0.190 0.000 srgb} bind def
+/col25 {0.630 0.250 0.000 srgb} bind def
+/col26 {0.750 0.380 0.000 srgb} bind def
+/col27 {1.000 0.500 0.500 srgb} bind def
+/col28 {1.000 0.630 0.630 srgb} bind def
+/col29 {1.000 0.750 0.750 srgb} bind def
+/col30 {1.000 0.880 0.880 srgb} bind def
+/col31 {1.000 0.840 0.000 srgb} bind def
+
+end
+
+/cp {closepath} bind def
+/ef {eofill} bind def
+/gr {grestore} bind def
+/gs {gsave} bind def
+/sa {save} bind def
+/rs {restore} bind def
+/l {lineto} bind def
+/m {moveto} bind def
+/rm {rmoveto} bind def
+/n {newpath} bind def
+/s {stroke} bind def
+/sh {show} bind def
+/slc {setlinecap} bind def
+/slj {setlinejoin} bind def
+/slw {setlinewidth} bind def
+/srgb {setrgbcolor} bind def
+/rot {rotate} bind def
+/sc {scale} bind def
+/sd {setdash} bind def
+/ff {findfont} bind def
+/sf {setfont} bind def
+/scf {scalefont} bind def
+/sw {stringwidth} bind def
+/tr {translate} bind def
+/tnt {dup dup currentrgbcolor
+  4 -2 roll dup 1 exch sub 3 -1 roll mul add
+  4 -2 roll dup 1 exch sub 3 -1 roll mul add
+  4 -2 roll dup 1 exch sub 3 -1 roll mul add srgb}
+  bind def
+/shd {dup dup currentrgbcolor 4 -2 roll mul 4 -2 roll mul
+  4 -2 roll mul srgb} bind def
+/$F2psBegin {$F2psDict begin /$F2psEnteredState save def} def
+/$F2psEnd {$F2psEnteredState restore end} def
+
+/pageheader {
+save
+newpath 0 579 moveto 0 0 lineto 367 0 lineto 367 579 lineto closepath clip newpath
+-214.7 613.3 translate
+1 -1 scale
+$F2psBegin
+10 setmiterlimit
+0 slj 0 slc
+ 0.06000 0.06000 sc
+} bind def
+/pagefooter {
+$F2psEnd
+restore
+} bind def
+%%EndProlog
+pageheader
+%
+% Fig objects follow
+%
+% 
+% here starts figure with depth 50
+% Polyline
+0 slj
+0 slc
+7.500 slw
+ [60] 0 sd
+n 6000 600 m
+ 6000 10200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 7200 600 m
+ 7200 10200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 7800 600 m
+ 7800 10200 l gs col1 s gr  [] 0 sd
+% Polyline
+15.000 slw
+n 4200 1200 m 9000 1200 l 9000 4200 l 4200 4200 l
+ cp gs col0 s gr 
+% Polyline
+7.500 slw
+ [60] 0 sd
+n 6600 600 m
+ 6600 10200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 9000 600 m
+ 9000 10200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 9600 600 m
+ 9600 10200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 3600 1200 m
+ 9600 1200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 3600 3600 m
+ 9600 3600 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 3600 4200 m
+ 9600 4200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 3600 4800 m
+ 9600 4800 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 3600 5400 m
+ 9600 5400 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 3600 6000 m
+ 9600 6000 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 3675 6600 m
+ 9675 6600 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 3600 7200 m
+ 9600 7200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 3600 7800 m
+ 9600 7800 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 3600 8400 m
+ 9600 8400 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 3600 9000 m
+ 9600 9000 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 3600 9600 m
+ 9600 9600 l gs col1 s gr  [] 0 sd
+% Polyline
+15.000 slw
+n 3600 600 m 9600 600 l 9600 10200 l 3600 10200 l
+ cp gs col0 s gr 
+% Polyline
+n 4800 4200 m 8400 4200 l 8400 7200 l 4800 7200 l
+ cp gs col0 s gr 
+% Polyline
+7.500 slw
+ [60] 0 sd
+n 4200 600 m
+ 4200 10200 l gs col1 s gr  [] 0 sd
+% Polyline
+15.000 slw
+n 5400 1500 m 7800 1500 l 7800 2700 l 5400 2700 l
+ cp gs col0 s gr 
+% Polyline
+7.500 slw
+ [60] 0 sd
+n 6300 1200 m
+ 6300 4200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 6900 1200 m
+ 6900 4200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 8700 1200 m
+ 8700 4200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 5100 4200 m
+ 5100 7200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 5700 4200 m
+ 5700 7200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 6300 4200 m
+ 6300 7200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 6900 4200 m
+ 6900 7200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 7500 4200 m
+ 7500 7200 l gs col1 s gr  [] 0 sd
+% Polyline
+n 8546 3900 m 8554 3900 l gs col0 s gr
+% Polyline
+ [60] 0 sd
+n 8100 1200 m
+ 8100 4200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 4200 3300 m
+ 9000 3300 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 4200 2100 m
+ 9000 2100 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 8700 1200 m
+ 8700 4200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 4200 1500 m
+ 9000 1500 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 8100 4200 m
+ 8100 7200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 4800 5100 m
+ 8400 5100 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 3600 5400 m
+ 9600 5400 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 4800 5700 m
+ 8400 5700 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 4800 6300 m
+ 8400 6300 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 4875 6900 m
+ 8475 6900 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 4800 4500 m
+ 8400 4500 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 5400 600 m
+ 5400 10200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 4800 600 m
+ 4800 10200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 7500 1200 m
+ 7500 4200 l gs col1 s gr  [] 0 sd
+% Polyline
+15.000 slw
+n 7800 1800 m 8400 1800 l 8400 3000 l 7800 3000 l
+ cp gs col0 s gr 
+% Polyline
+7.500 slw
+ [60] 0 sd
+n 8400 600 m
+ 8400 10200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 4800 1950 m
+ 8400 1950 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 4800 2250 m
+ 8400 2250 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 4800 2550 m
+ 8400 2550 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 5100 1200 m
+ 5100 4200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 7800 2850 m
+ 8400 2850 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 5400 1650 m
+ 7800 1650 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 3600 1800 m
+ 9600 1800 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 4200 2700 m
+ 9000 2700 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 5700 1200 m
+ 5700 4200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 4200 3900 m
+ 9000 3900 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 4500 1200 m
+ 4500 4200 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 3600 2400 m
+ 9600 2400 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 3600 3000 m
+ 9600 3000 l gs col1 s gr  [] 0 sd
+% Polyline
+15.000 slw
+n 5400 1800 m 4800 1800 l 4800 3000 l 5400 3000 l
+ cp gs col0 s gr 
+% Polyline
+7.500 slw
+ [60] 0 sd
+n 4800 2850 m
+ 5400 2850 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 4950 1800 m
+ 4950 3000 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 5250 1800 m
+ 5250 3000 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 7950 1800 m
+ 7950 3000 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 8250 1800 m
+ 8250 3000 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 7650 1500 m
+ 7650 2700 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 7050 1500 m
+ 7050 2700 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 6750 1500 m
+ 6750 2700 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 6450 1500 m
+ 6450 2700 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 6150 1500 m
+ 6150 2700 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 5850 1500 m
+ 5850 2700 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 5550 1500 m
+ 5550 2700 l gs col1 s gr  [] 0 sd
+% Polyline
+ [60] 0 sd
+n 7350 1500 m
+ 7350 2700 l gs col1 s gr  [] 0 sd
+/Times-Roman ff 300.00 scf sf
+3975 9975 m
+gs 1 -1 sc (a) col0 sh gr
+/Times-Roman ff 300.00 scf sf
+5100 6600 m
+gs 1 -1 sc (ab) col0 sh gr
+/Times-Roman ff 300.00 scf sf
+4500 3825 m
+gs 1 -1 sc (ac) col0 sh gr
+/Times-Roman ff 300.00 scf sf
+4875 2775 m
+gs 1 -1 sc (aca) col0 sh gr
+/Times-Roman ff 300.00 scf sf
+5700 2400 m
+gs 1 -1 sc (acb) col0 sh gr
+/Times-Roman ff 300.00 scf sf
+7850 2775 m
+gs 1 -1 sc (acc) col0 sh gr
+% here ends figure;
+pagefooter
+showpage
+%%Trailer
+%EOF