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 6 entries in the algebraic variable array.
C There are a total of 1 entries in each of the rate and state variable arrays.
C There are a total of 21 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is ANM in component aldosterone (dimensionless).
C CONSTS(2) is CKE in component aldosterone (monovalent_mEq_per_litre).
C CONSTS(16) is ANMAL in component angiotensin_control_of_aldosterone_secretion (dimensionless).
C CONSTS(3) is ANMALD in component parameter_values (dimensionless).
C CONSTS(17) is OSMAL in component osmotic_control_of_aldosterone_secretion (dimensionless).
C CONSTS(21) is AMR1 in component aldosterone_secretion (dimensionless).
C CONSTS(4) is AMKMUL in component parameter_values (dimensionless).
C CONSTS(5) is ALDINF in component parameter_values (dimensionless).
C CONSTS(6) is ALDKNS in component parameter_values (dimensionless).
C CONSTS(18) is AMRBSC in component aldosterone_secretion (dimensionless).
C CONSTS(19) is AMRT in component aldosterone_secretion (dimensionless).
C CONSTS(20) is AMR in component aldosterone_secretion (dimensionless).
C STATES(1) is AMC in component aldosterone_concentration (dimensionless).
C CONSTS(7) is AMT in component parameter_values (minute).
C ALGBRC(2) is AM in component general_aldosterone_multiplier (dimensionless).
C CONSTS(8) is AM1UL in component parameter_values (dimensionless).
C CONSTS(9) is AM1LL in component parameter_values (dimensionless).
C CONSTS(10) is AMCSNS in component parameter_values (dimensionless).
C CONSTS(11) is ALDMM in component parameter_values (dimensionless).
C ALGBRC(1) is AM1 in component general_aldosterone_multiplier (dimensionless).
C ALGBRC(5) is AMK in component aldosterone_effect_on_cell_membrane_K_transport (dimensionless).
C CONSTS(12) is AMKM in component parameter_values (dimensionless).
C ALGBRC(3) is AMKT in component aldosterone_effect_on_cell_membrane_K_transport (dimensionless).
C ALGBRC(6) is AMNA in component aldosterone_effect_on_cell_membrane_Na_transport (dimensionless).
C CONSTS(13) is AMNAM in component parameter_values (dimensionless).
C CONSTS(14) is AMNAUL in component parameter_values (dimensionless).
C CONSTS(15) is AMNALL in component parameter_values (dimensionless).
C ALGBRC(4) is AMNAT in component aldosterone_effect_on_cell_membrane_Na_transport (dimensionless).
C RATES(1) is d/dt AMC in component aldosterone_concentration (dimensionless).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 0.987545
CONSTS(2) = 4.44092
CONSTS(3) = 2.5
CONSTS(4) = 12
CONSTS(5) = 0
CONSTS(6) = 0
STATES(1) = 1.0
CONSTS(7) = 60
CONSTS(8) = 5
CONSTS(9) = 0
CONSTS(10) = 0.65
CONSTS(11) = 2.5
CONSTS(12) = 0.5
CONSTS(13) = 0.8
CONSTS(14) = 15
CONSTS(15) = 0.04
CONSTS(16) = (CONSTS(1) - 1.00000)*CONSTS(3)+1.00000
CONSTS(17) = (CONSTS(2) - 3.30000)/1.00000
CONSTS(18) = CONSTS(16)*0.909000*CONSTS(17)
CONSTS(19) = (CONSTS(18) - 1.00000)*CONSTS(4)+1.00000
CONSTS(20) = TERNRY(CONSTS(19).LT.0.00000, 0.00000, CONSTS(19))
CONSTS(21) = TERNRY(CONSTS(6).GT.0.00000, CONSTS(6), CONSTS(20)+CONSTS(5))
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(1) = (CONSTS(21) - STATES(1))/CONSTS(7)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = CONSTS(8) - (CONSTS(8) - 1.00000)/( ((CONSTS(9) - 1.00000)/(CONSTS(9) - CONSTS(8)))*(STATES(1) - 1.00000)*CONSTS(10)+1.00000)
ALGBRC(2) = (ALGBRC(1) - 1.00000)*CONSTS(11)+1.00000
ALGBRC(3) = (ALGBRC(2) - 1.00000)*CONSTS(12)+1.00000
ALGBRC(4) = (ALGBRC(2) - 1.00000)*CONSTS(13)+1.00000
ALGBRC(5) = TERNRY(ALGBRC(3).LT.0.200000, 0.200000, ALGBRC(3))
ALGBRC(6) = TERNRY(ALGBRC(4).LT.CONSTS(15), CONSTS(15), TERNRY(ALGBRC(4).GT.CONSTS(14), CONSTS(14), ALGBRC(4))
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END
