#! /bin/csh -f -e

/bin/rm -f caquinhos-2-*-mag.pgm
cat caquinhos-2.pgm | pnmcut 030 100 60 60 | pnmflip -rotate180  > caquinhos-2-a-mag.pgm
cat caquinhos-2.pgm | pnmcut 300 400 60 60 | pnmflip -rotate270  > caquinhos-2-b-mag.pgm
pnmcat -lr caquinhos-2-{a,b}-mag.pgm > caquinhos-2-ab-mag.pgm
xv caquinhos-2-ab-mag.pgm caquinhos-2-?-mag.pgm &

set grey = "112"

/bin/rm -f caquinhos-2-????.pgm
PZSplit \
  -inFile caquinhos-2 \
  -outFile caquinhos-2 \
  -black ${grey}
ls -l caquinhos-2-????.pgm
xv caquinhos-2-????.pgm caquinhos-2.pgm &

foreach f ( 0000 0001 )
  PZBoundary \
    -inFile  caquinhos-2-$f \
    -black ${grey} \
    -grey ${grey} \
    -outFile caquinhos-2-$f
  PZDraw \
    -window  caquinhos-2-$f \
    -outName caquinhos-2-$f \
    -grid 50 \
    -eps \
    < caquinhos-2-$f.flc
  ghostview -openwindows caquinhos-2-$f.eps
  PZFilter \
    -inFile caquinhos-2-$f \
    -outFile caquinhos-2-$f-f \
    -lMin 12 \
    -lMax 48 \
    -nIter 20 \
    -step 0.5

  if ( -r caquinhos-2-$f-d.wnd ) then
    foreach g ( $f $f-f )
      PZDraw \
        -window  caquinhos-2-$f-d \
        -outName caquinhos-2-$g-d \
        -eps \
        < caquinhos-2-$g.flc
      ghostview -openwindows caquinhos-2-$g-d.eps
    end
  endif

  PZInvar \
    -inFile caquinhos-2-$f-f \
    -outFile caquinhos-2-$f-c \
    -squeeze 100.0
  /bin/rm -f caquinhos-2-$f-c.plt
  cat caquinhos-2-$f-c.flc \
  | egrep -e '^[-0-9]' \
  > caquinhos-2-$f-c.plt
  show-curvature-plot caquinhos-2-$f-c
end

PZMatrix \
  -xChain caquinhos-2-0000-c \
  -yChain caquinhos-2-0001-c \
  -outFile caquinhos-2-0000x0001-${grey}

xv caquinhos-2-0000x0001-${grey}.pgm &