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 6 entries in each of the rate and state variable arrays.
C There are a total of 25 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is j_1 in component parameters (flux).
C CONSTS(2) is j_2 in component parameters (flux).
C CONSTS(3) is j_3 in component parameters (flux).
C CONSTS(4) is v_d1 in component parameters (flux).
C CONSTS(5) is v_d2 in component parameters (flux).
C CONSTS(6) is v_d3 in component parameters (flux).
C CONSTS(7) is k_d1 in component parameters (first_order_rate_constant).
C CONSTS(8) is k_d2 in component parameters (first_order_rate_constant).
C CONSTS(9) is k_d3 in component parameters (first_order_rate_constant).
C CONSTS(10) is k_c1 in component parameters (first_order_rate_constant).
C CONSTS(11) is k_c2 in component parameters (first_order_rate_constant).
C CONSTS(12) is k_c3 in component parameters (first_order_rate_constant).
C CONSTS(13) is k_m1 in component parameters (nanomolar).
C CONSTS(14) is k_m2 in component parameters (nanomolar).
C CONSTS(15) is k_m3 in component parameters (nanomolar).
C CONSTS(16) is v_12 in component parameters (flux).
C CONSTS(17) is v_11 in component parameters (flux).
C CONSTS(18) is v_10 in component parameters (flux).
C CONSTS(19) is k_120 in component parameters (nanomolar).
C CONSTS(20) is k_110 in component parameters (nanomolar).
C CONSTS(21) is k_100 in component parameters (nanomolar).
C CONSTS(22) is k_d4 in component parameters (first_order_rate_constant).
C CONSTS(23) is k_d5 in component parameters (first_order_rate_constant).
C CONSTS(24) is k_d6 in component parameters (first_order_rate_constant).
C CONSTS(25) is n in component parameters (dimensionless).
C STATES(1) is C_1 in component C_1 (nanomolar).
C STATES(2) is C_2 in component C_2 (nanomolar).
C STATES(3) is T_1 in component T_1 (nanomolar).
C STATES(4) is C_3 in component C_3 (nanomolar).
C STATES(5) is T_2 in component T_2 (nanomolar).
C STATES(6) is T_3 in component T_3 (nanomolar).
C RATES(1) is d/dt C_1 in component C_1 (nanomolar).
C RATES(2) is d/dt C_2 in component C_2 (nanomolar).
C RATES(4) is d/dt C_3 in component C_3 (nanomolar).
C RATES(3) is d/dt T_1 in component T_1 (nanomolar).
C RATES(5) is d/dt T_2 in component T_2 (nanomolar).
C RATES(6) is d/dt T_3 in component T_3 (nanomolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 0.9
      CONSTS(2) = 0.5
      CONSTS(3) = 0.6
      CONSTS(4) = 6
      CONSTS(5) = 1.052
      CONSTS(6) = 3
      CONSTS(7) = 0.8
      CONSTS(8) = 0.9
      CONSTS(9) = 0.8
      CONSTS(10) = 0.2
      CONSTS(11) = 0.22
      CONSTS(12) = 0.6
      CONSTS(13) = 5
      CONSTS(14) = 5
      CONSTS(15) = 5
      CONSTS(16) = 15
      CONSTS(17) = 15
      CONSTS(18) = 15
      CONSTS(19) = 10
      CONSTS(20) = 10
      CONSTS(21) = 10
      CONSTS(22) = 0.16
      CONSTS(23) = 0.16
      CONSTS(24) = 0.16
      CONSTS(25) = 2
      STATES(1) = 0
      STATES(2) = 0
      STATES(3) = 6
      STATES(4) = 0
      STATES(5) = 5
      STATES(6) = 1
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = ( CONSTS(16)*STATES(3) ** CONSTS(25))/(CONSTS(19) ** CONSTS(25)+STATES(3) ** CONSTS(25)+STATES(2) ** CONSTS(25)) -  CONSTS(22)*STATES(1)
      RATES(2) = ( CONSTS(17)*STATES(5) ** CONSTS(25))/(CONSTS(20) ** CONSTS(25)+STATES(5) ** CONSTS(25)+STATES(4) ** CONSTS(25)) -  CONSTS(23)*STATES(2)
      RATES(4) = ( CONSTS(18)*STATES(6) ** CONSTS(25))/(CONSTS(21) ** CONSTS(25)+STATES(6) ** CONSTS(25)+STATES(1) ** CONSTS(25)) -  CONSTS(24)*STATES(4)
      RATES(3) = (CONSTS(1)+( CONSTS(4)*STATES(6) ** CONSTS(25))/(CONSTS(13) ** CONSTS(25)+STATES(6) ** CONSTS(25))+ CONSTS(10)*STATES(1)) -  CONSTS(7)*STATES(3)
      RATES(5) = (CONSTS(2)+( CONSTS(5)*STATES(3) ** CONSTS(25))/(CONSTS(14) ** CONSTS(25)+STATES(3) ** CONSTS(25))+ CONSTS(11)*STATES(2)) -  CONSTS(8)*STATES(5)
      RATES(6) = (CONSTS(3)+( CONSTS(6)*STATES(5) ** CONSTS(25))/(CONSTS(15) ** CONSTS(25)+STATES(5) ** CONSTS(25))+ CONSTS(12)*STATES(4)) -  CONSTS(9)*STATES(6)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RETURN
      END