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