C
C There are a total of 0 entries in the algebraic variable array.
C There are a total of 0 entries in each of the rate and state variable arrays.
C There are a total of 30 entries in the constant variable array.
C
C
C CONSTS(29) is SHbO2 in component SHbO2 (dimensionless).
C CONSTS(27) is KHbO2 in component KHbO2 (per_molar).
C CONSTS(21) is O2 in component O2 (molar).
C CONSTS(30) is SHbCO2 in component SHbCO2 (dimensionless).
C CONSTS(28) is KHbCO2 in component KHbCO2 (per_molar).
C CONSTS(23) is CO2 in component CO2 (molar).
C CONSTS(25) is Hrbc in component model_parameters (molar).
C CONSTS(1) is K2 in component model_parameters (per_molar).
C CONSTS(2) is K2_ in component model_parameters (molar).
C CONSTS(3) is K3 in component model_parameters (per_molar).
C CONSTS(4) is K3_ in component model_parameters (molar).
C CONSTS(26) is K4 in component K4 (per_molar).
C CONSTS(5) is K5_ in component model_parameters (molar).
C CONSTS(6) is K6_ in component model_parameters (molar).
C CONSTS(7) is O2_S in component model_parameters (micromolar).
C CONSTS(8) is H_S in component model_parameters (nanomolar).
C CONSTS(9) is n1 in component model_parameters (dimensionless).
C CONSTS(10) is n2 in component model_parameters (dimensionless).
C CONSTS(11) is CO2_S in component model_parameters (millimolar).
C CONSTS(12) is K4_ in component model_parameters (per_molar).
C CONSTS(13) is n0 in component model_parameters (dimensionless).
C CONSTS(20) is alpha_O2 in component alpha_O2 (M_mmHg).
C CONSTS(14) is PO2 in component model_parameters (mmHg).
C CONSTS(22) is alpha_CO2 in component alpha_CO2 (M_mmHg).
C CONSTS(15) is PCO2 in component model_parameters (mmHg).
C CONSTS(16) is Wpl in component model_parameters (ml_ml).
C CONSTS(17) is T in component model_parameters (celsius).
C CONSTS(18) is Rrbc in component model_parameters (dimensionless).
C CONSTS(24) is Hpl in component model_parameters (molar).
C CONSTS(19) is pHpl in component model_parameters (pH).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 29.5
      CONSTS(2) = 1E-6
      CONSTS(3) = 25.1
      CONSTS(4) = 1E-6
      CONSTS(5) = 2.63E-8
      CONSTS(6) = 1.91E-8
      CONSTS(7) = 146.0
      CONSTS(8) = 57.5
      CONSTS(9) = 1.06
      CONSTS(10) = 0.12
      CONSTS(11) = 1.31
      CONSTS(12) = 202123.0
      CONSTS(13) = 1.7
      CONSTS(14) = 100.0
      CONSTS(15) = 40.0
      CONSTS(16) = 0.94
      CONSTS(17) = 37.0
      CONSTS(18) = 0.69
      CONSTS(19) = 7.24
      CONSTS(20) =  ((1.37000 -  0.0137000*(CONSTS(17) - 37.0000))+ 0.000580000*CONSTS(17) - 37.0000 ** 2.00000)*(1.00000e-06/CONSTS(16))
      CONSTS(21) =  CONSTS(20)*CONSTS(14)
      CONSTS(22) =  ((3.07000 -  0.0570000*(CONSTS(17) - 37.0000))+ 0.00200000*CONSTS(17) - 37.0000 ** 2.00000)*(1.00000e-05/CONSTS(16))
      CONSTS(23) =  CONSTS(22)*CONSTS(15)
      CONSTS(24) = 10.0000 ** - CONSTS(19)
      CONSTS(25) = CONSTS(24)/CONSTS(18)
      CONSTS(26) =  CONSTS(12)*CONSTS(21)/CONSTS(7) ** CONSTS(13)*CONSTS(25)/CONSTS(8) ** - CONSTS(9)*CONSTS(23)/CONSTS(11) ** - CONSTS(10)
      CONSTS(27) = ( CONSTS(26)*( CONSTS(3)*CONSTS(23)*(1.00000+CONSTS(4)/CONSTS(25))+(1.00000+CONSTS(25)/CONSTS(6))))/( CONSTS(1)*CONSTS(23)*(1.00000+CONSTS(2)/CONSTS(25))+(1.00000+CONSTS(25)/CONSTS(5)))
      CONSTS(28) = ( CONSTS(1)*(1.00000+CONSTS(2)/CONSTS(25))+ CONSTS(3)*CONSTS(26)*(1.00000+CONSTS(4)/CONSTS(25))*CONSTS(21))/(1.00000+CONSTS(25)/CONSTS(5)+ CONSTS(26)*(1.00000+CONSTS(25)/CONSTS(6))*CONSTS(21))
      CONSTS(29) = ( CONSTS(27)*CONSTS(21))/(1.00000+ CONSTS(27)*CONSTS(21))
      CONSTS(30) = ( CONSTS(28)*CONSTS(23))/(1.00000+ CONSTS(28)*CONSTS(23))
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RETURN
      END