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 2 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 (minute).
C STATES(1) is GSH in component GSH (millimolar).
C CONSTS(1) is V in component GSH (flux).
C CONSTS(2) is Kms in component GSH (millimolar).
C CONSTS(3) is Kmp in component GSH (millimolar).
C CONSTS(4) is Kmq in component GSH (millimolar).
C CONSTS(5) is Keq in component GSH (millimolar).
C STATES(2) is SDLGSH in component SDLGSH (millimolar).
C CONSTS(6) is D_lactate in component D_lactate (millimolar).
C CONSTS(7) is HTA in component HTA (millimolar).
C CONSTS(8) is V in component SDLGSH (flux).
C CONSTS(9) is Kms in component SDLGSH (millimolar).
C CONSTS(10) is Kmp in component SDLGSH (millimolar).
C CONSTS(11) is Keq in component SDLGSH (dimensionless).
C RATES(1) is d/dt GSH in component GSH (millimolar).
C RATES(2) is d/dt SDLGSH in component SDLGSH (millimolar).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = 1.0
CONSTS(1) = 3.44E-3
CONSTS(2) = 0.49
CONSTS(3) = 0.49
CONSTS(4) = 0.49
CONSTS(5) = 0.49
STATES(2) = 1.0
CONSTS(6) = 0.0
CONSTS(7) = 1.0
CONSTS(8) = 8.12E-2
CONSTS(9) = 0.61
CONSTS(10) = 0.61
CONSTS(11) = 0.61
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(1) = ( (CONSTS(1)/CONSTS(2))*(STATES(2) - ( STATES(1)*CONSTS(6))/CONSTS(5)))/(1.00000+STATES(2)/CONSTS(2)+STATES(1)/CONSTS(3)+CONSTS(6)/CONSTS(4))
RATES(2) = ( (CONSTS(8)/CONSTS(9))*(CONSTS(7) - STATES(2)/CONSTS(11)))/(1.00000+CONSTS(7)/CONSTS(9)+STATES(2)/CONSTS(10))
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RETURN
END
