honeybee-grasshopper-radiance 1.35.1__py3-none-any.whl

honeybee_grasshopper_radiance/__init__.py

@@ -0,0 +1,7 @@
+"""Honeybee Grasshopper Radiance Plugin.
+
+Note that this package is not intended to run with cPython and it only possesses
+the Grasshopper components. In order to run the plugin, the core libraries must
+be installed in a manner that they can be discovered by Rhino.
+The package includes both the userobjects (.ghuser) and the Python source (.py).
+"""

honeybee_grasshopper_radiance/src/HB Adjust HDR.py

@@ -0,0 +1,107 @@
+# This file is part of Honeybee.
+#
+# Copyright (c) 2025, Ladybug Tools.
+# You should have received a copy of the GNU Affero General Public License
+# along with Honeybee; If not, see <http://www.gnu.org/licenses/>.
+# 
+# @license AGPL-3.0-or-later <https://spdx.org/licenses/AGPL-3.0-or-later>
+
+"""
+Adjust and format a High Dynamic Range (HDR) image file.
+_
+Possible adjustments include chaging the exposure of the image to mimic what would
+be seen by a human eye and adding an optional text label to the image.
+-
+
+    Args:
+        _hdr: Path to a High Dynamic Range (HDR) image file.
+        adj_expos_: Boolean to note whether the exposure of the image should be adjusted to
+            mimic the human visual response in the output. The goal of this process
+            is to output an image that correlates more strongly with a person’s
+            subjective impression of a scene rather than the absolute birghtness
+            of the scene. (Default: False).
+        label_: Optional text label to be appended to the bottom of the image. This
+            is useful when one has several images and would like to easily
+            identify them while scrolling through them.
+        label_hgt_: An integer for the height of the label text in pixels. (Default: 32).
+
+    Returns:
+        hdr: Path to the resulting adjusted HDR image file.
+"""
+
+ghenv.Component.Name = 'HB Adjust HDR'
+ghenv.Component.NickName = 'AdjustHDR'
+ghenv.Component.Message = '1.9.0'
+ghenv.Component.Category = 'HB-Radiance'
+ghenv.Component.SubCategory = '4 :: Results'
+ghenv.Component.AdditionalHelpFromDocStrings = '3'
+
+import os
+
+try:  # import honeybee_radiance_command dependencies
+    from honeybee_radiance_command.pcond import Pcond
+    from honeybee_radiance_command.psign import Psign
+    from honeybee_radiance_command.pcompos import Pcompos
+except ImportError as e:
+    raise ImportError('\nFailed to import honeybee_radiance_command:\n\t{}'.format(e))
+
+try:  # import honeybee_radiance dependencies
+    from honeybee_radiance.config import folders as rad_folders
+except ImportError as e:
+    raise ImportError('\nFailed to import honeybee_radiance:\n\t{}'.format(e))
+
+try:  # import ladybug_rhino dependencies
+    from ladybug_rhino.grasshopper import all_required_inputs
+except ImportError as e:
+    raise ImportError('\nFailed to import ladybug_rhino:\n\t{}'.format(e))
+
+# check the Radiance date of the installed radiance
+try:  # import lbt_recipes dependencies
+    from lbt_recipes.version import check_radiance_date
+except ImportError as e:
+    raise ImportError('\nFailed to import lbt_recipes:\n\t{}'.format(e))
+check_radiance_date()
+
+
+if all_required_inputs(ghenv.Component):
+    # set up the paths for the various files used in translation
+    img_dir = os.path.dirname(_hdr)
+    input_image = os.path.basename(_hdr)
+    hdr = _hdr
+    commands = []
+
+    # add the command to adjust the exposure to reflect human vision (if requested)
+    if adj_expos_:
+        adj_image = input_image.lower().replace('.hdr', '_h.HDR')
+        pcond = Pcond(input=input_image, output=adj_image)
+        pcond.options.h = True
+        commands.append(pcond)
+        hdr = os.path.join(img_dir, adj_image)
+        input_image = adj_image
+
+    # add the command to for a text label (if requested)
+    if label_:
+        label_images = []
+        for i, l_tex in enumerate(reversed(label_.split('\n'))):
+            label_image = 'label{}.HDR'.format(i)
+            psign = Psign(text=l_tex, output=label_image)
+            psign.options.cb = (0, 0, 0)
+            psign.options.cf = (1, 1, 1)
+            psign.options.h = label_hgt_ if label_hgt_ is not None else 32
+            commands.append(psign)
+            label_images.append(label_image)
+
+        lbl_image = input_image.lower().replace('.hdr', '_label.HDR')
+        pcompos = Pcompos(input=label_images + [input_image], output=lbl_image)
+        pcompos.options.a = 1
+        commands.append(pcompos)
+        hdr = os.path.join(img_dir, lbl_image)
+        input_image = lbl_image
+
+    # run the commands in series and load the global horizontal irradiance
+    env = None
+    if rad_folders.env != {}:
+        env = rad_folders.env
+    env = dict(os.environ, **env) if env else None
+    for r_cmd in commands:
+        r_cmd.run(env, cwd=img_dir)

