Generated Code

The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
C
C There are a total of 17 entries in the algebraic variable array.
C There are a total of 5 entries in each of the rate and state variable arrays.
C There are a total of 20 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C CONSTS(1) is C_m in component model_parameters (uF_per_cm2).
C STATES(1) is V_s in component soma_compartment (mV).
C CONSTS(2) is V_Na in component soma_compartment (mV).
C CONSTS(3) is V_K in component soma_compartment (mV).
C STATES(2) is V_D in component dendritic_compartment (mV).
C ALGBRC(14) is I_Na in component soma_compartment (uA_per_cm2).
C ALGBRC(16) is I_soma in component soma_compartment (uA_per_cm2).
C ALGBRC(1) is I_h in component soma_compartment (uA_per_cm2).
C ALGBRC(7) is I_K_DR in component soma_compartment (uA_per_cm2).
C CONSTS(4) is g_K_DR in component soma_compartment (mS_per_cm2).
C CONSTS(5) is g_Na in component soma_compartment (mS_per_cm2).
C CONSTS(6) is g_c in component model_parameters (mS_per_cm2).
C CONSTS(7) is g_h in component soma_compartment (mS_per_cm2).
C CONSTS(8) is p in component model_parameters (dimensionless).
C STATES(3) is n in component gating_variables (dimensionless).
C STATES(4) is h in component gating_variables (dimensionless).
C ALGBRC(11) is m_infinity in component gating_variables (dimensionless).
C CONSTS(9) is V_L in component dendritic_compartment (mV).
C ALGBRC(17) is I_D in component dendritic_compartment (uA_per_cm2).
C ALGBRC(12) is I_L in component dendritic_compartment (uA_per_cm2).
C ALGBRC(8) is I_pump in component dendritic_compartment (uA_per_cm2).
C CONSTS(20) is I_pump_ss in component dendritic_compartment (uA_per_cm2).
C ALGBRC(15) is I_NMDA in component dendritic_compartment (uA_per_cm2).
C ALGBRC(13) is I_Na_NMDA in component dendritic_compartment (uA_per_cm2).
C CONSTS(15) is R_pump in component dendritic_compartment (uA_per_cm2).
C ALGBRC(2) is f_NMDA in component dendritic_compartment (dimensionless).
C CONSTS(10) is alpha in component dendritic_compartment (mMcm2_per_uAs).
C CONSTS(16) is g_NMDA in component dendritic_compartment (mS_per_cm2).
C CONSTS(17) is g_Na_NMDA in component dendritic_compartment (mS_per_cm2).
C CONSTS(18) is g_L in component dendritic_compartment (mS_per_cm2).
C STATES(5) is Na in component dendritic_compartment (mM).
C CONSTS(11) is Na_eq in component dendritic_compartment (mM).
C CONSTS(12) is K_p in component dendritic_compartment (mM).
C CONSTS(13) is K_Na in component dendritic_compartment (mM).
C CONSTS(14) is q in component dendritic_compartment (mV).
C ALGBRC(3) is n_infinity in component gating_variables (dimensionless).
C ALGBRC(4) is h_infinity in component gating_variables (dimensionless).
C ALGBRC(5) is r_infinity in component gating_variables (dimensionless).
C ALGBRC(9) is tau_h in component gating_variables (second).
C ALGBRC(10) is tau_n in component gating_variables (second).
C ALGBRC(6) is tau_mL in component gating_variables (second).
C CONSTS(19) is tau_r in component gating_variables (second).
C RATES(1) is d/dt V_s in component soma_compartment (mV).
C RATES(2) is d/dt V_D in component dendritic_compartment (mV).
C RATES(5) is d/dt Na in component dendritic_compartment (mM).
C RATES(4) is d/dt h in component gating_variables (dimensionless).
C RATES(3) is d/dt n in component gating_variables (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 1
      STATES(1) = -64
      CONSTS(2) = 55
      CONSTS(3) = -85
      STATES(2) = -77
      CONSTS(4) = 3.2
      CONSTS(5) = 3.2
      CONSTS(6) = 0.1
      CONSTS(7) = 0.1
      CONSTS(8) = 0.5
      STATES(3) = 0.002
      STATES(4) = 1
      CONSTS(9) = -50
      CONSTS(10) = 0.173
      STATES(5) = 5.09
      CONSTS(11) = 8
      CONSTS(12) = 15
      CONSTS(13) = 15
      CONSTS(14) = 12.5
      CONSTS(15) =  18.0000*(CONSTS(8)/(1.00000 - CONSTS(8)))
      CONSTS(16) =  1.25000*(CONSTS(8)/(1.00000 - CONSTS(8)))
      CONSTS(17) =  1.00000*(CONSTS(8)/(1.00000 - CONSTS(8)))
      CONSTS(18) =  0.180000*(CONSTS(8)/(1.00000 - CONSTS(8)))
      CONSTS(19) = 190.000
      CONSTS(20) = ( (( CONSTS(15)*CONSTS(8))/(1.00000 - CONSTS(8)))*CONSTS(11) ** 3.00000)/(CONSTS(13) ** 3.00000+CONSTS(11) ** 3.00000)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+30.0000)/8.30000))
      ALGBRC(9) = 0.430000+0.860000/(1.00000+EXP((STATES(1)+25.0000)/5.00000))
      RATES(4) = (ALGBRC(4) - STATES(4))/ALGBRC(9)
      ALGBRC(3) = 1.00000/(1.00000+EXP(- (STATES(1)+31.0000)/5.30000))
      ALGBRC(10) = (0.800000+1.60000/(1.00000+EXP( 0.100000*(STATES(1)+25.0000))))/(1.00000+EXP( - 0.100000*(STATES(1)+70.0000)))
      RATES(3) = (ALGBRC(3) - STATES(3))/ALGBRC(10)
      ALGBRC(8) = ( (( CONSTS(15)*CONSTS(8))/(1.00000 - CONSTS(8)))* STATES(2)*1.00000+STATES(5) ** 3.00000)/(CONSTS(12) ** 3.00000+ STATES(2)*1.00000+STATES(5) ** 3.00000)
      ALGBRC(2) = 1.00000/(1.00000+ 0.141000*EXP(- STATES(2)/CONSTS(14)))
      ALGBRC(13) =  CONSTS(17)*ALGBRC(2)*(STATES(2) - CONSTS(2))
      RATES(5) =  CONSTS(10)*(- ALGBRC(13) -  3.00000*(ALGBRC(8) - CONSTS(20)))
      ALGBRC(11) = 1.00000/(1.00000+EXP(- (STATES(1)+35.0000)/6.20000))
      ALGBRC(14) =  CONSTS(5)*STATES(4)*(STATES(1) - CONSTS(2))*ALGBRC(11) ** 3.00000
      ALGBRC(1) =  CONSTS(7)*STATES(4)*(STATES(1)+30.0000)
      ALGBRC(7) =  CONSTS(4)*(STATES(1) - CONSTS(3))*STATES(3) ** 2.00000
      ALGBRC(16) = ALGBRC(14)+ALGBRC(7)+ALGBRC(1)
      RATES(1) = (- ALGBRC(16) -  (CONSTS(6)/CONSTS(8))*(STATES(2) - STATES(1)))/CONSTS(1)
      ALGBRC(12) =  CONSTS(18)*(STATES(2) - CONSTS(9))
      ALGBRC(15) =  CONSTS(16)*ALGBRC(2)*STATES(2)
      ALGBRC(17) = ((ALGBRC(15)+ALGBRC(8)) - CONSTS(20))+ALGBRC(12)
      RATES(2) = (- ALGBRC(17)+ (CONSTS(6)/(1.00000 - CONSTS(8)))*(STATES(1) - STATES(2)))/CONSTS(1)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+30.0000)/8.30000))
      ALGBRC(9) = 0.430000+0.860000/(1.00000+EXP((STATES(1)+25.0000)/5.00000))
      ALGBRC(3) = 1.00000/(1.00000+EXP(- (STATES(1)+31.0000)/5.30000))
      ALGBRC(10) = (0.800000+1.60000/(1.00000+EXP( 0.100000*(STATES(1)+25.0000))))/(1.00000+EXP( - 0.100000*(STATES(1)+70.0000)))
      ALGBRC(8) = ( (( CONSTS(15)*CONSTS(8))/(1.00000 - CONSTS(8)))* STATES(2)*1.00000+STATES(5) ** 3.00000)/(CONSTS(12) ** 3.00000+ STATES(2)*1.00000+STATES(5) ** 3.00000)
      ALGBRC(2) = 1.00000/(1.00000+ 0.141000*EXP(- STATES(2)/CONSTS(14)))
      ALGBRC(13) =  CONSTS(17)*ALGBRC(2)*(STATES(2) - CONSTS(2))
      ALGBRC(11) = 1.00000/(1.00000+EXP(- (STATES(1)+35.0000)/6.20000))
      ALGBRC(14) =  CONSTS(5)*STATES(4)*(STATES(1) - CONSTS(2))*ALGBRC(11) ** 3.00000
      ALGBRC(1) =  CONSTS(7)*STATES(4)*(STATES(1)+30.0000)
      ALGBRC(7) =  CONSTS(4)*(STATES(1) - CONSTS(3))*STATES(3) ** 2.00000
      ALGBRC(16) = ALGBRC(14)+ALGBRC(7)+ALGBRC(1)
      ALGBRC(12) =  CONSTS(18)*(STATES(2) - CONSTS(9))
      ALGBRC(15) =  CONSTS(16)*ALGBRC(2)*STATES(2)
      ALGBRC(17) = ((ALGBRC(15)+ALGBRC(8)) - CONSTS(20))+ALGBRC(12)
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+80.0000)/8.00000))
      ALGBRC(6) = 0.400000/( 5.00000*EXP(- (STATES(2)+11.0000)/8.30000)+(- (STATES(2)+11.0000)/8.30000)/(EXP(- (STATES(2)+11.0000)/8.30000) - 1.00000))
      RETURN
      END