NRNUNIT_USE_LEGACY: consider using old NEURON constant values

[anilbey]

Problem

NEURON 8.0 introduces a breaking change in the constant values [1, 2].

The NEURON our supercomputer (the one used in generating the packages) does not take this change into consideration (therefore uses the old constants).

The NEURON on Pypi does take this change into consideration (uses the new constants).

Here are the changes in those values over years. R is some gas constant. FARADAY is the Faraday constant.

Year     R           FARADAY   
2014    8.3144598   96485.33289
2010    8.3144621   96485.3365
2002    8.314472    96485.3383
1998    8.314472    96485.3415
1986    8.314510    96485.309
1973    8.31441     96486.95

Solution

We probably need to set the NRNUNIT_USE_LEGACY flag in EmodelRunner so that users will use the old constants as we used them in our supercomputer.

We can set the parameter either with export NRNUNIT_USE_LEGACY=1 or with h.nrnunit_use_legacy(True) [3].

References

1. https://nrn.readthedocs.io/en/8.2.2/changelog.html#id23
2. https://www.neuron.yale.edu/phpBB/viewtopic.php?t=3783
3. https://nrn.readthedocs.io/en/8.2.2/python/envvariables.html#nrnunit-use-legacy