# Range of E/I synaptic conductances
    Nvals = 31  # Number of values for each dimension
    startG = 0.0  # nS
    endG = 6120.0  # nS

    bumpCurrentSlope = getBumpCurrentSlope(p['noise_sigma'],
                                           threshold=-np.infty,
                                           type='no_theta')
    extraIterparams = {'bumpCurrentSlope': bumpCurrentSlope}
    #extraIterparams['bumpCurrentSlope'] = [1.0]

    ###############################################################################
    submitParamSweep(p,
                     startG,
                     endG,
                     Nvals,
                     ENV,
                     simRootDir,
                     simLabel,
                     appName,
                     rtLimit,
                     numCPU,
                     blocking,
                     timePrefix,
                     numRepeat,
                     dry_run,
                     extraIterparams,
                     rc=parser.rowcol,
                     extra_qsub_params=o.extra_qsub_params)
    p['master_seed']      = 123456
    p['time']             = o.time or 600e3  # ms
    p['nthreads']         = 1
    p['ntrials']          = o.ntrials
    p['velON']            = 1
    p['constantPosition'] = 0
    p['verbosity']        = o.verbosity

    # No theta parameters
    p['Iext_e_const'] = 487.5   # pA
    p['Iext_i_const'] = 212.5   # pA
    p['Iext_e_theta'] = 0       # pA
    p['Iext_i_theta'] = 0       # pA

    # Range of E/I synaptic conductances
    Nvals  = 31      # Number of values for each dimension
    startG = 0.0     # nS
    endG   = 6120.0  # nS

    bumpCurrentSlope = getBumpCurrentSlope(p['noise_sigma'],
            threshold=-np.infty, type='no_theta')
    extraIterparams = {'bumpCurrentSlope' : bumpCurrentSlope}
    #extraIterparams['bumpCurrentSlope'] = [1.0]

    ###############################################################################
    submitParamSweep(p, startG, endG, Nvals, ENV, simRootDir, simLabel,
                     appName, rtLimit, numCPU, blocking, timePrefix, numRepeat,
                     dry_run, extraIterparams, rc=parser.rowcol,
                     extra_qsub_params=o.extra_qsub_params)
    simLabel    = '{0}pA'.format(int(p['noise_sigma']))
    appName     = 'simulation_velocity.py'
    rtLimit     = o.rtLimit or '12:00:00'
    numCPU      = 1
    blocking    = True
    timePrefix  = False
    numRepeat   = 1
    dry_run     = o.dry_run


    p['master_seed'] = 123456
    p['time']        = 10e3  # ms
    p['nthreads']    = 1
    p['ntrials']     = 10

    p['IvelMax']     = 100
    p['dIvel']       = 10

    p['verbosity']   = 'DEBUG'


    # Range of E/I synaptic conductances
    Nvals  = 31      # Number of values for each dimension
    startG = 0.0     # nS
    endG   = 6120.0  # nS

    ###############################################################################

    submitParamSweep(p, startG, endG, Nvals, ENV, simRootDir, simLabel,
            appName, rtLimit, numCPU, blocking, timePrefix, numRepeat, dry_run)
    p['ntrials'] = o.ntrials
    p['bumpCurrentSlope'] = -1
    p['velON'] = 0
    p['constantPosition'] = 1
    p['pcON'] = 1
    p['stateMonDur'] = 1e3  # ms
    p['verbosity'] = o.verbosity

    # Range of E/I synaptic conductances
    Nvals = 31  # Number of values for each dimension
    startG = 0.0  # nS
    endG = 6120.0  # nS

    ###############################################################################
    submitParamSweep(p,
                     startG,
                     endG,
                     Nvals,
                     ENV,
                     simRootDir,
                     simLabel,
                     appName,
                     rtLimit,
                     numCPU,
                     blocking,
                     timePrefix,
                     numRepeat,
                     dry_run,
                     rc=parser.rowcol,
                     printout=o.printout)