Changeset 27023 in project
 Timestamp:
 07/08/12 10:03:03 (9 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

release/4/nemo/trunk/examples/CGC.nemo
r27006 r27023 12 12 (ica from ioncurrents)) 13 13 14 (const Vrest = 68) 15 16 (const fix_celsius = 30) 17 18 19 20 (defun sigmoid (x y) (1 / (exp (x / y) + 1))) 21 14 15 (defun sigm (x y) (1 / (exp (x / y) + 1))) 16 22 17 23 18 (defun linoid (x y) … … 26 21 else (x / (exp (x / y)  1)) 27 22 )) 23 24 (const fix_celsius = 30) 28 25 29 26 (component (type decayingpool) (name ca) … … 40 37 (output ca) 41 38 ) 42 43 39 44 40 (component (type ioniccurrent) (name CaHVA ) … … 75 71 76 72 (defun alpha_u (v) 77 (Q10 * Aalpha_ s * exp((v  V0alpha_s) / Kalpha_s)))73 (Q10 * Aalpha_u * exp((v  V0alpha_u) / Kalpha_u))) 78 74 79 75 (defun beta_u (v) … … 96 92 ) 97 93 98 99 94 (component (type pore) 100 95 (const gbar = 0.00046) … … 104 99 (const e = 129.33) 105 100 (output e )) 101 102 103 106 104 107 105 ) ;; end CaHVA current 106 107 (component (type voltageclamp) (name CaHVA) 108 109 (const vchold = 71) 110 (const vcbase = 69) 111 (const vcinc = 10) 112 (const vcsteps = 8) 113 (const vchdur = 30) 114 (const vcbdur = 100) 115 116 (output vchold vcbase vcinc vcsteps vchdur vcbdur) 117 ) 108 118 109 119 … … 140 150 141 151 (defun alpha_a (v) 142 (Q10 * Aalpha_a * sigm oid((v  V0alpha_a) Kalpha_a)))152 (Q10 * Aalpha_a * sigm((v  V0alpha_a) Kalpha_a))) 143 153 144 154 (defun beta_a (v) 145 (Q10 * Abeta_a * exp((v  V0beta_a) / Kbeta_a)))155 (Q10 * (Abeta_a / exp((v  V0beta_a) / Kbeta_a)))) 146 156 147 157 (defun alpha_b (v) 148 (Q10 * Aalpha_b * sigm oid((v  V0alpha_b) Kalpha_b)))158 (Q10 * Aalpha_b * sigm((v  V0alpha_b) Kalpha_b))) 149 159 150 160 (defun beta_b (v) 151 (Q10 * Abeta_b * sigmoid((v  V0beta_b) Kbeta_b))) 161 (Q10 * Abeta_b * sigm((v  V0beta_b) Kbeta_b))) 162 163 152 164 153 165 (a_inf = (1 / (1 + exp ((v  V0_ainf) / K_ainf)))) … … 181 193 182 194 ) ;; end KA current 183 195 196 197 (component (type voltageclamp) (name KA) 198 199 (const vchold = 71) 200 (const vcbase = 69) 201 (const vcinc = 10) 202 (const vcsteps = 8) 203 (const vchdur = 30) 204 (const vcbdur = 100) 205 206 (output vchold vcbase vcinc vcsteps vchdur vcbdur) 207 ) 208 184 209 185 210 (component (type ioniccurrent) (name KCa ) … … 259 284 (mpower 1) 260 285 (hpower 0) 261 (m inf(alpha_d (v)))262 (m tau(beta_d (v)))286 (malpha (alpha_d (v))) 287 (mbeta (beta_d (v))) 263 288 )) 264 289 … … 276 301 277 302 ) ;; end Kir current 278 279 280 303 304 305 306 (component (type voltageclamp) (name Kir) 307 308 (const vchold = 71) 309 (const vcbase = 69) 310 (const vcinc = 10) 311 (const vcsteps = 8) 312 (const vchdur = 30) 313 (const vcbdur = 100) 314 315 (output vchold vcbase vcinc vcsteps vchdur vcbdur) 316 ) 281 317 282 318 … … 311 347 (mpower 1) 312 348 (hpower 0) 349 (mtau (1 / (alpha_n(v) + beta_n (v)) )) 313 350 (minf (1 / (1 + exp((neg (v  V0_ninf)) / B_ninf)))) 314 (mtau (1 / (alpha_n(v) + beta_n (v)) ))315 351 )) 316 352 ) … … 327 363 328 364 365 (component (type voltageclamp) (name KM) 366 367 (const vchold = 71) 368 (const vcbase = 69) 369 (const vcinc = 10) 370 (const vcsteps = 8) 371 (const vchdur = 30) 372 (const vcbdur = 100) 373 374 (output vchold vcbase vcinc vcsteps vchdur vcbdur) 375 ) 329 376 330 377 (component (type ioniccurrent) (name KV ) … … 348 395 349 396 (defun beta_n (v) 350 (Q10 * Abeta_n * linoid ((v  V0beta_n) Kbeta_n)))397 (Q10 * Abeta_n * exp((v  V0beta_n) / Kbeta_n) )) 351 398 352 399 (hhionicgate … … 369 416 370 417 ) ;; end KV current 371 372 373 (component (type ioniccurrent) (name pNa ) 374 375 (component (type gate) 376 377 ;; rate constants 378 (const Q10 = (pow (3 ((fix_celsius  30) / 10)))) 379 380 (const Aalpha_m = 0.091) 381 (const Kalpha_m = 5) 382 (const V0alpha_m = 42) 383 (const Abeta_m = 0.062) 384 (const Kbeta_m = 5) 385 (const V0beta_m = 42) 386 (const V0_minf = 42) 387 (const B_minf = 5) 388 389 ;; rate functions 390 (defun alpha_m (v) 391 (Q10 * Aalpha_m * linoid( (v  V0alpha_m) Kalpha_m))) 392 393 (defun beta_m (v) 394 (Q10 * Abeta_m * linoid ( (v  V0beta_m) Kbeta_m) )) 395 396 397 (hhionicgate 398 (pNa ;; ion name: exported variables will be of the form {ion}_{id} 399 (initialm ((alpha_d (v)) / (alpha_d (v) + beta_d (v)) )) 400 (mpower 1) 401 (hpower 0) 402 (minf (1 / (1 + exp((neg (v  V0_minf)) / B_minf)))) 403 (mtau (5 / (alpha_m (v) + beta_m (v)))) 404 )) 405 406 ) 407 408 409 (component (type pore) 410 (const gbar = 2e5) 411 (output gbar )) 412 413 414 (component (type permeatingion) (name na) 415 (const e = 87.39) 416 (output e )) 417 418 ) ;; end pNa current 419 418 419 420 (component (type voltageclamp) (name KV) 421 422 (const vchold = 71) 423 (const vcbase = 69) 424 (const vcinc = 10) 425 (const vcsteps = 8) 426 (const vchdur = 30) 427 (const vcbdur = 100) 428 429 (output vchold vcbase vcinc vcsteps vchdur vcbdur) 430 ) 431 420 432 433 (component (type ioniccurrent) (name Lkg1) 434 435 (component (type pore) 436 (const gbar = (5.68e5)) 437 (output gbar)) 438 439 (component (type permeatingion) (name nonspecific) 440 (const e = 58) 441 (output e )) 442 443 ) ;; end leak current 444 445 446 (component (type ioniccurrent) (name Lkg2) 447 448 (component (type pore) 449 (const ggaba = (2.17e5)) 450 (output ggaba)) 451 452 (component (type permeatingion) (name nonspecific) 453 (const egaba = 65) 454 (output egaba )) 455 456 ) ;; end leak current 457 421 458 422 459 (component (type ioniccurrent) (name Na ) … … 444 481 445 482 ;; rate functions 446 (defun alpha_m (v) 483 (defun alpha_m (v) 447 484 (Q10 * Aalpha_m * linoid((v  V0alpha_m) Kalpha_m) )) 448 485 … … 451 488 452 489 (defun alpha_h (v) 453 (Q10 * Aalpha_h * linoid((v  V0alpha_h) Kalpha_h) ))454 490 (Q10 * Aalpha_h * exp((v  V0alpha_h) / Kalpha_h) )) 491 455 492 (defun beta_h (v) 456 (Q10 * Abeta_h * exp((v  V0beta_h) / Kbeta_h)))493 (Q10 * Abeta_h / (1 + exp((v  V0beta_h) / Kbeta_h) ))) 457 494 458 495 … … 483 520 484 521 522 (component (type voltageclamp) (name Na) 523 524 (const vchold = 71) 525 (const vcbase = 60) 526 (const vcinc = 10) 527 (const vcsteps = 9) 528 (const vchdur = 30) 529 (const vcbdur = 100) 530 531 (output vchold vcbase vcinc vcsteps vchdur vcbdur)) 532 533 485 534 (component (type ioniccurrent) (name Nar ) 486 535 … … 543 592 544 593 ) ;; end Nar current 545 546 547 (component (type ioniccurrent) (name Lkg1) 548 549 (component (type pore) 550 (const gbar = (5.68e5)) 551 (output gbar)) 552 553 (component (type permeatingion) (name nonspecific) 554 (const e = 58) 555 (output e )) 556 557 ) ;; end leak current 558 559 (component (type ioniccurrent) (name Lkg2) 560 561 (component (type pore) 562 (const ggaba = (2.17e5)) 563 (output ggaba)) 564 565 (component (type permeatingion) (name nonspecific) 566 (const egaba = 65) 567 (output egaba )) 568 569 ) ;; end leak current 570 571 594 595 596 (component (type voltageclamp) (name Nar) 597 598 (const vchold = 71) 599 (const vcbase = 60) 600 (const vcinc = 10) 601 (const vcsteps = 9) 602 (const vchdur = 30) 603 (const vcbdur = 100) 604 605 (output vchold vcbase vcinc vcsteps vchdur vcbdur)) 606 607 (component (type ioniccurrent) (name pNa ) 608 609 (component (type gate) 610 611 ;; rate constants 612 (const Q10 = (pow (3 ((fix_celsius  30) / 10)))) 613 614 (const Aalpha_m = 0.091) 615 (const Kalpha_m = 5) 616 (const V0alpha_m = 42) 617 (const Abeta_m = 0.062) 618 (const Kbeta_m = 5) 619 (const V0beta_m = 42) 620 (const V0_minf = 42) 621 (const B_minf = 5) 622 623 ;; rate functions 624 (defun alpha_m (v) 625 (Q10 * Aalpha_m * linoid( (v  V0alpha_m) Kalpha_m))) 626 627 (defun beta_m (v) 628 (Q10 * Abeta_m * linoid ( (v  V0beta_m) Kbeta_m) )) 629 630 631 (hhionicgate 632 (pNa ;; ion name: exported variables will be of the form {ion}_{id} 633 (initialm ((alpha_m (v)) / (alpha_m (v) + beta_m (v)) )) 634 (mpower 1) 635 (hpower 0) 636 (minf (1 / (1 + exp((neg (v  V0_minf)) / B_minf)))) 637 (mtau (5 / (alpha_m (v) + beta_m (v)))) 638 )) 639 640 ) 641 642 643 (component (type pore) 644 (const gbar = 2e5) 645 (output gbar )) 646 647 648 (component (type permeatingion) (name na) 649 (const e = 87.39) 650 (output e )) 651 652 ) ;; end pNa current 653 654 (component (type voltageclamp) (name pNa) 655 656 (const vchold = 71) 657 (const vcbase = 60) 658 (const vcinc = 10) 659 (const vcsteps = 9) 660 (const vchdur = 30) 661 (const vcbdur = 100) 662 663 (output vchold vcbase vcinc vcsteps vchdur vcbdur)) 664 572 665 ))
Note: See TracChangeset
for help on using the changeset viewer.