Changeset 27547 in project
 Timestamp:
 10/05/12 03:56:12 (7 years ago)
 Location:
 release/4/nemo/trunk/examples/CGC
 Files:

 73 added
 4 deleted
 2 edited
Legend:
 Unmodified
 Added
 Removed

release/4/nemo/trunk/examples/CGC/DeSouzaCGC.nemo
r27532 r27547 147 147 ;; rate functions 148 148 149 (defun alpha_a (v )149 (defun alpha_a (v Q10) 150 150 (Q10 * Aalpha_a * sigm((v  V0alpha_a) Kalpha_a))) 151 151 152 (defun beta_a (v )152 (defun beta_a (v Q10) 153 153 (Q10 * (Abeta_a / exp((v  V0beta_a) / Kbeta_a)))) 154 154 155 (defun alpha_b (v )155 (defun alpha_b (v Q10) 156 156 (Q10 * Aalpha_b * sigm((v  V0alpha_b) Kalpha_b))) 157 157 158 (defun beta_b (v )158 (defun beta_b (v Q10) 159 159 (Q10 * Abeta_b * sigm((v  V0beta_b) Kbeta_b))) 160 160 … … 162 162 163 163 (a_inf = (1 / (1 + exp ((v  V0_ainf) / K_ainf)))) 164 (tau_a = (1 / (alpha_a (v ) + beta_a (v)) ))164 (tau_a = (1 / (alpha_a (v Q10) + beta_a (v Q10)) )) 165 165 (b_inf = (1 / (1 + exp ((v  V0_binf) / K_binf)))) 166 (tau_b = (1 / (alpha_b (v ) + beta_b (v)) ))166 (tau_b = (1 / (alpha_b (v Q10) + beta_b (v Q10)) )) 167 167 168 168 (hhionicgate … … 280 280 281 281 ;; rate constants 282 (Q10 = (pow (3 ((celsius  30) / 10))))282 (Q10 = (pow (3 ((celsius  20) / 10)))) 283 283 284 284 (const Aalpha_d = 0.13289) … … 349 349 (const Kbeta_n = 20) 350 350 (const V0beta_n = 30) 351 (const V0_ninf = 3 0)351 (const V0_ninf = 35) 352 352 (const B_ninf = 6) 353 353 … … 363 363 (hhionicgate 364 364 (KM ;; ion name: exported variables will be of the form {ion}_{id} 365 (initialm ( alpha_n (v Q10)) / (alpha_n (v Q10) + beta_n (v Q10)))365 (initialm (1 / (1 + exp((neg (v  V0_ninf)) / B_ninf)))) 366 366 (mpower 1) 367 367 (hpower 0) … … 509 509 (Ooff = (Q10 * AOoff)) 510 510 (a = (pow ((Oon / Con) (1.0 / 4.0)))) 511 (b = (pow ((Ooff /Coff) (1.0 / 4.0))))511 (b = (pow ((Ooff / Coff) (1.0 / 4.0)))) 512 512 513 513 ;; rate functions … … 516 516 517 517 (defun beta (v Q10) 518 (Q10 * Abeta * exp (((neg v) / Vbeta))))518 (Q10 * Abeta * exp (neg (v) / Vbeta))) 519 519 520 520 (defun teta (v Q10) 521 (Q10 * Ateta * exp (((neg v) / Vteta))))521 (Q10 * Ateta * exp (neg (v) / Vteta))) 522 522 523 523 … … 540 540 (<> C1 I1 Con Coff) 541 541 (<> C2 I2 (Con * a) (Coff * b)) 542 (<> C3 I3 (Con * pow (a , 2)) (Coff * pow (b,2)))543 (<> C4 I4 (Con * pow (a , 3)) (Coff * pow (b,3)))544 (<> C5 I5 (Con * pow (a , 4)) (Coff * pow (b,4)))545 546 (O <> I6 Oon Ooff) 542 (<> C3 I3 (Con * pow (a 2)) (Coff * pow (b 2))) 543 (<> C4 I4 (Con * pow (a 3)) (Coff * pow (b 3))) 544 (<> C5 I5 (Con * pow (a 4)) (Coff * pow (b 4))) 545 546 (O <> I6 Oon Ooff)) 547 547 548 548 (conserve (1 = (C1 + C2 + C3 + C4 + C5 + O + B + I1 + I2 + I3 + I4 + I5 + I6))) … … 564 564 565 565 ) ;; end Na current 566 567 566 568 567 (component (type voltageclamp) (name Na) 
release/4/nemo/trunk/examples/CGC/plot_original_neuron_vclamp.m
r27445 r27547 2 2 graphics_toolkit("gnuplot"); 3 3 4 CaHVA0 = load (" parts/original_vclamp/CaHVA.dat");5 CaHVA1 = load (" parts/neuron_vclamp/CaHVA.dat");4 CaHVA0 = load ("vclamp/original/CaHVA.dat"); 5 CaHVA1 = load ("vclamp/nemo_nmodl_generated/CaHVA.dat"); 6 6 7 KA0 = load (" parts/original_vclamp/KA.dat");8 KA1 = load (" parts/neuron_vclamp/KA.dat");7 KA0 = load ("vclamp/original/KA.dat"); 8 KA1 = load ("vclamp/nemo_nmodl_generated/KA.dat"); 9 9 10 K ir0 = load ("parts/original_vclamp/Kir.dat");11 K ir1 = load ("parts/neuron_vclamp/Kir.dat");10 KCa0 = load ("vclamp/original/KCa.dat"); 11 KCa1 = load ("vclamp/nemo_nmodl_generated/KCa.dat"); 12 12 13 K M0 = load ("parts/original_vclamp/KM.dat");14 K M1 = load ("parts/neuron_vclamp/KM.dat");13 Kir0 = load ("vclamp/original/Kir.dat"); 14 Kir1 = load ("vclamp/nemo_nmodl_generated/Kir.dat"); 15 15 16 K V0 = load ("parts/original_vclamp/KV.dat");17 K V1 = load ("parts/neuron_vclamp/KV.dat");16 KM0 = load ("vclamp/original/KM.dat"); 17 KM1 = load ("vclamp/nemo_nmodl_generated/KM.dat"); 18 18 19 Na0 = load ("parts/original_vclamp/Na.dat");20 Na1 = load ("parts/neuron_vclamp/Na.dat");19 KV0 = load ("vclamp/original/KV.dat"); 20 KV1 = load ("vclamp/nemo_nmodl_generated/KV.dat"); 21 21 22 NaR0 = load ("parts/original_vclamp/NaR.dat"); 23 NaR1 = load ("parts/neuron_vclamp/NaR.dat"); 24 25 pNa0 = load ("parts/original_vclamp/pNa.dat"); 26 pNa1 = load ("parts/neuron_vclamp/pNa.dat"); 22 Na0 = load ("vclamp/original/Na.dat"); 23 Na1 = load ("vclamp/nemo_nmodl_generated/Na.dat"); 27 24 28 25 … … 36 33 37 34 subplot(3,3,3); 35 plot(KCa0(:,1),KCa0(:,2),KCa1(:,1),KCa1(:,2),'linewidth',2); 36 title ("KCa current"); 37 38 subplot(3,3,4); 38 39 plot(Kir0(:,1),Kir0(:,2),Kir1(:,1),Kir1(:,2),'linewidth',2); 39 40 title ("Kir current"); 40 41 41 subplot(3,3, 4);42 subplot(3,3,5); 42 43 plot(KM0(:,1),KM0(:,2),KM1(:,1),KM1(:,2),'linewidth',2); 43 44 title ("KM current"); 44 45 45 subplot(3,3, 5);46 subplot(3,3,6); 46 47 plot(KV0(:,1),KV0(:,2),KV1(:,1),KV1(:,2),'linewidth',2); 47 48 title ("KV current"); 48 49 49 subplot(3,3, 6);50 subplot(3,3,7); 50 51 plot(Na0(:,1),Na0(:,2),Na1(:,1),Na1(:,2),'linewidth',2); 51 52 title ("Na current"); 52 53 53 subplot(3,3,7);54 plot(NaR0(:,1),NaR0(:,2),NaR1(:,1),NaR1(:,2),'linewidth',2);55 title ("NaR current");56 57 subplot(3,3,8);58 plot(pNa0(:,1),pNa0(:,2),pNa1(:,1),pNa1(:,2),'linewidth',2);59 title ("pNa current");60 61 54 print ("NEURON_Original_Vclamp.eps", "depsc"); 62 55
Note: See TracChangeset
for help on using the changeset viewer.