Changeset 11996 in project
- Timestamp:
- 09/25/08 09:56:46 (12 years ago)
- Location:
- release/3/nemo/trunk
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
release/3/nemo/trunk/examples/AKP06/PotIhCa.scm
r11995 r11996 116 116 ) ;; end Kv4 current 117 117 118 118 119 119 120 (component (type ion-channel) (name Ih) … … 141 142 (h-power 0) 142 143 (m-inf (Ih_inf v)) 143 (m-tau (/ temp_adj (Ih_tau v)))144 (m-tau (/ (Ih_tau v) temp_adj)) 144 145 )) 145 146 … … 147 148 148 149 (component (type pore) 149 (const gbar_Ih = 0. 2)150 (const gbar_Ih = 0.0002) 150 151 (output gbar_Ih )) 151 152 … … 156 157 ) ;; end Ih current 157 158 158 (component (type ion-channel) (name CaP)159 160 (component (type gate)161 162 ;; rate functions163 (defun CaP_inf (v)164 (let ((cv 19)165 (ck 5.5))166 (/ 1.0 (+ 1.0 (exp (neg (/ (+ v cv) ck)))))))167 168 (defun CaP_tau (v temp_adj)169 (if (> v -50)170 (* 1e3 (+ 0.000191 (* 0.00376 (pow (exp (neg (/ (+ v 41.9) 27.8))) 2))))171 (* 1e3 (+ 0.00026367 (* 0.1278 (exp (* 0.10327 v)))))))172 173 (hh-ionic-conductance174 (CaP ;; ion name: exported variables will be of the form {ion}_{id}175 (initial-m (CaP_inf Vrest))176 (m-power 1)177 (h-power 0)178 (m-inf (CaP_inf v))179 (m-tau (/ (CaP_tau v) temp_adj))))180 181 )182 183 (component (type pore)184 (const gmax_CaP = 0.01667)185 (gbar_CaP = (* gmax_CaP (ghk v celsius cai cao)))186 (output gbar_CaP ))187 188 (component (type accumulating-substance) (name ca) )189 190 191 ) ;; end CaP current159 ;; (component (type ion-channel) (name CaP) 160 161 ;; (component (type gate) 162 163 ;; ;; rate functions 164 ;; (defun CaP_inf (v) 165 ;; (let ((cv 19) 166 ;; (ck 5.5)) 167 ;; (/ 1.0 (+ 1.0 (exp (neg (/ (+ v cv) ck))))))) 168 169 ;; (defun CaP_tau (v) 170 ;; (if (> v -50) 171 ;; (* 1e3 (+ 0.000191 (* 0.00376 (pow (exp (neg (/ (+ v 41.9) 27.8))) 2)))) 172 ;; (* 1e3 (+ 0.00026367 (* 0.1278 (exp (* 0.10327 v))))))) 173 174 ;; (hh-ionic-conductance 175 ;; (CaP ;; ion name: exported variables will be of the form {ion}_{id} 176 ;; (initial-m (CaP_inf Vrest)) 177 ;; (m-power 1) 178 ;; (h-power 0) 179 ;; (m-inf (CaP_inf v)) 180 ;; (m-tau (/ (CaP_tau v) temp_adj)))) 181 182 ;; ) 183 184 ;; (component (type pore) 185 ;; (const gmax_CaP = 0.01667) 186 ;; (gbar_CaP = (* gmax_CaP (ghk v celsius cai cao))) 187 ;; (output gbar_CaP )) 188 189 ;; (component (type accumulating-substance) (name ca) ) 190 191 192 ;; ) ;; end CaP current 192 193 193 194 … … 268 269 269 270 270 )) ;; end model271 )) ;; end model -
release/3/nemo/trunk/examples/AKP06/morphology_mechanisms.hoc
r11848 r11996 16 16 //insert CaBK 17 17 insert Caint 18 //insert CaP18 insert CaP 19 19 //insert Ih 20 20 insert leak 21 insert PotIhCa 21 22 22 insert PotIhCa23 23 soma.ena = 60 // (mV) 24 24 soma.ek = -88 // (mV) 25 //soma.eh_Ih = -30 26 soma.e_Ih_PotIhCa = -30 25 27 soma.e_leak = -61 26 28 soma.cao = 2 // (mM) 27 soma.e_Ih_PotIhCa = -30 28 //soma.e_Ih = -30 29 29 30 30 soma.gmax_Narsg = 0.016 // (mho/cm2) 31 31 soma.gmax_Na = 0.014 32 //soma.gbar_Kv1 = 0.011 33 soma.gbar_Kv1_PotIhCa = 0.011 34 //soma.gbar_Kv4 = 0.0039 35 soma.gbar_Kv4_PotIhCa = 0.0039 32 36 soma.gmax_Kbin = 0.0016 37 //soma.gkbar_CaBK = 0.014 38 soma.gbar_CaBK_PotIhCa = 0.014 39 soma.pcabar_CaP = 0.00006 40 //soma.gmax_CaP_PotIhCa = 0.06 41 //soma.ghbar_Ih = 0.0002 42 soma.gbar_Ih_PotIhCa = 0.0002 33 43 soma.gbar_leak = 9e-5 34 35 soma.gbar_Kv1_PotIhCa = 0.01136 soma.gbar_Kv4_PotIhCa = 0.003937 //soma.gmax_Kv4 = 0.003938 soma.gbar_Ih_PotIhCa = 0.000239 //soma.gmax_Ih = 0.000240 soma.gmax_CaP_PotIhCa = 0.0641 //soma.gmax_CaP = 0.0642 43 soma.gbar_CaBK_PotIhCa = 0.01444 //soma.gmax_CaBK = 0.014 -
release/3/nemo/trunk/extensions/nemo-hh.scm
r11995 r11996 95 95 (m-inf (lookup-field 'm-inf alst)) 96 96 (m-tau (lookup-field 'm-tau alst)) 97 (m-alpha (or (lookup-field 'm-alpha alst) `(/ ,m-inf ,m-tau))) 98 (m-beta (or (lookup-field 'm-beta alst) `(/ (- 1 ,m-inf) ,m-tau))) 97 (m-inf-sym (p$ ion 'm-inf)) 98 (m-tau-sym (p$ ion 'm-tau)) 99 (m-alpha (or (lookup-field 'm-alpha alst) `(/ ,m-inf-sym ,m-tau-sym))) 100 (m-beta (or (lookup-field 'm-beta alst) `(/ (- 1 ,m-inf-sym) ,m-tau-sym))) 99 101 (open (p$ ion 'mO)) 100 102 (closed (p$ ion 'mC)) … … 102 104 (transitions (-> ,closed ,open ,m-alpha) 103 105 (-> ,open ,closed ,m-beta))))) 106 (if m-inf (env-extend! m-inf-sym '(asgn) 'none `(rhs ,m-inf))) 107 (if m-tau (env-extend! m-tau-sym '(asgn) 'none `(rhs ,m-tau))) 104 108 (apply env-extend! (cons* (p$ ion 'm) '(tscomp) initial-m mst)) 105 109 (add-external! (p$ ion 'm) 'output) -
release/3/nemo/trunk/nemo.scm
r11995 r11996 168 168 169 169 (define (model->nmodl options model) 170 (nemo:nmodl-translator model (lookup-def ' nmodl-method options) (lookup-def 'table options) -150 150 1))170 (nemo:nmodl-translator model (lookup-def 'method options) (lookup-def 'table options) -150 150 1)) 171 171 172 172 -
release/3/nemo/trunk/nmodl.scm
r11995 r11996 178 178 (define (bind ks vs e) `(let ,(zip ks vs) ,e)) 179 179 180 (define (name-normalize expr) 181 (match expr 182 (('if c t e) `(if ,(name-normalize c) ,(name-normalize t) ,(name-normalize e))) 183 (('let bs e) 184 `(let ,(map (lambda (b) `(,(car b) ,(name-normalize (cadr b)))) bs) ,(name-normalize e))) 185 ((f . es) 186 (cons f (map name-normalize es))) 187 ((? symbol? ) (nmodl-name expr)) 188 ((? atom? ) expr))) 180 189 181 190 (define (if-convert expr) … … 230 239 (let* ((expr1 (if-convert expr)) 231 240 (expr2 (subst-convert expr1 subst-empty)) 232 (expr3 (let-lift expr2))) 233 expr3))) 241 (expr3 (let-lift expr2)) 242 (expr4 (name-normalize expr3))) 243 expr4))) 234 244 235 245 (define (format-expr/NMODL indent expr . rest) … … 607 617 (map second acc-ions) (map fourth acc-ions))) 608 618 (asgns1 (delete-duplicates asgns0))) 609 (for-each (lambda (x) (pp indent+ , x)) asgns1)619 (for-each (lambda (x) (pp indent+ ,(nmodl-name x))) asgns1) 610 620 (pp indent "}") 611 621 … … 616 626 (if (not (null? locals)) (pp indent+ (LOCAL ,(sl\ ", " locals))))) 617 627 (for-each (lambda (def) 618 (let ((n ( first def))628 (let ((n (nmodl-name (first def)) ) 619 629 (b (second def))) 620 630 (pp indent+ ,(expr->string/NMODL b n)))) asgn-eq-defs)
Note: See TracChangeset
for help on using the changeset viewer.