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 76 entries in the algebraic variable array.
C There are a total of 29 entries in each of the rate and state variable arrays.
C There are a total of 86 entries in the constant variable array.
C
C
C VOI is time in component environment (millisecond).
C STATES(1) is Vm in component Vm (millivolt).
C CONSTS(1) is Cm in component model_parameters (picoF).
C ALGBRC(32) is i_CaT in component i_CaT (picoA).
C ALGBRC(15) is i_CaL in component i_CaL (picoA).
C ALGBRC(50) is i_f in component i_f (picoA).
C ALGBRC(51) is i_st in component i_st (picoA).
C ALGBRC(38) is i_Kr in component i_Kr (picoA).
C ALGBRC(43) is i_Ks in component i_Ks (picoA).
C ALGBRC(45) is i_to in component AP_sensitive_currents (picoA).
C ALGBRC(47) is i_sus in component AP_sensitive_currents (picoA).
C ALGBRC(53) is i_NaK in component i_NaK (picoA).
C ALGBRC(67) is i_NaCa in component i_NaCa (picoA).
C ALGBRC(54) is i_b_Ca in component i_b_Ca (picoA).
C ALGBRC(52) is i_b_Na in component i_b_Na (picoA).
C CONSTS(76) is E_Na in component electric_potentials (millivolt).
C CONSTS(77) is E_K in component electric_potentials (millivolt).
C CONSTS(78) is E_Ks in component electric_potentials (millivolt).
C CONSTS(2) is R in component model_parameters (joule_per_kilomole_kelvin).
C CONSTS(3) is T in component model_parameters (kelvin).
C CONSTS(4) is F in component model_parameters (coulomb_per_mole).
C CONSTS(5) is Nai in component model_parameters (millimolar).
C CONSTS(6) is Nao in component model_parameters (millimolar).
C CONSTS(7) is Ki in component model_parameters (millimolar).
C CONSTS(8) is Ko in component model_parameters (millimolar).
C CONSTS(9) is E_CaL in component i_CaL (millivolt).
C CONSTS(10) is g_CaL in component i_CaL (nanoS_per_picoF).
C STATES(2) is Ca_sub in component calcium_dynamics (millimolar).
C STATES(3) is dL in component i_CaL_dL_gate (dimensionless).
C STATES(4) is fL in component i_CaL_fL_gate (dimensionless).
C STATES(5) is fCa in component i_CaL_fCa_gate (dimensionless).
C ALGBRC(1) is dL_infinity in component i_CaL_dL_gate (dimensionless).
C ALGBRC(46) is tau_dL in component i_CaL_dL_gate (millisecond).
C ALGBRC(33) is alpha_dL in component i_CaL_dL_gate (per_millisecond).
C ALGBRC(44) is beta_dL in component i_CaL_dL_gate (per_millisecond).
C ALGBRC(17) is adVm in component i_CaL_dL_gate (millivolt).
C ALGBRC(39) is bdVm in component i_CaL_dL_gate (millivolt).
C ALGBRC(2) is fL_infinity in component i_CaL_fL_gate (dimensionless).
C ALGBRC(18) is tau_fL in component i_CaL_fL_gate (millisecond).
C CONSTS(11) is alpha_fCa in component i_CaL_fCa_gate (per_millisecond).
C ALGBRC(3) is fCa_infinity in component i_CaL_fCa_gate (dimensionless).
C ALGBRC(19) is tau_fCa in component i_CaL_fCa_gate (millisecond).
C CONSTS(12) is Km_fCa in component i_CaL_fCa_gate (millimolar).
C CONSTS(13) is g_CaT in component i_CaT (nanoS_per_picoF).
C CONSTS(14) is E_CaT in component i_CaT (millivolt).
C STATES(6) is dT in component i_CaT_dT_gate (dimensionless).
C STATES(7) is fT in component i_CaT_fT_gate (dimensionless).
C ALGBRC(4) is dT_infinity in component i_CaT_dT_gate (dimensionless).
C ALGBRC(20) is tau_dT in component i_CaT_dT_gate (millisecond).
C ALGBRC(5) is fT_infinity in component i_CaT_fT_gate (dimensionless).
C ALGBRC(21) is tau_fT in component i_CaT_fT_gate (millisecond).
C CONSTS(15) is g_Kr in component i_Kr (nanoS_per_picoF).
C STATES(8) is paS in component i_Kr_pa_gate (dimensionless).
C STATES(9) is paF in component i_Kr_pa_gate (dimensionless).
C STATES(10) is pi_ in component i_Kr_pi_gate (dimensionless).
C ALGBRC(6) is pa_infinity in component i_Kr_pa_gate (dimensionless).
C ALGBRC(22) is tau_paS in component i_Kr_pa_gate (millisecond).
C ALGBRC(23) is tau_paF in component i_Kr_pa_gate (millisecond).
C ALGBRC(7) is pi_infinity in component i_Kr_pi_gate (dimensionless).
C ALGBRC(24) is tau_pi in component i_Kr_pi_gate (millisecond).
C CONSTS(16) is g_Ks in component i_Ks (nanoS_per_picoF).
C STATES(11) is n in component i_Ks_n_gate (dimensionless).
C ALGBRC(34) is n_infinity in component i_Ks_n_gate (dimensionless).
C ALGBRC(40) is tau_n in component i_Ks_n_gate (millisecond).
C ALGBRC(8) is alpha_n in component i_Ks_n_gate (per_millisecond).
C ALGBRC(25) is beta_n in component i_Ks_n_gate (per_millisecond).
C CONSTS(17) is g_to in component AP_sensitive_currents (nanoS_per_picoF).
C CONSTS(18) is g_sus in component AP_sensitive_currents (nanoS_per_picoF).
C STATES(12) is q in component AP_sensitive_currents_q_gate (dimensionless).
C STATES(13) is r in component AP_sensitive_currents_r_gate (dimensionless).
C ALGBRC(9) is q_infinity in component AP_sensitive_currents_q_gate (dimensionless).
C ALGBRC(26) is tau_q in component AP_sensitive_currents_q_gate (millisecond).
C ALGBRC(10) is r_infinity in component AP_sensitive_currents_r_gate (dimensionless).
C ALGBRC(27) is tau_r in component AP_sensitive_currents_r_gate (millisecond).
C ALGBRC(48) is i_f_Na in component i_f (picoA).
C ALGBRC(49) is i_f_K in component i_f (picoA).
C CONSTS(19) is g_if in component i_f (nanoS_per_picoF).
C STATES(14) is y in component i_f_y_gate (dimensionless).
C ALGBRC(11) is y_infinity in component i_f_y_gate (dimensionless).
C ALGBRC(28) is tau_y in component i_f_y_gate (millisecond).
C CONSTS(20) is VIf_half in component i_f_y_gate (millivolt).
C CONSTS(21) is g_st in component i_st (nanoS_per_picoF).
C CONSTS(22) is E_st in component i_st (millivolt).
C STATES(15) is qa in component i_st_qa_gate (dimensionless).
C STATES(16) is qi in component i_st_qi_gate (dimensionless).
C ALGBRC(12) is qa_infinity in component i_st_qa_gate (dimensionless).
C ALGBRC(41) is tau_qa in component i_st_qa_gate (millisecond).
C ALGBRC(29) is alpha_qa in component i_st_qa_gate (per_millisecond).
C ALGBRC(35) is beta_qa in component i_st_qa_gate (per_millisecond).
C ALGBRC(36) is qi_infinity in component i_st_qi_gate (dimensionless).
C ALGBRC(42) is tau_qi in component i_st_qi_gate (millisecond).
C ALGBRC(13) is alpha_qi in component i_st_qi_gate (per_millisecond).
C ALGBRC(30) is beta_qi in component i_st_qi_gate (per_millisecond).
C CONSTS(23) is g_b_Na in component i_b_Na (nanoS_per_picoF).
C CONSTS(24) is Km_Kp in component i_NaK (millimolar).
C CONSTS(25) is Km_Nap in component i_NaK (millimolar).
C CONSTS(26) is i_NaK_max in component i_NaK (picoA_per_picoF).
C CONSTS(27) is g_b_Ca in component i_b_Ca (nanoS_per_picoF).
C CONSTS(28) is kNaCa in component i_NaCa (picoA_per_picoF).
C ALGBRC(64) is x1 in component i_NaCa (dimensionless).
C ALGBRC(60) is x2 in component i_NaCa (dimensionless).
C ALGBRC(65) is x3 in component i_NaCa (dimensionless).
C ALGBRC(66) is x4 in component i_NaCa (dimensionless).
C ALGBRC(55) is k41 in component i_NaCa (dimensionless).
C CONSTS(81) is k34 in component i_NaCa (dimensionless).
C ALGBRC(63) is k23 in component i_NaCa (dimensionless).
C ALGBRC(62) is k21 in component i_NaCa (dimensionless).
C ALGBRC(59) is k32 in component i_NaCa (dimensionless).
C CONSTS(80) is k43 in component i_NaCa (dimensionless).
C ALGBRC(57) is k12 in component i_NaCa (dimensionless).
C ALGBRC(58) is k14 in component i_NaCa (dimensionless).
C CONSTS(29) is Qci in component i_NaCa (dimensionless).
C CONSTS(30) is Qn in component i_NaCa (dimensionless).
C CONSTS(31) is Qco in component i_NaCa (dimensionless).
C CONSTS(32) is K3ni in component i_NaCa (millimolar).
C CONSTS(33) is Kci in component i_NaCa (millimolar).
C CONSTS(34) is K1ni in component i_NaCa (millimolar).
C CONSTS(35) is K2ni in component i_NaCa (millimolar).
C CONSTS(36) is Kcni in component i_NaCa (millimolar).
C CONSTS(37) is K3no in component i_NaCa (millimolar).
C CONSTS(38) is K1no in component i_NaCa (millimolar).
C CONSTS(39) is K2no in component i_NaCa (millimolar).
C CONSTS(40) is Kco in component i_NaCa (millimolar).
C CONSTS(79) is RTOnF in component i_NaCa (millivolt).
C ALGBRC(61) is do in component i_NaCa (dimensionless).
C ALGBRC(56) is di in component i_NaCa (dimensionless).
C CONSTS(41) is Cao in component model_parameters (millimolar).
C ALGBRC(68) is j_SRCarel in component j_SRCarel (millimolar_per_millisecond).
C STATES(17) is R in component j_SRCarel (dimensionless).
C STATES(18) is O in component j_SRCarel (dimensionless).
C STATES(19) is I in component j_SRCarel (dimensionless).
C STATES(20) is RI in component j_SRCarel (dimensionless).
C CONSTS(42) is ks in component j_SRCarel (per_millisecond).
C CONSTS(43) is MaxSR in component j_SRCarel (dimensionless).
C CONSTS(44) is MinSR in component j_SRCarel (dimensionless).
C CONSTS(45) is EC50_SR in component j_SRCarel (millimolar).
C CONSTS(46) is HSR in component j_SRCarel (dimensionless).
C ALGBRC(31) is koSRCa in component j_SRCarel (per_millimolar2_millisecond).
C ALGBRC(37) is kiSRCa in component j_SRCarel (per_millimolar_millisecond).
C CONSTS(47) is koCa in component j_SRCarel (per_millimolar2_millisecond).
C CONSTS(48) is kiCa in component j_SRCarel (per_millimolar_millisecond).
C ALGBRC(14) is kCaSR in component j_SRCarel (dimensionless).
C CONSTS(49) is kim in component j_SRCarel (per_millisecond).
C CONSTS(50) is kom in component j_SRCarel (per_millisecond).
C STATES(21) is Ca_jsr in component calcium_dynamics (millimolar).
C ALGBRC(69) is j_Ca_dif in component intracellular_calcium_fluxes (millimolar_per_millisecond).
C ALGBRC(71) is j_up in component intracellular_calcium_fluxes (millimolar_per_millisecond).
C ALGBRC(73) is j_tr in component intracellular_calcium_fluxes (millimolar_per_millisecond).
C CONSTS(51) is tau_dif_Ca in component intracellular_calcium_fluxes (millisecond).
C CONSTS(52) is tau_tr in component intracellular_calcium_fluxes (millisecond).
C CONSTS(53) is P_up in component intracellular_calcium_fluxes (millimolar_per_millisecond).
C CONSTS(54) is K_up in component intracellular_calcium_fluxes (millimolar).
C STATES(22) is Ca_nsr in component calcium_dynamics (millimolar).
C STATES(23) is Cai in component calcium_dynamics (millimolar).
C CONSTS(55) is TC_tot in component calcium_buffering (millimolar).
C CONSTS(56) is TMC_tot in component calcium_buffering (millimolar).
C CONSTS(57) is CM_tot in component calcium_buffering (millimolar).
C CONSTS(58) is CQ_tot in component calcium_buffering (millimolar).
C ALGBRC(72) is delta_fTC in component calcium_buffering (per_millisecond).
C ALGBRC(74) is delta_fTMC in component calcium_buffering (per_millisecond).
C ALGBRC(70) is delta_fCMs in component calcium_buffering (per_millisecond).
C ALGBRC(76) is delta_fCMi in component calcium_buffering (per_millisecond).
C ALGBRC(75) is delta_fCQ in component calcium_buffering (per_millisecond).
C ALGBRC(16) is delta_fTMM in component calcium_buffering (per_millisecond).
C STATES(24) is fTMM in component calcium_buffering (dimensionless).
C STATES(25) is fCMi in component calcium_buffering (dimensionless).
C STATES(26) is fCMs in component calcium_buffering (dimensionless).
C STATES(27) is fTC in component calcium_buffering (dimensionless).
C STATES(28) is fTMC in component calcium_buffering (dimensionless).
C STATES(29) is fCQ in component calcium_buffering (dimensionless).
C CONSTS(59) is kf_TC in component calcium_buffering (per_millimolar_millisecond).
C CONSTS(60) is kf_TMM in component calcium_buffering (per_millimolar_millisecond).
C CONSTS(61) is kf_TMC in component calcium_buffering (per_millimolar_millisecond).
C CONSTS(62) is kf_CM in component calcium_buffering (per_millimolar_millisecond).
C CONSTS(63) is kf_CQ in component calcium_buffering (per_millimolar_millisecond).
C CONSTS(64) is kb_TC in component calcium_buffering (per_millisecond).
C CONSTS(65) is kb_TMC in component calcium_buffering (per_millisecond).
C CONSTS(66) is kb_TMM in component calcium_buffering (per_millisecond).
C CONSTS(67) is kb_CM in component calcium_buffering (per_millisecond).
C CONSTS(68) is kb_CQ in component calcium_buffering (per_millisecond).
C CONSTS(69) is Mgi in component model_parameters (millimolar).
C CONSTS(85) is V_i in component model_parameters (picolitre).
C CONSTS(84) is V_jsr in component model_parameters (picolitre).
C CONSTS(86) is V_nsr in component model_parameters (picolitre).
C CONSTS(83) is V_sub in component model_parameters (picolitre).
C CONSTS(82) is V_cell in component model_parameters (picolitre).
C CONSTS(70) is V_jsr_part in component model_parameters (dimensionless).
C CONSTS(71) is V_i_part in component model_parameters (dimensionless).
C CONSTS(72) is V_nsr_part in component model_parameters (dimensionless).
C CONSTS(73) is R_cell in component model_parameters (micrometre).
C CONSTS(74) is L_cell in component model_parameters (micrometre).
C CONSTS(75) is L_sub in component model_parameters (micrometre).
C RATES(1) is d/dt Vm in component Vm (millivolt).
C RATES(3) is d/dt dL in component i_CaL_dL_gate (dimensionless).
C RATES(4) is d/dt fL in component i_CaL_fL_gate (dimensionless).
C RATES(5) is d/dt fCa in component i_CaL_fCa_gate (dimensionless).
C RATES(6) is d/dt dT in component i_CaT_dT_gate (dimensionless).
C RATES(7) is d/dt fT in component i_CaT_fT_gate (dimensionless).
C RATES(8) is d/dt paS in component i_Kr_pa_gate (dimensionless).
C RATES(9) is d/dt paF in component i_Kr_pa_gate (dimensionless).
C RATES(10) is d/dt pi_ in component i_Kr_pi_gate (dimensionless).
C RATES(11) is d/dt n in component i_Ks_n_gate (dimensionless).
C RATES(12) is d/dt q in component AP_sensitive_currents_q_gate (dimensionless).
C RATES(13) is d/dt r in component AP_sensitive_currents_r_gate (dimensionless).
C RATES(14) is d/dt y in component i_f_y_gate (dimensionless).
C RATES(15) is d/dt qa in component i_st_qa_gate (dimensionless).
C RATES(16) is d/dt qi in component i_st_qi_gate (dimensionless).
C RATES(17) is d/dt R in component j_SRCarel (dimensionless).
C RATES(18) is d/dt O in component j_SRCarel (dimensionless).
C RATES(19) is d/dt I in component j_SRCarel (dimensionless).
C RATES(20) is d/dt RI in component j_SRCarel (dimensionless).
C RATES(27) is d/dt fTC in component calcium_buffering (dimensionless).
C RATES(28) is d/dt fTMC in component calcium_buffering (dimensionless).
C RATES(24) is d/dt fTMM in component calcium_buffering (dimensionless).
C RATES(25) is d/dt fCMi in component calcium_buffering (dimensionless).
C RATES(26) is d/dt fCMs in component calcium_buffering (dimensionless).
C RATES(29) is d/dt fCQ in component calcium_buffering (dimensionless).
C RATES(23) is d/dt Cai in component calcium_dynamics (millimolar).
C RATES(2) is d/dt Ca_sub in component calcium_dynamics (millimolar).
C RATES(22) is d/dt Ca_nsr in component calcium_dynamics (millimolar).
C RATES(21) is d/dt Ca_jsr in component calcium_dynamics (millimolar).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = -65
CONSTS(1) = 32
CONSTS(2) = 8314.4
CONSTS(3) = 310.15
CONSTS(4) = 96485
CONSTS(5) = 10
CONSTS(6) = 140
CONSTS(7) = 140
CONSTS(8) = 5.4
CONSTS(9) = 45
CONSTS(10) = 0.464
STATES(2) = 0.000223
STATES(3) = 0
STATES(4) = 1
STATES(5) = 1
CONSTS(11) = 0.021
CONSTS(12) = 0.00035
CONSTS(13) = 0.1832
CONSTS(14) = 45
STATES(6) = 0
STATES(7) = 1
CONSTS(15) = 0.08113973
STATES(8) = 0
STATES(9) = 0
STATES(10) = 1
CONSTS(16) = 0.0259
STATES(11) = 0
CONSTS(17) = 0.252
CONSTS(18) = 0.02
STATES(12) = 1
STATES(13) = 0
CONSTS(19) = 0.15
STATES(14) = 1
CONSTS(20) = -64
CONSTS(21) = 0.003
CONSTS(22) = 37.4
STATES(15) = 0
STATES(16) = 1
CONSTS(23) = 0.00486
CONSTS(24) = 1.4
CONSTS(25) = 14
CONSTS(26) = 2.88
CONSTS(27) = 0.0006
CONSTS(28) = 187.5
CONSTS(29) = 0.1369
CONSTS(30) = 0.4315
CONSTS(31) = 0
CONSTS(32) = 26.44
CONSTS(33) = 0.0207
CONSTS(34) = 395.3
CONSTS(35) = 2.289
CONSTS(36) = 26.44
CONSTS(37) = 4.663
CONSTS(38) = 1628
CONSTS(39) = 561.4
CONSTS(40) = 3.663
CONSTS(41) = 2
STATES(17) = 0.7499955
STATES(18) = 3.4e-6
STATES(19) = 1.1e-6
STATES(20) = 0.25
CONSTS(42) = 250000
CONSTS(43) = 15
CONSTS(44) = 1
CONSTS(45) = 0.45
CONSTS(46) = 2.5
CONSTS(47) = 10
CONSTS(48) = 0.5
CONSTS(49) = 0.005
CONSTS(50) = 0.06
STATES(21) = 0.029
CONSTS(51) = 0.04
CONSTS(52) = 40
CONSTS(53) = 0.012
CONSTS(54) = 0.0006
STATES(22) = 1.35
STATES(23) = 0.0001
CONSTS(55) = 0.031
CONSTS(56) = 0.062
CONSTS(57) = 0.045
CONSTS(58) = 10
STATES(24) = 0.69
STATES(25) = 0.042
STATES(26) = 0.089
STATES(27) = 0.02
STATES(28) = 0.22
STATES(29) = 0.032
CONSTS(59) = 88.8
CONSTS(60) = 2.277
CONSTS(61) = 227.7
CONSTS(62) = 227.7
CONSTS(63) = 0.534
CONSTS(64) = 0.446
CONSTS(65) = 0.00751
CONSTS(66) = 0.751
CONSTS(67) = 0.542
CONSTS(68) = 0.445
CONSTS(69) = 2.5
CONSTS(70) = 0.0012
CONSTS(71) = 0.46
CONSTS(72) = 0.0116
CONSTS(73) = 4
CONSTS(74) = 70
CONSTS(75) = 0.02
CONSTS(76) = (( CONSTS(2)*CONSTS(3))/CONSTS(4))*log(CONSTS(6)/CONSTS(5))
CONSTS(77) = (( CONSTS(2)*CONSTS(3))/CONSTS(4))*log(CONSTS(8)/CONSTS(7))
CONSTS(78) = (( CONSTS(2)*CONSTS(3))/CONSTS(4))*log((CONSTS(8)+ 0.120000*CONSTS(6))/(CONSTS(7)+ 0.120000*CONSTS(5)))
CONSTS(79) = ( CONSTS(2)*CONSTS(3))/CONSTS(4)
CONSTS(80) = CONSTS(5)/(CONSTS(32)+CONSTS(5))
CONSTS(81) = CONSTS(6)/(CONSTS(37)+CONSTS(6))
CONSTS(82) = 0.00100000* 3.14159265358979*CONSTS(73) ** 2.00000*CONSTS(74)
CONSTS(83) = 0.00100000*2.00000* 3.14159265358979*CONSTS(75)*(CONSTS(73) - CONSTS(75)/2.00000)*CONSTS(74)
CONSTS(84) = CONSTS(70)*CONSTS(82)
CONSTS(85) = CONSTS(71)*CONSTS(82) - CONSTS(83)
CONSTS(86) = CONSTS(72)*CONSTS(82)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(16) = CONSTS(60)*CONSTS(69)*(1.00000 - (STATES(28)+STATES(24))) - CONSTS(66)*STATES(24)
RATES(24) = ALGBRC(16)
ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+35.0000)/7.30000))
ALGBRC(18) = 44.3000+ 257.100*EXP(- (STATES(1)+32.5000)/13.9000 ** 2.00000)
RATES(4) = (ALGBRC(2) - STATES(4))/ALGBRC(18)
ALGBRC(3) = CONSTS(12)/(CONSTS(12)+STATES(2))
ALGBRC(19) = ALGBRC(3)/CONSTS(11)
RATES(5) = (ALGBRC(3) - STATES(5))/ALGBRC(19)
ALGBRC(4) = 1.00000/(1.00000+EXP(- (STATES(1)+26.3000)/6.00000))
ALGBRC(20) = 1.00000/( 1.06800*EXP((STATES(1)+26.3000)/30.0000)+ 1.06800*EXP(- (STATES(1)+26.3000)/30.0000))
RATES(6) = (ALGBRC(4) - STATES(6))/ALGBRC(20)
ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+61.7000)/5.60000))
ALGBRC(21) = 1.00000/( 0.0153000*EXP(- (STATES(1)+61.7000)/83.3000)+ 0.0150000*EXP((STATES(1)+61.7000)/15.3800))
RATES(7) = (ALGBRC(5) - STATES(7))/ALGBRC(21)
ALGBRC(6) = 1.00000/(1.00000+EXP(- (STATES(1)+23.2000)/10.6000))
ALGBRC(22) = 0.846554/( 0.00420000*EXP(STATES(1)/17.0000)+ 0.000150000*EXP(- STATES(1)/21.6000))
RATES(8) = (ALGBRC(6) - STATES(8))/ALGBRC(22)
ALGBRC(23) = 0.846554/( 0.0372000*EXP(STATES(1)/15.9000)+ 0.000960000*EXP(- STATES(1)/22.5000))
RATES(9) = (ALGBRC(6) - STATES(9))/ALGBRC(23)
ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+28.6000)/17.1000))
ALGBRC(24) = 1.00000/( 0.100000*EXP(- STATES(1)/54.6450)+ 0.656000*EXP(STATES(1)/106.157))
RATES(10) = (ALGBRC(7) - STATES(10))/ALGBRC(24)
ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+49.0000)/13.0000))
ALGBRC(26) = 6.06000+39.1020/( 0.570000*EXP( - 0.0800000*(STATES(1)+44.0000))+ 0.0650000*EXP( 0.100000*(STATES(1)+45.9300)))
RATES(12) = (ALGBRC(9) - STATES(12))/ALGBRC(26)
ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(1) - 19.3000)/15.0000))
ALGBRC(27) = 2.75352+14.4052/( 1.03700*EXP( 0.0900000*(STATES(1)+30.6100))+ 0.369000*EXP( - 0.120000*(STATES(1)+23.8400)))
RATES(13) = (ALGBRC(10) - STATES(13))/ALGBRC(27)
ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1) - CONSTS(20))/13.5000))
ALGBRC(28) = 0.716653/(EXP(- (STATES(1)+386.900)/45.3020)+EXP((STATES(1) - 73.0800)/19.2310))
RATES(14) = (ALGBRC(11) - STATES(14))/ALGBRC(28)
ALGBRC(14) = CONSTS(43) - (CONSTS(43) - CONSTS(44))/(1.00000+CONSTS(45)/STATES(21) ** CONSTS(46))
ALGBRC(31) = CONSTS(47)/ALGBRC(14)
ALGBRC(37) = CONSTS(48)*ALGBRC(14)
RATES(17) = ( CONSTS(49)*STATES(20) - ALGBRC(37)*STATES(2)*STATES(17)) - ( ALGBRC(31)*STATES(2) ** 2.00000*STATES(17) - CONSTS(50)*STATES(18))
RATES(18) = ( ALGBRC(31)*STATES(2) ** 2.00000*STATES(17) - CONSTS(50)*STATES(18)) - ( ALGBRC(37)*STATES(2)*STATES(18) - CONSTS(49)*STATES(19))
RATES(19) = ( ALGBRC(37)*STATES(2)*STATES(18) - CONSTS(49)*STATES(19)) - ( CONSTS(50)*STATES(19) - ALGBRC(31)*STATES(2) ** 2.00000*STATES(20))
RATES(20) = ( CONSTS(50)*STATES(19) - ALGBRC(31)*STATES(2) ** 2.00000*STATES(20)) - ( CONSTS(49)*STATES(20) - ALGBRC(37)*STATES(2)*STATES(17))
ALGBRC(8) = 0.0140000/(1.00000+EXP(- (STATES(1) - 40.0000)/9.00000))
ALGBRC(25) = 0.00100000*EXP(- STATES(1)/45.0000)
ALGBRC(34) = ALGBRC(8)/(ALGBRC(8)+ALGBRC(25))
ALGBRC(40) = 1.00000/(ALGBRC(8)+ALGBRC(25))
RATES(11) = (ALGBRC(34) - STATES(11))/ALGBRC(40)
ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+57.0000)/5.00000))
ALGBRC(29) = 1.00000/( 0.150000*EXP(- STATES(1)/11.0000)+ 0.200000*EXP(- STATES(1)/700.000))
ALGBRC(35) = 1.00000/( 16.0000*EXP(STATES(1)/8.00000)+ 15.0000*EXP(STATES(1)/50.0000))
ALGBRC(41) = 1.00000/(ALGBRC(29)+ALGBRC(35))
RATES(15) = (ALGBRC(12) - STATES(15))/ALGBRC(41)
ALGBRC(13) = 1.00000/( 3100.00*EXP(STATES(1)/13.0000)+ 700.000*EXP(STATES(1)/70.0000))
ALGBRC(30) = 1.00000/( 95.0000*EXP(- STATES(1)/10.0000)+ 50.0000*EXP(- STATES(1)/700.000))+0.000229000/(1.00000+EXP(- STATES(1)/5.00000))
ALGBRC(36) = ALGBRC(13)/(ALGBRC(13)+ALGBRC(30))
ALGBRC(42) = 6.65000/(ALGBRC(13)+ALGBRC(30))
RATES(16) = (ALGBRC(36) - STATES(16))/ALGBRC(42)
ALGBRC(1) = 1.00000/(1.00000+EXP(- (STATES(1)+13.5000)/6.00000))
ALGBRC(17) = TERNRY(STATES(1).EQ.- 35.0000, - 35.0000, TERNRY(STATES(1).EQ.0.00000, 1.00000e-05, STATES(1))
ALGBRC(33) = ( - 0.0283900*(ALGBRC(17)+35.0000))/(EXP(- (ALGBRC(17)+35.0000)/2.50000) - 1.00000) - ( 0.0849000*ALGBRC(17))/(EXP(- ALGBRC(17)/4.80000) - 1.00000)
ALGBRC(39) = TERNRY(STATES(1).EQ.5.00000, 5.00001, STATES(1))
ALGBRC(44) = ( 0.0114300*(ALGBRC(39) - 5.00000))/(EXP((ALGBRC(39) - 5.00000)/2.50000) - 1.00000)
ALGBRC(46) = 1.00000/(ALGBRC(33)+ALGBRC(44))
RATES(3) = (ALGBRC(1) - STATES(3))/ALGBRC(46)
ALGBRC(32) = CONSTS(1)*CONSTS(13)*(STATES(1) - CONSTS(14))*STATES(6)*STATES(7)
ALGBRC(15) = CONSTS(1)*CONSTS(10)*(STATES(1) - CONSTS(9))*STATES(3)*STATES(4)*STATES(5)
ALGBRC(48) = CONSTS(1)*0.383300*CONSTS(19)*(STATES(1) - CONSTS(76))*STATES(14) ** 2.00000
ALGBRC(49) = CONSTS(1)*0.616700*CONSTS(19)*(STATES(1) - CONSTS(77))*STATES(14) ** 2.00000
ALGBRC(50) = ALGBRC(48)+ALGBRC(49)
ALGBRC(51) = CONSTS(1)*CONSTS(21)*(STATES(1) - CONSTS(22))*STATES(15)*STATES(16)
ALGBRC(38) = CONSTS(1)*CONSTS(15)*(STATES(1) - CONSTS(77))*( 0.600000*STATES(9)+ 0.400000*STATES(8))*STATES(10)
ALGBRC(43) = CONSTS(1)*CONSTS(16)*(STATES(1) - CONSTS(78))*STATES(11) ** 2.00000
ALGBRC(45) = CONSTS(1)*CONSTS(17)*(STATES(1) - CONSTS(77))*STATES(12)*STATES(13)
ALGBRC(47) = CONSTS(1)*CONSTS(18)*(STATES(1) - CONSTS(77))*STATES(13)
ALGBRC(53) = ( CONSTS(1)*CONSTS(26))/( (1.00000+CONSTS(24)/CONSTS(8) ** 1.20000)*(1.00000+CONSTS(25)/CONSTS(5) ** 1.30000)*(1.00000+EXP(- ((STATES(1) - CONSTS(76))+120.000)/30.0000)))
ALGBRC(55) = EXP(( - CONSTS(30)*STATES(1))/( 2.00000*CONSTS(79)))
ALGBRC(61) = 1.00000+ (CONSTS(41)/CONSTS(40))*(1.00000+EXP(( CONSTS(31)*STATES(1))/CONSTS(79)))+ (CONSTS(6)/CONSTS(38))*(1.00000+ (CONSTS(6)/CONSTS(39))*(1.00000+CONSTS(6)/CONSTS(37)))
ALGBRC(63) = ( (( (CONSTS(6)/CONSTS(38))*CONSTS(6))/CONSTS(39))*(1.00000+CONSTS(6)/CONSTS(37))*EXP(( - CONSTS(30)*STATES(1))/( 2.00000*CONSTS(79))))/ALGBRC(61)
ALGBRC(62) = ( (CONSTS(41)/CONSTS(40))*EXP(( CONSTS(31)*STATES(1))/CONSTS(79)))/ALGBRC(61)
ALGBRC(59) = EXP(( CONSTS(30)*STATES(1))/( 2.00000*CONSTS(79)))
ALGBRC(64) = ALGBRC(55)*CONSTS(81)*(ALGBRC(63)+ALGBRC(62))+ ALGBRC(62)*ALGBRC(59)*(CONSTS(80)+ALGBRC(55))
ALGBRC(56) = 1.00000+ (STATES(2)/CONSTS(33))*(1.00000+EXP(( - CONSTS(29)*STATES(1))/CONSTS(79))+CONSTS(5)/CONSTS(36))+ (CONSTS(5)/CONSTS(34))*(1.00000+ (CONSTS(5)/CONSTS(35))*(1.00000+CONSTS(5)/CONSTS(32)))
ALGBRC(57) = ( (STATES(2)/CONSTS(33))*EXP(( - CONSTS(29)*STATES(1))/CONSTS(79)))/ALGBRC(56)
ALGBRC(58) = ( (( (CONSTS(5)/CONSTS(34))*CONSTS(5))/CONSTS(35))*(1.00000+CONSTS(5)/CONSTS(32))*EXP(( CONSTS(30)*STATES(1))/( 2.00000*CONSTS(79))))/ALGBRC(56)
ALGBRC(60) = ALGBRC(59)*CONSTS(80)*(ALGBRC(58)+ALGBRC(57))+ ALGBRC(55)*ALGBRC(57)*(CONSTS(81)+ALGBRC(59))
ALGBRC(65) = ALGBRC(58)*CONSTS(80)*(ALGBRC(63)+ALGBRC(62))+ ALGBRC(57)*ALGBRC(63)*(CONSTS(80)+ALGBRC(55))
ALGBRC(66) = ALGBRC(63)*CONSTS(81)*(ALGBRC(58)+ALGBRC(57))+ ALGBRC(58)*ALGBRC(62)*(CONSTS(81)+ALGBRC(59))
ALGBRC(67) = ( CONSTS(1)*CONSTS(28)*( ALGBRC(60)*ALGBRC(62) - ALGBRC(64)*ALGBRC(57)))/(ALGBRC(64)+ALGBRC(60)+ALGBRC(65)+ALGBRC(66))
ALGBRC(54) = CONSTS(1)*CONSTS(27)*(STATES(1) - CONSTS(9))
ALGBRC(52) = CONSTS(1)*CONSTS(23)*(STATES(1) - CONSTS(76))
RATES(1) = - (ALGBRC(15)+ALGBRC(32)+ALGBRC(50)+ALGBRC(51)+ALGBRC(38)+ALGBRC(43)+ALGBRC(45)+ALGBRC(47)+ALGBRC(53)+ALGBRC(67)+ALGBRC(54)+ALGBRC(52))/CONSTS(1)
ALGBRC(70) = CONSTS(62)*STATES(2)*(1.00000 - STATES(26)) - CONSTS(67)*STATES(26)
RATES(26) = ALGBRC(70)
ALGBRC(68) = CONSTS(42)*STATES(18)*(STATES(21) - STATES(2))
ALGBRC(69) = (STATES(2) - STATES(23))/CONSTS(51)
RATES(2) = ( ALGBRC(68)*CONSTS(84))/CONSTS(83) - (((ALGBRC(15)+ALGBRC(32)+ALGBRC(54)) - 2.00000*ALGBRC(67))/( 2.00000*CONSTS(4)*CONSTS(83))+ALGBRC(69)+ CONSTS(57)*ALGBRC(70))
ALGBRC(72) = CONSTS(59)*STATES(23)*(1.00000 - STATES(27)) - CONSTS(64)*STATES(27)
RATES(27) = ALGBRC(72)
ALGBRC(71) = CONSTS(53)/(1.00000+CONSTS(54)/STATES(23))
ALGBRC(73) = (STATES(22) - STATES(21))/CONSTS(52)
RATES(22) = ALGBRC(71) - ( ALGBRC(73)*CONSTS(84))/CONSTS(86)
ALGBRC(74) = CONSTS(61)*STATES(23)*(1.00000 - (STATES(28)+STATES(24))) - CONSTS(65)*STATES(28)
RATES(28) = ALGBRC(74)
ALGBRC(75) = CONSTS(63)*STATES(21)*(1.00000 - STATES(29)) - CONSTS(68)*STATES(29)
RATES(29) = ALGBRC(75)
RATES(21) = ALGBRC(73) - (ALGBRC(68)+ CONSTS(58)*ALGBRC(75))
ALGBRC(76) = CONSTS(62)*STATES(23)*(1.00000 - STATES(25)) - CONSTS(67)*STATES(25)
RATES(25) = ALGBRC(76)
RATES(23) = ( ALGBRC(69)*CONSTS(83) - ALGBRC(71)*CONSTS(86))/CONSTS(85) - ( CONSTS(57)*ALGBRC(76)+ CONSTS(55)*ALGBRC(72)+ CONSTS(56)*ALGBRC(74))
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(16) = CONSTS(60)*CONSTS(69)*(1.00000 - (STATES(28)+STATES(24))) - CONSTS(66)*STATES(24)
ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+35.0000)/7.30000))
ALGBRC(18) = 44.3000+ 257.100*EXP(- (STATES(1)+32.5000)/13.9000 ** 2.00000)
ALGBRC(3) = CONSTS(12)/(CONSTS(12)+STATES(2))
ALGBRC(19) = ALGBRC(3)/CONSTS(11)
ALGBRC(4) = 1.00000/(1.00000+EXP(- (STATES(1)+26.3000)/6.00000))
ALGBRC(20) = 1.00000/( 1.06800*EXP((STATES(1)+26.3000)/30.0000)+ 1.06800*EXP(- (STATES(1)+26.3000)/30.0000))
ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+61.7000)/5.60000))
ALGBRC(21) = 1.00000/( 0.0153000*EXP(- (STATES(1)+61.7000)/83.3000)+ 0.0150000*EXP((STATES(1)+61.7000)/15.3800))
ALGBRC(6) = 1.00000/(1.00000+EXP(- (STATES(1)+23.2000)/10.6000))
ALGBRC(22) = 0.846554/( 0.00420000*EXP(STATES(1)/17.0000)+ 0.000150000*EXP(- STATES(1)/21.6000))
ALGBRC(23) = 0.846554/( 0.0372000*EXP(STATES(1)/15.9000)+ 0.000960000*EXP(- STATES(1)/22.5000))
ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+28.6000)/17.1000))
ALGBRC(24) = 1.00000/( 0.100000*EXP(- STATES(1)/54.6450)+ 0.656000*EXP(STATES(1)/106.157))
ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+49.0000)/13.0000))
ALGBRC(26) = 6.06000+39.1020/( 0.570000*EXP( - 0.0800000*(STATES(1)+44.0000))+ 0.0650000*EXP( 0.100000*(STATES(1)+45.9300)))
ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(1) - 19.3000)/15.0000))
ALGBRC(27) = 2.75352+14.4052/( 1.03700*EXP( 0.0900000*(STATES(1)+30.6100))+ 0.369000*EXP( - 0.120000*(STATES(1)+23.8400)))
ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1) - CONSTS(20))/13.5000))
ALGBRC(28) = 0.716653/(EXP(- (STATES(1)+386.900)/45.3020)+EXP((STATES(1) - 73.0800)/19.2310))
ALGBRC(14) = CONSTS(43) - (CONSTS(43) - CONSTS(44))/(1.00000+CONSTS(45)/STATES(21) ** CONSTS(46))
ALGBRC(31) = CONSTS(47)/ALGBRC(14)
ALGBRC(37) = CONSTS(48)*ALGBRC(14)
ALGBRC(8) = 0.0140000/(1.00000+EXP(- (STATES(1) - 40.0000)/9.00000))
ALGBRC(25) = 0.00100000*EXP(- STATES(1)/45.0000)
ALGBRC(34) = ALGBRC(8)/(ALGBRC(8)+ALGBRC(25))
ALGBRC(40) = 1.00000/(ALGBRC(8)+ALGBRC(25))
ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+57.0000)/5.00000))
ALGBRC(29) = 1.00000/( 0.150000*EXP(- STATES(1)/11.0000)+ 0.200000*EXP(- STATES(1)/700.000))
ALGBRC(35) = 1.00000/( 16.0000*EXP(STATES(1)/8.00000)+ 15.0000*EXP(STATES(1)/50.0000))
ALGBRC(41) = 1.00000/(ALGBRC(29)+ALGBRC(35))
ALGBRC(13) = 1.00000/( 3100.00*EXP(STATES(1)/13.0000)+ 700.000*EXP(STATES(1)/70.0000))
ALGBRC(30) = 1.00000/( 95.0000*EXP(- STATES(1)/10.0000)+ 50.0000*EXP(- STATES(1)/700.000))+0.000229000/(1.00000+EXP(- STATES(1)/5.00000))
ALGBRC(36) = ALGBRC(13)/(ALGBRC(13)+ALGBRC(30))
ALGBRC(42) = 6.65000/(ALGBRC(13)+ALGBRC(30))
ALGBRC(1) = 1.00000/(1.00000+EXP(- (STATES(1)+13.5000)/6.00000))
ALGBRC(17) = TERNRY(STATES(1).EQ.- 35.0000, - 35.0000, TERNRY(STATES(1).EQ.0.00000, 1.00000e-05, STATES(1))
ALGBRC(33) = ( - 0.0283900*(ALGBRC(17)+35.0000))/(EXP(- (ALGBRC(17)+35.0000)/2.50000) - 1.00000) - ( 0.0849000*ALGBRC(17))/(EXP(- ALGBRC(17)/4.80000) - 1.00000)
ALGBRC(39) = TERNRY(STATES(1).EQ.5.00000, 5.00001, STATES(1))
ALGBRC(44) = ( 0.0114300*(ALGBRC(39) - 5.00000))/(EXP((ALGBRC(39) - 5.00000)/2.50000) - 1.00000)
ALGBRC(46) = 1.00000/(ALGBRC(33)+ALGBRC(44))
ALGBRC(32) = CONSTS(1)*CONSTS(13)*(STATES(1) - CONSTS(14))*STATES(6)*STATES(7)
ALGBRC(15) = CONSTS(1)*CONSTS(10)*(STATES(1) - CONSTS(9))*STATES(3)*STATES(4)*STATES(5)
ALGBRC(48) = CONSTS(1)*0.383300*CONSTS(19)*(STATES(1) - CONSTS(76))*STATES(14) ** 2.00000
ALGBRC(49) = CONSTS(1)*0.616700*CONSTS(19)*(STATES(1) - CONSTS(77))*STATES(14) ** 2.00000
ALGBRC(50) = ALGBRC(48)+ALGBRC(49)
ALGBRC(51) = CONSTS(1)*CONSTS(21)*(STATES(1) - CONSTS(22))*STATES(15)*STATES(16)
ALGBRC(38) = CONSTS(1)*CONSTS(15)*(STATES(1) - CONSTS(77))*( 0.600000*STATES(9)+ 0.400000*STATES(8))*STATES(10)
ALGBRC(43) = CONSTS(1)*CONSTS(16)*(STATES(1) - CONSTS(78))*STATES(11) ** 2.00000
ALGBRC(45) = CONSTS(1)*CONSTS(17)*(STATES(1) - CONSTS(77))*STATES(12)*STATES(13)
ALGBRC(47) = CONSTS(1)*CONSTS(18)*(STATES(1) - CONSTS(77))*STATES(13)
ALGBRC(53) = ( CONSTS(1)*CONSTS(26))/( (1.00000+CONSTS(24)/CONSTS(8) ** 1.20000)*(1.00000+CONSTS(25)/CONSTS(5) ** 1.30000)*(1.00000+EXP(- ((STATES(1) - CONSTS(76))+120.000)/30.0000)))
ALGBRC(55) = EXP(( - CONSTS(30)*STATES(1))/( 2.00000*CONSTS(79)))
ALGBRC(61) = 1.00000+ (CONSTS(41)/CONSTS(40))*(1.00000+EXP(( CONSTS(31)*STATES(1))/CONSTS(79)))+ (CONSTS(6)/CONSTS(38))*(1.00000+ (CONSTS(6)/CONSTS(39))*(1.00000+CONSTS(6)/CONSTS(37)))
ALGBRC(63) = ( (( (CONSTS(6)/CONSTS(38))*CONSTS(6))/CONSTS(39))*(1.00000+CONSTS(6)/CONSTS(37))*EXP(( - CONSTS(30)*STATES(1))/( 2.00000*CONSTS(79))))/ALGBRC(61)
ALGBRC(62) = ( (CONSTS(41)/CONSTS(40))*EXP(( CONSTS(31)*STATES(1))/CONSTS(79)))/ALGBRC(61)
ALGBRC(59) = EXP(( CONSTS(30)*STATES(1))/( 2.00000*CONSTS(79)))
ALGBRC(64) = ALGBRC(55)*CONSTS(81)*(ALGBRC(63)+ALGBRC(62))+ ALGBRC(62)*ALGBRC(59)*(CONSTS(80)+ALGBRC(55))
ALGBRC(56) = 1.00000+ (STATES(2)/CONSTS(33))*(1.00000+EXP(( - CONSTS(29)*STATES(1))/CONSTS(79))+CONSTS(5)/CONSTS(36))+ (CONSTS(5)/CONSTS(34))*(1.00000+ (CONSTS(5)/CONSTS(35))*(1.00000+CONSTS(5)/CONSTS(32)))
ALGBRC(57) = ( (STATES(2)/CONSTS(33))*EXP(( - CONSTS(29)*STATES(1))/CONSTS(79)))/ALGBRC(56)
ALGBRC(58) = ( (( (CONSTS(5)/CONSTS(34))*CONSTS(5))/CONSTS(35))*(1.00000+CONSTS(5)/CONSTS(32))*EXP(( CONSTS(30)*STATES(1))/( 2.00000*CONSTS(79))))/ALGBRC(56)
ALGBRC(60) = ALGBRC(59)*CONSTS(80)*(ALGBRC(58)+ALGBRC(57))+ ALGBRC(55)*ALGBRC(57)*(CONSTS(81)+ALGBRC(59))
ALGBRC(65) = ALGBRC(58)*CONSTS(80)*(ALGBRC(63)+ALGBRC(62))+ ALGBRC(57)*ALGBRC(63)*(CONSTS(80)+ALGBRC(55))
ALGBRC(66) = ALGBRC(63)*CONSTS(81)*(ALGBRC(58)+ALGBRC(57))+ ALGBRC(58)*ALGBRC(62)*(CONSTS(81)+ALGBRC(59))
ALGBRC(67) = ( CONSTS(1)*CONSTS(28)*( ALGBRC(60)*ALGBRC(62) - ALGBRC(64)*ALGBRC(57)))/(ALGBRC(64)+ALGBRC(60)+ALGBRC(65)+ALGBRC(66))
ALGBRC(54) = CONSTS(1)*CONSTS(27)*(STATES(1) - CONSTS(9))
ALGBRC(52) = CONSTS(1)*CONSTS(23)*(STATES(1) - CONSTS(76))
ALGBRC(70) = CONSTS(62)*STATES(2)*(1.00000 - STATES(26)) - CONSTS(67)*STATES(26)
ALGBRC(68) = CONSTS(42)*STATES(18)*(STATES(21) - STATES(2))
ALGBRC(69) = (STATES(2) - STATES(23))/CONSTS(51)
ALGBRC(72) = CONSTS(59)*STATES(23)*(1.00000 - STATES(27)) - CONSTS(64)*STATES(27)
ALGBRC(71) = CONSTS(53)/(1.00000+CONSTS(54)/STATES(23))
ALGBRC(73) = (STATES(22) - STATES(21))/CONSTS(52)
ALGBRC(74) = CONSTS(61)*STATES(23)*(1.00000 - (STATES(28)+STATES(24))) - CONSTS(65)*STATES(28)
ALGBRC(75) = CONSTS(63)*STATES(21)*(1.00000 - STATES(29)) - CONSTS(68)*STATES(29)
ALGBRC(76) = CONSTS(62)*STATES(23)*(1.00000 - STATES(25)) - CONSTS(67)*STATES(25)
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END