C
C There are a total of 1 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 5 entries in the constant variable array.
C
C
C VOI is time in component environment (hour).
C ALGBRC(1) is GnRH in component GnRH (nanomolar).
C STATES(1) is F in component F (dimensionless).
C CONSTS(1) is kfb in component model_parameters (second_order_rate_constant).
C CONSTS(2) is kbf in component model_parameters (first_order_rate_constant).
C STATES(2) is B in component B (dimensionless).
C STATES(3) is R in component R (dimensionless).
C CONSTS(3) is s in component model_parameters (first_order_rate_constant).
C CONSTS(4) is a1 in component model_parameters (first_order_rate_constant).
C CONSTS(5) is a2 in component model_parameters (first_order_rate_constant).
C STATES(4) is C in component C (dimensionless).
C RATES(1) is d/dt F in component F (dimensionless).
C RATES(2) is d/dt B in component B (dimensionless).
C RATES(3) is d/dt R in component R (dimensionless).
C RATES(4) is d/dt C in component C (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 1.0
      CONSTS(1) = 19.35
      CONSTS(2) = 9.91
      STATES(2) = 0.0
      STATES(3) = 2115.0
      CONSTS(3) = 6.80
      CONSTS(4) = 0.0328
      CONSTS(5) = 0.0830
      STATES(4) = 0.0
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(3) = CONSTS(3) -  (CONSTS(4)+ CONSTS(5)*STATES(2))*STATES(3)
      RATES(4) =  (CONSTS(4)+ CONSTS(5)*STATES(2))*STATES(3)
      ALGBRC(1) = TERNRY(VOI.GE.0.00000.AND.VOI.LT.0.0666667, 0.500000, TERNRY(VOI.GE.0.0666667.AND.VOI.LT.0.400000, 0.00000, TERNRY(VOI.GE.0.400000.AND.VOI.LT.0.466667, 0.500000, TERNRY(VOI.GE.0.466667.AND.VOI.LT.2.46667, 0.00000, TERNRY(VOI.GE.2.46667.AND.VOI.LT.2.53333, 0.500000, TERNRY(VOI.GE.2.53333.AND.VOI.LT.2.61667, 0.00000, TERNRY(VOI.GE.2.61667.AND.VOI.LT.2.68333, 0.500000, TERNRY(VOI.GE.2.68333.AND.VOI.LT.4.68333, 0.00000, TERNRY(VOI.GE.4.68333.AND.VOI.LT.4.75000, 0.500000, TERNRY(VOI.GE.4.75000.AND.VOI.LT.4.91667, 0.00000, TERNRY(VOI.GE.4.91667.AND.VOI.LT.4.98333, 0.500000, TERNRY(VOI.GE.4.98333.AND.VOI.LT.6.98333, 0.00000, TERNRY(VOI.GE.6.98333.AND.VOI.LT.7.06667, 0.500000, TERNRY(VOI.GE.7.06667.AND.VOI.LT.7.73333, 0.00000, TERNRY(VOI.GE.7.73333.AND.VOI.LT.7.80000, 0.500000, TERNRY(VOI.GE.7.80000.AND.VOI.LT.9.80000, 0.00000, 0.0/0.0)
      RATES(1) =  CONSTS(2)*STATES(2) -  CONSTS(1)*STATES(1)*ALGBRC(1)
      RATES(2) =  CONSTS(1)*STATES(1)*ALGBRC(1) -  CONSTS(2)*STATES(2)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = TERNRY(VOI.GE.0.00000.AND.VOI.LT.0.0666667, 0.500000, TERNRY(VOI.GE.0.0666667.AND.VOI.LT.0.400000, 0.00000, TERNRY(VOI.GE.0.400000.AND.VOI.LT.0.466667, 0.500000, TERNRY(VOI.GE.0.466667.AND.VOI.LT.2.46667, 0.00000, TERNRY(VOI.GE.2.46667.AND.VOI.LT.2.53333, 0.500000, TERNRY(VOI.GE.2.53333.AND.VOI.LT.2.61667, 0.00000, TERNRY(VOI.GE.2.61667.AND.VOI.LT.2.68333, 0.500000, TERNRY(VOI.GE.2.68333.AND.VOI.LT.4.68333, 0.00000, TERNRY(VOI.GE.4.68333.AND.VOI.LT.4.75000, 0.500000, TERNRY(VOI.GE.4.75000.AND.VOI.LT.4.91667, 0.00000, TERNRY(VOI.GE.4.91667.AND.VOI.LT.4.98333, 0.500000, TERNRY(VOI.GE.4.98333.AND.VOI.LT.6.98333, 0.00000, TERNRY(VOI.GE.6.98333.AND.VOI.LT.7.06667, 0.500000, TERNRY(VOI.GE.7.06667.AND.VOI.LT.7.73333, 0.00000, TERNRY(VOI.GE.7.73333.AND.VOI.LT.7.80000, 0.500000, TERNRY(VOI.GE.7.80000.AND.VOI.LT.9.80000, 0.00000, 0.0/0.0)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END