Cut a cylinder out of the current configuration.

<p>Usage: <tt> input = [ <i>zind yind sx0 sy0 rrem flag szmin szmax</i> ]  makecylinder</tt>

<p>
<tt><i>zind</i></tt> = 1, 2 or 3 specifies the direction of the cylinder axis (to be
parallel to <tt><b>c</b><sub>1</sub></tt>, <tt><b>c</b><sub>2</sub></tt>, 
or <tt><b>c</b><sub>3</sub></tt>).
<br>
<tt><i>yind</i></tt> = 1, 2 or 3 specifies a direction perpendicular to the cylinder axis.
This input simplies the implementation of this function --- the code does not have to guess
the coordinate system.
<br>
<tt><i>sx0, sy0</i></tt> specifies the location of the cylinder axis in scaled coordinates.
<br>
<tt><i>rrem</i></tt> specifies the radius (in &Aring;) of the cylinder.
<br>
If <tt><i>flag</i> = 0</tt>, then all atoms outside the cylinder is removed.
If <tt><i>flag</i> = 1</tt>, then all atoms outside the cylinder is fixed.
<br>
The operation is only applied for atoms whose <tt><i>zind</i></tt>-th scaled coordinate 
lies between <tt>szmin</tt> and <tt>szmax</tt>.

<p>
Example: <tt>MD++/scripts/CSD-book/chap3/sect3.1/ta-screw.script</tt>

<p>
For details, see <a href="http://micro.stanford.edu/~caiwei/CSD-Book/">Computer Simulations
of Dislocations</a>, Chapter 3, Section 3.1.
