  COMMAND PROTFUNC
-----------------------------------------------------------------
    Purpose: modify particle-angles algorithmically
      Usage: PROTFUNC(<particleID>,<expression>,<parameter t0>,<"X", "Y" or "Z">)
   Examples: protfunc(1,"cos(t0*(x0+y0))",0,x)
