Ignore:
Timestamp:
09/25/08 09:56:46 (13 years ago)
Author:
Ivan Raikov
Message:

Added some useful shortcuts for tau/inf to alpha/beta conversion.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • release/3/nemo/trunk/examples/AKP06/PotIhCa.scm

    r11995 r11996  
    116116                        ) ;; end Kv4 current
    117117             
     118
    118119             
    119120             (component (type ion-channel) (name Ih)
     
    141142                                     (h-power   0)
    142143                                     (m-inf     (Ih_inf v))
    143                                      (m-tau     (/ temp_adj (Ih_tau v)))
     144                                     (m-tau     (/ (Ih_tau v) temp_adj))
    144145                                     ))
    145146                                   
     
    147148
    148149                        (component (type pore)
    149                                    (const  gbar_Ih  = 0.2)
     150                                   (const  gbar_Ih  = 0.0002)
    150151                                   (output gbar_Ih ))
    151152                       
     
    156157                        ) ;; end Ih current
    157158
    158              (component (type ion-channel) (name CaP)
    159                        
    160                         (component (type gate)
    161                                    
    162                                    ;; rate functions
    163                                    (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-conductance
    174                                     (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 current
     159;;           (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
    192193
    193194             
     
    268269             
    269270             
    270              )) ;; end model
     271             )) ;; end model
Note: See TracChangeset for help on using the changeset viewer.