procedure sb2srf

This subroutine plots a 3-d surface given a 2-d unit-cell of density. All (x,y,z) values are taken to be given in world coordinates. The z-component of the eye-poisition should be positive and that of all the lattice-vertices should be negative; the viewing-screen is fixed at z=0.

Syntax
sb2srf EYE,LATICE,DENS,N1,N2,DLOW,DHIGH,DVERT,IC1,IC2,NCBAND,LIGHT,LSHINE

Arguments
ARGUMENT    TYPE I/O DIMENSION DESCRIPTION
EYE:   R*4 I 3 (x,y,z) coordinate of eye-position.
LATICE:   R*4 I 3 x 3 (x,y,z) coordinates of the origin
   and the a and b lattice-vertices.
DENS:   R*4 I (N1+1) The density at regular points within
   x (N2+1) the unit cell, wrapped around so
   that DENS(0,J)=DENS(N1,J) etc..
N1,N2:   I*4 I - The dimensions of the unit-cell grid.
DLOW:   R*4 I - Lowest density to be plotted.
DHIGH:   R*4 I - Highest density to be plotted.
DVERT:   R*4 I - "Vertical" world-coordinate length
   corresponding to density-range.
IC1,IC2:   I*4 I - Lowest and highest colour-index to
   be used for the rendering.
NCBAND:   I*4 I - Number of colour-bands for the
   height, so that the number of shades
   per band = (IC2-IC1+1)/NCBAND.
LIGHT:   R*4 I 3 (x,y,z) direction of flood-light.
LSHINE:   L*1 I - Shiny surface if TRUE, else diffuse.