honeybee_grasshopper_radiance/src/HB Ambient Resolution.py

@@ -0,0 +1,63 @@
+# Honeybee: A Plugin for Environmental Analysis
+# This file is part of Honeybee.
+#
+# Copyright (c) 2025, Ladybug Tools.
+# You should have received a copy of the GNU Affero General Public License
+# along with Honeybee; If not, see <http://www.gnu.org/licenses/>.
+# 
+# @license AGPL-3.0-or-later <https://spdx.org/licenses/AGPL-3.0-or-later>
+
+"""
+Get the recommended ambient resoluation (-ar) needed to resolve details with a
+given dimension in model units.
+_
+This recommendation is derived from the overall dimensions of the Radince scene
+being simulated as well as the ambient accuracy (-aa) being used in the simulation.
+_
+The result from this component can be plugged directly into the additional_par_
+of the "HB Radiance Parameter" component or into the radiance_par of any
+recipe components.
+
+-
+    Args:
+        _model: The Honeybee Model being used for Radiance simulation.
+        _detail_dim: A number in model units that represents the dimension of the
+            smallest detail that must be resolved in the Radiance simulation.
+        _aa_: An number for ambient accuracy (-aa) being used in the Radiance smiulation.
+            This value should be matched between this component and the component
+            into which the ouput ar is being input. (Default: 0.25 for
+            low-resolution Radiance studies).
+
+    Returns:
+        ar: The abmient resolution needed to resolve the _detail_dim as a text string.
+            These can be plugged into the additional_par_ of the "HB Radiance
+            Parameter" component or the radiance_par_ input of the recipes.
+"""
+
+ghenv.Component.Name = 'HB Ambient Resolution'
+ghenv.Component.NickName = 'AR'
+ghenv.Component.Message = '1.9.0'
+ghenv.Component.Category = 'HB-Radiance'
+ghenv.Component.SubCategory = '3 :: Recipes'
+ghenv.Component.AdditionalHelpFromDocStrings = '0'
+
+try:
+    from ladybug_rhino.grasshopper import all_required_inputs, recipe_result
+except ImportError as e:
+    raise ImportError('\nFailed to import ladybug_rhino:\n\t{}'.format(e))
+
+
+if all_required_inputs(ghenv.Component):
+    # set the default -aa
+    aa = 0.25 if _aa_ is None else _aa_
+
+    # get the longest dimension
+    min_pt, max_pt = _model.min, _model.max
+    x_dim = max_pt.x - min_pt.x
+    y_dim = max_pt.y - min_pt.y
+    z_dim = max_pt.z - min_pt.z
+    longest_dim = max((x_dim, y_dim, z_dim))
+    
+    # calculate the ambient resolution.
+    a_res = int((longest_dim * aa) / _detail_dim)
+    ar = '-ar {}'.format(a_res)

honeybee_grasshopper_radiance/src/HB Annual Average Values.py

