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 2 entries in the algebraic variable array.
C There are a total of 34 entries in each of the rate and state variable arrays.
C There are a total of 55 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is v_c in component v_c (nanomolar).
C ALGBRC(1) is dv_c_dt in component dv_c_dt (flux).
C CONSTS(1) is k_growth in component rate_constants (first_order_rate_constant).
C STATES(2) is v_n in component v_n (nanomolar).
C ALGBRC(2) is k_volume in component k_volume (dimensionless).
C STATES(3) is mcln2_c in component mcln2_c (nanomolar).
C CONSTS(2) is k_50 in component rate_constants (first_order_rate_constant).
C CONSTS(3) is k_10 in component rate_constants (first_order_rate_constant).
C STATES(4) is mcln2_n in component mcln2_n (nanomolar).
C STATES(5) is mclb5_c in component mclb5_c (nanomolar).
C CONSTS(4) is k_51 in component rate_constants (first_order_rate_constant).
C CONSTS(5) is k_11 in component rate_constants (first_order_rate_constant).
C STATES(6) is mclb5_n in component mclb5_n (nanomolar).
C CONSTS(6) is k_1 in component rate_constants (first_order_rate_constant).
C STATES(7) is sbf_n in component sbf_n (nanomolar).
C CONSTS(7) is k_2 in component rate_constants (first_order_rate_constant).
C CONSTS(8) is k_39 in component rate_constants (first_order_rate_constant).
C CONSTS(9) is k_34 in component rate_constants (second_order_rate_constant).
C CONSTS(10) is k_35 in component rate_constants (first_order_rate_constant).
C STATES(8) is sbfwhi5p_n in component sbfwhi5p_n (nanomolar).
C STATES(9) is whi5_n in component whi5_n (nanomolar).
C STATES(10) is sbfwhi5_n in component sbfwhi5_n (nanomolar).
C STATES(11) is cln3_c in component cln3_c (nanomolar).
C CONSTS(11) is k_6 in component rate_constants (flux).
C CONSTS(12) is k_43 in component rate_constants (first_order_rate_constant).
C CONSTS(13) is k_15 in component rate_constants (first_order_rate_constant).
C STATES(12) is cln2_c in component cln2_c (nanomolar).
C CONSTS(14) is k_3 in component rate_constants (first_order_rate_constant).
C CONSTS(15) is k_26 in component rate_constants (second_order_rate_constant).
C CONSTS(16) is k_27 in component rate_constants (first_order_rate_constant).
C CONSTS(17) is k_12 in component rate_constants (first_order_rate_constant).
C STATES(13) is cdk1_c in component cdk1_c (nanomolar).
C STATES(14) is cdk1cln2_c in component cdk1cln2_c (nanomolar).
C STATES(15) is clb5_c in component clb5_c (nanomolar).
C CONSTS(18) is k_4 in component rate_constants (first_order_rate_constant).
C CONSTS(19) is k_28 in component rate_constants (second_order_rate_constant).
C CONSTS(20) is k_29 in component rate_constants (first_order_rate_constant).
C CONSTS(21) is k_13 in component rate_constants (first_order_rate_constant).
C STATES(16) is cdk1clb5_c in component cdk1clb5_c (nanomolar).
C CONSTS(22) is k_7 in component rate_constants (flux).
C CONSTS(23) is k_44 in component rate_constants (first_order_rate_constant).
C CONSTS(24) is k_49 in component rate_constants (first_order_rate_constant).
C CONSTS(25) is k_16 in component rate_constants (first_order_rate_constant).
C STATES(17) is cdk1_n in component cdk1_n (nanomolar).
C STATES(18) is cln3_n in component cln3_n (nanomolar).
C CONSTS(26) is k_24 in component rate_constants (second_order_rate_constant).
C CONSTS(27) is k_25 in component rate_constants (first_order_rate_constant).
C CONSTS(28) is k_20 in component rate_constants (first_order_rate_constant).
C STATES(19) is cdk1cln3_n in component cdk1cln3_n (nanomolar).
C CONSTS(29) is k_21 in component rate_constants (first_order_rate_constant).
C CONSTS(30) is k_46 in component rate_constants (first_order_rate_constant).
C CONSTS(31) is k_53 in component rate_constants (first_order_rate_constant).
C STATES(20) is cdk1cln2_n in component cdk1cln2_n (nanomolar).
C CONSTS(32) is k_33 in component rate_constants (first_order_rate_constant).
C CONSTS(33) is k_32 in component rate_constants (second_order_rate_constant).
C CONSTS(34) is k_48 in component rate_constants (first_order_rate_constant).
C STATES(21) is cdk1clb5sic1_c in component cdk1clb5sic1_c (nanomolar).
C STATES(22) is sic1_c in component sic1_c (nanomolar).
C STATES(23) is cdk1clb5_n in component cdk1clb5_n (nanomolar).
C CONSTS(35) is k_41 in component rate_constants (first_order_rate_constant).
C STATES(24) is cdk1clb5sic1p_n in component cdk1clb5sic1p_n (nanomolar).
C CONSTS(36) is k_31 in component rate_constants (first_order_rate_constant).
C CONSTS(37) is k_30 in component rate_constants (second_order_rate_constant).
C CONSTS(38) is k_40 in component rate_constants (first_order_rate_constant).
C STATES(25) is cdk1cln3far1_n in component cdk1cln3far1_n (nanomolar).
C STATES(26) is far1_n in component far1_n (nanomolar).
C STATES(27) is cdk1cln3far1p_n in component cdk1cln3far1p_n (nanomolar).
C CONSTS(39) is k_36 in component rate_constants (second_order_rate_constant).
C CONSTS(40) is k_37 in component rate_constants (second_order_rate_constant).
C STATES(28) is cdk1clb5sic1_n in component cdk1clb5sic1_n (nanomolar).
C CONSTS(41) is k_47 in component rate_constants (first_order_rate_constant).
C CONSTS(42) is k_38 in component rate_constants (second_order_rate_constant).
C CONSTS(43) is k_9 in component rate_constants (flux).
C CONSTS(44) is k_18 in component rate_constants (first_order_rate_constant).
C STATES(29) is whi5_c in component whi5_c (nanomolar).
C CONSTS(45) is k_8 in component rate_constants (flux).
C CONSTS(46) is k_45 in component rate_constants (first_order_rate_constant).
C CONSTS(47) is k_17 in component rate_constants (first_order_rate_constant).
C STATES(30) is far1_c in component far1_c (nanomolar).
C CONSTS(48) is k_5 in component rate_constants (flux).
C CONSTS(49) is k_42 in component rate_constants (first_order_rate_constant).
C CONSTS(50) is k_14 in component rate_constants (first_order_rate_constant).
C STATES(31) is whi5p_c in component whi5p_c (nanomolar).
C CONSTS(51) is k_52 in component rate_constants (first_order_rate_constant).
C STATES(32) is whi5p_n in component whi5p_n (nanomolar).
C STATES(33) is sic1P_n in component sic1P_n (nanomolar).
C STATES(34) is far1P_n in component far1P_n (nanomolar).
C CONSTS(52) is k_22 in component rate_constants (first_order_rate_constant).
C CONSTS(53) is k_19 in component rate_constants (first_order_rate_constant).
C CONSTS(54) is k_23 in component rate_constants (first_order_rate_constant).
C RATES(1) is d/dt v_c in component v_c (nanomolar).
C RATES(2) is d/dt v_n in component v_n (nanomolar).
C RATES(3) is d/dt mcln2_c in component mcln2_c (nanomolar).
C RATES(5) is d/dt mclb5_c in component mclb5_c (nanomolar).
C RATES(4) is d/dt mcln2_n in component mcln2_n (nanomolar).
C RATES(6) is d/dt mclb5_n in component mclb5_n (nanomolar).
C RATES(7) is d/dt sbf_n in component sbf_n (nanomolar).
C RATES(11) is d/dt cln3_c in component cln3_c (nanomolar).
C RATES(12) is d/dt cln2_c in component cln2_c (nanomolar).
C RATES(15) is d/dt clb5_c in component clb5_c (nanomolar).
C RATES(13) is d/dt cdk1_c in component cdk1_c (nanomolar).
C RATES(18) is d/dt cln3_n in component cln3_n (nanomolar).
C RATES(17) is d/dt cdk1_n in component cdk1_n (nanomolar).
C RATES(14) is d/dt cdk1cln2_c in component cdk1cln2_c (nanomolar).
C RATES(16) is d/dt cdk1clb5_c in component cdk1clb5_c (nanomolar).
C RATES(20) is d/dt cdk1cln2_n in component cdk1cln2_n (nanomolar).
C RATES(23) is d/dt cdk1clb5_n in component cdk1clb5_n (nanomolar).
C RATES(19) is d/dt cdk1cln3_n in component cdk1cln3_n (nanomolar).
C RATES(10) is d/dt sbfwhi5_n in component sbfwhi5_n (nanomolar).
C RATES(8) is d/dt sbfwhi5p_n in component sbfwhi5p_n (nanomolar).
C RATES(25) is d/dt cdk1cln3far1_n in component cdk1cln3far1_n (nanomolar).
C RATES(28) is d/dt cdk1clb5sic1_n in component cdk1clb5sic1_n (nanomolar).
C RATES(24) is d/dt cdk1clb5sic1p_n in component cdk1clb5sic1p_n (nanomolar).
C RATES(27) is d/dt cdk1cln3far1p_n in component cdk1cln3far1p_n (nanomolar).
C RATES(21) is d/dt cdk1clb5sic1_c in component cdk1clb5sic1_c (nanomolar).
C RATES(22) is d/dt sic1_c in component sic1_c (nanomolar).
C RATES(29) is d/dt whi5_c in component whi5_c (nanomolar).
C RATES(30) is d/dt far1_c in component far1_c (nanomolar).
C RATES(31) is d/dt whi5p_c in component whi5p_c (nanomolar).
C RATES(33) is d/dt sic1P_n in component sic1P_n (nanomolar).
C RATES(34) is d/dt far1P_n in component far1P_n (nanomolar).
C RATES(9) is d/dt whi5_n in component whi5_n (nanomolar).
C RATES(26) is d/dt far1_n in component far1_n (nanomolar).
C RATES(32) is d/dt whi5p_n in component whi5p_n (nanomolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0.5
      CONSTS(1) = 0.0051
      STATES(2) = 0.5
      STATES(3) = 0
      CONSTS(2) = 0.6
      CONSTS(3) = 0.12
      STATES(4) = 0
      STATES(5) = 0
      CONSTS(4) = 0.6
      CONSTS(5) = 0.12
      STATES(6) = 0
      CONSTS(6) = 0.03523
      STATES(7) = 0
      CONSTS(7) = 0.03523
      CONSTS(8) = 1
      CONSTS(9) = 8.46
      CONSTS(10) = 0.0005
      STATES(8) = 0
      STATES(9) = 0
      STATES(10) = 0.025544
      STATES(11) = 0.000485
      CONSTS(11) = 0.00001
      CONSTS(12) = 0.005
      CONSTS(13) = 0.01
      STATES(12) = 0
      CONSTS(14) = 0.32
      CONSTS(15) = 2.82
      CONSTS(16) = 0.55
      CONSTS(17) = 0.1
      STATES(13) = 0.333333
      STATES(14) = 0
      STATES(15) = 0
      CONSTS(18) = 0.32
      CONSTS(19) = 2.82
      CONSTS(20) = 0.55
      CONSTS(21) = 0.35
      STATES(16) = 0
      CONSTS(22) = 0.01
      CONSTS(23) = 0.005
      CONSTS(24) = 0.001
      CONSTS(25) = 0.03
      STATES(17) = 0.0074127
      STATES(18) = 0
      CONSTS(26) = 2.82
      CONSTS(27) = 0.55
      CONSTS(28) = 0.01
      STATES(19) = 0
      CONSTS(29) = 0
      CONSTS(30) = 0.1
      CONSTS(31) = 0.001
      STATES(20) = 0
      CONSTS(32) = 0.55
      CONSTS(33) = 84.6
      CONSTS(34) = 0.012
      STATES(21) = 0
      STATES(22) = 0.039234
      STATES(23) = 0
      CONSTS(35) = 1
      STATES(24) = 0
      CONSTS(36) = 0.55
      CONSTS(37) = 42300
      CONSTS(38) = 1
      STATES(25) = 0
      STATES(26) = 0
      STATES(27) = 0
      CONSTS(39) = 4363.6
      CONSTS(40) = 4363.6
      STATES(28) = 0
      CONSTS(41) = 1
      CONSTS(42) = 4363.6
      CONSTS(43) = 0.00005
      CONSTS(44) = 0.0008
      STATES(29) = 0.073564
      CONSTS(45) = 0.00004
      CONSTS(46) = 0.005
      CONSTS(47) = 0.01
      STATES(30) = 0.0037926
      CONSTS(48) = 0.000042
      CONSTS(49) = 0.005
      CONSTS(50) = 0.01
      STATES(31) = 0
      CONSTS(51) = 0.005
      STATES(32) = 0
      STATES(33) = 0
      STATES(34) = 0
      CONSTS(52) = 0.01
      CONSTS(53) = 0.01
      CONSTS(54) = 0.01
      CONSTS(54) = 0.00000
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(2) = CONSTS(54)
      RATES(4) =  CONSTS(6)*STATES(7) -  CONSTS(2)*STATES(4)
      RATES(6) =  CONSTS(7)*STATES(7) -  CONSTS(4)*STATES(6)
      RATES(7) = ( CONSTS(8)*STATES(8) -  CONSTS(9)*STATES(7)*STATES(9))+ CONSTS(10)*STATES(10)
      RATES(19) = ((( CONSTS(26)*STATES(18)*STATES(17) -  CONSTS(27)*STATES(19))+ CONSTS(36)*STATES(25)) -  CONSTS(37)*STATES(26)*STATES(19))+ CONSTS(38)*STATES(27)
      RATES(10) = ( CONSTS(9)*STATES(7)*STATES(9) -  CONSTS(10)*STATES(10)) -  CONSTS(39)*STATES(10)*STATES(19)
      RATES(8) =  CONSTS(39)*STATES(10)*STATES(19) -  CONSTS(8)*STATES(8)
      RATES(25) = ( CONSTS(37)*STATES(26)*STATES(19) -  CONSTS(36)*STATES(25)) -  CONSTS(40)*STATES(25)*STATES(20)
      RATES(24) =  CONSTS(42)*STATES(28)*STATES(20) -  CONSTS(35)*STATES(24)
      RATES(27) =  CONSTS(40)*STATES(25)*STATES(20) -  CONSTS(38)*STATES(27)
      RATES(33) =  CONSTS(35)*STATES(24)
      RATES(34) =  CONSTS(38)*STATES(27)
      RATES(32) = ( CONSTS(8)*STATES(8) -  CONSTS(51)*STATES(32)) -  CONSTS(54)*STATES(32)
      ALGBRC(1) =  CONSTS(1)*STATES(1)
      RATES(1) = ALGBRC(1)
      RATES(11) = ((CONSTS(11) -  CONSTS(12)*STATES(11)) -  CONSTS(13)*STATES(11)) -  (ALGBRC(1)/STATES(1))*STATES(11)
      RATES(12) = ((( CONSTS(14)*STATES(3) -  CONSTS(15)*STATES(13)*STATES(12))+ CONSTS(16)*STATES(14)) -  CONSTS(17)*STATES(12)) -  (ALGBRC(1)/STATES(1))*STATES(12)
      RATES(15) = ((( CONSTS(18)*STATES(5) -  CONSTS(19)*STATES(15)*STATES(13))+ CONSTS(20)*STATES(16)) -  CONSTS(21)*STATES(15)) -  (ALGBRC(1)/STATES(1))*STATES(15)
      RATES(13) = ((((((CONSTS(22) -  CONSTS(23)*STATES(13))+ CONSTS(24)*STATES(17)+ CONSTS(16)*STATES(14)) -  CONSTS(15)*STATES(13)*STATES(12)) -  CONSTS(19)*STATES(13)*STATES(15))+ CONSTS(20)*STATES(16)) -  CONSTS(25)*STATES(13)) -  (ALGBRC(1)/STATES(1))*STATES(13)
      RATES(16) = (((( CONSTS(19)*STATES(13)*STATES(15) -  CONSTS(20)*STATES(16))+ CONSTS(32)*STATES(21)) -  CONSTS(33)*STATES(22)*STATES(16)) -  CONSTS(34)*STATES(16)) -  (ALGBRC(1)/STATES(1))*STATES(16)
      RATES(21) = (( CONSTS(33)*STATES(22)*STATES(16) -  CONSTS(32)*STATES(21)) -  CONSTS(41)*STATES(21)) -  (ALGBRC(1)/STATES(1))*STATES(21)
      RATES(22) = (((CONSTS(43) -  CONSTS(33)*STATES(22)*STATES(16))+ CONSTS(32)*STATES(21)) -  CONSTS(44)*STATES(22)) -  (ALGBRC(1)/STATES(1))*STATES(22)
      RATES(29) = ((CONSTS(45) -  CONSTS(46)*STATES(29)) -  CONSTS(47)*STATES(29)) -  (ALGBRC(1)/STATES(1))*STATES(29)
      RATES(30) = ((CONSTS(48) -  CONSTS(49)*STATES(30)) -  CONSTS(50)*STATES(30)) -  (ALGBRC(1)/STATES(1))*STATES(30)
      ALGBRC(2) = STATES(2)/STATES(1)
      RATES(3) = ( CONSTS(2)*STATES(4)*ALGBRC(2) -  CONSTS(3)*STATES(3)) -  (ALGBRC(1)/STATES(1))*STATES(3)
      RATES(5) = ( CONSTS(4)*STATES(6)*ALGBRC(2) -  CONSTS(5)*STATES(5)) -  (ALGBRC(1)/STATES(1))*STATES(5)
      RATES(18) = ((( CONSTS(12)*STATES(11))/ALGBRC(2) -  CONSTS(26)*STATES(18)*STATES(17))+ CONSTS(27)*STATES(19)) -  CONSTS(28)*STATES(18)
      RATES(17) = (((( CONSTS(23)*STATES(13))/ALGBRC(2) -  CONSTS(24)*STATES(17)) -  CONSTS(26)*STATES(18)*STATES(17))+ CONSTS(27)*STATES(19)) -  CONSTS(29)*STATES(17)
      RATES(14) = ((( CONSTS(15)*STATES(13)*STATES(12) -  CONSTS(16)*STATES(14)) -  CONSTS(30)*STATES(14))+ CONSTS(31)*STATES(20)*ALGBRC(2)) -  (ALGBRC(1)/STATES(1))*STATES(14)
      RATES(20) = ( CONSTS(30)*STATES(14))/ALGBRC(2) -  CONSTS(31)*STATES(20)
      RATES(23) =  CONSTS(35)*STATES(24)+( CONSTS(34)*STATES(16))/ALGBRC(2)
      RATES(28) = ( CONSTS(41)*STATES(21))/ALGBRC(2) -  CONSTS(42)*STATES(28)*STATES(20)
      RATES(31) =  CONSTS(51)*STATES(32)*ALGBRC(2)
      RATES(9) = (( CONSTS(46)*STATES(29))/ALGBRC(2) -  CONSTS(9)*STATES(9)*STATES(7)) -  CONSTS(52)*STATES(9)
      RATES(26) = ((( CONSTS(49)*STATES(30))/ALGBRC(2) -  CONSTS(37)*STATES(26)*STATES(19))+ CONSTS(36)*STATES(25)) -  CONSTS(53)*STATES(26)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) =  CONSTS(1)*STATES(1)
      ALGBRC(2) = STATES(2)/STATES(1)
      RETURN
      END