echo "-75 75   -75 75   -75 75   -75 75" \
 | unu reshape -s 2 2 2 \
 | unu resample -s 151 151 151 -k tent -c node \
 | unu 2op pow - 2 -o xx

unu swap -i xx -a 0 1 -o yy
unu swap -i xx -a 0 2 -o zz

unu 3op + xx yy zz \
 | unu axinfo -a 0 -c cell \
 | unu axinfo -a 1 -c cell \
 | unu axinfo -a 2 -c cell -o tmp.nhdr

echo "space dimension: 3" >> tmp.nhdr
echo "space origin: (-75,-75,-75)" >> tmp.nhdr
echo "space directions: (1,0,0) (0,1,0) (0,0,1)" >> tmp.nhdr
echo "kinds: space space space" >> tmp.nhdr

../../bin/dnorm -i tmp.nhdr -o parab3-150.nrrd

unu resample -i parab3-150.nrrd -s 30 =  =  -o parab3-x30.nrrd
unu resample -i parab3-150.nrrd -s =  30 =  -o parab3-y30.nrrd
unu resample -i parab3-150.nrrd -s =  =  30 -o parab3-z30.nrrd




echo "-74.5 74.5   -74.5 74.5" \
 | unu reshape -s 2 2 \
 | unu resample -s 150 150 -k tent -c node \
 | unu 2op pow - 2 -o xx

unu swap -i xx -a 0 1 -o yy

unu 2op + xx yy \
 | unu axinfo -a 0 -c cell \
 | unu axinfo -a 1 -c cell -o tmp.nhdr

echo "space dimension: 2" >> tmp.nhdr
echo "space origin: (-74.5,-74.5)" >> tmp.nhdr
echo "space directions: (1,0) (0,1)" >> tmp.nhdr
echo "kinds: space space" >> tmp.nhdr

../../bin/dnorm -i tmp.nhdr -o parab2-150.nrrd

unu resample -i parab2-150.nrrd -s 30 =  -o parab2-x30.nrrd
unu resample -i parab2-150.nrrd -s =  30 -o parab2-y30.nrrd

rm -f xx yy zz tmp.nhdr tmp.raw