@@ -0,0 +1,205 @@
+# This file is part of Honeybee.
+#
+# Copyright (c) 2025, Ladybug Tools.
+# You should have received a copy of the GNU Affero General Public License
+# along with Honeybee; If not, see <http://www.gnu.org/licenses/>.
+# 
+# @license AGPL-3.0-or-later <https://spdx.org/licenses/AGPL-3.0-or-later>
+
+"""
+Get average illuminance or irradiance values over an annual daylight or irradiance
+simulation.
+_
+The _hoys_ input can also be used to filter the data for a particular time period or
+hour/timestep of the simulation.
+
+-
+    Args:
+        _results: An list of annual Radiance result files from either the "HB Annual Daylight"
+            or the "HB Annual Irradiance" component (containing the .ill files and
+            the sun-up-hours.txt). This can also be just the path to the folder
+            containing these result files.
+        dyn_sch_: Optional dynamic Aperture Group Schedules from the "HB Aperture Group
+            Schedule" component, which will be used to customize the behavior
+            of any dyanmic aperture geometry in the output metrics. If unsupplied,
+            all dynamic aperture groups will be in their default state in for
+            the output metrics.
+        _hoys_: An optional numbers or list of numbers to select the hours of the year (HOYs)
+            for which results will be computed. These HOYs can be obtained from the
+            "LB Calculate HOY" or the "LB Analysis Period" components. If None, all
+            hours of the results will be used.
+        median_: Set to True to get the median values instead of the average. The median
+            values can only be calculated for a results folder from the
+            "HB Annual Daylight" component. (Default: False).
+        grid_filter_: The name of a grid or a pattern to filter the grids. For instance,
+            first_floor_* will simulate only the sensor grids that have an
+            identifier that starts with first_floor_. By default all the grids
+            will be processed.
+
+    Returns:
+        report: Reports, errors, warnings, etc.
+        values: Average illuminance or irradiance valules for each sensor in lux or W/m2.
+            Each value is for a different sensor of the grid. These can be plugged
+            into the "LB Spatial Heatmap" component along with meshes of the sensor
+            grids to visualize results.
+"""
+
+ghenv.Component.Name = 'HB Annual Average Values'
+ghenv.Component.NickName = 'AvgValues'
+ghenv.Component.Message = '1.9.0'
+ghenv.Component.Category = 'HB-Radiance'
+ghenv.Component.SubCategory = '4 :: Results'
+ghenv.Component.AdditionalHelpFromDocStrings = '2'
+
+import os
+import subprocess
+
+try:
+    from ladybug.futil import write_to_file
+except ImportError as e:
+    raise ImportError('\nFailed to import ladybug:\n\t{}'.format(e))
+
+try:
+    from honeybee.config import folders
+except ImportError as e:
+    raise ImportError('\nFailed to import honeybee:\n\t{}'.format(e))
+
+try:
+    from honeybee_radiance.postprocess.annualdaylight import _process_input_folder
+except ImportError as e:
+    raise ImportError('\nFailed to import honeybee_radiance:\n\t{}'.format(e))
+
+try:
+    from honeybee_radiance_postprocess.dynamic import DynamicSchedule
+except ImportError as e:
+    raise ImportError('\nFailed to import honeybee_radiance:\n\t{}'.format(e))
+
+try:
+    from pollination_handlers.outputs.helper import read_sensor_grid_result
+except ImportError as e:
+    raise ImportError('\nFailed to import pollination_handlers:\n\t{}'.format(e))
+
+try:
+    from ladybug_rhino.grasshopper import all_required_inputs, list_to_data_tree, \
+        give_warning
+except ImportError as e:
+    raise ImportError('\nFailed to import ladybug_rhino:\n\t{}'.format(e))
+
+
+def parse_sun_up_hours(sun_up_hours, hoys, timestep):
+    """Parse the sun-up hours from the result file .txt file.
+
+    Args:
+        sun_up_hours: A list of integers for the sun-up hours.
+        hoys: A list of 8760 * timestep values for the hoys to select. If an empty
+            list is passed, None will be returned.
+        timestep: Integer for the timestep of the analysis.
+    """
+    if len(hoys) != 0:
+        schedule = [False] * (8760 * timestep)
+        for hr in hoys:
+            schedule[int(hr * timestep)] = True
+        su_pattern = [schedule[int(h * timestep)] for h in sun_up_hours]
+        return su_pattern
+
+
+def average_values(ill_file, su_pattern, full_len):
+    """Compute average values for a given result file."""
+    avg_vals = []
+    with open(ill_file) as results:
+        if su_pattern is None:  # no HOY filter on results
+            for pt_res in results:
+                values = [float(r) for r in pt_res.split()]
+                total_val = sum(values)
+                avg_vals.append(total_val / full_len)
+        else: 
+            for pt_res in results:
+                values = [float(r) for r, is_hoy in zip(pt_res.split(), su_pattern) if is_hoy]
+                total_val = sum(values)
+                try:
+                    avg_vals.append(total_val / full_len)
+                except ZeroDivisionError:
+                    avg_vals.append(0)
+    return avg_vals
+
+
+if all_required_inputs(ghenv.Component):
+    # set up the default values
+    median_ = False if median_ is None else median_
+    grid_filter_ = '*' if grid_filter_ is None else grid_filter_
+    res_folder = os.path.dirname(_results[0]) if os.path.isfile(_results[0]) \
+        else _results[0]
+
+    # check to see if results use the newer numpy arrays
+    if os.path.isdir(os.path.join(res_folder, '__static_apertures__')) or \
+        os.path.isfile(os.path.join(res_folder, 'grid_states.json')):
+        res_type = 'average' if median_ is False else 'median'
+        cmds = [folders.python_exe_path, '-m', 'honeybee_radiance_postprocess',
+                'post-process', '{}-values'.format(res_type), res_folder, '-sf',
+                'metrics']
+        if len(_hoys_) != 0:
+            hoys_str = '\n'.join(str(h) for h in _hoys_)
+            hoys_file = os.path.join(res_folder, 'hoys.txt')
+            write_to_file(hoys_file, hoys_str)
+            cmds.extend(['--hoys-file', hoys_file])
+        if grid_filter_ != '*':
+            cmds.extend(['--grids-filter', grid_filter_])
+        if len(dyn_sch_) != 0:
+            if os.path.isfile(os.path.join(res_folder, 'grid_states.json')):
+                dyn_sch = dyn_sch_[0] if isinstance(dyn_sch_[0], DynamicSchedule) else \
+                    DynamicSchedule.from_group_schedules(dyn_sch_)
+                dyn_sch_file = dyn_sch.to_json(folder=res_folder)
+                cmds.extend(['--states', dyn_sch_file])
+            else:
+                msg = 'No dynamic aperture groups were found in the Model.\n' \
+                    'The input dynamic schedules will be ignored.'
+                print(msg)
+                give_warning(ghenv.Component, msg)
+        use_shell = True if os.name == 'nt' else False
+        custom_env = os.environ.copy()
+        custom_env['PYTHONHOME'] = ''
+        process = subprocess.Popen(
+            cmds, cwd=res_folder, shell=use_shell, env=custom_env,
+            stdout=subprocess.PIPE, stderr=subprocess.PIPE)
+        stdout = process.communicate()  # wait for the process to finish
+        if stdout[-1] != '':
+            print(stdout[-1])
+            raise ValueError('Failed to compute {} values.'.format(res_type))
+        res_dir = os.path.join(res_folder, 'metrics', '{}_values'.format(res_type))
+        if os.path.isdir(res_dir):
+            values = read_sensor_grid_result(res_dir, res_type,'full_id', False)
+            values = list_to_data_tree(values)
+
+    else:
+        if len(dyn_sch_) != 0:
+            msg = 'Dynamic Schedules are currently only supported for Annual Daylight ' \
+                'simulations.\nThe input schedules will be ignored.'
+            print(msg)
+            give_warning(ghenv.Component, msg)
+        if median_:
+            raise ValueError('The median values can only be calculated for a '
+                             'results folder from the "HB Annual Daylight" '
+                             'component.')
+        # extract the timestep if it exists
+        timestep = 1
+        tstep_file = os.path.join(res_folder, 'timestep.txt')
+        if os.path.isfile(tstep_file):
+            with open(tstep_file) as tf:
+                timestep = int(tf.readline())
+        full_len = 8760 * timestep if len(_hoys_) == 0 else len(_hoys_)
+
+        # parse the sun-up-hours
+        grids, sun_up_hours = _process_input_folder(res_folder, grid_filter_)
+        su_pattern = parse_sun_up_hours(sun_up_hours, _hoys_, timestep)
+
+        # compute the average values
+        values = []
+        for grid_info in grids:
+            ill_file = os.path.join(res_folder, '%s.ill' % grid_info['full_id'])
+            dgp_file = os.path.join(res_folder, '%s.dgp' % grid_info['full_id'])
+            if os.path.isfile(dgp_file):
+                avg = average_values(dgp_file, su_pattern, full_len)
+            else:
+                avg = average_values(ill_file, su_pattern, full_len)
+            values.append(avg)
+        values = list_to_data_tree(values)

