The reference model code needed is available on ModelDB Accession Number 182797:
https://senselab.med.yale.edu/ModelDB/ShowModel.cshtml?model=182797

Then, merge the files in this directory with the reference model code.


Example usage
-------------

Simulations with Poisson synaptic inputs inserted into the soma, no modulation:

nrngui -c 'rank=326' -c 'blh=0' -c 'scen_type=4' -c 'atol_=1e-06' -c 'ts=10000' syn_inputs.hoc

Simulations with Poisson synaptic inputs inserted into the soma and dendrites, no modulation:

nrngui -c 'rank=326' -c 'blh=0' -c 'scen_type=5' -c 'atol_=1e-06' -c 'ts=10000' syn_inputs_dend.hoc

Simulations with Poisson synaptic inputs inserted into the soma, with modulation:

nrngui -c 'rank=326' -c 'blh=0' -c 'scen_type=4' -c 'atol_=1e-06' -c 'ts=10000' -c 'mo=0.4' -c 'mf=8' syn_inputs_mod.hoc

Simulations with Poisson synaptic inputs inserted into the soma and dendrites, with modulation:

nrngui -c 'rank=326' -c 'blh=0' -c 'scen_type=5' -c 'atol_=1e-06' -c 'ts=10000' -c 'mo=0.4' -c 'mf=8' syn_inputs_mod_dend.hoc



Command-line arguments
-----------------------

Explanation for above and additional model invocation parameters are as follows.

'rank' - Model rank; see Table 1 in paper. Each rank comes with a set of fitted passive
         properties and active conductance densities.
'blh' - set to 1 to block Ih; 0 otherwise
'blks' - set to 1 to block Kdrs; 0 otherwise
'blka' - set to 1 to block KA; 0 otherwise
'atol' - absolute tolerance for adaptive timestep integration (CVODE) in NEURON
'ts' - duration of simulation (ms)
'excw' - set specific excitatory peak conductance (uS)
'inhw' - set specific excitatory peak conductance (uS)
'ia' - additional somatic bias current (nA)
'mo' - modulation amount as percentage (real number [0..1])
'mf' - modulation frequency (Hz)
