Protein (PDB) viewer written in AWK

plot-protein.awk

#!/usr/bin/awk -f

BEGIN {
  gnuplot = "gnuplot -p";

  print "set terminal wxt noraise" | gnuplot;
  print "set view equal xyz" | gnuplot;
  print "set linetype 1 linecolor palette z linewidth 5" | gnuplot;
  print "unset tics" | gnuplot;
  print "unset border" | gnuplot;
  print "unset colorbox" | gnuplot;
  print "splot '-' notitle with lines" | gnuplot;
}

$1 == "ATOM" && $3 == "CA" {
  printf "%f %f %f\n", $7, $8, $9 | gnuplot;
}

END {
  print "e\n" | gnuplot;
}