#!/bin/sh -x

cd "`dirname $0`" &&
mkdir -p tmp &&

dataset=data/NED_55723387/ned_55723387/w001001.adf
gdal_translate -of PNM -ot UInt16 -scale 0 400 0 65535 "$dataset" tmp/sf.pgm &&
echo 'P2 1 1 1 1' > tmp/trivial.pgm &&
out/arrange \
  tmp/sf.pgm@ysize=3in,zsize=1in \
  tmp/trivial.pgm@xsize=2in,ysize=0.5in,y=0.5in,x=0.5in,z=0.6in \
  tmp/sf.pgm@xsize=1in,ysize=1in,zsize=0.5in,x=-0.5in,y=2.5in,z=0.5in \
  > tmp/sfbig.pgm &&
out/contour 0.01in tmp/sfbig.pgm > tmp/sfcontour.pgm &&
eog tmp/sfcontour.pgm