C
C There are a total of 7 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 17 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is MDFLW in component angiotensin (L_per_minute).
C CONSTS(16) is ANGSCR in component instantaneous_angiotensin_formation (dimensionless).
C CONSTS(15) is MDFLW3 in component instantaneous_angiotensin_formation (L_per_minute).
C STATES(1) is ANX1 in component time_delayed_angiotensin_formation (dimensionless).
C CONSTS(2) is ANXM in component parameter_values (dimensionless).
C CONSTS(3) is ANV in component parameter_values (minute).
C CONSTS(17) is ANX in component time_delayed_angiotensin_formation (dimensionless).
C ALGBRC(3) is ANPR in component total_angiotensin_formation (dimensionless).
C CONSTS(4) is REK in component parameter_values (dimensionless).
C ALGBRC(1) is ANPRT in component total_angiotensin_formation (dimensionless).
C ALGBRC(5) is ANPR1 in component artificial_angiotensin_formation (dimensionless).
C CONSTS(5) is ANGKNS in component parameter_values (dimensionless).
C CONSTS(6) is ANGINF in component parameter_values (dimensionless).
C STATES(2) is ANC in component angiotensin_concentration (dimensionless).
C CONSTS(7) is ANT in component parameter_values (minute).
C ALGBRC(2) is ANM in component general_angiotensin_multiplier (dimensionless).
C CONSTS(8) is ANMUL in component parameter_values (dimensionless).
C CONSTS(9) is ANMLL in component parameter_values (dimensionless).
C CONSTS(10) is ANCSNS in component parameter_values (dimensionless).
C ALGBRC(6) is ANU in component angiotensin_effect_on_circulation (dimensionless).
C CONSTS(11) is ANUM in component parameter_values (dimensionless).
C CONSTS(12) is ANULL in component parameter_values (dimensionless).
C ALGBRC(4) is ANU1 in component angiotensin_effect_on_circulation (dimensionless).
C ALGBRC(7) is ANUVN in component angiotensin_effect_on_venous_constriction (dimensionless).
C CONSTS(13) is ANUVM in component parameter_values (dimensionless).
C CONSTS(14) is Z12 in component parameter_values (dimensionless).
C RATES(1) is d/dt ANX1 in component time_delayed_angiotensin_formation (dimensionless).
C RATES(2) is d/dt ANC in component angiotensin_concentration (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 1.00051
      STATES(1) = 0.0
      CONSTS(2) = 0
      CONSTS(3) = 5000
      CONSTS(4) = 1
      CONSTS(5) = 0
      CONSTS(6) = 0
      STATES(2) = 0.859476
      CONSTS(7) = 12
      CONSTS(8) = 1.8
      CONSTS(9) = 0.7
      CONSTS(10) = 0.4
      CONSTS(11) = 6
      CONSTS(12) = 0.8
      CONSTS(13) = 0
      CONSTS(14) = 5
      CONSTS(15) = CONSTS(1)
      CONSTS(16) = TERNRY(CONSTS(15).GT.1.00000, 1.00000/(1.00000+ (CONSTS(15) - 1.00000)*72.0000), 10.0000 - 9.00000/(1.00000+ (1.00000 - CONSTS(15))*8.00000))
      CONSTS(17) =  (CONSTS(16) - 1.00000)*CONSTS(2)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = (CONSTS(17) - STATES(1))/CONSTS(3)
      ALGBRC(1) =  (CONSTS(16)+STATES(1))*CONSTS(4)
      ALGBRC(3) = TERNRY(ALGBRC(1).LT.1.00000e-05, 1.00000e-05, ALGBRC(1))
      ALGBRC(5) = TERNRY(CONSTS(5).GT.0.00000, CONSTS(5), ALGBRC(3)+CONSTS(6))
      RATES(2) = (ALGBRC(5) - STATES(2))/CONSTS(7)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) =  (CONSTS(16)+STATES(1))*CONSTS(4)
      ALGBRC(3) = TERNRY(ALGBRC(1).LT.1.00000e-05, 1.00000e-05, ALGBRC(1))
      ALGBRC(5) = TERNRY(CONSTS(5).GT.0.00000, CONSTS(5), ALGBRC(3)+CONSTS(6))
      ALGBRC(2) = CONSTS(8) - (CONSTS(8) - 1.00000)/( ((CONSTS(9) - 1.00000)/(CONSTS(9) - CONSTS(8)))*(STATES(2) - 1.00000)*CONSTS(10)+1.00000)
      ALGBRC(4) =  (ALGBRC(2) - 1.00000)*CONSTS(11)+1.00000
      ALGBRC(6) = TERNRY(ALGBRC(4).LT.CONSTS(12), CONSTS(12), ALGBRC(4))
      ALGBRC(7) =  (ALGBRC(6) - 1.00000)*CONSTS(13)+1.00000
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END