set S="38 55 76"

echo "-1 1   -1 1   -1 1   -1 1" \
 | unu reshape -s 2 2 2 \
 | unu resample -s $S -k tent -c node \
 | unu 2op pow - 2 -o xx.nrrd

echo "-1 1   -1 1   -1 1   -1 1" \
 | unu reshape -s 2 2 2 \
 | unu swap -a 0 1 \
 | unu resample -s $S -k tent -c node \
 | unu 2op pow - 2 -o yy.nrrd

echo "-1 1   -1 1   -1 1   -1 1" \
 | unu reshape -s 2 2 2 \
 | unu swap -a 0 2 \
 | unu resample -s $S -k tent -c node \
 | unu 2op pow - 2 -o zz.nrrd

unu 3op + {xx,yy,zz}.nrrd \
 | unu 1op sqrt \
 | unu 3op gauss - 0.5 0.15 \
 | unu axinfo -a 0 -mm -1 1 \
 | unu axinfo -a 1 -mm -1 1 \
 | unu axinfo -a 2 -mm -1 1 \
 | ../../bin/dnorm -i - -o - \
 | unu pad -min -2 -2 -2 -max M+2 M+2 M+2 -o sphere.nrrd

rm -f {xx,yy,zz}.nrrd