CHARMM c28a3 graphx.doc

File: Graphx ]-[ Node: Top
Up: (commands.doc) -=- Next:


Graphics is a subparser of charmm, invoked by via the GRAPH command.
All of the miscellaneous commands (miscom.doc), coordinate commands
(corman.doc), and internal coordinate commands (intcor.doc) are
available from the GRAPHX> prompt.  Only the 1st three characters are
used for primary graphics commands, but many of the options require
the 1st four characters.

     The graphics facility has been extended to provide general X11
support, and the original Apollo GPR screen display has been dropped; 
a NODISPLAY version can also be built, which will generate all of the
derived files.  The other major enhancement is the production of
PostScript output files, in either color or grayscale; both X11 and
PostScript use the Apollo imaging model.  Additional information on X11
usage tips and compiling for X11 are given at the end of this document.
Finally, a recent addition is the production of input files for POV-Ray,
an excellent freeware ray tracing package for making high quality 
molecular images.  See

     Option keywords are indicated by the use of upper case; lower
case terms are variable values, generally real numbers, but decimal
points are not required.  Triplets ( x y z ) are position dependent;
omitted values are assumed to be zero.  Items enclosed in square 
brackets are [optional] but their absence often implies a default 
choice.  Default choices are indicated with an asterisk (*) in syntax 
listings where apropriate.

* Menu:

* Summary::       Syntax and Command Summary
* Description::   Detailed Command Description with Examples
* Output::        PostScript, FDAT, LIGHT, and POV-Ray file formats
* Addendum::      X11 Usage and Compiling Tips, Other Useful Programs

File: Graphx ]-[ Node: Summary
Up: Top -=- Previous: Top -=- Next: Description

                       GRAPHX COMMAND SUMMARY

	Initializing, exiting:

GRAPhx  [XXSZ iwidth] [XYSZ iheight]  [ NVERtex integer ]

OFF                 :: disable all graphics and exit (see END, below).
END                 :: exit from command parser only; preserve window
RESet               :: revert to default view, scale
DEFault-colors      :: revert to default colors

HELp                :: provides a command listing

	Available within the graphics subcommand parser:

miscellaneous-command-spec         ! see *note miscom:(miscom.doc).
IC        ic-subcommand-spec       ! see *note intcor:(intcor.doc).
COOR      coor-subcommand-spec     ! see *note corman:(corman.doc).

	These commands affect what is viewed:

DRAw [atom-selection]   :: change displayed atoms, force a redraw

DISplay [ON]* [MAIN]* [COMP] [BOND]* [VECT] [ATOM] [TEXT] -
                      [HBONds] [LABEls] [AXES]    :: change displayed objects

COLor color-name [brightfactor] [COMP] atom-selection
                 [ LEVel izcue-level ]

  color-name =   NONE   RED  BLACk   YELLow GREEn  WHITe
                 BLUE   CYAN MAGEnta GRAY   ORANge BROWn
                 PURPle TURQuoise    CHARtreuse  DKBLue
                    (or an integer from 0 to 15)

  izcue-level = an integer from 0 to 11; 0 == brightest, 11 == dimmest

HBStyle  [COLOR color-name]  [WIDTH iwidth]  [DASH idash]  :: HBOND style
LBL  label-type  label-atoms  SIZE label-size COLOr label-color
                 USER user-label
   label-atoms = FIRSt and/or atom-selection
   label-size  = VSMAll SMALl MEDIum LARGe                   default: SMALL
   label-color = color-name ( see COLOR command above ... )  default: YELLOW
   user-label  = up to 8 characters

LINe iwidth         :: bonds or vectors; pixels
RADii [DEFaults] [PARam] atom-scale [bond-scale] atom-sel
AXEs  [XLIM xmin xmax] [YLIM ymin ymax] [ZLIM zmin zmax]
      [DEFAult]  (all limits set to -25, +25 A)
STEreo [ON/OFF] [dist] [angle]

TEXt [text-body]    :: set title text, null entry clears
FONT [ VSMAll | SMALl | MEDIum | LARGe ]  :: text font; default MEDIUM

	These commands change the view only:

SCAle  factor [MOL/LAB] [REP int]
BOXsize size  [MOL/LAB]
CENter [atom-selection]
POInt  x y z
ROTate rx ry rz  [MOL/LAB] [REP int]
TRAnslate x y z  [MOL/LAB] [REP int]
ZCLip [low] high
ZCUe [ [low] high ]/[AUTO]

	These commands do not affect the display:

AUTo [ON*/OFF]      :: redraw after every command
ERAse [ON*/OFF]     :: screen clear prior to next drawing
EXEcute pathname    :: execute a program; no arguments are passed

	Output Files: the UNIT must be OPENed first...

PLUto UNIT n                                 :: PLUTO FDAT file
MAKE  UNIT n                                 :: LIGHT .atm file

PSC   UNIT n [COLOr] [BWREv] [PORTrait]      :: PostScript file

POV   UNIT n  [ INIT ]                       :: POV-Ray file
              [ TERM ]
              [ INCLude n ]  [ UOBJ n ]  
              [ OBJect  obj-spec ]

	(see the description for full syntax details)

File: Graphx ]-[ Node: Description
Up: Top -=- Previous: Summary -=- Next: Output



GRAPhx  [XXSZ iwidth]   [XYSZ iheight]     [NVERtex integer]

Invoked from the main CHARMM command parser; if already initialized
(i.e. GRAPHX ... END) the previous graphic states are retained.

The XXSZ and XYSZ options set the X11 window width and height for
the duration of the graphics session; a window resize is NOT passed
back to graphics code at this time.

The NVERtex option increases the allocated storage for the vertex array
used by the POV-Ray object code, especially for the MESH and RIBBON
objects; the default is set to NATOM. but may need to be increased for
scripts which create a lot of POV objects, especially RIBBON objects.

The NOWIndow option suppresses the X11 window for batch mode or remote
usage; ideal for making PostScript files, etc.



List the available commands and syntax.



Disable all graphics and exit the graphics subcommand parser; see
also the END command.



Exit from command parser only; allows other CHARMM commands to be
performed w/o losing the current display.  Re-invoking graphics
does not re-initialize the graphics settings.  Useful for trajectory
movies using MERGE DRAW, e.g.

open unit 12 read file name dyn.trj.10
merge draw firstu 12 nunit 1

Overlays can be created (on screen only) by preceding the END command 
in the above example with ERAse OFF.



Restores view settings to the program defaults ( scale, translation
and rotations).


DISplay [ON*]  [MAIN*] [COMP] [VECT] [ATOM] [BOND*] [TEXT] [HBONds] -
               [LABEls] [AXES]

Turns the display of various graphic features on or off; the default
is DISPLAY ON MAIN BOND TEXT which will show the connectivity of
the atoms in main coordinate set, using the default CHARMM title.
The options are:

ON*/OFF        enable (default) or disable the display of a feature
MAIN*          the main coordinate set; displayed by default
COMP           the comparison coordinate set; both may be displayed
BOND*          atom connectivity as atom-colored half-bonds; default
ATOM           filled circles using current radius (see RADII)
TEXT*          title display
HBONDS         current HBOND list, using double width lines
LABELS         residue names, atom types, user labels ...
AXES           lab frame axes; + end solid, labeled; - end dashed


display atom on          ! enable atom display; MAIN assumed
display text             ! toggles title display
display hbonds off       ! disable H-bond display


DRAw [atom-selection]

Forces a redraw when AUTO mode is off; also used to to change which
set of atoms is currently being displayed.  All display modes and
output files from PLOT, PSC, PLUTO, XMOLE, and MAKE commands use this
atom selection.  The initial selection is all atoms and may be restored 

draw sele all end


COLor color-name [brightfactor] [COMP] atom-selection
                 [LEVel i]

Sets the color of individual atoms according to the atom
selection, using one of the color names below:

      PURPle TURQuoise    CHARtreuse  DKBLue

The color applies to the main coordinate set unless COMP is
specified; brightfactor is a relative intensity, 0.0 -- 1.0
Alternatively, the zcue levels may be set directly to one of the
indices, an integer in the range 0:11 with 0 being fully bright.

Setting all atoms to WHITE, and using the brightfactor gives the
best control over how PostScript grayscale output will look on the
printed page; note that NONE is background, BLACK is the basic carbon
color is really gray, and WHITE will be black in PostScript grayscale
output or in color output with black/white reversal.


     All carbons for segment s are colored cyan:

color cyan sele type c* .and. segid s end

     Colored based on weighting array:

color green lev 11 sele prop wmain .gt. 1.0 end
color green lev 10 sele prop wmain .gt. 2.0 end



Restore the default color assignments, based on element type.




LINe iwidth  (bonds or vectors; pixels)

Set the line width for bonds & vectors, in pixels (integer).


line 2


RADii [DEFaults] [PARam] scale [bond] atom-sel

Sets the radius for displaying atoms, and for output files produced
by the PLOT, PSC, XMOLE, and MAKE commands.  The options are:

     scale          required if no other options are specified, and
                    assumed to be 1.0 if omitted; performs a relative
                    scale if used by itself, or scales the radii set
                    by the DEFAULT or PARAM options

     DEFAULTS       set radii to a convenient size for display, based
                    on atom type ( C, N, O, ... )

     PARAM          use VDW radii

     bond           value for bond radii for LIGHT program (see MAKE)

     atom-sel       atom selection to apply the radii command to


rad 0.8                       ! reduces radii to 80% of current value
radii param .5                ! set radii to 50% of VDW
radii param .5  .15           ! bond radii to 0.15 A
rad 1.5 sele type H* end      ! enlarge all H atoms by 50%


LBL  label-type  label-atoms  SIZE label-size  COLOR label-color
   label-atoms = FIRST* and/or atom-selection
   label-color = color-name ( see COLOR command above; default: YELLOW )
   user-label  = up to 8 characters

The LBL command identifies which atoms are to be labeled, what
atom attributes are to be included in the label, and the relative
size of the labels; the defaults are marked with an asterisk [*].

The INIT option clears all labels, and any other options are ignored.

One or more the following attributes may be included in the label,
by simply including the keyword(s) in the LBL statement:

      SEGId       segment name (from GENErate; A4)
      RESN        residue name (from the RTF; A4)
      RESId       residue ID, a numeral (A4)
      TYPE        atom type, e.g. N, CA, CB ...  (A4)
      CHEM        atom parameter type code (A4)
      CHARge      atomic charge (G12.4)
      WEIGht      value stored in the weight vector (G12.4)
      USER        arbitrary user-specified text (A8)

The label length (24 bytes) is such that all attributes may NOT be
displayed simulaneously; in particular, CHARge and WEIGht may not be
displayed at the same time for the same atom.

SIZE is specified by one the keywords VSMALL, SMALL, MEDIUM, LARGE,
with a default of SMALL.  The COLOR keyword allows setting the label
color, using the same color names as the COLor command; the default
label color is yellow.  Each use of the LBL command can create a group
of labels of a different size and color, for atoms which don't overlap
with any previous label atom selections.

The blank delimited word following the keyword USER, up to 8 characters, allows
the use of any text string as a label for the selected atoms.

The default is to label the first atom of each residue; an atom
selection overrides this, unless the FIRSt keyword is present;
in this case, the first atom of each selected residue is labeled.


! the first atom of each residue is labeled by name with normal text

! the first atom of each residue in the segment MAIN is labeled
! by name and number with very small text

! all oxygen atoms are labeled by charge with small text

! all alpha-carbons are labeled by the weight vector with medium text

! enter a null label; may be used to selectively "blank" labels
! in this case, all alpha-carbon labels are set to a string of blanks
! for display efficiency, LBL INIT is preferable

! show the location of formal charges on amino acid side chains


HBStyle  [COLOR color-name]  [WIDTH iwidth]  [DASH idash] 

Set the style for representing HBONDS; color-name is as for the
COLOR command, and iwidth and idash are integers in pixel units.
Specifying HBSTYLE alone resets to the default style, which is
equivalent to


If at least one option (COLOR, WIDTH, or DASH) is specified, the
remaining options are unchanged; thus HBSTYLE COLOR WHITE will
not reset the WIDTH to 4 pixels, but leave it as it was.


AXEs  [XLIM xmin xmax] [YLIM ymin ymax] [ZLIM zmin zmax]

Changes the length of the displayed axes; the default, which is
restored via the DEFAult keyword, is -25 to +25 A for each dimension.
Only the axes lines specified are changed, e.g. AXES ZLIM -5 5 only
changes the endpoints of the Z axes line.


STEreo [ON/OFF] [dist] [angle]

Invoke side-by-side stereo mode for screen display and for output
files produced by the PLOT and MAKE commands, when the ON keyword
is used, or when stereo is off.  The  dist  option controls the
separation between the two images; the  angle  option specifies
the parallax angle for left and right eye views.  The OFF keyword
is used to return to mono mode, and is assumed if the command is
used while in stereo mode.  The default angle is 7 degrees, but the
dist parameter defaults to an interval derived from the window width;
wider windows are best for stereo.


stereo on 16.0 7 


TEXT [ text-body | "text-body" ]

Supply the text for the title display; quotes override the conversion
to all upper case, and a null entry clears the current title.  The
quotes are removed, and are not displayed as part of the title; the
current CHARMM title is the default graphics title.


FONt [ VSMAll | SMALl | MEdium | LArge ]

Change the title font to one of four sizes; the initial setting
is MEDIUM, as is the default if no size is specified.


SCAle  factor [MOL/LAB] [REP int]

Change the Angstrom/pixel scale factor; initially, 1 A = 32 pixels.
The repeat factor can provide an ersatz zoom effect.  The default is
the LAB frame (i.e. scale is done relative to screen center rather
then the system center).


BOXsize size  [MOL/LAB]

Set the viewing 'box' to (size)x(size), in Angstroms.


CENter [atom-selection]

Move the selected atoms to the center of display space.



Scales the molecule to fit in the display window.


POInt  x y z

The point specified by  x y z  becomes the center of display space.


ROTate rx ry rz  [MOL/LAB] [REP int]

Apply a rotation to the viewing transform; does not affect the
coordinates.  The default is the LAB frame (i.e. a Z rotation
always spins the screen).  If the MOL (molecule) frame is used,
then the rotation will be about the origin of the molecular coordinate
system (i.e. does not depend on the current view matrix).
(NOTE1: the ROTate command uses a left handed system, multiply rotation
angles by -1 as necessary).
(NOTE2: if more than one rotation angle is given as nonzero, then the
rotations will occur sequentially, first x, then y, and then z).


rot 0 90       ! rotate by 90 deg around the y axis
rotate 180     ! rotate by 180 deg around the x axis
rot 90 0 90    ! rotate by 90 deg around x, then 90 deg around z


TRAnslate x y z  [MOL/LAB] [REP int]

Apply a translation to the viewing transform; does not affect the


tran 2 9.5     ! translate +2 A along x axis, +9.5 along y axis
tra 0 0 4      ! translate +4 axis along the z axis


ZCLip [low] high

Set hither and yon clip limits; atoms outside the limits are not
displayed, whether selected or not.


zclip 10       ! atoms outside z = ( -10 .. +10 ) are not displayed
zclip -5 10


ZCUe [[low] high ]/[AUTO]

Controls the z coordinate range over which depth cueing will be
applied.  AUTO use the displayed atom coordinates to set the limits.


zcue 10        ! zcue from -10 to +10
zcue -4 8
zcue auto



Enables or disables automatic redraw after every command; the
initial setting is ON, and the command functions as a toggle.
AUTO OFF is useful for making multiple changes without the
time required for a redraw.



Enables or disables a screen clear prior to the next drawing; the
initial setting is ON, and the command functions as a toggle.
ERASE OFF is useful for overlaying trajectory frames or other
related collections of structures.  Possibly best used with AUTO OFF,
and using DRAW for each structure.  Also applies to POV and PSC
commands; after using the INIT keyword, the file is not closed until 
one of the following:

  (1) PSC/POV command with TERMinate keyword; only correct method
  (2) the UNIT is CLOSEd; no title or page eject, etc.
  (3) the program terminates; no title or page eject, etc.

Creates an overlay of 50 consecutive trajectory frames (coord traj 
open on unit 20; .ps file open on unit 12):

erase off
psc unit 12 init
set n 1
label frmlp
coor read file unit 20 ifile -1
psc unit 12
let n += 1
if n .lt. 50 goto frmlp
psc unit 12 term


EXEcute pathname

Execute a program w/o arguments.


exe /bin/ls         ! list the current directory  (unix)

Equivalent to the SYSTem command (miscom.doc)

File: Graphx ]-[ Node: Output
Up: Top -=- Previous: Description -=- Next: Addendum

                        Output Facility

PSC UNIT n [COLOr] [PORTrait] [BWREverse] [INIT | TERM]

Writes out a PostScript display file using the current viewing transform
and selected atoms; the default is grayscale in landscape mode, and may
be changed using the COLOr and/or PORTrait keywords.  Colors are direct
translations of the RGB color map used for Apollo and X11 displays, or
arbitrary conversions to grayscale.  Direct control over grayscale may 
be achieved by setting all atoms to COLOR WHITE, and using the 
'brightfactor' option of the COLOR command to set the various gray 
levels; note that full WHITE is black, since the background (paper) is 
white in grayscale mode.  The default background is black in COLOR mode, but may changed to white
with BWREverse keyword; atoms or text colored white on the screen will 
be printed in black.  As always, the UNIT must be opened first 
(recommended file extension .ps).

The INIT keyword forces device initialization in ERASE OFF mode, and 
should only be used in that context; it is not normally required. 
Likewise, the TERM keyword writes the final part of the PostScript file in ERASE OFF mode,
and does not draw the structure; only the graphics title and the final 
few PostScript commands are written to the file.  See the ERAse command 
for additional details and a usage example.



Writes out atom coordinates and connectivity based on the current atom
selection and view transform, in CSD FDAT format.  
WARNING: 999 atom limit!

Stereo mode settings are ignored, as are radii and color.  Also, atoms
should be renamed to their element types to get proper radii, etc within
pluto (e.g. CA is not carbon alpha, it's calcium).  As with PLOT, the
UNIT must be opened first ( recommended file extension  .fdat ).


rename atom C1 sele type C  end
rename atom C2 sele type CA end
rename atom N1 sele type N  end
rename atom O1 sele type O  end
rename atom H1 sele type H1 end
   ! etc., etc., etc.
scale .5
rot 0 90
open unit 50 write card name molecule.fdat
pluto unit 50



Writes atom coordinates, radii, and color to a file formatted for
the LIGHT program (available from the N.I.H.) which produces nice
ray-traced images using current stereo settings and view transform.

A 1280x1024 window is the optimum size for best scaling and centering
of the image produced by LIGHT.

The UNIT must be opened first; the LIGHT program requires that an
extension of .atm be use for the file name.


POV  UNIT int  [ INIT ]  [ INCLude int ]  [ UOBJ int ]
               [ TERM ]

        Write atom-based objects; bonds, atoms, H-bonds, labels

  UNIT int    :: final POV file; user must open; also open int+1 for stereo
  INCLude int :: read file w. user override of default camera, lighting, ...
  UOBJ int    :: read file w. OBJect definitions, i.e. CYLI, SLAB, ...

  INIT        :: intialize POV file in ERASE OFF mode (overlays)
  TERM        :: close POV file in ERASE OFF mode; other options ignored

All units must be opened; in STEreo mode, UNIT numbers are used in pairs
as N, N+1 for the left- and right-eye views.  The object file must be 
created first with POV OBJ commands (described below).  The atom, bonds,
H-bonds, and labels currently displayed are written to the file; just 
as for the X11 display and PostScript output, this is determined by the 
last atom selection used with a DRAw command, and graphics elements 
enabled via the DISplay command.  The current viewpoint is used for 
transforming coordinates for the POV output file, and the graphics 
SCAle determines the camera distance.  The atom radii are determined 
from the current graphics RADii, while the bond cylinders are scaled 
based on the graphics LINe width.

The POV file produced (recommended extension .pov) is a flat text file, 
and may be manually edited to further customize the image that will be 
ultimately created by processing the POV file.  Note that CHARMM only 
creates an input file for POV-Ray; it is assumed that the user has 
access to version 3.x of the POV-Ray program.  When using the povray3 
program on Unix systems, the PPM bitmap output format is recommended.  
For details on obtaining and using the POV-Ray program, see the web 
site at

The INCLude file provides an easy mechanism to add user customization, 
e.g. changing the light sources or camera position, or adding other 
POV-Ray objects to be displayed.  Comments in the POV input file 
created indicate the section which is completely replaced by the 
userfile; using the CHARMM supplied default setup as a template for 
INCLude files is strongly encouraged.

The INIT keyword forces device initialization in ERASE OFF mode, and 
should only be used in that context; it is not normally required.  
Likewise, the TERM keyword writes the final part of the POV file in 
ERASE OFF mode, and does not output a structure; only the graphics 
title is written to the file.  See the ERAse command for additional 
details and a PSC (PostScript) usage example.

	Example: write a single POV file with some added objects

open unit 11 read card name
open unit 21 write card name protein.pov
pov unit 21 uobj 11

	Example: write two POV files for a stereo image pair

open unit 21 write card name protein-l.pov
open unit 22 write card name protein-r.pov
pov unit 21 


        Object definition command; add cylinders, slabs, spheres, etc.

POV  UNIT int  OBJect  obj-spec

This command requires the user to both OPEN and CLOSe the file; this 
allows multiple usages of the POV OBJ command, using the open UNIT 

  UNIT int :: object include file; user must open AND close; see UOBJ above

  obj-spec :: CYLInder    cyl-spec
              SLAB        slab-spec
              SPHEre      sphere-spec
              MESH        mesh-spec
              RIBBon      ribbon-spec
              BACKground  bkgrnd-spec
              RGB         rgb-spec

  cyl-spec :: PTS  x1 y1 z1  x2 y2 z2  [ RAD real ]  [ LEN real ]  [ obj-mat ]
              AXIS                          "             "             "
              HELIx     atom-sele           "             "             "
              PAX       atom-sele           "             "             "

  slab-spec :: PTS   x1 y1 z1  x2 y2 z2               [ obj-mat ]
               PAX   atom-sele  [ THICkness real ]         "

  sphere-spec :: x y z RAD real   [ SCALe real ]  [ obj-mat ]

  mesh-spec   :: ATOM atom-sele [ obj-mat ]
                 PTS  num-tri  [ PUNIt integer ]   [ obj-mat ]  
		                xyz1    xyz2    xyz3
				 :       :       :
				xyzN-2  xyzN-1  xyzN

  num-tri     :: number of triangles, i.e. data lines
  xyz1, xyzN  :: an ordered triple of real numbers; N = 3 * num-tri

  ribbon-spec :: ATOM atom-sele [ obj-mat ]

  bkgrnd-spec :: SIMPle   clr-spec
                 PLANe    [ ZPOS real ] [ SCALe real ] [ obj-mat ]

  rgb-spec :: r g b NAME "pov-color"      (see clr-spec below)

  atom-sele :: standard atom selection

  obj-mat  :: GLASs     [ clr-spec ]
              TRANsp    [ clr-spec ] 
              GRADient  [ clr-spec clr-spec ]  [ finish-spec ]
              STRIpe    [ clr-spec clr-spec ]  [ finish-spec ]
              CHECker   [ clr-spec clr-spec ]  [ finish-spec ]
              HEXAgon   [ clr-spec clr-spec clr-spec ]  [ finish-spec ]
              MARBle    [ clr-spec clr-spec clr-spec ]  [ finish-spec ]
              GRANite   [ clr-spec clr-spec clr-spec ]  [ finish-spec ]
              WOOD      [ wood-spec ]  [ finish-spec ]
              AGATe     [ finish-spec ]
              JADE      [ finish-spec ]

  finish-spec  :: SHINy | METAl | DULL | GLOSsy | MIRRor | LUMInous

  wood-spec :: CHERry | PINE | SPRUce | OAK | ROSEwood | MAPLe

  clr-spec     :: COLOr charmm-color | "pov-color"
  charmm-color :: standard CHARMM graphics color usage
  pov-color    :: POV color names; no err check, case sensitive (need ")

The syntax and options for the POV OBJect commands are a bit complex, 
in order to provide flexibility and a reasonable sampling of POV 
objects, textures, finishes, and colors.  Currently, one may define 
spheres, cylinders, meshes, ribbons, slabs, custom backgrounds, and 
custom colors.  This is an initial implementation, and additional 
object types are planned, such as toruses and other shapes, and 
possibly some surfaces.  A variety of POV textures are supported as 
the object material, with the standard POV finishes.  

Colors used may be standard CHARMM colors, standard POV colors, or 
user defined POV colors (via POV OBJ RGB command).  Since POV color 
names are case sensitive, they must always be enclosed within double 
quotes (") in a CHARMM command.  Note that the textures have default 
colors, which may be changed for some textures; in many cases two or 
three colors are needed to define the pattern, and in those cases all 
(not 1 or 2) colors must be specified.  For example, one must always 
either use the default colors or supply 3 color names for the HEXAgon 

The CYLInder object may be positioned in four different ways: (a) via 
the endpoints of its principal axis [PTS], (b) via the vector defined 
by the last COOR AXIS command [AXIS], (c) via the vector defined by 
the last COOR HELIX command [HELIX], or (d) via the principal axis of 
SELEcted atoms [PAX].  Each method has its own means of determining a 
default length and radius, which may be changed by using the LEN and 
RAD options.  Examples:

coor axis sele type CA .and. resid 23:48 end
pov unit 11 obj cyli axis rad 2.0 marble shiny

define back1 sele resid 23:48 .and. ( type N .or. type CA .or. type C ) end
pov unit 11 obj cyli pax sele back1 end glass "OldGold"

The SLAB object may be positioned either via the (left,bottom,front) 
and (right,top,back) corners [PTS], or via the principal axes of 
SELEcted atoms [PAX].  For PAX, the coordinate distribution determines 
the shape, but the smallest dimension may be changed via the THICkness 
option.  Example:

define beta1 sele type CA .and. ( resid 62:75 or resid 80:93 ) end
pov unit 11 obj slab pax sele beta1 end thick 0.25 wood oak glossy

The SPHEre object is placed by its center, and a RADius is required;
the SCALe option applies to the pattern, i.e. coarser or finer detail
in the texture (marble veins, hexagons, etc.).  Example:

coor stat sele atom A 1 N end
pov unit 11 obj sphere ?XAVE ?YAVE ?ZAVE rad 5.0 transp "NeonBlue"

The MESH and RIBBON objects define a set of triangles based on an 
atom selection or list of points to create a surface; the RIBBON 
uses "smoothed" triangles, and the triangle edges are not seen in 
the resulting displayed surface, as they are for the MESH.  The 
triangles are defined based on the atom coordinates in the sequence 
they are defined in the PSF, which in turn depends on the order of 
declaration of atoms in the RESIdue definitions in the original RTF 
used.  Protein backbone chrome ribbon usage example:

pov unit 11 obj ribbon chrome  sele type N .or. type CA .or. type C end

For the PTS option (MESH only), 3 xyz ordered triples are required for 
each triangle, meaning each data line has 9 real numbers.  For example, 
a parallelpiped could be defined by using two triangles for each face, 
or twelve data lines: 

pov unit 11 obj mesh glass "OldGold" pts 12
  x1  y1  z1    x2  y2  z2    x3  y3  z3
  :   :   :     :   :   :     :   :   :
  x34 y34 z34   x35 y35 z35   x36 y36 z36  

Note that 108 (12*9) real numbers are expected here, with the 3 xyz
triples that define a triangle on each subsequent line.  Alternatively,
the PUNIT option reads the same data (CARD data, 9 reals per line)
from the previously OPENed unit number.

The BACKground may be set to any single color [SIMPle], or a plane 
with any of the object material textures mapped onto it [PLANe].  The 
default PLANe position is Z=100 (hopefully behind the POV objects and
molecules), but may be changed via the ZPOS option.  As for the SPHEre,
the SCALe option applies to the texture; the pattern (not the object) 
is enlarged or shrunk by changing the SCALe value.  Examples:

pov unit 11 obj back simple "Gray60"

pov unit 11 obj back plane zpos 50. gradient "MidnightBlue" "SummerSky"

The RGB option may be used to add custom POV colors for objects defined 
later in the input stream; a color must be defined before it can be 
used.  Example:

pov unit 11 obj rgb 0.7 0.7 0.2 name "RvAmber"
pov unit 11 obj sphere  0.0  0.0  0.0 rad 20. transp "RvAmber"


	A final complete POV example, illustrating the basic points:

open unit 11 write card name
pov unit 11 obj back plane gradient "DarkPurple" "LightPurple"
pov unit 11 obj cyli pax sele resid 25:52 .and. type CA end aluminum
define beta1 sele type CA .and. ( resid 62:75 or resid 80:93 ) end
pov unit 11 obj slab pax sele beta1 end agate shiny
close unit 11

open unit 11 read card name
open unit 21 write card name protein.pov
pov unit 21 uobj 11

File: Graphx ]-[ Node: Addendum
Up: Top -=- Previous: Output -=- Next: Top

                          X11 Usage Tips

On some workstations, the graphics display colors will only be
correct when that window has "focus".

The recommended "focus" policy is pointer-has-focus, i.e. the active
window is whichever one the mouse pointer is in; this permits typing
graphics commands in a mostly obscured terminal window (xterm, winterm,
hpterm, aixterm, etc.).  The alternative is to "lower" the text command
window via the window manager.  Automatic "topping" of the window with
focus is not recommended.

Although the graphics display window can be resized, the size change is
not currently passed back to the drawing routines.

                       X11 Compiling Problems

     The default assumptions are:

(1) required include files (Xlib.h, Xutil.h) reside in /usr/include/X11

(2) the library libX11.a resides in /usr/lib/X11

These files *must* exist for the X11 graphics to be compiled, and of
course the keyword XDISPLAY must be in the pref.dat file.  If the files
are NOT in the standard places (such as under HP-UX), the compile option
-I may be needed to point to the directory with the include files, and
the corresponding -L option may be needed to point to the directory 
which contains libX11.a for the final link step.

                   Other Useful Programs (all "freeware")

(1) ghostscript -- besides allowing on screen preview of PostScript files
prior to printing, this useful utility can convert the .ps files output by
the PSC command to other formats (e.g. GIF, PBM) and can also be used to make
an "encapsulated" PostScript (EPS) file.  Available via FTP from almost any
software archive that distributes GNU software (Free Software Foundation).

(2) xfig -- this X11-based technical drawing package can read EPS files, and
provides an ideal way to do "pasteup" of PostScript files from diverse sources
into a composite figure.  Available via FTP from sites carrying X11 software;
hardware specific sites may have versions pre-configured for SGI, HP, etc.

(3) pbmplus, netpbm -- this suite of bitmap manipulation programs offers the
best looking conversion to GIF files, using a PBM file from ghostscript as an
intermediate file type.  The names pbmplus and netpbm indicate minor variants
of the same suite, with netpbm being more recent, according to UseNet lore.
Available at the MIT X11 site, and similar FTP sites.

(4) light -- a ray-tracing program written by BR Brooks at NIH, which uses
the file format ouput by the graphics MAK command; used to produce movies
and slick-looking single images for presentations.  Contact the author.

(5) gnuplot -- 2D and 3D graphics, with X11 and PostScript support (and about
50 other output devices as well); useful for time series plots, etc., which
can be combined with molecular graphics via xfig (see above).  Available at
most GNU sites and from the "home" site at Dartmouth.

(6) ImageMagick -- a powerful suite of image manipulation programs, including
animation capabilities; available at

(7) MPEG -- mpegencode and mpegplay produce and replay animations; multiple
images in PBM format can be combined into a movie

CHARMM .doc Homepage

Information and HTML Formatting Courtesy of:

NIH/DCRT/Laboratory for Structural Biology
FDA/CBER/OVRR Biophysics Laboratory
Modified, updated and generalized by C.L. Brooks, III
The Scripps Research Institute