CRootBox/python/example2b.py at master · Plant-Root-Soil-Interactions-Modelling/CRootBox · GitHub

33 lines (27 loc) · 920 Bytes

"""multiple root systems"""

import py_rootbox as rb

name = "Zea_mays_4_Leitner_2014"

simtime = 120

N = 3 # number of columns and rows

dist = 40 # distance between the root systems [cm]

# Initializes N*N root systems

allRS = []

for i in range(0, N):

for j in range(0, N):

rs = rb.RootSystem()

rs.readParameters("modelparameter/" + name + ".xml")

rs.getRootSystemParameter().seedPos = rb.Vector3d(dist * i, dist * j, -3.) # cm

rs.initialize()

allRS.append(rs)

# Simulate

for rs in allRS:

rs.simulate(simtime, True)

# Export results as single vtp files (as polylines)

ana = rb.SegmentAnalyser() # see example 3b

for i, rs in enumerate(allRS):

vtpname = "../results/example_2b_" + str(i) + ".vtp"

rs.write(vtpname)

ana.addSegments(rs) # collect all

# Write all into single file (segments)

ana.write("../results/example_2b_all.vtp")

print("done.")

Expand file tree