from parameters import * from subsystem import * from geometry import * # ---- instantiation ---- # ---- release sites ---- rel_A = m.ReleaseSite(name='rel_A', complex=m.Complex('A', orientation=m.Orientation.UP), region=up, number_to_release=100) rel_B = m.ReleaseSite(name='rel_B', complex=m.Complex('B', orientation=m.Orientation.UP), region=bottom, number_to_release=100) # ---- surface classes assignment ---- # ---- compartments assignment ---- # ---- create instantiation object and add components ---- instantiation = m.Instantiation() instantiation.add_geometry_object(up) instantiation.add_geometry_object(bottom) instantiation.add_release_site(rel_A) instantiation.add_release_site(rel_B)
import os MCELL_PATH = os.environ.get('MCELL_PATH', '') if MCELL_PATH: sys.path.append(os.path.join(MCELL_PATH, 'lib')) else: print("Error: variable MCELL_PATH that is used to find the mcell library was not set.") sys.exit(1) import mcell as m # single property box1 = m.geometry_utils.create_box('b1', 1) box2 = m.geometry_utils.create_box('b1', 1.1) i1 = m.Instantiation() i1.add_geometry_object(box1) # error try: i1.add_geometry_object(box2) assert False except ValueError as err: print(err) box3 = m.geometry_utils.create_box('b2', 1) box4 = m.geometry_utils.create_box('b2', 1) i2 = m.Instantiation() i2.add_geometry_object(box3) # warning