honeybee_grasshopper_radiance/src/HB Annual Cumulative Values.py

@@ -0,0 +1,191 @@
+# This file is part of Honeybee.
+#
+# Copyright (c) 2025, Ladybug Tools.
+# You should have received a copy of the GNU Affero General Public License
+# along with Honeybee; If not, see <http://www.gnu.org/licenses/>.
+# 
+# @license AGPL-3.0-or-later <https://spdx.org/licenses/AGPL-3.0-or-later>
+
+"""
+Get cumulative radiation (or sum of illuminance) values over an annual irradiance
+or daylight simulation.
+_
+The _hoys_ input can also be used to filter the data for a particular time period or
+hour/timestep of the simulation.
+
+-
+    Args:
+        _results: An list of annual Radiance result files from either the "HB Annual Daylight"
+            or the "HB Annual Irradiance" component (containing the .ill files and
+            the sun-up-hours.txt). This can also be just the path to the folder
+            containing these result files.
+        dyn_sch_: Optional dynamic Aperture Group Schedules from the "HB Aperture Group
+            Schedule" component, which will be used to customize the behavior
+            of any dyanmic aperture geometry in the output metrics. If unsupplied,
+            all dynamic aperture groups will be in their default state in for
+            the output metrics.
+        _hoys_: An optional numbers or list of numbers to select the hours of the year (HOYs)
+            for which results will be computed. These HOYs can be obtained from the
+            "LB Calculate HOY" or the "LB Analysis Period" components. If None, all
+            hours of the results will be used.
+        grid_filter_: The name of a grid or a pattern to filter the grids. For instance,
+            first_floor_* will simulate only the sensor grids that have an
+            identifier that starts with first_floor_. By default all the grids
+            will be processed.
+
+    Returns:
+        report: Reports, errors, warnings, etc.
+        values: In the case of an annual irradaince simulation, this is the cumulative
+            radiation valules for each sensor in Wh/m2. For annual daylight, it is
+            cumulative illuminance (lux-hours). These can be plugged into the "LB
+            Spatial Heatmap" component along with meshes of the sensor
+            grids to visualize results.
+"""
+
+ghenv.Component.Name = 'HB Annual Cumulative Values'
+ghenv.Component.NickName = 'CumulValues'
+ghenv.Component.Message = '1.9.0'
+ghenv.Component.Category = 'HB-Radiance'
+ghenv.Component.SubCategory = '4 :: Results'
+ghenv.Component.AdditionalHelpFromDocStrings = '2'
+
+import os
+import subprocess
+
+try:
+    from ladybug.futil import write_to_file
+except ImportError as e:
+    raise ImportError('\nFailed to import ladybug:\n\t{}'.format(e))
+
+try:
+    from honeybee.config import folders
+except ImportError as e:
+    raise ImportError('\nFailed to import honeybee:\n\t{}'.format(e))
+
+try:
+    from honeybee_radiance.postprocess.annualdaylight import _process_input_folder
+except ImportError as e:
+    raise ImportError('\nFailed to import honeybee_radiance:\n\t{}'.format(e))
+
+try:
+    from honeybee_radiance_postprocess.dynamic import DynamicSchedule
+except ImportError as e:
+    raise ImportError('\nFailed to import honeybee_radiance:\n\t{}'.format(e))
+
+try:
+    from pollination_handlers.outputs.helper import read_sensor_grid_result
+except ImportError as e:
+    raise ImportError('\nFailed to import pollination_handlers:\n\t{}'.format(e))
+
+try:
+    from ladybug_rhino.grasshopper import all_required_inputs, list_to_data_tree, \
+        give_warning
+except ImportError as e:
+    raise ImportError('\nFailed to import ladybug_rhino:\n\t{}'.format(e))
+
+
+def parse_sun_up_hours(sun_up_hours, hoys, timestep):
+    """Parse the sun-up hours from the result file .txt file.
+
+    Args:
+        sun_up_hours: A list of integers for the sun-up hours.
+        hoys: A list of 8760 * timestep values for the hoys to select. If an empty
+            list is passed, None will be returned.
+        timestep: Integer for the timestep of the analysis.
+    """
+    if len(hoys) != 0:
+        schedule = [False] * (8760 * timestep)
+        for hr in hoys:
+            schedule[int(hr * timestep)] = True
+        su_pattern = [schedule[int(h * timestep)] for h in sun_up_hours]
+        return su_pattern
+
+
+def cumulative_values(ill_file, su_pattern, timestep):
+    """Compute average values for a given result file."""
+    cumul_vals = []
+    with open(ill_file) as results:
+        if su_pattern is None:  # no HOY filter on results
+            for pt_res in results:
+                values = [float(r) for r in pt_res.split()]
+                cumul_vals.append(sum(values) / timestep)
+        else: 
+            for pt_res in results:
+                values = [float(r) for r, is_hoy in zip(pt_res.split(), su_pattern) if is_hoy]
+                cumul_vals.append(sum(values) / timestep)
+    return cumul_vals
+
+
+if all_required_inputs(ghenv.Component):
+    # set up the default values
+    grid_filter_ = '*' if grid_filter_ is None else grid_filter_
+    res_folder = os.path.dirname(_results[0]) if os.path.isfile(_results[0]) \
+        else _results[0]
+
+    # check to see if results use the newer numpy arrays
+    if os.path.isdir(os.path.join(res_folder, '__static_apertures__'))  or \
+            os.path.isfile(os.path.join(res_folder, 'grid_states.json')):
+        cmds = [folders.python_exe_path, '-m', 'honeybee_radiance_postprocess',
+                'post-process', 'cumulative-values', res_folder, '-sf', 'metrics']
+        if len(_hoys_) != 0:
+            hoys_str = '\n'.join(str(h) for h in _hoys_)
+            hoys_file = os.path.join(res_folder, 'hoys.txt')
+            write_to_file(hoys_file, hoys_str)
+            cmds.extend(['--hoys-file', hoys_file])
+        if grid_filter_ != '*':
+            cmds.extend(['--grids-filter', grid_filter_])
+        if len(dyn_sch_) != 0:
+            if os.path.isfile(os.path.join(res_folder, 'grid_states.json')):
+                dyn_sch = dyn_sch_[0] if isinstance(dyn_sch_[0], DynamicSchedule) else \
+                    DynamicSchedule.from_group_schedules(dyn_sch_)
+                dyn_sch_file = dyn_sch.to_json(folder=res_folder)
+                cmds.extend(['--states', dyn_sch_file])
+            else:
+                msg = 'No dynamic aperture groups were found in the Model.\n' \
+                    'The input dynamic schedules will be ignored.'
+                print(msg)
+                give_warning(ghenv.Component, msg)
+        use_shell = True if os.name == 'nt' else False
+        custom_env = os.environ.copy()
+        custom_env['PYTHONHOME'] = ''
+        process = subprocess.Popen(
+            cmds, cwd=res_folder, shell=use_shell, env=custom_env,
+            stdout=subprocess.PIPE, stderr=subprocess.PIPE)
+        stdout = process.communicate()  # wait for the process to finish
+        if stdout[-1] != '':
+            print(stdout[-1])
+            raise ValueError('Failed to compute cumulative values.')
+        avg_dir = os.path.join(res_folder, 'metrics', 'cumulative_values')
+        if os.path.isdir(avg_dir):
+            values = read_sensor_grid_result(avg_dir, 'cumulative','full_id', False)
+            values = list_to_data_tree(values)
+
+    else:
+        if len(dyn_sch_) != 0:
+            msg = 'Dynamic Schedules are currently only supported for Annual Daylight ' \
+                'simulations.\nThe input schedules will be ignored.'
+            print(msg)
+            give_warning(ghenv.Component, msg)
+
+        # extract the timestep if it exists
+        timestep = 1
+        tstep_file = os.path.join(res_folder, 'timestep.txt')
+        if os.path.isfile(tstep_file):
+            with open(tstep_file) as tf:
+                timestep = int(tf.readline())
+
+        # parse the sun-up-hours
+        grids, sun_up_hours = _process_input_folder(res_folder, grid_filter_)
+        su_pattern = parse_sun_up_hours(sun_up_hours, _hoys_, timestep)
+
+        # compute the average values
+        values = []
+        for grid_info in grids:
+            ill_file = os.path.join(res_folder, '%s.ill' % grid_info['full_id'])
+            dgp_file = os.path.join(res_folder, '%s.dgp' % grid_info['full_id'])
+            if os.path.isfile(dgp_file):
+                cumul = cumulative_values(dgp_file, su_pattern, timestep)
+            else:
+                cumul = cumulative_values(ill_file, su_pattern, timestep)
+            values.append(cumul)
+        values = list_to_data_tree(values)