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 0 entries in the algebraic variable array.
C There are a total of 11 entries in each of the rate and state variable arrays.
C There are a total of 18 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C CONSTS(1) is k11 in component kinetic_constants (second_order_rate_constant).
C CONSTS(2) is k11_ in component kinetic_constants (first_order_rate_constant).
C CONSTS(3) is k71 in component kinetic_constants (first_order_rate_constant).
C CONSTS(4) is k72 in component kinetic_constants (first_order_rate_constant).
C CONSTS(5) is k12 in component kinetic_constants (second_order_rate_constant).
C CONSTS(6) is k12_ in component kinetic_constants (first_order_rate_constant).
C CONSTS(7) is k51 in component kinetic_constants (first_order_rate_constant).
C CONSTS(8) is k52 in component kinetic_constants (first_order_rate_constant).
C CONSTS(9) is k41 in component kinetic_constants (first_order_rate_constant).
C CONSTS(10) is k31 in component kinetic_constants (first_order_rate_constant).
C CONSTS(11) is k32 in component kinetic_constants (first_order_rate_constant).
C CONSTS(12) is k42 in component kinetic_constants (first_order_rate_constant).
C CONSTS(13) is k22 in component kinetic_constants (first_order_rate_constant).
C CONSTS(14) is k22_ in component kinetic_constants (first_order_rate_constant).
C CONSTS(15) is k21 in component kinetic_constants (first_order_rate_constant).
C CONSTS(16) is k21_ in component kinetic_constants (first_order_rate_constant).
C STATES(1) is PB in component PB (micromolar).
C CONSTS(17) is kreset1 in component kreset1 (first_order_rate_constant).
C STATES(2) is A1 in component A1 (micromolar).
C STATES(3) is A2 in component A2 (micromolar).
C STATES(4) is PC1 in component PC1 (micromolar).
C STATES(5) is PF1 in component PF1 (micromolar).
C STATES(6) is PG1 in component PG1 (micromolar).
C STATES(7) is PC2 in component PC2 (micromolar).
C STATES(8) is PF2 in component PF2 (micromolar).
C STATES(9) is PD1 in component PD1 (micromolar).
C STATES(10) is PE1 in component PE1 (micromolar).
C STATES(11) is PD2 in component PD2 (micromolar).
C CONSTS(18) is PE2 in component PE2 (micromolar).
C RATES(1) is d/dt PB in component PB (micromolar).
C RATES(2) is d/dt A1 in component A1 (micromolar).
C RATES(3) is d/dt A2 in component A2 (micromolar).
C RATES(4) is d/dt PC1 in component PC1 (micromolar).
C RATES(9) is d/dt PD1 in component PD1 (micromolar).
C RATES(10) is d/dt PE1 in component PE1 (micromolar).
C RATES(5) is d/dt PF1 in component PF1 (micromolar).
C RATES(6) is d/dt PG1 in component PG1 (micromolar).
C RATES(7) is d/dt PC2 in component PC2 (micromolar).
C RATES(11) is d/dt PD2 in component PD2 (micromolar).
C RATES(8) is d/dt PF2 in component PF2 (micromolar).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 110.0
CONSTS(2) = 25.0
CONSTS(3) = 0.1
CONSTS(4) = 60.0
CONSTS(5) = 110.0
CONSTS(6) = 0.2
CONSTS(7) = 7.0
CONSTS(8) = 0.0
CONSTS(9) = 60.0
CONSTS(10) = 500.0
CONSTS(11) = 50.0
CONSTS(12) = 70.0
CONSTS(13) = 100.0
CONSTS(14) = 17.0
CONSTS(15) = 100.0
CONSTS(16) = 0.2
STATES(1) = 0.1
CONSTS(17) = 35.0
STATES(2) = 0.1
STATES(3) = 0.1
STATES(4) = 0.1
STATES(5) = 0.1
STATES(6) = 0.1
STATES(7) = 0.1
STATES(8) = 0.1
STATES(9) = 0.1
STATES(10) = 0.1
STATES(11) = 0.1
CONSTS(18) = 0.1
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(1) = ( CONSTS(2)*STATES(4)+ CONSTS(3)*STATES(5)+ CONSTS(17)*STATES(6)+ CONSTS(6)*STATES(7)+ CONSTS(4)*STATES(8)) - ( CONSTS(1)*STATES(2)*STATES(1)+ CONSTS(5)*STATES(3)*STATES(1))
RATES(2) = CONSTS(2)*STATES(4) - CONSTS(1)*STATES(2)*STATES(1)
RATES(3) = CONSTS(6)*STATES(7) - CONSTS(5)*STATES(3)*STATES(1)
RATES(4) = ( CONSTS(1)*STATES(2)*STATES(1) - (CONSTS(2)+CONSTS(15))*STATES(4))+ CONSTS(16)*STATES(9)
RATES(9) = CONSTS(15)*STATES(4) - (CONSTS(16)+CONSTS(10))*STATES(9)
RATES(10) = CONSTS(10)*STATES(9) - CONSTS(9)*STATES(10)
RATES(5) = CONSTS(9)*STATES(10) - (CONSTS(7)+CONSTS(3))*STATES(5)
RATES(6) = CONSTS(7)*STATES(5) - CONSTS(17)*STATES(6)
RATES(7) = ( CONSTS(5)*STATES(3)*STATES(1) - (CONSTS(6)+CONSTS(13))*STATES(7))+ CONSTS(14)*STATES(11)
RATES(11) = CONSTS(13)*STATES(7) - (CONSTS(14)+CONSTS(11))*STATES(11)
RATES(8) = CONSTS(12)*CONSTS(18) - (CONSTS(8)+CONSTS(4))*STATES(8)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RETURN
END
