www.neuron.yale.edu

Thank you for the quick response!

Hi, I was trying to make this mod file work with variable time step. The mod file uses dt to calculate the mexp and hexp PROCEDURE trates(v) { LOCAL tinc TABLE minf, mexp, hinf, hexp DEPEND dt FROM vmin TO vmax WITH 199 rates(v): not consistently executed from here if usetable == 1 tinc = -dt mexp =...

Thank you, Ted!

Hi, I am trying to plot voltage clamp currents of multiple channels for different mod files. For this I use Vector.record(). As there as multiple channels, I am passing a string to Vector.record(). As a test case, I tried it with hh channels. 'hh2' is same as hh with ik made a RANGE variable. channe...

Hi, If you are interested in clamping concentration to a particular value using a mod mechanism, see this code: https://github.com/darshanmandge/Tools/tree/master/ConcClamp. Screenshot here: https://github.com/darshanmandge/Tools/blob/master/README.md (2. Ionic Concentration Clamp) Cheers, Darshan

Thank you for clearing these doubts!

Hi, I have a doubt in mA/cm2 to nA conversion: if I want to convert the current density INA [mA/cm2] in nA I should divide the current in mA/cm2 by the factor (area*(0.01)): INA/(area*(0.01)) True, assuming that area is in units of um2. Shouldn't the current in mA/cm2 be multiplied by (area*0.01) to...

Thanks for the prompt reply and example of CONDUCTANCE keyword.

Thanks again for the elaborate explanation. So, if those statements (minf, ninf) are in BREAKPOINT, are they undergoing unnecessary recalculations? I guess by moving them to the DERIVATIVE block will save some simulation time which might be more evident in long or large-scale simulations, Also, I co...

Thank you Michael. Using the recursive multiplication approach works for these expressions when using derivimplicit. I also tried using function and LOCAL. It gave me the same results both of them (though not desirable). I have to fix the parameter values. I will use functions for greater power valu...

Thanks. This shortcut for copying text works perfectly fine. I guess the key assigned for copy was changed to 'return' in newer NEURON versions.

I will also post here if I find a solution/workaround to the compilation error.

Best Wishes,
Darshan

Thank you Ted for your prompt reply! 1. I expanded the terms for adpi' and atpi' expressions: atpi' = ((nu - atpi* (1+(k2^2)*(adpi^2)+2*k2*adpi))/tauc) adpi' = ((atpi*(1+(k2^2)*(adpi^2)+2*k2*adpi) - eta*adpi)/tauc) and also changed the "k" which I declared as constant earlier to "k2&q...

Hi, I was trying to translate this code given here: https://www.math.fsu.edu/~bertram/software/islet/JTB_04.ode for pancreatic beta-cells, published by Wierschem and Bertram, J. Theor. Biol., 228:513-521, 2004. Here is the nmodl code I made: TITLE weirshem and Bertram 2004 UNITS { (molar) = (1/liter...

Hi, I get the following error when I run the working example given in the first post: C:\nrn\bin\nrniv.exe: cai not a section variable in D:/NEURON_Tools/rxd_hoc/rxd_hoc.hoc near line 32 g.addvar("ca", &soma.cai(0.5), 1, 1) NEURON -- VERSION 7.6.2 master (f5a1284) 2018-08-15 Duke, Yale...

Thanks for your help, Ted!

Darshan