SURFPTS


Generates points on boundary surfaces previously defined by the SURFACE or REGION command.

The points are generated by shooting rays through a user specified set of points from an origin point, line or plane and finding the surface intersection of each ray.

SYNTAX

surfpts/ itype/iname/ inside outside both/ifirst,ilast,istride/geom/ray_origin

itype is the name of the surface or region.

ifirst,ilast,istride or pset,get,psetname define the set of points to shoot rays through.

inside outside or both indicates the point location for a region relative to the surface.

Options for Ray Geometry

The parameters geom / ray_origin are the geometry and ray origin with parameters based on the geometry.

xyz/ x1,y1,z1/x2,y2,z2/x3,y3,z3 where points 1, 2, 3 define the plane to shoot rays from that are normal to the plane.

rtz/ x1,y1,z1/x2,y2,z2 where points 1, 2, define the line to shoot rays from that are perpendicular to the line.

rtp/ xcen,ycen,zcen where point is the center.

points / iffirst,iflast,ifstride define a set of points to shoot rays from.

EXAMPLES

* coordinates and elevations for watertable plane
define x1 499029.5961
define x2 500174.3616
define y1 538579.7712
define y2 539378.5612
define uleft 1779.9725
define uright 1778.7472
define lleft 1779.7481
define lright 1778.5228

* point distribution for ray shooting
define / XMIN / 498400. 
define / YMIN / 538900.  
define / XMAX / 500200. 
define / YMAX / 540400. 
define / NZ / 1 
define / NX / 61 
define / NY / 51 

cmo/create/cmos
surface/swtr/intrface/plane/x2 y1 lright/x2 y2 uright/x1 y2 uleft

createpts/xyz/NX NY NZ / XMIN YMIN ZMAX / XMAX YMAX ZMAX 
pset/rayend/seq/1 0 0

* SURFPTS/CTYPE/INAME/IREGPT/IPFIRST,IPLAST,ISTRIDE/GEOM/RAY ORIGIN
* ctype = surface or region
* iname = name of surface or region
* iregpt = inside, outside, or both (ignored for surface)
* geom = center pt, line, or plane
*  

surfpts/surface/swtr/inside/pset,get,rayend/ &
    xyz /0. 0. 0./1. 0. 0./1. 1. 1./

* remove ray points and duplicate points
rmpoint/pset get rayend
filter/1 0 0
rmpoint/compress

Generate a set of points on a plane