HydroShoot Tutorial - Grapevine#
Authors : R.Albasha, C.Fournier, C.Pradal, M.Chelle, J.A Prieto, G.Louarn, T.Simonneau, E.Lebon#
Publication : https://doi.org/10.1093/insilicoplants/diz007#
HydroShoot is a functional-structural plant model that simulates the interactions between shoot’s hydraulic structure, gas-exchange and energy-budget, at the organ level.
In this tutorial, we will :
reconstruct in 3D the shoot architecture of a vine from plant digitalization data
simulate hydraulic structure, gas and energy exchange dynamics
3D reconstruction of vine architecture#
Let’s import all the packages we need for this tutorial.
[2]:
from os import getcwd
from pathlib import Path
from openalea.widgets.plantgl import *
from openalea.widgets.mtg import *
from openalea.mtg import traversal
from openalea.plantgl.all import Scene
from openalea.hydroshoot import architecture, display, model
We load the shoot architecture of a vine from plant digitalization data with HydroShoot which will return a Multi-scale Tree Graph. Then, we display it.
[3]:
# Path for plant digitalization data.
g = architecture.vine_mtg('grapevine_pot.csv')
plot(g)
[3]:
[4]:
# Local Coordinates Correction
for v in traversal.iter_mtg2(g, g.root):
n = g.node(g.Trunk(v, Scale=1)[0])
theta = 180 if int(n.index()) < 200 else -90 if int(n.index()) < 300 else 0
architecture.vine_orientation(g, v, theta, local_rotation=True)
# Scene rotation
for v in traversal.iter_mtg2(g, g.root):
architecture.vine_orientation(g, v, 90., local_rotation=False)
for v in traversal.iter_mtg2(g, g.root):
architecture.vine_phyto_modular(g, v)
architecture.vine_mtg_properties(g, v)
architecture.vine_mtg_geometry(g, v)
architecture.vine_transform(g, v)
Now, we create a scene from the MTG and we display it with the PantGL() function.
[5]:
scene = display.visu(g, def_elmnt_color_dict=True, scene=Scene(), view_result=False)
PlantGL(scene)
[5]:
With the MTG() function, we can visualize a particuliar property of the MTG through a different colormap. For instance, we display the vine considering the length of every part of the shoot.
[6]:
MTG(g, 'Length')
[6]:
Simulation of hydraulic structure, gas and energy exchange dynamics#
Once we tested some basic features, we run the model and simulate hydraulic structure, gas and energy exchange dynamics of the shoot.
[16]:
model.run(g, Path(getcwd()), Path(str(getcwd()) + '/meteo.input'), Scene = scene, psi_soil=-0.2, gdd_since_budbreak=100.)
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+ Project: /home/bauget/Documents/Dev/oawidgets/doc/examples/hydroshoot_tutorial
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
------------------------------------------------------------------------
GDD since budbreak: 100.0 °Cd
Soil class: Loamy_Sand
Site location: (43.61 N, 3.87 E, 44.00 m amsl)
Calculate Energy budget: True
Ignore shoot resistance: False
Calculate Hydraulic structure: True
Add rhyzoshpere resistance: False
Use user form factors: False
Use user irradiance data: False
Use user nitrogen data: False
------------------------------------------------------------------------
========================================================================
Date: 2012-08-01 00:00:00
psi_error = 0.015 :: Nb_iter = 6 ipsi_step = 0.500000
t_error = 6.029 counter = 1 t_iter = 9 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 3.117 counter = 2 t_iter = 8 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 1.732 counter = 3 t_iter = 7 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 1.064 counter = 4 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.653 counter = 5 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.401 counter = 6 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.246 counter = 7 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.151 counter = 8 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.093 counter = 9 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.057 counter = 10 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.035 counter = 11 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.022 counter = 12 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.017 counter = 13 t_iter = 0 it_step = 0.25
---------------------------
psi_soil -0.2000
psi_collar -0.2008
psi_leaf -0.2071
gs: 0.0200
flux H2O 2.6863
flux C2O -0.3422689259835871
Tleaf 15.58 Tair 20.74
========================================================================
========================================================================
Date: 2012-08-01 01:00:00
psi_error = 0.015 :: Nb_iter = 6 ipsi_step = 0.500000
t_error = 4.479 counter = 1 t_iter = 8 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 2.286 counter = 2 t_iter = 7 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 1.194 counter = 3 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.698 counter = 4 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.408 counter = 5 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.239 counter = 6 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.14 counter = 7 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.082 counter = 8 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.048 counter = 9 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.028 counter = 10 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.022 counter = 11 t_iter = 1 it_step = 0.25
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.02 counter = 12 t_iter = 0 it_step = 0.125
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.019 counter = 13 t_iter = 0 it_step = 0.0625
---------------------------
psi_soil -0.2000
psi_collar -0.2009
psi_leaf -0.2072
gs: 0.0200
flux H2O 3.7178
flux C2O -0.3718426133976098
Tleaf 16.79 Tair 20.47
========================================================================
========================================================================
Date: 2012-08-01 02:00:00
psi_error = 0.015 :: Nb_iter = 6 ipsi_step = 0.500000
t_error = 2.968 counter = 1 t_iter = 8 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 1.502 counter = 2 t_iter = 7 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.761 counter = 3 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.41 counter = 4 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.228 counter = 5 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.127 counter = 6 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.071 counter = 7 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.039 counter = 8 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.022 counter = 9 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.017 counter = 10 t_iter = 0 it_step = 0.25
---------------------------
psi_soil -0.2000
psi_collar -0.2010
psi_leaf -0.2073
gs: 0.0200
flux H2O 4.6508
flux C2O -0.39588780934058665
Tleaf 17.71 Tair 20.07
========================================================================
========================================================================
Date: 2012-08-01 03:00:00
psi_error = 0.015 :: Nb_iter = 6 ipsi_step = 0.500000
t_error = 2.752 counter = 1 t_iter = 8 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 1.392 counter = 2 t_iter = 7 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.705 counter = 3 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.376 counter = 4 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.209 counter = 5 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.116 counter = 6 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.064 counter = 7 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.036 counter = 8 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.02 counter = 9 t_iter = 0 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.015 counter = 10 t_iter = 0 it_step = 0.25
---------------------------
psi_soil -0.2000
psi_collar -0.2010
psi_leaf -0.2073
gs: 0.0200
flux H2O 4.1292
flux C2O -0.38843073450752547
Tleaf 17.42 Tair 19.60
========================================================================
========================================================================
Date: 2012-08-01 04:00:00
psi_error = 0.014 :: Nb_iter = 6 ipsi_step = 0.500000
t_error = 3.268 counter = 1 t_iter = 8 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 1.662 counter = 2 t_iter = 7 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.847 counter = 3 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.477 counter = 4 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.271 counter = 5 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.154 counter = 6 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.088 counter = 7 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.05 counter = 8 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.028 counter = 9 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.016 counter = 10 t_iter = 0 it_step = 0.5
---------------------------
psi_soil -0.2000
psi_collar -0.2008
psi_leaf -0.2071
gs: 0.0200
flux H2O 2.7616
flux C2O -0.36040269169941547
Tleaf 16.30 Tair 18.93
========================================================================
========================================================================
Date: 2012-08-01 05:00:00
psi_error = 0.015 :: Nb_iter = 6 ipsi_step = 0.500000
t_error = 3.055 counter = 1 t_iter = 8 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 1.541 counter = 2 t_iter = 7 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.779 counter = 3 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.429 counter = 4 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.241 counter = 5 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.135 counter = 6 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.076 counter = 7 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.042 counter = 8 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.024 counter = 9 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.019 counter = 10 t_iter = 0 it_step = 0.25
---------------------------
psi_soil -0.2000
psi_collar -0.2009
psi_leaf -0.2072
gs: 0.0245
flux H2O 3.6076
flux C2O -0.22411278542749147
Tleaf 16.23 Tair 18.66
========================================================================
========================================================================
Date: 2012-08-01 06:00:00
psi_error = 0.019 :: Nb_iter = 7 ipsi_step = 0.500000
t_error = 2.675 counter = 1 t_iter = 8 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 2 ipsi_step = 0.500000
t_error = 1.202 counter = 2 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 2 ipsi_step = 0.500000
t_error = 0.552 counter = 3 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.257 counter = 4 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.125 counter = 5 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.062 counter = 6 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.031 counter = 7 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.015 counter = 8 t_iter = 0 it_step = 0.5
---------------------------
psi_soil -0.2000
psi_collar -0.2034
psi_leaf -0.2100
gs: 0.1161
flux H2O 22.9546
flux C2O 1.993459944737447
Tleaf 18.46 Tair 20.42
========================================================================
========================================================================
Date: 2012-08-01 07:00:00
psi_error = 0.029 :: Nb_iter = 7 ipsi_step = 0.500000
t_error = 3.97 counter = 1 t_iter = 8 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 1.352 counter = 2 t_iter = 7 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.455 counter = 3 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.195 counter = 4 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.084 counter = 5 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.036 counter = 6 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.015 counter = 7 t_iter = 0 it_step = 0.5
---------------------------
psi_soil -0.2000
psi_collar -0.2133
psi_leaf -0.2210
gs: 0.3865
flux H2O 94.7064
flux C2O 7.4354996417591845
Tleaf 22.66 Tair 22.91
========================================================================
========================================================================
Date: 2012-08-01 08:00:00
psi_error = 0.029 :: Nb_iter = 7 ipsi_step = 0.500000
t_error = 2.299 counter = 1 t_iter = 7 it_step = 0.5
psi_error = 0.001 :: Nb_iter = 4 ipsi_step = 0.500000
t_error = 0.966 counter = 2 t_iter = 6 it_step = 0.5
psi_error = 0.001 :: Nb_iter = 3 ipsi_step = 0.500000
t_error = 0.411 counter = 3 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.176 counter = 4 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.075 counter = 5 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.032 counter = 6 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.014 counter = 7 t_iter = 0 it_step = 0.5
---------------------------
psi_soil -0.2000
psi_collar -0.2110
psi_leaf -0.2185
gs: 0.3471
flux H2O 78.4185
flux C2O 6.67335026955239
Tleaf 23.25 Tair 24.51
========================================================================
========================================================================
Date: 2012-08-01 09:00:00
psi_error = 0.036 :: Nb_iter = 8 ipsi_step = 0.500000
t_error = 2.7 counter = 1 t_iter = 8 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 2 ipsi_step = 0.500000
t_error = 0.982 counter = 2 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.347 counter = 3 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.121 counter = 4 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.043 counter = 5 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.017 counter = 6 t_iter = 0 it_step = 0.5
---------------------------
psi_soil -0.2000
psi_collar -0.2193
psi_leaf -0.2277
gs: 0.4262
flux H2O 134.7067
flux C2O 8.150045008913985
Tleaf 25.44 Tair 25.59
========================================================================
========================================================================
Date: 2012-08-01 10:00:00
psi_error = 0.045 :: Nb_iter = 8 ipsi_step = 0.500000
t_error = 2.721 counter = 1 t_iter = 8 it_step = 0.5
psi_error = 0.001 :: Nb_iter = 3 ipsi_step = 0.500000
t_error = 0.959 counter = 2 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 2 ipsi_step = 0.500000
t_error = 0.344 counter = 3 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.14 counter = 4 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.057 counter = 5 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.023 counter = 6 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.009 counter = 7 t_iter = 0 it_step = 0.5
---------------------------
psi_soil -0.2000
psi_collar -0.2262
psi_leaf -0.2354
gs: 0.4510
flux H2O 178.4525
flux C2O 8.453443081095463
Tleaf 27.14 Tair 27.53
========================================================================
========================================================================
Date: 2012-08-01 11:00:00
psi_error = 0.047 :: Nb_iter = 8 ipsi_step = 0.500000
t_error = 2.555 counter = 1 t_iter = 7 it_step = 0.5
psi_error = 0.001 :: Nb_iter = 3 ipsi_step = 0.500000
t_error = 0.867 counter = 2 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 2 ipsi_step = 0.500000
t_error = 0.337 counter = 3 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.136 counter = 4 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.055 counter = 5 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.022 counter = 6 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.009 counter = 7 t_iter = 0 it_step = 0.5
---------------------------
psi_soil -0.2000
psi_collar -0.2278
psi_leaf -0.2371
gs: 0.4394
flux H2O 187.7694
flux C2O 8.327728105660174
Tleaf 26.34 Tair 26.84
========================================================================
========================================================================
Date: 2012-08-01 12:00:00
psi_error = 0.041 :: Nb_iter = 8 ipsi_step = 0.500000
t_error = 2.506 counter = 1 t_iter = 7 it_step = 0.5
psi_error = 0.001 :: Nb_iter = 3 ipsi_step = 0.500000
t_error = 0.867 counter = 2 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 2 ipsi_step = 0.500000
t_error = 0.308 counter = 3 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.125 counter = 4 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.052 counter = 5 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.021 counter = 6 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.009 counter = 7 t_iter = 0 it_step = 0.5
---------------------------
psi_soil -0.2000
psi_collar -0.2231
psi_leaf -0.2319
gs: 0.4368
flux H2O 159.0264
flux C2O 8.148085201258707
Tleaf 26.07 Tair 26.42
========================================================================
========================================================================
Date: 2012-08-01 13:00:00
psi_error = 0.049 :: Nb_iter = 8 ipsi_step = 0.500000
t_error = 2.436 counter = 1 t_iter = 7 it_step = 0.5
psi_error = 0.002 :: Nb_iter = 4 ipsi_step = 0.500000
t_error = 0.888 counter = 2 t_iter = 6 it_step = 0.5
psi_error = 0.001 :: Nb_iter = 3 ipsi_step = 0.500000
t_error = 0.358 counter = 3 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 2 ipsi_step = 0.500000
t_error = 0.145 counter = 4 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.059 counter = 5 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.024 counter = 6 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.01 counter = 7 t_iter = 0 it_step = 0.5
---------------------------
psi_soil -0.2000
psi_collar -0.2278
psi_leaf -0.2371
gs: 0.4346
flux H2O 187.6562
flux C2O 8.111484153914883
Tleaf 26.17 Tair 27.15
========================================================================
========================================================================
Date: 2012-08-01 14:00:00
psi_error = 0.047 :: Nb_iter = 8 ipsi_step = 0.500000
t_error = 2.249 counter = 1 t_iter = 7 it_step = 0.5
psi_error = 0.002 :: Nb_iter = 4 ipsi_step = 0.500000
t_error = 0.793 counter = 2 t_iter = 6 it_step = 0.5
psi_error = 0.001 :: Nb_iter = 3 ipsi_step = 0.500000
t_error = 0.286 counter = 3 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.104 counter = 4 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.04 counter = 5 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.017 counter = 6 t_iter = 0 it_step = 0.5
---------------------------
psi_soil -0.2000
psi_collar -0.2266
psi_leaf -0.2358
gs: 0.4526
flux H2O 180.1850
flux C2O 8.440244346879666
Tleaf 27.24 Tair 28.17
========================================================================
========================================================================
Date: 2012-08-01 15:00:00
psi_error = 0.046 :: Nb_iter = 8 ipsi_step = 0.500000
t_error = 2.552 counter = 1 t_iter = 7 it_step = 0.5
psi_error = 0.003 :: Nb_iter = 5 ipsi_step = 0.500000
t_error = 0.866 counter = 2 t_iter = 6 it_step = 0.5
psi_error = 0.001 :: Nb_iter = 3 ipsi_step = 0.500000
t_error = 0.316 counter = 3 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 2 ipsi_step = 0.500000
t_error = 0.12 counter = 4 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.045 counter = 5 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.019 counter = 6 t_iter = 0 it_step = 0.5
---------------------------
psi_soil -0.2000
psi_collar -0.2246
psi_leaf -0.2337
gs: 0.4490
flux H2O 168.2701
flux C2O 8.359262513747856
Tleaf 27.35 Tair 28.63
========================================================================
========================================================================
Date: 2012-08-01 16:00:00
psi_error = 0.038 :: Nb_iter = 8 ipsi_step = 0.500000
t_error = 2.037 counter = 1 t_iter = 7 it_step = 0.5
psi_error = 0.002 :: Nb_iter = 4 ipsi_step = 0.500000
t_error = 0.786 counter = 2 t_iter = 6 it_step = 0.5
psi_error = 0.001 :: Nb_iter = 3 ipsi_step = 0.500000
t_error = 0.307 counter = 3 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 2 ipsi_step = 0.500000
t_error = 0.12 counter = 4 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.047 counter = 5 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.019 counter = 6 t_iter = 0 it_step = 0.5
---------------------------
psi_soil -0.2000
psi_collar -0.2185
psi_leaf -0.2269
gs: 0.4259
flux H2O 129.2807
flux C2O 8.105457749145756
Tleaf 25.56 Tair 26.67
========================================================================
========================================================================
Date: 2012-08-01 17:00:00
psi_error = 0.032 :: Nb_iter = 8 ipsi_step = 0.500000
t_error = 1.641 counter = 1 t_iter = 7 it_step = 0.5
psi_error = 0.002 :: Nb_iter = 4 ipsi_step = 0.500000
t_error = 0.66 counter = 2 t_iter = 6 it_step = 0.5
psi_error = 0.001 :: Nb_iter = 3 ipsi_step = 0.500000
t_error = 0.268 counter = 3 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 2 ipsi_step = 0.500000
t_error = 0.109 counter = 4 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.045 counter = 5 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.018 counter = 6 t_iter = 0 it_step = 0.5
---------------------------
psi_soil -0.2000
psi_collar -0.2136
psi_leaf -0.2215
gs: 0.3735
flux H2O 96.7823
flux C2O 7.270791714292491
Tleaf 24.34 Tair 25.39
========================================================================
========================================================================
Date: 2012-08-01 18:00:00
psi_error = 0.023 :: Nb_iter = 7 ipsi_step = 0.500000
t_error = 1.671 counter = 1 t_iter = 7 it_step = 0.5
psi_error = 0.001 :: Nb_iter = 3 ipsi_step = 0.500000
t_error = 0.696 counter = 2 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 2 ipsi_step = 0.500000
t_error = 0.292 counter = 3 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 2 ipsi_step = 0.500000
t_error = 0.123 counter = 4 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.052 counter = 5 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.022 counter = 6 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.01 counter = 7 t_iter = 0 it_step = 0.5
---------------------------
psi_soil -0.2000
psi_collar -0.2068
psi_leaf -0.2138
gs: 0.2286
flux H2O 47.8853
flux C2O 4.350952108728501
Tleaf 23.22 Tair 24.32
========================================================================
========================================================================
Date: 2012-08-01 19:00:00
psi_error = 0.016 :: Nb_iter = 6 ipsi_step = 0.500000
t_error = 2.149 counter = 1 t_iter = 7 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 1.015 counter = 2 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 2 ipsi_step = 0.500000
t_error = 0.481 counter = 3 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.23 counter = 4 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.11 counter = 5 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.054 counter = 6 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.027 counter = 7 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.014 counter = 8 t_iter = 0 it_step = 0.5
---------------------------
psi_soil -0.2000
psi_collar -0.2021
psi_leaf -0.2085
gs: 0.0686
flux H2O 12.8998
flux C2O 0.7577001469349423
Tleaf 21.62 Tair 23.20
========================================================================
========================================================================
Date: 2012-08-01 20:00:00
psi_error = 0.015 :: Nb_iter = 6 ipsi_step = 0.500000
t_error = 2.507 counter = 1 t_iter = 7 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 1.259 counter = 2 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.633 counter = 3 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.319 counter = 4 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.169 counter = 5 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.091 counter = 6 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.049 counter = 7 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.027 counter = 8 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.014 counter = 9 t_iter = 0 it_step = 0.5
---------------------------
psi_soil -0.2000
psi_collar -0.2009
psi_leaf -0.2072
gs: 0.0200
flux H2O 3.5279
flux C2O -0.4988267577765756
Tleaf 21.29 Tair 23.22
========================================================================
========================================================================
Date: 2012-08-01 21:00:00
psi_error = 0.015 :: Nb_iter = 6 ipsi_step = 0.500000
t_error = 2.991 counter = 1 t_iter = 8 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 1.505 counter = 2 t_iter = 7 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.759 counter = 3 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.388 counter = 4 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.213 counter = 5 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.117 counter = 6 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.065 counter = 7 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.035 counter = 8 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.02 counter = 9 t_iter = 0 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.015 counter = 10 t_iter = 0 it_step = 0.25
---------------------------
psi_soil -0.2000
psi_collar -0.2008
psi_leaf -0.2071
gs: 0.0200
flux H2O 2.9644
flux C2O -0.47801247774871825
Tleaf 20.62 Tair 22.95
========================================================================
========================================================================
Date: 2012-08-01 22:00:00
psi_error = 0.014 :: Nb_iter = 6 ipsi_step = 0.500000
t_error = 4.027 counter = 1 t_iter = 8 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 2.037 counter = 2 t_iter = 7 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 1.033 counter = 3 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.569 counter = 4 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.323 counter = 5 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.184 counter = 6 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.105 counter = 7 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.059 counter = 8 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.034 counter = 9 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.019 counter = 10 t_iter = 0 it_step = 0.5
---------------------------
psi_soil -0.2000
psi_collar -0.2007
psi_leaf -0.2069
gs: 0.0200
flux H2O 1.9235
flux C2O -0.446092964897299
Tleaf 19.57 Tair 22.76
========================================================================
========================================================================
Date: 2012-08-01 23:00:00
psi_error = 0.014 :: Nb_iter = 6 ipsi_step = 0.500000
t_error = 4.085 counter = 1 t_iter = 8 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 2.068 counter = 2 t_iter = 7 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 1.049 counter = 3 t_iter = 6 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.58 counter = 4 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.331 counter = 5 t_iter = 5 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.189 counter = 6 t_iter = 4 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.108 counter = 7 t_iter = 3 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.062 counter = 8 t_iter = 2 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.035 counter = 9 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.02 counter = 10 t_iter = 1 it_step = 0.5
psi_error = 0.0 :: Nb_iter = 1 ipsi_step = 0.500000
t_error = 0.016 counter = 11 t_iter = 0 it_step = 0.25
---------------------------
psi_soil -0.2000
psi_collar -0.2005
psi_leaf -0.2067
gs: 0.0200
flux H2O 0.4022
flux C2O -0.44303677153944687
Tleaf 19.47 Tair 22.70
========================================================================
beg time 2025-12-08 15:55:42.874062
end time 2025-12-08 15:56:57.032907
--- Total runtime: 74 sec ---
[16]:
| Rg | An | E | Tleaf | psi_soil | psi_collar | psi_leaf | theta_soil | |
|---|---|---|---|---|---|---|---|---|
| 2012-08-01 00:00:00 | 0.000000 | -0.342269 | 2.686274 | 15.579503 | -0.2 | -0.200797 | -0.207051 | 0.057313 |
| 2012-08-01 01:00:00 | 0.000000 | -0.371843 | 3.717761 | 16.789593 | -0.2 | -0.200928 | -0.207199 | 0.057313 |
| 2012-08-01 02:00:00 | 0.000000 | -0.395888 | 4.650826 | 17.714749 | -0.2 | -0.201048 | -0.207335 | 0.057313 |
| 2012-08-01 03:00:00 | 0.000000 | -0.388431 | 4.129240 | 17.421603 | -0.2 | -0.200982 | -0.207260 | 0.057313 |
| 2012-08-01 04:00:00 | 0.000000 | -0.360403 | 2.761584 | 16.301959 | -0.2 | -0.200808 | -0.207065 | 0.057313 |
| 2012-08-01 05:00:00 | 0.389796 | -0.224113 | 3.607560 | 16.228794 | -0.2 | -0.200916 | -0.207186 | 0.057313 |
| 2012-08-01 06:00:00 | 7.839554 | 1.993460 | 22.954559 | 18.461949 | -0.2 | -0.203430 | -0.210016 | 0.057313 |
| 2012-08-01 07:00:00 | 57.518970 | 7.435500 | 94.706417 | 22.659515 | -0.2 | -0.213285 | -0.221050 | 0.057313 |
| 2012-08-01 08:00:00 | 34.072722 | 6.673350 | 78.418502 | 23.254350 | -0.2 | -0.210995 | -0.218484 | 0.057313 |
| 2012-08-01 09:00:00 | 63.225379 | 8.150045 | 134.706708 | 25.438236 | -0.2 | -0.219289 | -0.227703 | 0.057313 |
| 2012-08-01 10:00:00 | 71.579147 | 8.453443 | 178.452534 | 27.135245 | -0.2 | -0.226247 | -0.235405 | 0.057313 |
| 2012-08-01 11:00:00 | 69.847101 | 8.327728 | 187.769425 | 26.339609 | -0.2 | -0.227808 | -0.237130 | 0.057313 |
| 2012-08-01 12:00:00 | 62.987194 | 8.148085 | 159.026378 | 26.067356 | -0.2 | -0.223097 | -0.231926 | 0.057313 |
| 2012-08-01 13:00:00 | 64.213154 | 8.111484 | 187.656202 | 26.169649 | -0.2 | -0.227789 | -0.237129 | 0.057313 |
| 2012-08-01 14:00:00 | 63.129949 | 8.440244 | 180.184993 | 27.243201 | -0.2 | -0.226581 | -0.235814 | 0.057313 |
| 2012-08-01 15:00:00 | 57.564103 | 8.359263 | 168.270090 | 27.353405 | -0.2 | -0.224633 | -0.233671 | 0.057313 |
| 2012-08-01 16:00:00 | 51.941885 | 8.105458 | 129.280656 | 25.563003 | -0.2 | -0.218493 | -0.226852 | 0.057313 |
| 2012-08-01 17:00:00 | 38.415012 | 7.270792 | 96.782262 | 24.339307 | -0.2 | -0.213647 | -0.221454 | 0.057313 |
| 2012-08-01 18:00:00 | 18.804525 | 4.350952 | 47.885339 | 23.217439 | -0.2 | -0.206760 | -0.213755 | 0.057313 |
| 2012-08-01 19:00:00 | 4.076470 | 0.757700 | 12.899775 | 21.624065 | -0.2 | -0.202120 | -0.208542 | 0.057313 |
| 2012-08-01 20:00:00 | 0.000000 | -0.498827 | 3.527888 | 21.290430 | -0.2 | -0.200907 | -0.207176 | 0.057313 |
| 2012-08-01 21:00:00 | 0.000000 | -0.478012 | 2.964441 | 20.622700 | -0.2 | -0.200834 | -0.207093 | 0.057313 |
| 2012-08-01 22:00:00 | 0.000000 | -0.446093 | 1.923494 | 19.570359 | -0.2 | -0.200703 | -0.206946 | 0.057313 |
| 2012-08-01 23:00:00 | 0.000000 | -0.443037 | 0.402219 | 19.465893 | -0.2 | -0.200508 | -0.206726 | 0.057313 |
The same way we visualized the “Length” property, we can visualize now more specific property once the model is run. For instance, we display the shoot through the stomatal conductance property.
[17]:
MTG(g, 'gs')
[17]:
[ ]: