# fracture.geomp21

Syntax

fp21 = fracture.geomp21(gset<,dfn>) (3D only)

Get the P21 of polygons. The P21 is the cumulative sum of intersection lengths with the geometry polygons per unit area of the geometry polygons. If dfn is not given the calculation is undertaken over all fractures of all DFNs.

Returns: fp21 - P21 value gset - geometry set pointer dfn - optional DFN pointer

Usage Example

The following example illustrates how fracture.geomp21 can be used to get the P21 on an outcrop for various DFNs, or for all fractures.

model new
model domain extent -5 5

;Create the outcrop geometry
geometry set 'outcrop'
geometry polygon create by-position (-5,-5,0) (-5,5,0) (5,5,0) (5,-5,0)

;Generate fractures using the default template
fracture generate dfn 'p1' fracture-count 100
fracture generate dfn 'p2' fracture-count 200
fracture generate dfn 'p3' fracture-count 300
fracture generate dfn 'p4' fracture-count 400

fish define access_p21
local geom = geom.set.find('outcrop')
local dfn1 = dfn.find('p1')
local dfn2 = dfn.find('p2')
local dfn3 = dfn.find('p3')
local dfn4 = dfn.find('p4')
io.out('P21 of dfn p1 = '+string(fracture.geomp21(geom,dfn1)))
io.out('P21 of dfn p2 = '+string(fracture.geomp21(geom,dfn2)))
io.out('P21 of dfn p3 = '+string(fracture.geomp21(geom,dfn3)))
io.out('P21 of dfn p4 = '+string(fracture.geomp21(geom,dfn4)))
io.out('P21 of all dfns = '+string(fracture.geomp21(geom)))
end

[access_p21]

program return