back

Equation 2: c_l, c_p, c_s

setenv FR "2 4 3"

echo "1  1  1  1    3  -4  -3" \
 | emap -i - -fr $FR -amb 1 1 1 -rh -o emap.nrrd

alias SOID soid -sphere -res 16 -r 0.01 -o -

alias OFF2EPS off2eps -i - -fr $FR -or -e emap.nrrd -rh \
  -ur -0.95 0.95 -vr -0.95 0.95 -wd 0.0 0.0 1.5 0.7 0.05

SOID -sc 0.85 0.20 0.20 | OFF2EPS -o c_l.eps 
SOID -sc 0.85 0.85 0.20 | OFF2EPS -o c_p.eps 
SOID -sc 0.75 0.75 0.75 | OFF2EPS -o c_s.eps 

rm -f emap.nrrd

foreach S ( l p s )
  eps2ppm c_${S}.eps 200 \
   | unu project -a 0 -m mean \
   | unu resample -s 256 256 -o c_${S}.png
  unu resample -i c_${S}.png -s 90 90 -o _c_${S}.png
end
unu join -i c_{l,p,s}.png -a 0 -o c_lps.png
unu join -i _c_{l,p,s}.png -a 0 -o _c_lps.png