This entry contains the scripts for simulating the model of our article "Excitatory and inhibitory effects of HCN channel modulation on excitability of layer V pyramidal cells".
Two folders are included in this ModelDB entry:
`modulhcn_almog`
`modulhcn_hay`
Before running any simulations, run "nrnivmodl" in both folders to compile the .mod mechanisms.
To run the simulations for Fig. 1 and draw the results, run the following commands:

    cd modulhcn_almog
    sh runfig1.sh
    cd ../modulhcn_hay
    sh runfig1.sh
    cd ..
    python3 drawfig1.py

In the same way, to draw Fig. 2-6, do the same procedure, replacing '1' by numbers 2-6. Note that Fig. 3 and Fig. 5 re-use data from Fig. 2, so run runfig2.sh in addition to runfig3.sh and runfig5.sh to be able to draw the figures (and similarly, Fig. 6 re-uses data from Fig. 4).
Approximate simulation times:

    modulhcn_almog/runfig1.sh: 1-2 days
    modulhcn_almog/runfig2.sh: 2-3 hours
    modulhcn_almog/runfig3.sh: 2 hours
    modulhcn_almog/runfig4.sh: 4-7 days

    modulhcn_hay/runfig1.sh: 3-6 hours
    modulhcn_hay/runfig2.sh: 20 min
    modulhcn_hay/runfig3.sh: 30 min
    modulhcn_hay/runfig4.sh: 1-2 days
    modulhcn_hay/runfig5.sh: 5-10 days
    modulhcn_hay/runfig6.sh: 2-3 days

### Changelog

2022-12: MOD-C-2-CPP: Compatibility fixed for NEURON 9.0, 8.2 and 8.1 Fix for 9.0.0: curr used as both variable and function in file epsp.mod