IMAS consistency check

OMAS can enforce consistency with IMAS data structure. There are three possible settings:

• consistency_check = False: no consistency check. False can be used to use OMAS for other purposes other than IMAS

• consistency_check = ‘warn’: will print a wanrning if entry is outside of IMAS scope. ‘warn’ can be useful when wanting to work with entries that are not yet supported by IMAS. Remember to open a JIRA issue (https://jira.iter.org) to start a conversation about making these entries officially part of IMAS.

• consistency_check = True: will raise an error if entry is outside of IMAS scope. Recommended when working with IMAS. The error raised will clearly say where the error is, and provide suggestions based on valid IMAS structure.

import os

from omas import *

print('*' * 20)
print('consistency_check = False')
print('*' * 20)
ods = ODS().sample()
ods.consistency_check = False
# can add entry with wrong dimensions
ods['equilibrium.time_slice[0].profiles_2d[0].b_tor'] = 1
# can add entry that is obsolescent
ods['equilibrium.time_slice[0].profiles_2d[0].b_tor'] = [[1, 1], [1, 1]]
# can add entry that does not exist in IMAS
ods['equilibrium.time_slice[0].does_not_exist.global_quantities.ip'] = 1

print('*' * 20)
print("consistency_check = 'warn'")
print('*' * 20)
ods = ODS().sample()
ods.consistency_check = 'warn'
# can add entry with wrong dimensions but warning message is printed
ods['equilibrium.time_slice[0].profiles_2d[0].b_tor'] = 1
# can add entry that is obsolescent but warning message is printed
ods['equilibrium.time_slice[0].profiles_2d[0].b_tor'] = [[1, 1], [1, 1]]
# can add entry that does not exist in IMAS but warning message is printed
ods['equilibrium.time_slice[0].does_not_exist.global_quantities.ip'] = 1

In the above example consistency_check = warn or True result in the following warnings:

equilibrium.time_slice.:.profiles_2d.:.b_tor must be an array with dimensions: ['equilibrium.time_slice[:].profiles_2d[:].grid.dim1', 'equilibrium.time_slice[:].profiles_2d[:].grid.dim2']
equilibrium.time_slice.:.profiles_2d.:.b_tor is in OBSOLESCENT state

print('*' * 20)
print("consistency_check = True")
print('*' * 20)
ods = ODS().sample()
ods.consistency_check = True  # this is the default
# add entry with wrong dimensions
try:
    ods['equilibrium.time_slice[0].profiles_2d[0].b_tor'] = 1
    raise Exception('Was able to add an entry with wrong dimensions!')
except ValueError:
    pass
# add entry that is obsolescent
ods['equilibrium.time_slice[0].profiles_2d[0].grid.dim1'] = [1, 2, 3]
ods['equilibrium.time_slice[0].profiles_2d[0].grid.dim2'] = [1, 2, 3, 4]
ods['equilibrium.time_slice[0].profiles_2d[0].b_tor'] = [[1, 1, 1] * 4]
# add entry that does not exist in IMAS
try:
    ods['equilibrium.time_slice[0].does_not_exist.global_quantities.ip'] = 1
    raise Exception('Was able to add an entry with that is not in IMAS even if consistency_check = True')
except LookupError:
    pass

In the above example consistency_check = True result in the following error:

LookupError: Not a valid IMAS 3.18.0 location: `equilibrium.time_slice.:.does_not_exist`
                                                                         ^^^^^^^^^^^^^^
Did you mean: ['profiles_2d', 'ggd', 'boundary', 'profiles_1d', 'constraints', 'global_quantities', 'coordinate_system', 'boundary_separatrix', 'time', 'convergence']