C
C There are a total of 13 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 41 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is PA in component autonomics (mmHg).
C CONSTS(2) is PO2ART in component autonomics (mmHg).
C CONSTS(3) is PLA in component autonomics (mmHg).
C CONSTS(4) is PRA in component autonomics (mmHg).
C CONSTS(5) is PPA in component autonomics (mmHg).
C CONSTS(32) is PA1 in component pressure_driving_autonomic_receptors (mmHg).
C CONSTS(6) is CRRFLX in component parameter_values (mmHg).
C CONSTS(7) is EXE in component parameter_values (mmHg).
C CONSTS(33) is AUC in component chemoreceptors_effect_of_PA (dimensionless).
C CONSTS(8) is AUC1 in component parameter_values (dimensionless).
C CONSTS(38) is AUC3 in component chemoreceptors_effect_of_art_PO2 (dimensionless).
C CONSTS(9) is O2CHMO in component parameter_values (per_mmHg).
C CONSTS(35) is AUC2 in component chemoreceptors_effect_of_art_PO2 (dimensionless).
C ALGBRC(1) is AU6C in component arterial_baroreceptor_reflex (dimensionless).
C CONSTS(10) is AUX in component parameter_values (dimensionless).
C CONSTS(11) is AUK in component parameter_values (per_minute).
C CONSTS(12) is BAROTC in component parameter_values (minute).
C CONSTS(34) is AUB in component arterial_baroreceptor_reflex (dimensionless).
C CONSTS(36) is A1B in component arterial_baroreceptor_reflex (dimensionless).
C CONSTS(37) is AU6A in component arterial_baroreceptor_reflex (dimensionless).
C CONSTS(13) is AU4 in component arterial_baroreceptor_reflex (dimensionless).
C STATES(1) is AU6 in component arterial_baroreceptor_reflex (dimensionless).
C CONSTS(39) is AUN in component CNS_ischemic_reflex (dimensionless).
C CONSTS(14) is AUN1 in component parameter_values (dimensionless).
C CONSTS(40) is AULP in component autonomic_response_to_vasculature_pressure (dimensionless).
C CONSTS(15) is AULPM in component parameter_values (dimensionless).
C CONSTS(41) is AUEX in component autonomic_response_to_exercise (dimensionless).
C CONSTS(16) is EXC in component parameter_values (dimensionless).
C CONSTS(17) is EXCXP in component parameter_values (dimensionless).
C ALGBRC(5) is AUTTL in component total_autonomic_stimulation (dimensionless).
C ALGBRC(3) is AUTTL1 in component total_autonomic_stimulation (dimensionless).
C CONSTS(18) is EXCML in component parameter_values (dimensionless).
C ALGBRC(4) is AU in component actual_autonomic_stimulation (dimensionless).
C CONSTS(19) is AUDMP in component parameter_values (minute).
C CONSTS(20) is AUMAX in component parameter_values (dimensionless).
C CONSTS(21) is AUMIN in component parameter_values (dimensionless).
C CONSTS(22) is AUSLP in component parameter_values (dimensionless).
C ALGBRC(8) is DAU in component actual_autonomic_stimulation (per_minute).
C STATES(2) is AU1 in component actual_autonomic_stimulation (dimensionless).
C ALGBRC(2) is AUT in component actual_autonomic_stimulation (dimensionless).
C ALGBRC(6) is VVR in component autonomic_drive_on_target_organs_and_tissues (litre).
C ALGBRC(9) is AUH in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
C ALGBRC(10) is AUR in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
C ALGBRC(11) is AOM in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
C ALGBRC(12) is AUM in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
C ALGBRC(13) is AVE in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
C CONSTS(23) is VV9 in component parameter_values (litre).
C CONSTS(24) is AUL in component parameter_values (litre).
C CONSTS(25) is AUV in component parameter_values (dimensionless).
C CONSTS(26) is AUS in component parameter_values (dimensionless).
C CONSTS(27) is O2A in component parameter_values (dimensionless).
C CONSTS(28) is AUM1 in component parameter_values (dimensionless).
C CONSTS(29) is AUM2 in component parameter_values (dimensionless).
C CONSTS(30) is AUY in component parameter_values (dimensionless).
C ALGBRC(7) is AUO in component autonomic_drive_on_target_organs_and_tissues (dimensionless).
C CONSTS(31) is MDMP in component parameter_values (dimensionless).
C RATES(1) is d/dt AU6 in component arterial_baroreceptor_reflex (dimensionless).
C RATES(2) is d/dt AU1 in component actual_autonomic_stimulation (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 103.525
      CONSTS(2) = 97.0439
      CONSTS(3) = 2
      CONSTS(4) = 0.00852183
      CONSTS(5) = 15.6376
      CONSTS(6) = 0
      CONSTS(7) = 0
      CONSTS(8) = 0.3
      CONSTS(9) = 0.01
      CONSTS(10) = 1
      CONSTS(11) = 0.004
      CONSTS(12) = 0.16
      CONSTS(13) = -0.060024
      STATES(1) = 1.00132
      CONSTS(14) = 0.5
      CONSTS(15) = 0
      CONSTS(16) = 1
      CONSTS(17) = 1.0
      CONSTS(18) = 0.01
      CONSTS(19) = 0.3
      CONSTS(20) = 5.0
      CONSTS(21) = 0.4
      CONSTS(22) = 1.5
      STATES(2) = 1.00007
      CONSTS(23) = 2.51
      CONSTS(24) = 1.5
      CONSTS(25) = 0.55
      CONSTS(26) = 1
      CONSTS(27) = 0.1
      CONSTS(28) = 3
      CONSTS(29) = 1
      CONSTS(30) = 0
      CONSTS(31) = 0
      CONSTS(32) = TERNRY(CONSTS(6).GT.1.00000e-07, CONSTS(6), CONSTS(1) - CONSTS(7))
      CONSTS(33) = TERNRY(CONSTS(32).LT.80.0000.AND.CONSTS(32).GE.40.0000,  0.00500000*(80.0000 - CONSTS(32))*CONSTS(8), TERNRY(CONSTS(32).LT.40.0000,  0.200000*CONSTS(8), 0.00000)
      CONSTS(34) = TERNRY(CONSTS(32).LT.160.000.AND.CONSTS(32).GE.80.0000,  0.0166670*(160.000 - CONSTS(32)), TERNRY(CONSTS(32).LT.80.0000, 1.33360, 0.00000)
      CONSTS(35) = TERNRY(CONSTS(2).LT.80.0000.AND.CONSTS(2).GE.40.0000,  CONSTS(9)*(80.0000 - CONSTS(2)), TERNRY(CONSTS(2).LT.40.0000,  CONSTS(9)*40.0000, 0.00000)
      CONSTS(36) =  (CONSTS(34) - 1.00000)*CONSTS(10)+1.00000
      CONSTS(37) = CONSTS(36) - CONSTS(13)
      CONSTS(38) = CONSTS(33)+CONSTS(35)
      CONSTS(39) = TERNRY(CONSTS(32).LT.40.0000,  0.0400000*(40.0000 - CONSTS(32))*CONSTS(14), 0.00000)
      CONSTS(40) =  (15.0000/(CONSTS(3)+CONSTS(4)+CONSTS(5)) - 1.00000)*CONSTS(15)+1.00000
      CONSTS(41) = CONSTS(16) ** CONSTS(17)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = (CONSTS(37) - STATES(1))/CONSTS(12)
      ALGBRC(1) = STATES(1)
      ALGBRC(3) =  ( CONSTS(41)*CONSTS(40)*(CONSTS(38)+ALGBRC(1)+CONSTS(39)) - 1.00000)*CONSTS(18)+1.00000
      ALGBRC(5) = TERNRY(ALGBRC(3).LT.0.00000, 0.00000, ALGBRC(3))
      ALGBRC(8) = (ALGBRC(5) - STATES(2))/CONSTS(19)
      RATES(2) = ALGBRC(8)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = STATES(1)
      ALGBRC(3) =  ( CONSTS(41)*CONSTS(40)*(CONSTS(38)+ALGBRC(1)+CONSTS(39)) - 1.00000)*CONSTS(18)+1.00000
      ALGBRC(5) = TERNRY(ALGBRC(3).LT.0.00000, 0.00000, ALGBRC(3))
      ALGBRC(8) = (ALGBRC(5) - STATES(2))/CONSTS(19)
      ALGBRC(2) = CONSTS(20) - (CONSTS(20) - 1.00000)/EXP( CONSTS(22)*(STATES(2) - 1.00000))
      ALGBRC(4) = TERNRY(ALGBRC(2).LT.CONSTS(21), CONSTS(21), ALGBRC(2))
      ALGBRC(6) = (CONSTS(23) -  ALGBRC(4)*CONSTS(24))+CONSTS(24)
      ALGBRC(7) = ALGBRC(4) - 1.00000
      ALGBRC(9) =  ALGBRC(7)*CONSTS(25)+1.00000
      ALGBRC(10) =  ALGBRC(7)*CONSTS(26)+1.00000
      ALGBRC(11) =  ALGBRC(7)*CONSTS(27)+1.00000
      ALGBRC(12) =  ALGBRC(7)*CONSTS(28)+1.00000 ** CONSTS(29)
      ALGBRC(13) =  ALGBRC(7)*CONSTS(30)+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