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 7 entries in each of the rate and state variable arrays.
C There are a total of 11 entries in the constant variable array.
C
C
C VOI is time in component environment (hour).
C STATES(1) is T in component T (cell_ml).
C CONSTS(1) is k1 in component model_parameters (cell_per_ml_molecule2_hour).
C STATES(2) is Lm in component Lm (molecule_cell).
C STATES(3) is RT in component RT (molecule_cell).
C STATES(4) is m in component m (cell_ml).
C STATES(5) is LT in component LT (molecule_cell).
C CONSTS(2) is k2 in component model_parameters (molecule_ml2_per_cell3_hour).
C CONSTS(3) is k3 in component model_parameters (per_hour).
C CONSTS(4) is k4 in component model_parameters (per_hour).
C CONSTS(5) is k5 in component model_parameters (ml2_per_cell_molecule_hour).
C STATES(6) is Rm in component Rm (molecule_cell).
C CONSTS(6) is k6 in component model_parameters (molecule_per_cell_hour).
C CONSTS(7) is k7 in component model_parameters (per_hour).
C CONSTS(8) is k8 in component model_parameters (per_molar_per_minute).
C STATES(7) is SL in component SL (molecule_ml).
C CONSTS(9) is k10 in component model_parameters (molecule_per_cell_hour).
C CONSTS(10) is k11 in component model_parameters (per_hour).
C CONSTS(11) is k9 in component model_parameters (per_molar_per_minute).
C RATES(1) is d/dt T in component T (cell_ml).
C RATES(5) is d/dt LT in component LT (molecule_cell).
C RATES(3) is d/dt RT in component RT (molecule_cell).
C RATES(4) is d/dt m in component m (cell_ml).
C RATES(2) is d/dt Lm in component Lm (molecule_cell).
C RATES(6) is d/dt Rm in component Rm (molecule_cell).
C RATES(7) is d/dt SL in component SL (molecule_ml).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 500.0
      CONSTS(1) = 8.38E-10
      STATES(2) = 1E3
      STATES(3) = 1E4
      STATES(4) = 500.0
      STATES(5) = 0.0
      CONSTS(2) = 6E-3
      CONSTS(3) = 5.9413
      CONSTS(4) = 0.35
      CONSTS(5) = 2.52E-9
      STATES(6) = 1E3
      CONSTS(6) = 2.244E3
      CONSTS(7) = 0.35
      CONSTS(8) = 1.92E10
      STATES(7) = 0.0
      CONSTS(9) = 3.11E3
      CONSTS(10) = 13.9
      CONSTS(11) = 87.3E8
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = - ( 1.00000*CONSTS(1)*STATES(4)*STATES(1)*STATES(2)*STATES(3))
      RATES(5) =  CONSTS(2)*STATES(1)*STATES(4) - ( CONSTS(3)*STATES(5)+ CONSTS(4)*STATES(5)+ CONSTS(5)*STATES(4)*STATES(1)*STATES(5)*STATES(6))
      RATES(3) = CONSTS(6) - ( CONSTS(7)*STATES(3)+ 1.00000*CONSTS(8)*STATES(3)*STATES(7))
      RATES(4) = - ( 1.00000*CONSTS(1)*STATES(4)*STATES(1)*STATES(5)*STATES(6))
      RATES(2) = CONSTS(9) - ( CONSTS(3)*STATES(2)+ CONSTS(4)*STATES(2)+ CONSTS(5)*STATES(4)*STATES(1)*STATES(2)*STATES(3))
      RATES(6) = CONSTS(6) - ( CONSTS(7)*STATES(6)+ 1.00000*CONSTS(8)*STATES(6)*STATES(7))
      RATES(7) =  CONSTS(3)*( STATES(2)*STATES(4)+ STATES(5)*STATES(1)) - ( CONSTS(10)*STATES(7)+ 1.00000*CONSTS(11)*( STATES(6)*STATES(4)+ STATES(3)*STATES(1))*STATES(7))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RETURN
      END