In [2]:

    

%matplotlib inline
import moose

import pylab
import rdesigneur as rd


    

In [3]:

    

rdes = rd.rdesigneur(
    turnOffElec = True,
    chemProto = [['makeChemOscillator()', 'osc']],
    chemDistrib = [['osc', 'soma', 'install', '1' ]],
    plotList = [['soma', '1', 'dend/a', 'conc', 'Concentration of a'],
        ['soma', '1', 'dend/b', 'conc', 'Concentration of b']],
    moogList = [['soma', '1', 'dend/a', 'conc', 'a Conc', 0, 360 ]]
)


    

In [4]:

    

rdes.buildModel()
bv = moose.vec( '/model/chem/dend/b' )
bv[0].concInit *= 2
bv[-1].concInit *= 2
moose.reinit()


    

    




Rdesigneur: Elec model has 1 compartments and 0 spines on 1 compartments.
Chem part of model has  250 dendrite voxels X 3 pools,
    

In [5]:

    

c=bv[0]
print c.vec[1].valueFields


    

    




('this', 'name', 'me', 'parent', 'children', 'path', 'className', 'numData', 'numField', 'idValue', 'index', 'fieldIndex', 'tick', 'dt', 'valueFields', 'sourceFields', 'destFields', 'msgOut', 'msgIn', 'n', 'nInit', 'diffConst', 'motorConst', 'conc', 'concInit', 'volume', 'speciesId')

In [6]:

    

rdes.displayMoogli( 1, 400, 0.001 ) # dt, runtime, rotation