INIT3D - Draws an isosurface when used in conjunction with entry DANDR.
CALL INIT3D (EYE,NU,NV,NW,ST1,LX,NY,IS2,IU,S)
This routine is part of the obsolete Isosrfhr utility. Isosrfhr has been replaced by the Isosurface utility. To see the Isosurface utility overview man page, type "man isosurface".
- (an input array of type REAL, dimensioned 3) is the position of the eye in the UVW coordinate system. Objects are considered to be in a box with 2 extreme corners at (1,1,1) and (NU,NV,NW). The eye position must have positive coordinates away from the coordinate planes U = 0, V = 0, and W = 0. While gaining experience with the package, use EYE(1) = 5.*NU, EYE(2) = 4.*NV, and EYE(3) = 3.*NW.
- (an input expression of type INTEGER) defines the range to be used for the U dimension of the 3-D UVW grid.
- (an input expression of type INTEGER) defines the range to be used for the V dimension of the 3-D UVW grid.
- (an input expression of type INTEGER) defines the range to be used for the W dimension of the 3-D UVW grid.
- (an input array of type REAL) used as scratch work space. It must be dimensioned at least of size NV*NW*2 words.
- (an input array of type INTEGER) which is the first dimension of array IS2. It must be of a length to hold NX bits. On a 64 bit machine it would be LX=(NX-1)/64+1.
- (an input expression of type INTEGER) defines the range of cells in the Y-direction of the XY pixel model of the image plane. NY is also the second dimension of array IS2.
- (an input array of type INTEGER) which holds the XY pixel model of the
image plane. It is dimensioned LX by NY. The pixel model contains NX by NY
elements. The sizes of NX and NY vary with the choice of options. For
example, if hidden lines are drawn, NX and NY (and LX) can be smaller than
when hidden lines are not drawn. If some visible lines are not drawn, NX
and NY (and LX) must be made larger. As a guide, some examples which have
been adequate are:
NU NV NW NX NY --------------------------------
100 100 60 200 200
60 60 60 110 110
40 40 40 75 75
- (an input expression of type INTEGER) defines the unit number of a scratch file for the package. ST1 will be written NU times on this file.
- (an input array of type REAL, dimensioned 4) contains the world
coordinates of the area where the picture is to be drawn. That is, all
plotting coordinates generated will be bounded as follows: X coordinates
will be between S(1) and S(2), Y coordinates will be between S(3) and
S(4). To prevent distortion, have S(2)-S(1) = S(4)-S(3).
All coordinates should be in a range between 10. and 1010.
Entries INIT3D and DANDR must be used in tandem to create a plot. First call entry INIT3D to create a set of 2-D perspective planes which are output to a scratch file. Then, in the calling program you must loop over the U dimension of UVW space starting at the last slab, NU, and moving forward to slab 1. Initially set all elements of IOBJS to zeros. Then, at all locations (J,K) where the surface of an object would appear, set IOBJS(J,K) = 1. When IOBJS(J,K) has been defined for that VW slab (U = NU - I + 1), call the entry DANDR, which draws any object contours based on the 2-D perspective defined for that slab.
An Isosrfhr example of the logical structure of a program follows:
C Entry INIT3D writes NU 2-D perspective planes, ST1(NV,NW,2),
C of the UVW grid, viewed from the position EYE, to a scratch file.
C NU is the extent of the U dimension.
DO 1 IBKWDS = 1,NU
I = NU+1-IBKWDS
C You must define the presence of the objects at each VW grid point for
C all slabs in the U direction. Ones are used to indicate the presence
C of an object. Zeros indicate no object.
C IF (OBJ.EQ.YES) IOBJS(J,K) = 1
C IF (OBJ.EQ.NO) IOBJS(J,K) = 0
C DANDR reads ST1(NV,NW,2) for slab I.
1 CALL DANDR(NV,NW,ST1,LX,NX,NY,IS2,IU,S,IOBJS,MV)
Use the ncargex command to see the following relevant examples: tisohr
To use INIT3D load the NCAR Graphics libraries ncarg, ncarg_gks, and ncarg_c, preferably in that order.
Online: isosrfhr, dandr, isosurface, isosurface_params, ezisos, isgeti, isgetr, isseti, issetr, pwrzi, ncarg_cbind
Copyright (C) 1987-2009
University Corporation for Atmospheric Research
The use of this Software is governed by a License Agreement.