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 4 entries in the algebraic variable array.
C There are a total of 17 entries in each of the rate and state variable arrays.
C There are a total of 43 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C STATES(1) is MPF_n in component MPF_n (dimensionless).
C CONSTS(1) is k_in in component parameters (first_order_rate_constant).
C CONSTS(2) is k_out in component parameters (first_order_rate_constant).
C CONSTS(3) is k_dn_ in component parameters (first_order_rate_constant).
C CONSTS(4) is k_dn__ in component parameters (first_order_rate_constant).
C CONSTS(5) is k_wee_ in component parameters (first_order_rate_constant).
C CONSTS(6) is k_wee__ in component parameters (first_order_rate_constant).
C CONSTS(7) is k_stg_ in component parameters (first_order_rate_constant).
C CONSTS(8) is k_stg__ in component parameters (first_order_rate_constant).
C STATES(2) is MPF_c in component MPF_c (dimensionless).
C STATES(3) is FZY in component FZY (dimensionless).
C STATES(4) is Wee1_n in component Wee1_n (dimensionless).
C STATES(5) is StgP_n in component StgP_n (dimensionless).
C STATES(6) is preMPF_n in component preMPF_n (dimensionless).
C STATES(7) is preMPF_c in component preMPF_c (dimensionless).
C CONSTS(9) is k_sc in component parameters (first_order_rate_constant).
C CONSTS(10) is k_dc_ in component parameters (first_order_rate_constant).
C CONSTS(11) is epsilon in component parameters (dimensionless).
C STATES(8) is N in component N (dimensionless).
C STATES(9) is StgP_c in component StgP_c (dimensionless).
C STATES(10) is Wee1_c in component Wee1_c (dimensionless).
C STATES(11) is IE in component IE (dimensionless).
C CONSTS(12) is j_aie in component parameters (dimensionless).
C CONSTS(13) is j_iie in component parameters (dimensionless).
C CONSTS(14) is k_aie in component parameters (first_order_rate_constant).
C CONSTS(15) is k_iie in component parameters (first_order_rate_constant).
C CONSTS(16) is j_afz in component parameters (dimensionless).
C CONSTS(17) is j_ifz in component parameters (dimensionless).
C CONSTS(18) is k_afz in component parameters (first_order_rate_constant).
C CONSTS(19) is k_ifz in component parameters (first_order_rate_constant).
C STATES(12) is Stg_m in component Stg_m (dimensionless).
C CONSTS(20) is k_dm_ in component parameters (first_order_rate_constant).
C CONSTS(21) is k_dm__ in component parameters (first_order_rate_constant).
C CONSTS(22) is j_m in component parameters (dimensionless).
C STATES(13) is Xp in component Xp (dimensionless).
C STATES(14) is Xm in component Xm (dimensionless).
C CONSTS(23) is k_sxm in component parameters (first_order_rate_constant).
C CONSTS(24) is k_sxp in component parameters (first_order_rate_constant).
C CONSTS(25) is k_ins in component parameters (first_order_rate_constant).
C CONSTS(26) is k_outs in component parameters (first_order_rate_constant).
C CONSTS(27) is k_astg_ in component parameters (first_order_rate_constant).
C CONSTS(28) is k_astg__ in component parameters (first_order_rate_constant).
C CONSTS(29) is k_istg in component parameters (first_order_rate_constant).
C CONSTS(30) is k_dstg in component parameters (first_order_rate_constant).
C CONSTS(31) is j_astg in component parameters (dimensionless).
C CONSTS(32) is j_istg in component parameters (dimensionless).
C STATES(15) is Stg_n in component Stg_n (dimensionless).
C STATES(16) is Stg_c in component Stg_c (dimensionless).
C CONSTS(33) is k_sstg in component parameters (first_order_rate_constant).
C CONSTS(34) is k_inw in component parameters (first_order_rate_constant).
C CONSTS(35) is k_outw in component parameters (first_order_rate_constant).
C CONSTS(36) is k_awee in component parameters (first_order_rate_constant).
C CONSTS(37) is k_iwee_ in component parameters (first_order_rate_constant).
C CONSTS(38) is k_iwee__ in component parameters (first_order_rate_constant).
C CONSTS(39) is j_awee in component parameters (dimensionless).
C CONSTS(40) is j_iwee in component parameters (dimensionless).
C STATES(17) is Wee1P_n in component Wee1P_n (dimensionless).
C ALGBRC(1) is Wee1P_c in component Wee1P_c (dimensionless).
C CONSTS(41) is Wee1_T in component Wee1P_c (dimensionless).
C ALGBRC(2) is CycB_T in component CycB_T (dimensionless).
C ALGBRC(3) is Stg_T in component Stg_T (dimensionless).
C ALGBRC(4) is StgP_T in component StgP_T (dimensionless).
C CONSTS(42) is k_ez in component parameters (first_order_rate_constant).
C RATES(1) is d/dt MPF_n in component MPF_n (dimensionless).
C RATES(6) is d/dt preMPF_n in component preMPF_n (dimensionless).
C RATES(2) is d/dt MPF_c in component MPF_c (dimensionless).
C RATES(7) is d/dt preMPF_c in component preMPF_c (dimensionless).
C RATES(11) is d/dt IE in component IE (dimensionless).
C RATES(3) is d/dt FZY in component FZY (dimensionless).
C RATES(12) is d/dt Stg_m in component Stg_m (dimensionless).
C RATES(14) is d/dt Xm in component Xm (dimensionless).
C RATES(13) is d/dt Xp in component Xp (dimensionless).
C RATES(5) is d/dt StgP_n in component StgP_n (dimensionless).
C RATES(15) is d/dt Stg_n in component Stg_n (dimensionless).
C RATES(9) is d/dt StgP_c in component StgP_c (dimensionless).
C RATES(16) is d/dt Stg_c in component Stg_c (dimensionless).
C RATES(4) is d/dt Wee1_n in component Wee1_n (dimensionless).
C RATES(10) is d/dt Wee1_c in component Wee1_c (dimensionless).
C RATES(17) is d/dt Wee1P_n in component Wee1P_n (dimensionless).
C RATES(8) is d/dt N in component N (dimensionless).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = 0
CONSTS(1) = 0.15
CONSTS(2) = 0
CONSTS(3) = 0.01
CONSTS(4) = 1.5
CONSTS(5) = 0.005
CONSTS(6) = 1
CONSTS(7) = 0.2
CONSTS(8) = 2
STATES(2) = 1
STATES(3) = 0
STATES(4) = 0
STATES(5) = 0
STATES(6) = 0
STATES(7) = 0
CONSTS(9) = 0.01
CONSTS(10) = 0.01
CONSTS(11) = 0.00007
STATES(8) = 1
STATES(9) = 0
STATES(10) = 0
STATES(11) = 0
CONSTS(12) = 0.01
CONSTS(13) = 0.01
CONSTS(14) = 1
CONSTS(15) = 0.4
CONSTS(16) = 0.01
CONSTS(17) = 0.01
CONSTS(18) = 1
CONSTS(19) = 0.2
STATES(12) = 1
CONSTS(20) = 0.002
CONSTS(21) = 0.2
CONSTS(22) = 0.05
STATES(13) = 0
STATES(14) = 0
CONSTS(23) = 0.0005
CONSTS(24) = 0.001
CONSTS(25) = 0.08
CONSTS(26) = 0.02
CONSTS(27) = 0
CONSTS(28) = 1
CONSTS(29) = 0.3
CONSTS(30) = 0.015
CONSTS(31) = 0.05
CONSTS(32) = 0.05
STATES(15) = 0
STATES(16) = 1
CONSTS(33) = 0.02
CONSTS(34) = 0.04
CONSTS(35) = 0.01
CONSTS(36) = 0.3
CONSTS(37) = 0.01
CONSTS(38) = 1
CONSTS(39) = 0.05
CONSTS(40) = 0.05
STATES(17) = 0
CONSTS(41) = 0.8
CONSTS(42) = 0.5
CONSTS(42) = 0.00000
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(8) = CONSTS(42)
RATES(1) = ((( CONSTS(1)*STATES(2) - CONSTS(2)*STATES(1)) - (CONSTS(3)+ CONSTS(4)*STATES(3))*STATES(1)) - (CONSTS(5)+ CONSTS(6)*STATES(4))*STATES(1))+ (CONSTS(7)+ CONSTS(8)*STATES(5))*STATES(6)
RATES(6) = ((( CONSTS(1)*STATES(7) - CONSTS(2)*STATES(6)) - (CONSTS(3)+ CONSTS(4)*STATES(3))*STATES(6))+ (CONSTS(5)+ CONSTS(6)*STATES(4))*STATES(1)) - (CONSTS(5)+ CONSTS(8)*STATES(5))*STATES(6)
RATES(2) = (((CONSTS(9) - ( CONSTS(11)*STATES(8)*( CONSTS(1)*STATES(2) - CONSTS(2)*STATES(1)))/(1.00000 - STATES(8)*CONSTS(11))) - CONSTS(10)*STATES(2))+ (CONSTS(7)+ CONSTS(8)*STATES(9))*STATES(7)) - (CONSTS(5)+ CONSTS(6)*STATES(10))*STATES(2)
RATES(7) = ((( - CONSTS(11)*STATES(8)*( CONSTS(1)*STATES(7) - CONSTS(2)*STATES(6)))/(1.00000 - STATES(8)*CONSTS(11)) - CONSTS(10)*STATES(7)) - (CONSTS(7)+ CONSTS(8)*STATES(9))*STATES(7))+ (CONSTS(5)+ CONSTS(6)*STATES(10))*STATES(2)
RATES(11) = ( CONSTS(14)*(1.00000 - STATES(11))*STATES(1))/((CONSTS(12)+1.00000) - STATES(11)) - ( CONSTS(15)*STATES(11))/(CONSTS(13)+STATES(11))
RATES(3) = ( CONSTS(18)*STATES(11)*(1.00000 - STATES(3)))/((CONSTS(16)+1.00000) - STATES(3)) - ( CONSTS(19)*STATES(3))/(CONSTS(17)+STATES(3))
RATES(12) = - (CONSTS(20)/(CONSTS(22)+STATES(12))+ CONSTS(21)*STATES(13))*STATES(12)
RATES(14) = CONSTS(23)*STATES(8)
RATES(13) = CONSTS(24)*STATES(14)
RATES(5) = ((( CONSTS(25)*STATES(9) - CONSTS(26)*STATES(5))+( (CONSTS(27)+ CONSTS(28)*STATES(1))*STATES(15))/(CONSTS(31)+STATES(15))) - ( CONSTS(29)*STATES(5))/(CONSTS(32)+STATES(5))) - CONSTS(30)*STATES(5)
RATES(15) = ((( CONSTS(25)*STATES(16) - CONSTS(26)*STATES(15)) - ( (CONSTS(27)+ CONSTS(28)*STATES(1))*STATES(15))/(CONSTS(31)+STATES(15)))+( CONSTS(29)*STATES(5))/(CONSTS(32)+STATES(5))) - CONSTS(30)*STATES(15)
RATES(9) = (( - CONSTS(30)*STATES(9) - ( CONSTS(11)*STATES(8)*( CONSTS(25)*STATES(9) - CONSTS(26)*STATES(5)))/(1.00000 - STATES(8)*CONSTS(11)))+( (CONSTS(27)+ CONSTS(28)*STATES(2))*STATES(16))/(CONSTS(31)+STATES(16))) - ( CONSTS(29)*STATES(9))/(CONSTS(32)+STATES(9))
RATES(16) = ((( CONSTS(33)*STATES(12) - CONSTS(30)*STATES(16)) - ( CONSTS(11)*STATES(8)*( CONSTS(25)*STATES(16) - CONSTS(26)*STATES(15)))/(1.00000 - STATES(8)*CONSTS(11))) - ( (CONSTS(27)+ CONSTS(28)*STATES(2))*STATES(16))/(CONSTS(31)+STATES(16)))+( CONSTS(29)*STATES(9))/(CONSTS(32)+STATES(9))
RATES(4) = (( CONSTS(34)*STATES(10) - CONSTS(35)*STATES(4))+( CONSTS(36)*STATES(17))/(CONSTS(39)+STATES(17))) - ( (CONSTS(37)+ CONSTS(38)*STATES(1))*STATES(4))/(CONSTS(40)+STATES(4))
ALGBRC(1) = (CONSTS(41) - STATES(8)*CONSTS(11)*(STATES(4)+STATES(17)))/(1.00000 - STATES(8)*CONSTS(11)) - STATES(10)
RATES(10) = (( - ( CONSTS(34)*STATES(10) - CONSTS(35)*STATES(4))*STATES(8)*CONSTS(11))/(1.00000 - STATES(8)*CONSTS(11))+( CONSTS(36)*ALGBRC(1))/(CONSTS(39)+ALGBRC(1))) - ( (CONSTS(37)+ CONSTS(38)*STATES(2))*STATES(10))/(CONSTS(40)+STATES(10))
RATES(17) = (( CONSTS(34)*ALGBRC(1) - CONSTS(35)*STATES(17)) - ( CONSTS(36)*STATES(17))/(CONSTS(39)+STATES(17)))+( (CONSTS(37)+ CONSTS(38)*STATES(1))*STATES(4))/(CONSTS(40)+STATES(4))
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = (CONSTS(41) - STATES(8)*CONSTS(11)*(STATES(4)+STATES(17)))/(1.00000 - STATES(8)*CONSTS(11)) - STATES(10)
ALGBRC(2) = (1.00000 - STATES(8)*CONSTS(11))*(STATES(2)+STATES(7))+ STATES(8)*CONSTS(11)*(STATES(1)+STATES(6))
ALGBRC(3) = (1.00000 - STATES(8)*CONSTS(11))*(STATES(16)+STATES(9))+ STATES(8)*CONSTS(11)*(STATES(15)+STATES(5))
ALGBRC(4) = (1.00000 - STATES(8)*CONSTS(11))*STATES(9)+ STATES(8)*CONSTS(11)*STATES(5)
RETURN
END