gdr = gd.getResultFile('FiberGeo.gdr')

# get result file and

# assign it to variable

# gdr

svp_reached = gdr['ResultMap']['SolidVolumePercentage']

# get reached solid volume

# percentage from result map and

# assign it to variable svp_reached

svp_target  = gdr['InputMap']['SolidVolumePercentage']

# get target solid volume

# percentage from input map and

# assign it to variable svp_target

svp_reached,svp_target = float(svp_reached),float(svp_target)

# transform string values to float

svp_error   = abs(svp_reached - svp_target)/svp_target

# compute relative error

# between reached and

# target svp

print(svp_error)

# print the computed error

VolumeInfo = gd.getVolumeFieldsInfo()

# get list of dictionaries of loaded

# Volume Fields and assign it to

# variable VolumeInfo

print(VolumeInfo)

# print all data contained in VolumeInfo

# to console / logfile

nx, ny, nz = gd.getVolDimensions()

# get the number of voxels in

# X-, Y- and Z-direction and

# assign the numbers to variables

Structure = gd.getStructure()

# assign 3D numpy array describing

# the loaded structure to the

# variable Structure

Name = VolumeInfo[0]['name']

# assign the name of the first

# volume field to

# the variable Name

VolumeField = gd.getVolumeField(Name)

# assign the numpy array describing

# the volume field

# to the variable VolumeField

Statistic = []

# Create empty list to store

# the statistical values

for k in range(nz):

# loop over all Z-layers

 value_sum = 0

# creating variable value_sum to

# sum up the result values

 value_count = 0

# creating variable value_count

# to count the summands

for j in range(ny):

# loop over all Y-layers

  for i in range(nx):

# loop over all X-layers

    if Structure[i][j][k] == 0:

# condition: if the kth Z-value in

# the jth Y-column of the ith X-layer

# is pore space (ID 0), the following

# indented section is executed

       value_sum = value_sum + VolumeField[i][j][k]

# add all pore space result values

# of the kth Z-layer to the sum value_sum

       value_count = value_count + 1

# count the summands of value_sum

 meanVal = value_sum / value_count

# compute mean value of all pore

# space result values in the kth

# Z-layer and assign it to the

# variable meanVal

 Statistic.append(meanVal)

# append the mean value of the

# Z-layer to the Statistic list

gDlg = gd.makeGraphDialog()

# create a graph dialog object

graph = gDlg.addGraph(Name, "Z layers", Name)

# add a graph object with the name

# of the volume field as title and

# Y-axis legend and Z-layers as

# X-axis legend  

Z_layers = list(range(1, nz + 1))

# writes the Z-layer numbers

# 1, 2, …, nz-1, nz into a

# list named Z-layers

graph.addData(Z_layers, Statistic, Name)

# add a single dataset with the

# Z-layers as X-values, the mean

# result values as Y-values and

# the name of the volume field

# as legend to this graph

gDlg.run()

# show graph dialog

Number = gd.getNumberOfTriangles()

# assigns the number of triangles

# to the variable Number

gd.msgBox ( f"The number of triangles is {Number}")

# show message dialog

Box = gd.getTriangulationBoundingBox()

# assigns the host name

# to the variable Host_name

X = Box[1][0]*10**6

# get the first entry of the

# second dictionary, i.e.

# the X-dimension in m,

# transform it to µm and assign

# it to the variable X

gd.msgBox ( f"The bounding box has {X} µm in X-direction.")

# show the result in

# a message dialog