function key = LPplot 
% CHANGE NAME OF FUNCTION FROM LPplot ABOVE TO YOUR FIRST NAME OR PSEUDONYM
% CHANGE THE NAME OF THIS FILE TO YOUR NAME.m or YOUR PSEUDONYM.m

% PREFIX TO START FILE NAMES
key.prefix = 'LPplot';
% CHANGE THE ABOVE STATEMENT TO REFLECT THE NEW NAME

% TEMPLATES TO CONVERT TO FAST RUNNING CODE 
key.steptemplate = 'LP_step_template.m';
key.meastemplate = 'LP_meas_template.m';

% PARAMETERS FOR GENERATING DATA
% MOST IMPORTANT PARAMETERS
key.gen.gNa = 120;
key.gen.gd  = 36;
key.gen.gl  = 0.3;
key.gen.ENa = 115;   % mV
key.gen.EK  = -12;   % mV
key.gen.El  = 10.6;  % mV
key.gen.cm  = 1.7;  % nF
% ALTERNATIVE VALUES FOR THE ABOVE PARAMETERS
%key.gen.gNa = 2300;
%key.gen.gd  = 3.5; % 0.35;
%key.gen.gl  = 0.1;
%key.gen.ENa = 50;   % mV
%key.gen.EK  = -80;  % mV
%key.gen.El  = -50;  % mV
%key.gen.cm = 1;     % uF/cm^2
% OTHER PARAMETERS
key.gen.goCa= 3.2;
key.gen.gA  = 2.2;
key.gen.gCa1= 0.21;
key.gen.gCa2= 0.047;
key.gen.gH  = 0.037;
key.gen.Eh  = -10;  % mV
key.gen.cn   = 180;   % /s
key.gen.koa  = 600;   % /s
key.gen.kob  = 35;    % /s
key.gen.kCa  = 360;   % /s
key.gen.kA1  = 140;   % /s
key.gen.kB  = 50;    % /s
key.gen.cA2  = 3.6;   % /s
key.gen.kaCa1= 50;    % /s
key.gen.kbCa1= 16;    % /s
key.gen.kaCa2= 10;    % /s
key.gen.cr   = 0.33;  % /s
key.gen.km   = 10000; % /s
key.gen.kh   = 500;   % /s
key.gen.hVn   = -25;  % mV
key.gen.hVkn  = 10;   % mV
key.gen.hVao1 = 0;    % mV
key.gen.hVao2 = -16;  % mV
key.gen.hVA1  = -12;  % mV
key.gen.hVB   = -62;  % mV
key.gen.hVA2  = -40;  % mV
key.gen.hVX   = 7;    % mV
key.gen.hVaCa1= -11;  % mV
key.gen.hVbCa1= -50;  % mV
key.gen.hVaCa2= 22;   % mV
key.gen.hVr   = -70;  % mV
key.gen.hVkr  = -110; % mV
key.gen.hVam  = -6;   % mV
key.gen.hVbm  = -34;  % mV
key.gen.hVah  = -39;  % mV
key.gen.hVbh  = -40;  % mV
key.gen.sn   = -17; % mV
key.gen.skn  = -22; % mV
key.gen.sao1 = -23; % mV
key.gen.sao2 = -5;  % mV
key.gen.sA1  = -26; % mV
key.gen.sB   = 6;   % mV
key.gen.sA2  = -12; % mV
key.gen.sX   = -15; % mV
key.gen.saCa1= -7;  % mV
key.gen.sbCa1= 8;   % mV
key.gen.saCa2= -7;  % mV
key.gen.sr   = 7;   % mV
key.gen.skr  = -13; % mV
key.gen.sam  = -20; % mV
key.gen.sbm  = -13; % mV
key.gen.sah  = -8;  % mV
key.gen.sbh  = -5;  % mV
key.gen.cf  = 0.6;  % mV/uM
key.gen.c1  = 2.5;  % uM
key.gen.c2  = 0.7;  % uM
key.gen.c3  = 0.6;  % uM
key.gen.CaO = 0.05; % uM
key.gen.cICa= 300;  % uM/nC
key.gen.cam = 0.11; % /mV
key.gen.cbm = 15; % /mV
key.gen.cah = 0.08; % /mV
key.gen.PNoiseLevel = 0.01;
key.gen.MNoiseLevel = 0.001;

% CHANGE CHOICE OF TWO PARAMETERS BELOW FOR THE PLOT AXES AND RANGES BELOW
% PLOT RANGES
key.plotrange.PARAMETER1 = LOW:STEP:HIGH; % no more than 30 steps please
key.plotrange.PARAMETER2 = LOW:STEP:HIGH; % no more than 30 steps please

% WHICH PARAMETERS TO FIT
key.fit = {'gNa','goCa'};

% ADD, REMOVE, AND OR CHANGE THE ASSUMPTONS BELOW
% ASSUMPTIONS ABOUT PARAMETERS & STATES, IF NOT SPECIFIED ASSUMED CORRECT
key.assume.gNa = 2300;
key.assume.gl = 0.1;

% INITIAL CONDITIONS FOR STATE VARIABLES
key.state.stateV = -65;
key.state.stateN = 0;
key.state.stateCa = .1;
key.state.stateAo= 0;
key.state.stateBo= 0;
key.state.stateAA= 0;
key.state.stateBA1= 0; 
key.state.stateBA2= 0;
key.state.stateACa1= 0; 
key.state.stateBCa1= 0;
key.state.stateACa2= 0;
key.state.stateR = 0;
key.state.stateM = 0;
key.state.stateH = 0;

% INPUT VARIABES AND CODE FOR CREATING INPUT
key.input = {'injectedCurrent'}; 
key.inputcreate = {@stepCurrent}; 
key.inputparams(1).holding = 0;
key.inputparams(1).step = 10;
key.inputparams(1).start = 100;
key.inputparams(1).stop = 200;

% NOISE VARIABLES 
key.noise.process = {'PNoise'};  
key.noise.measurement = {'MNoise'}; 

% MEASUREMENTS
key.meas = {'measuredVoltage'};

% TRAJECTORY TIME PARAMETERS
key.time.oneDT = 0.001;    % msec 
key.time.sqrtDT = sqrt(0.001); % msec^(1/2) 
key.time.TMAX = 300;       % msec
key.time.fitStart = 9;    % msec

% RANDOM SEED
key.seed = 0;

% DEBUGGING
key.debug = {'steady_n','relax_n','steady_ao','relax_ao','steady_bo',...
    'relax_bo','steady_aA','relax_aA','steady_bA1','relax_bA1', ...
    'steady_bA2','relax_bA2','weightX','steady_aCa1','relax_aCa1', ...
    'steady_bCa1','relax_bCa1','steady_aCa2','relax_aCa2','ECa', ...
    'steady_r','relax_r','alpha_m','beta_m','steady_m','relax_m', ...
    'alpha_h','beta_h','steady_h','relax_h','id','ioCa','iA','iCa', ...
    'ih','iNa','il','ionicCurrent'};
