storemode

Sometimes you will be interested in the contribution of a single channel type to the behavior of a complex multi-compartmental model, say you are interested in the overall $\mathrm{Ca}^{2+}$ influx. Computing this variable would imply getting the current or conductance from each channel of the specified type in each compartment and then adding them up. From viewpoint of single elements that act as a model description for hsolve, such values are simply not accessible, nevertheless since hsolve has access to all computed variables in chanmode 4 and 5, it will compute such an integrated current or conductance when the storemode field is set. You will not be able to find these variables with the findsolvefield command, but with a small trick, hsolve will report in your terminal where to find these variables.

These variables are defined as the summation of the field Ik or Gk for each channel of the same type in the model. For this to work correctly, each channel of the same type must have the same name in each compartment. Then, before the SETUP action, the storemode field should be set to one of the following values:

0: the default value, no sums are stored.

1: total currents are stored.

2: total conductances are stored.

Afterwards, you use hsolve as usual (set other fields, call SETUP etc.).

Now if you want hsolve to report where it will store the totals of currents or conductances for each channel type, set the silent flag of Genesis to a negative value, call the RESET action on the hsolve element and reset the silent flag to its original value. Hsolve will then report where to find the computed fields for each channel type (or better for each unique channel name that is encountered in the model):

    genesis > silent -1
    genesis > call solver RESET
    storing leak in itotal[0]
    storing CaP in itotal[1]
    storing KC in itotal[2]
    storing K2 in itotal[3]
    storing Ca_pool in itotal[4]
    storing Ca_nernst in itotal[5]
    storing stellate in itotal[6]
    storing Kdr in itotal[7]
    storing NaF in itotal[8]
    transferring element field values into solve arrays
    genesis > silent 0

Now you can use the getfield and showfield commands to examine the values of the variables.

    genesis > step
    genesis > showfield solve itotal[1]
    [ /Purkinje/solve ]
    itotal[3]            = 8.777549691e-10

The storemode field does not influence any other functionality of hsolve. So you are still able to use findsolvefield to inspect other computed values of interest.