C
C There are a total of 13 entries in the algebraic variable array.
C There are a total of 4 entries in each of the rate and state variable arrays.
C There are a total of 14 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C ALGBRC(5) is phi3_c1 in component phi (per_second).
C STATES(1) is h1 in component h1 (dimensionless).
C ALGBRC(1) is phi1_c1 in component phi (second_order_rate).
C ALGBRC(3) is phi2_c1 in component phi (per_second).
C CONSTS(1) is p in component model_parameters (micro_molar).
C ALGBRC(2) is phi_1_c1 in component phi (per_second).
C ALGBRC(9) is phi3_c2 in component phi (per_second).
C STATES(2) is h2 in component h2 (dimensionless).
C ALGBRC(4) is phi1_c2 in component phi (second_order_rate).
C ALGBRC(7) is phi2_c2 in component phi (per_second).
C ALGBRC(6) is phi_1_c2 in component phi (per_second).
C CONSTS(2) is r2 in component model_parameters (second_order_rate).
C CONSTS(3) is R1 in component model_parameters (micro_molar).
C CONSTS(4) is k1 in component model_parameters (micro_molar_per_second).
C CONSTS(5) is R3 in component model_parameters (micro_molar).
C CONSTS(6) is k2 in component model_parameters (micro_molar_per_second).
C CONSTS(7) is r4 in component model_parameters (per_second).
C CONSTS(8) is k3 in component model_parameters (micro_molar_per_second).
C CONSTS(9) is R5 in component model_parameters (micro_molar).
C STATES(3) is c1 in component c1 (micro_molar).
C STATES(4) is c2 in component c2 (micro_molar).
C CONSTS(10) is Vp in component model_parameters (micro_molar_per_second).
C CONSTS(11) is Kp in component model_parameters (micro_molar).
C ALGBRC(8) is j_pump_c1 in component j_pump (micro_molar_per_second).
C ALGBRC(10) is j_pump_c2 in component j_pump (micro_molar_per_second).
C CONSTS(12) is kf in component model_parameters (micro_molar_per_second).
C ALGBRC(11) is j_receptor_c1 in component j_receptor (micro_molar_per_second).
C ALGBRC(12) is j_receptor_c2 in component j_receptor (micro_molar_per_second).
C ALGBRC(13) is j_diffusion in component j_diffusion (micro_molar_per_second).
C CONSTS(13) is D in component model_parameters (per_second).
C CONSTS(14) is j_leak in component model_parameters (micro_molar_per_second).
C RATES(1) is d/dt h1 in component h1 (dimensionless).
C RATES(2) is d/dt h2 in component h2 (dimensionless).
C RATES(3) is d/dt c1 in component c1 (micro_molar).
C RATES(4) is d/dt c2 in component c2 (micro_molar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 0.8
      CONSTS(1) = 0.2778
      STATES(2) = 0.1
      CONSTS(2) = 100
      CONSTS(3) = 6
      CONSTS(4) = 44
      CONSTS(5) = 50
      CONSTS(6) = 26.5
      CONSTS(7) = 20
      CONSTS(8) = 1.6
      CONSTS(9) = 1.6
      STATES(3) = 0.3
      STATES(4) = 0.1
      CONSTS(10) = 1.2
      CONSTS(11) = 0.18
      CONSTS(12) = 28
      CONSTS(13) = 0.01
      CONSTS(14) = 0.2
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(5) = CONSTS(8)/(CONSTS(9)+STATES(3))
      ALGBRC(1) = ( CONSTS(2)*STATES(3))/(CONSTS(3)+STATES(3))
      ALGBRC(3) = (CONSTS(6)+ CONSTS(7)*STATES(3))/(CONSTS(5)+STATES(3))
      ALGBRC(2) = CONSTS(4)/(CONSTS(5)+STATES(3))
      RATES(1) =  ALGBRC(5)*(1.00000 - STATES(1)) - ( ALGBRC(1)*ALGBRC(3)*STATES(1)*CONSTS(1))/( ALGBRC(1)*CONSTS(1)+ALGBRC(2))
      ALGBRC(9) = CONSTS(8)/(CONSTS(9)+STATES(4))
      ALGBRC(4) = ( CONSTS(2)*STATES(4))/(CONSTS(3)+STATES(4))
      ALGBRC(7) = (CONSTS(6)+ CONSTS(7)*STATES(4))/(CONSTS(5)+STATES(4))
      ALGBRC(6) = CONSTS(4)/(CONSTS(5)+STATES(4))
      RATES(2) =  ALGBRC(9)*(1.00000 - STATES(2)) - ( ALGBRC(4)*ALGBRC(7)*STATES(2)*CONSTS(1))/( ALGBRC(4)*CONSTS(1)+ALGBRC(6))
      ALGBRC(8) = ( CONSTS(10)*STATES(3) ** 2.00000)/(CONSTS(11) ** 2.00000+STATES(3) ** 2.00000)
      ALGBRC(11) =  CONSTS(12)*( CONSTS(1)*STATES(1)*ALGBRC(1))/( ALGBRC(1)*CONSTS(1)+ALGBRC(2)) ** 4.00000
      ALGBRC(13) =  CONSTS(13)*(STATES(4) - STATES(3))
      RATES(3) = (ALGBRC(11) - ALGBRC(8))+CONSTS(14)+ALGBRC(13)
      ALGBRC(10) = ( CONSTS(10)*STATES(4) ** 2.00000)/(CONSTS(11) ** 2.00000+STATES(4) ** 2.00000)
      ALGBRC(12) =  CONSTS(12)*( CONSTS(1)*STATES(2)*ALGBRC(4))/( ALGBRC(4)*CONSTS(1)+ALGBRC(6)) ** 4.00000
      RATES(4) = (ALGBRC(12) - ALGBRC(10))+CONSTS(14)+ALGBRC(13)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(5) = CONSTS(8)/(CONSTS(9)+STATES(3))
      ALGBRC(1) = ( CONSTS(2)*STATES(3))/(CONSTS(3)+STATES(3))
      ALGBRC(3) = (CONSTS(6)+ CONSTS(7)*STATES(3))/(CONSTS(5)+STATES(3))
      ALGBRC(2) = CONSTS(4)/(CONSTS(5)+STATES(3))
      ALGBRC(9) = CONSTS(8)/(CONSTS(9)+STATES(4))
      ALGBRC(4) = ( CONSTS(2)*STATES(4))/(CONSTS(3)+STATES(4))
      ALGBRC(7) = (CONSTS(6)+ CONSTS(7)*STATES(4))/(CONSTS(5)+STATES(4))
      ALGBRC(6) = CONSTS(4)/(CONSTS(5)+STATES(4))
      ALGBRC(8) = ( CONSTS(10)*STATES(3) ** 2.00000)/(CONSTS(11) ** 2.00000+STATES(3) ** 2.00000)
      ALGBRC(11) =  CONSTS(12)*( CONSTS(1)*STATES(1)*ALGBRC(1))/( ALGBRC(1)*CONSTS(1)+ALGBRC(2)) ** 4.00000
      ALGBRC(13) =  CONSTS(13)*(STATES(4) - STATES(3))
      ALGBRC(10) = ( CONSTS(10)*STATES(4) ** 2.00000)/(CONSTS(11) ** 2.00000+STATES(4) ** 2.00000)
      ALGBRC(12) =  CONSTS(12)*( CONSTS(1)*STATES(2)*ALGBRC(4))/( ALGBRC(4)*CONSTS(1)+ALGBRC(6)) ** 4.00000
      RETURN
      END