// genesis


include fullLTD77flash.g
include saveLTD.g

float normdt  = 0.0001
float finedt  = 5e-06 
float roughdt = 0.001
float plotdt  = 0.1
float Cabasal = 0.056902  /* default value */

function do_ramp (peakCa, duration)
	float peakCa
	float duration
	str filename

        /* Equilibrium at the basal state for 1 sec */
        /* initialze */
	setfield /kinetics/CaReg/Ca2+/proto/CaSpine_inc CoInit 0
	setfield /kinetics/CaReg/Ca2+/proto/CaSpine_inc Co     0
	setfield /kinetics/CaReg/ramp kf 0

        reset 
        setclock 0 {finedt}
        setclock 1 {finedt}
        setclock 3 {plotdt}
        MAXTIME = 1
        step {MAXTIME} -time   /* t = 0 -> 0.1 */

	/* Apply ramp Ca */
        setclock 0 {finedt}
        setclock 1 {finedt}
        setclock 3 {plotdt}
	setfield /kinetics/CaReg/ramp kf { {peakCa - Cabasal}/duration }
	MAXTIME = duration
	step {MAXTIME} -time

	/* Quit ramp Ca */
	setclock 0 {finedt}
        setclock 1 {finedt}
        setclock 3 {plotdt}
	setfield /kinetics/CaReg/Ca2+/proto/CaSpine_inc CoInit 0
	setfield /kinetics/CaReg/Ca2+/proto/CaSpine_inc Co     0
	setfield /kinetics/CaReg/ramp kf 0
	MAXTIME = 1
	step {MAXTIME} -time

	/* After Ca ramp */
	setclock 0 {normdt}
        setclock 1 {normdt}
        setclock 3 {plotdt}
	MAXTIME = {1000 - duration}
	step {MAXTIME} -time

	setclock 0 {roughdt}
        setclock 1 {roughdt}
        setclock 3 {plotdt}
	MAXTIME = 5000.1
	step {MAXTIME} -time

	saveplots f77 {"ramppeakCa" @ {peakCa} @ "for" @ {duration} @ "sec"}
        echo "saved as peakCa = " {peakCa} "and duration = " {duration}

end


do_ramp 0.5 0.5
do_ramp 1   0.5
do_ramp 1.5 0.5
do_ramp 2   0.5
do_ramp 2.5 0.5
do_ramp 3   0.5
do_ramp 3.5 0.5
do_ramp 4   0.5

do_ramp 0.5 1
do_ramp 1   1
do_ramp 1.5 1
do_ramp 2   1
do_ramp 2.5 1
do_ramp 3   1
do_ramp 3.5 1
do_ramp 4   1

do_ramp 0.5 15
do_ramp 1   15
do_ramp 1.5 15
do_ramp 2   15
do_ramp 2.5 15
do_ramp 3   15
do_ramp 3.5 15
do_ramp 4   15

do_ramp 0.5 30
do_ramp 1   30
do_ramp 1.5 30
do_ramp 2   30
do_ramp 2.5 30
do_ramp 3   30
do_ramp 3.5 30
do_ramp 4   30