C
C There are a total of 126 entries in the algebraic variable array.
C There are a total of 46 entries in each of the rate and state variable arrays.
C There are a total of 175 entries in the constant variable array.
C
C
C VOI is t in component main (second).
C CONSTS(1) is rho in component main (Js2_per_m5).
C CONSTS(2) is g in component main (m_per_s2).
C CONSTS(3) is theta_deg in component main (dim).
C CONSTS(175) is theta_rad in component main (radian).
C ALGBRC(1) is z in component main (dim).
C CONSTS(4) is T in component main (second).
C ALGBRC(3) is mt in component main (dim).
C ALGBRC(2) is mta in component main (dim).
C CONSTS(5) is delay in component main (dim).
C CONSTS(6) is To_pk in component main (dim).
C STATES(1) is q_LV in component main (litre).
C CONSTS(7) is q_LV_us in component main (litre).
C ALGBRC(7) is u_LV in component main (kPa).
C ALGBRC(5) is E_LV in component main (kPa_per_L).
C CONSTS(8) is E_LVmax in component main (kPa_per_L).
C CONSTS(9) is E_LVmin in component main (kPa_per_L).
C STATES(2) is q_RV in component main (litre).
C CONSTS(10) is q_RV_us in component main (litre).
C ALGBRC(11) is u_RV in component main (kPa).
C ALGBRC(9) is E_RV in component main (kPa_per_L).
C CONSTS(11) is E_RVmax in component main (kPa_per_L).
C CONSTS(12) is E_RVmin in component main (kPa_per_L).
C STATES(3) is q_LA in component main (litre).
C CONSTS(13) is q_LA_us in component main (litre).
C ALGBRC(13) is u_LA in component main (kPa).
C CONSTS(14) is E_LA in component main (kPa_per_L).
C STATES(4) is q_RA in component main (litre).
C CONSTS(15) is q_RA_us in component main (litre).
C ALGBRC(15) is u_RA in component main (kPa).
C CONSTS(16) is E_RA in component main (kPa_per_L).
C ALGBRC(25) is v_AV in component main (L_per_s).
C CONSTS(17) is R_AV in component main (kPa_s_per_L).
C ALGBRC(24) is v_PV in component main (L_per_s).
C CONSTS(18) is R_PV in component main (kPa_s_per_L).
C ALGBRC(17) is v_TV in component main (L_per_s).
C CONSTS(19) is R_TV in component main (kPa_s_per_L).
C ALGBRC(19) is v_MV in component main (L_per_s).
C CONSTS(20) is R_MV in component main (kPa_s_per_L).
C STATES(5) is q_pulmA in component main (litre).
C CONSTS(21) is q_pulmA_us in component main (litre).
C ALGBRC(21) is u_pulmA in component main (kPa).
C CONSTS(22) is E_pulmA in component main (kPa_per_L).
C ALGBRC(28) is v_lungA in component main (L_per_s).
C CONSTS(23) is R_lungA in component main (kPa_s_per_L).
C STATES(6) is q_lungC in component main (litre).
C CONSTS(24) is q_lungC_us in component main (litre).
C ALGBRC(26) is u_lungC in component main (kPa).
C CONSTS(25) is E_lungC in component main (kPa_per_L).
C ALGBRC(30) is v_pulmV in component main (L_per_s).
C CONSTS(26) is R_pulmV in component main (kPa_s_per_L).
C STATES(7) is q_AA in component main (litre).
C CONSTS(27) is q_AA_us in component main (litre).
C ALGBRC(22) is u_AA in component main (kPa).
C CONSTS(28) is E_AA in component main (kPa_per_L).
C ALGBRC(29) is v_brachiocephalicA in component main (L_per_s).
C CONSTS(29) is R_brachiocephalicA in component main (kPa_s_per_L).
C STATES(8) is q_brachiocephalicA in component main (litre).
C CONSTS(30) is q_brachiocephalicA_us in component main (litre).
C ALGBRC(27) is u_brachiocephalicA in component main (kPa).
C CONSTS(31) is E_brachiocephalicA in component main (kPa_per_L).
C ALGBRC(32) is v_rSubclavianA in component main (L_per_s).
C CONSTS(32) is R_rSubclavianA in component main (kPa_s_per_L).
C ALGBRC(36) is v_rCCA in component main (L_per_s).
C CONSTS(33) is R_rCCA in component main (kPa_s_per_L).
C STATES(9) is q_rSubclavianA in component main (litre).
C CONSTS(34) is q_rSubclavianA_us in component main (litre).
C ALGBRC(31) is u_rSubclavianA in component main (kPa).
C CONSTS(35) is E_rSubclavianA in component main (kPa_per_L).
C ALGBRC(45) is v_rBrachialA in component main (L_per_s).
C CONSTS(36) is R_rBrachialA in component main (kPa_s_per_L).
C ALGBRC(37) is v_rVertebralA in component main (L_per_s).
C CONSTS(37) is R_rVertebralA in component main (kPa_s_per_L).
C ALGBRC(38) is v_lCCA in component main (L_per_s).
C CONSTS(38) is R_lCCA in component main (kPa_s_per_L).
C ALGBRC(34) is v_lSubclavianA in component main (L_per_s).
C CONSTS(39) is R_lSubclavianA in component main (kPa_s_per_L).
C STATES(10) is q_lSubclavianA in component main (litre).
C CONSTS(40) is q_lSubclavianA_us in component main (litre).
C ALGBRC(33) is u_lSubclavianA in component main (kPa).
C CONSTS(41) is E_lSubclavianA in component main (kPa_per_L).
C ALGBRC(42) is v_lBrachialA in component main (L_per_s).
C CONSTS(42) is R_lBrachialA in component main (kPa_s_per_L).
C ALGBRC(39) is v_lVertebralA in component main (L_per_s).
C CONSTS(43) is R_lVertebralA in component main (kPa_s_per_L).
C STATES(11) is q_headC in component main (litre).
C CONSTS(44) is q_headC_us in component main (litre).
C ALGBRC(35) is u_headC in component main (kPa).
C CONSTS(45) is E_headC in component main (kPa_per_L).
C ALGBRC(40) is v_headV in component main (L_per_s).
C CONSTS(46) is R_headV in component main (kPa_s_per_L).
C STATES(12) is q_lArmC in component main (litre).
C CONSTS(47) is q_lArmC_us in component main (litre).
C ALGBRC(41) is u_lArmC in component main (kPa).
C CONSTS(48) is E_lArmC in component main (kPa_per_L).
C ALGBRC(43) is v_lArmV in component main (L_per_s).
C CONSTS(49) is R_lArmV in component main (kPa_s_per_L).
C STATES(13) is q_rArmC in component main (litre).
C CONSTS(50) is q_rArmC_us in component main (litre).
C ALGBRC(44) is u_rArmC in component main (kPa).
C CONSTS(51) is E_rArmC in component main (kPa_per_L).
C ALGBRC(46) is v_rArmV in component main (L_per_s).
C CONSTS(52) is R_rArmV in component main (kPa_s_per_L).
C STATES(14) is q_A1 in component main (litre).
C CONSTS(53) is q_A1_us in component main (litre).
C STATES(15) is q_A2 in component main (litre).
C CONSTS(54) is q_A2_us in component main (litre).
C STATES(16) is q_A3 in component main (litre).
C CONSTS(55) is q_A3_us in component main (litre).
C STATES(17) is q_A4 in component main (litre).
C CONSTS(56) is q_A4_us in component main (litre).
C STATES(18) is q_A5 in component main (litre).
C CONSTS(57) is q_A5_us in component main (litre).
C STATES(19) is q_A6 in component main (litre).
C CONSTS(58) is q_A6_us in component main (litre).
C STATES(20) is q_A7 in component main (litre).
C CONSTS(59) is q_A7_us in component main (litre).
C STATES(21) is q_A8 in component main (litre).
C CONSTS(60) is q_A8_us in component main (litre).
C ALGBRC(55) is v_A1 in component main (L_per_s).
C ALGBRC(56) is v_A2 in component main (L_per_s).
C ALGBRC(57) is v_A3 in component main (L_per_s).
C ALGBRC(58) is v_A4 in component main (L_per_s).
C ALGBRC(59) is v_A5 in component main (L_per_s).
C ALGBRC(60) is v_A6 in component main (L_per_s).
C ALGBRC(61) is v_A7 in component main (L_per_s).
C ALGBRC(62) is v_A8 in component main (L_per_s).
C ALGBRC(47) is u_A1 in component main (kPa).
C ALGBRC(48) is u_A2 in component main (kPa).
C ALGBRC(49) is u_A3 in component main (kPa).
C ALGBRC(50) is u_A4 in component main (kPa).
C ALGBRC(51) is u_A5 in component main (kPa).
C ALGBRC(52) is u_A6 in component main (kPa).
C ALGBRC(53) is u_A7 in component main (kPa).
C ALGBRC(54) is u_A8 in component main (kPa).
C CONSTS(61) is E_A1 in component main (kPa_per_L).
C CONSTS(62) is E_A2 in component main (kPa_per_L).
C CONSTS(63) is E_A3 in component main (kPa_per_L).
C CONSTS(64) is E_A4 in component main (kPa_per_L).
C CONSTS(65) is E_A5 in component main (kPa_per_L).
C CONSTS(66) is E_A6 in component main (kPa_per_L).
C CONSTS(67) is E_A7 in component main (kPa_per_L).
C CONSTS(68) is E_A8 in component main (kPa_per_L).
C CONSTS(69) is R_A1 in component main (kPa_s_per_L).
C CONSTS(70) is R_A2 in component main (kPa_s_per_L).
C CONSTS(71) is R_A3 in component main (kPa_s_per_L).
C CONSTS(72) is R_A4 in component main (kPa_s_per_L).
C CONSTS(73) is R_A5 in component main (kPa_s_per_L).
C CONSTS(74) is R_A6 in component main (kPa_s_per_L).
C CONSTS(75) is R_A7 in component main (kPa_s_per_L).
C CONSTS(76) is R_A8 in component main (kPa_s_per_L).
C ALGBRC(64) is v_lInfPhrenicA in component main (L_per_s).
C CONSTS(77) is R_lInfPhrenicA in component main (kPa_s_per_L).
C STATES(22) is q_lAdrenalGlandC in component main (litre).
C CONSTS(78) is q_lAdrenalGlandC_us in component main (litre).
C ALGBRC(63) is u_lAdrenalGlandC in component main (kPa).
C CONSTS(79) is E_lAdrenalGlandC in component main (kPa_per_L).
C ALGBRC(66) is v_lAdrenalV in component main (L_per_s).
C CONSTS(80) is R_lAdrenalV in component main (kPa_s_per_L).
C STATES(23) is q_lInfPhrenicV in component main (litre).
C CONSTS(81) is q_lInfPhrenicV_us in component main (litre).
C ALGBRC(65) is u_lInfPhrenicV in component main (kPa).
C CONSTS(82) is E_lInfPhrenicV in component main (kPa_per_L).
C ALGBRC(105) is v_lInfPhrenicV in component main (L_per_s).
C CONSTS(83) is R_lInfPhrenicV in component main (kPa_s_per_L).
C ALGBRC(68) is v_rInfPhrenicA in component main (L_per_s).
C CONSTS(84) is R_rInfPhrenicA in component main (kPa_s_per_L).
C STATES(24) is q_rAdrenalGlandC in component main (litre).
C CONSTS(85) is q_rAdrenalGlandC_us in component main (litre).
C ALGBRC(67) is u_rAdrenalGlandC in component main (kPa).
C CONSTS(86) is E_rAdrenalGlandC in component main (kPa_per_L).
C ALGBRC(70) is v_rAdrenalV in component main (L_per_s).
C CONSTS(87) is R_rAdrenalV in component main (kPa_s_per_L).
C STATES(25) is q_rInfPhrenicV in component main (litre).
C CONSTS(88) is q_rInfPhrenicV_us in component main (litre).
C ALGBRC(69) is u_rInfPhrenicV in component main (kPa).
C CONSTS(89) is E_rInfPhrenicV in component main (kPa_per_L).
C ALGBRC(107) is v_rInfPhrenicV in component main (L_per_s).
C CONSTS(90) is R_rInfPhrenicV in component main (kPa_s_per_L).
C ALGBRC(72) is v_celiacA in component main (L_per_s).
C CONSTS(91) is R_celiacA in component main (kPa_s_per_L).
C STATES(26) is q_celiacA in component main (litre).
C CONSTS(92) is q_celiacA_us in component main (litre).
C ALGBRC(71) is u_celiacA in component main (kPa).
C CONSTS(93) is E_celiacA in component main (kPa_per_L).
C CONSTS(94) is z_celiacA in component main (meter).
C ALGBRC(74) is v_supMesA in component main (L_per_s).
C CONSTS(95) is R_supMesA in component main (kPa_s_per_L).
C STATES(27) is q_supMesA in component main (litre).
C CONSTS(96) is q_supMesA_us in component main (litre).
C ALGBRC(73) is u_supMesA in component main (kPa).
C CONSTS(97) is E_supMesA in component main (kPa_per_L).
C CONSTS(98) is z_supMesA in component main (meter).
C ALGBRC(76) is v_lRenalA in component main (L_per_s).
C CONSTS(99) is R_lRenalA in component main (kPa_s_per_L).
C STATES(28) is q_lGlomerulusC in component main (litre).
C CONSTS(100) is q_lGlomerulusC_us in component main (litre).
C ALGBRC(75) is u_lGlomerulusC in component main (kPa).
C CONSTS(101) is E_lGlomerulusC in component main (kPa_per_L).
C ALGBRC(78) is v_lEffArteriol in component main (L_per_s).
C CONSTS(102) is R_lEffArteriol in component main (kPa_s_per_L).
C STATES(29) is q_lPeritubularC in component main (litre).
C CONSTS(103) is q_lPeritubularC_us in component main (litre).
C ALGBRC(77) is u_lPeritubularC in component main (kPa).
C CONSTS(104) is E_lPeritubularC in component main (kPa_per_L).
C ALGBRC(111) is v_lRenalV in component main (L_per_s).
C CONSTS(105) is R_lRenalV in component main (kPa_s_per_L).
C ALGBRC(80) is v_rRenalA in component main (L_per_s).
C CONSTS(106) is R_rRenalA in component main (kPa_s_per_L).
C STATES(30) is q_rGlomerulusC in component main (litre).
C CONSTS(107) is q_rGlomerulusC_us in component main (litre).
C ALGBRC(79) is u_rGlomerulusC in component main (kPa).
C CONSTS(108) is E_rGlomerulusC in component main (kPa_per_L).
C ALGBRC(82) is v_rEffArteriol in component main (L_per_s).
C CONSTS(109) is R_rEffArteriol in component main (kPa_s_per_L).
C STATES(31) is q_rPeritubularC in component main (litre).
C CONSTS(110) is q_rPeritubularC_us in component main (litre).
C ALGBRC(81) is u_rPeritubularC in component main (kPa).
C CONSTS(111) is E_rPeritubularC in component main (kPa_per_L).
C ALGBRC(113) is v_rRenalV in component main (L_per_s).
C CONSTS(112) is R_rRenalV in component main (kPa_s_per_L).
C ALGBRC(93) is v_infMesA in component main (L_per_s).
C CONSTS(113) is R_infMesA in component main (kPa_s_per_L).
C ALGBRC(84) is v_stomachA in component main (L_per_s).
C CONSTS(114) is R_stomachA in component main (kPa_s_per_L).
C STATES(32) is q_stomachC in component main (litre).
C CONSTS(115) is q_stomachC_us in component main (litre).
C ALGBRC(83) is u_stomachC in component main (kPa).
C CONSTS(116) is E_stomachC in component main (kPa_per_L).
C CONSTS(117) is z_stomachC in component main (meter).
C ALGBRC(96) is v_stomachV in component main (L_per_s).
C CONSTS(118) is R_stomachV in component main (kPa_s_per_L).
C ALGBRC(86) is v_spleenA in component main (L_per_s).
C CONSTS(119) is R_spleenA in component main (kPa_s_per_L).
C STATES(33) is q_spleenC in component main (litre).
C CONSTS(120) is q_spleenC_us in component main (litre).
C ALGBRC(85) is u_spleenC in component main (kPa).
C CONSTS(121) is E_spleenC in component main (kPa_per_L).
C CONSTS(122) is z_spleenC in component main (meter).
C ALGBRC(97) is v_spleenV in component main (L_per_s).
C CONSTS(123) is R_spleenV in component main (kPa_s_per_L).
C ALGBRC(88) is v_pancreas1A in component main (L_per_s).
C ALGBRC(89) is v_pancreas2A in component main (L_per_s).
C CONSTS(124) is R_pancreas1A in component main (kPa_s_per_L).
C CONSTS(125) is R_pancreas2A in component main (kPa_s_per_L).
C STATES(34) is q_pancreasC in component main (litre).
C CONSTS(126) is q_pancreasC_us in component main (litre).
C ALGBRC(87) is u_pancreasC in component main (kPa).
C CONSTS(127) is E_pancreasC in component main (kPa_per_L).
C CONSTS(128) is z_pancreasC in component main (meter).
C ALGBRC(98) is v_pancreasV in component main (L_per_s).
C CONSTS(129) is R_pancreasV in component main (kPa_s_per_L).
C ALGBRC(91) is v_intestineA in component main (L_per_s).
C CONSTS(130) is R_intestineA in component main (kPa_s_per_L).
C STATES(35) is q_intestineC in component main (litre).
C CONSTS(131) is q_intestineC_us in component main (litre).
C ALGBRC(90) is u_intestineC in component main (kPa).
C CONSTS(132) is E_intestineC in component main (kPa_per_L).
C CONSTS(133) is z_ in component main (meter).
C ALGBRC(99) is v_intestineV in component main (L_per_s).
C CONSTS(134) is R_intestineV in component main (kPa_s_per_L).
C ALGBRC(94) is v_colonA in component main (L_per_s).
C CONSTS(135) is R_colonA in component main (kPa_s_per_L).
C STATES(36) is q_colonC in component main (litre).
C CONSTS(136) is q_colonC_us in component main (litre).
C ALGBRC(92) is u_colonC in component main (kPa).
C CONSTS(137) is E_colonC in component main (kPa_per_L).
C CONSTS(138) is z_colonC in component main (meter).
C ALGBRC(100) is v_colonV in component main (L_per_s).
C CONSTS(139) is R_colonV in component main (kPa_s_per_L).
C STATES(37) is q_portalV in component main (litre).
C CONSTS(140) is q_portalV_us in component main (litre).
C ALGBRC(95) is u_portalV in component main (kPa).
C CONSTS(141) is E_portalV in component main (kPa_per_L).
C CONSTS(142) is z_portalV in component main (meter).
C ALGBRC(102) is v_portalV in component main (L_per_s).
C CONSTS(143) is R_portalV in component main (kPa_s_per_L).
C ALGBRC(103) is v_liverA in component main (L_per_s).
C CONSTS(144) is R_liverA in component main (kPa_s_per_L).
C STATES(38) is q_liverC in component main (litre).
C CONSTS(145) is q_liverC_us in component main (litre).
C ALGBRC(101) is u_liverC in component main (kPa).
C CONSTS(146) is E_liverC in component main (kPa_per_L).
C CONSTS(147) is z_liverC in component main (meter).
C ALGBRC(109) is v_liverV in component main (L_per_s).
C CONSTS(148) is R_liverV in component main (kPa_s_per_L).
C STATES(39) is q_V1 in component main (litre).
C CONSTS(149) is q_V1_us in component main (litre).
C STATES(40) is q_V2 in component main (litre).
C CONSTS(150) is q_V2_us in component main (litre).
C STATES(41) is q_V3 in component main (litre).
C CONSTS(151) is q_V3_us in component main (litre).
C STATES(42) is q_V4 in component main (litre).
C CONSTS(152) is q_V4_us in component main (litre).
C STATES(43) is q_V5 in component main (litre).
C CONSTS(153) is q_V5_us in component main (litre).
C STATES(44) is q_V6 in component main (litre).
C CONSTS(154) is q_V6_us in component main (litre).
C ALGBRC(115) is v_V1 in component main (L_per_s).
C ALGBRC(116) is v_V2 in component main (L_per_s).
C ALGBRC(117) is v_V3 in component main (L_per_s).
C ALGBRC(118) is v_V4 in component main (L_per_s).
C ALGBRC(119) is v_V5 in component main (L_per_s).
C ALGBRC(120) is v_V6 in component main (L_per_s).
C ALGBRC(104) is u_V1 in component main (kPa).
C ALGBRC(106) is u_V2 in component main (kPa).
C ALGBRC(108) is u_V3 in component main (kPa).
C ALGBRC(110) is u_V4 in component main (kPa).
C ALGBRC(112) is u_V5 in component main (kPa).
C ALGBRC(114) is u_V6 in component main (kPa).
C CONSTS(155) is E_V1 in component main (kPa_per_L).
C CONSTS(156) is E_V2 in component main (kPa_per_L).
C CONSTS(157) is E_V3 in component main (kPa_per_L).
C CONSTS(158) is E_V4 in component main (kPa_per_L).
C CONSTS(159) is E_V5 in component main (kPa_per_L).
C CONSTS(160) is E_V6 in component main (kPa_per_L).
C CONSTS(161) is R_V1 in component main (kPa_s_per_L).
C CONSTS(162) is R_V2 in component main (kPa_s_per_L).
C CONSTS(163) is R_V3 in component main (kPa_s_per_L).
C CONSTS(164) is R_V4 in component main (kPa_s_per_L).
C CONSTS(165) is R_V5 in component main (kPa_s_per_L).
C CONSTS(166) is R_V6 in component main (kPa_s_per_L).
C ALGBRC(122) is v_lCommonIliacA in component main (L_per_s).
C CONSTS(167) is R_lCommonIliacA in component main (kPa_s_per_L).
C STATES(45) is q_lLegC in component main (litre).
C CONSTS(168) is q_lLegC_us in component main (litre).
C ALGBRC(121) is u_lLegC in component main (kPa).
C CONSTS(169) is E_lLegC in component main (kPa_per_L).
C ALGBRC(123) is v_lCommonIliacV in component main (L_per_s).
C CONSTS(170) is R_lCommonIliacV in component main (kPa_s_per_L).
C ALGBRC(125) is v_rCommonIliacA in component main (L_per_s).
C CONSTS(171) is R_rCommonIliacA in component main (kPa_s_per_L).
C STATES(46) is q_rLegC in component main (litre).
C CONSTS(172) is q_rLegC_us in component main (litre).
C ALGBRC(124) is u_rLegC in component main (kPa).
C CONSTS(173) is E_rLegC in component main (kPa_per_L).
C ALGBRC(126) is v_rCommonIliacV in component main (L_per_s).
C CONSTS(174) is R_rCommonIliacV in component main (kPa_s_per_L).
C ALGBRC(4) is q_TotHeart in component main (litre).
C ALGBRC(6) is q_TotLung in component main (litre).
C ALGBRC(8) is q_TotHead in component main (litre).
C ALGBRC(10) is q_TotThoracic in component main (litre).
C ALGBRC(12) is q_TotAorta in component main (litre).
C ALGBRC(14) is q_TotRenal in component main (litre).
C ALGBRC(16) is q_TotGut in component main (litre).
C ALGBRC(18) is q_TotVenaCava in component main (litre).
C ALGBRC(20) is q_TotLimbs in component main (litre).
C ALGBRC(23) is q_TotBlood in component main (litre).
C RATES(1) is d/dt q_LV in component main (litre).
C RATES(2) is d/dt q_RV in component main (litre).
C RATES(3) is d/dt q_LA in component main (litre).
C RATES(4) is d/dt q_RA in component main (litre).
C RATES(5) is d/dt q_pulmA in component main (litre).
C RATES(6) is d/dt q_lungC in component main (litre).
C RATES(7) is d/dt q_AA in component main (litre).
C RATES(8) is d/dt q_brachiocephalicA in component main (litre).
C RATES(9) is d/dt q_rSubclavianA in component main (litre).
C RATES(10) is d/dt q_lSubclavianA in component main (litre).
C RATES(11) is d/dt q_headC in component main (litre).
C RATES(12) is d/dt q_lArmC in component main (litre).
C RATES(13) is d/dt q_rArmC in component main (litre).
C RATES(14) is d/dt q_A1 in component main (litre).
C RATES(15) is d/dt q_A2 in component main (litre).
C RATES(16) is d/dt q_A3 in component main (litre).
C RATES(17) is d/dt q_A4 in component main (litre).
C RATES(18) is d/dt q_A5 in component main (litre).
C RATES(19) is d/dt q_A6 in component main (litre).
C RATES(20) is d/dt q_A7 in component main (litre).
C RATES(21) is d/dt q_A8 in component main (litre).
C RATES(22) is d/dt q_lAdrenalGlandC in component main (litre).
C RATES(23) is d/dt q_lInfPhrenicV in component main (litre).
C RATES(24) is d/dt q_rAdrenalGlandC in component main (litre).
C RATES(25) is d/dt q_rInfPhrenicV in component main (litre).
C RATES(26) is d/dt q_celiacA in component main (litre).
C RATES(27) is d/dt q_supMesA in component main (litre).
C RATES(28) is d/dt q_lGlomerulusC in component main (litre).
C RATES(29) is d/dt q_lPeritubularC in component main (litre).
C RATES(30) is d/dt q_rGlomerulusC in component main (litre).
C RATES(31) is d/dt q_rPeritubularC in component main (litre).
C RATES(32) is d/dt q_stomachC in component main (litre).
C RATES(33) is d/dt q_spleenC in component main (litre).
C RATES(34) is d/dt q_pancreasC in component main (litre).
C RATES(35) is d/dt q_intestineC in component main (litre).
C RATES(36) is d/dt q_colonC in component main (litre).
C RATES(37) is d/dt q_portalV in component main (litre).
C RATES(38) is d/dt q_liverC in component main (litre).
C RATES(39) is d/dt q_V1 in component main (litre).
C RATES(40) is d/dt q_V2 in component main (litre).
C RATES(41) is d/dt q_V3 in component main (litre).
C RATES(42) is d/dt q_V4 in component main (litre).
C RATES(43) is d/dt q_V5 in component main (litre).
C RATES(44) is d/dt q_V6 in component main (litre).
C RATES(45) is d/dt q_lLegC in component main (litre).
C RATES(46) is d/dt q_rLegC in component main (litre).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 1000
      CONSTS(2) = 9.81
      CONSTS(3) = 90
      CONSTS(4) = 1
      CONSTS(5) = 0.1
      CONSTS(6) = 0.1
      STATES(1) = 2.24894867772612
      CONSTS(7) = 0.005
      CONSTS(8) = 320
      CONSTS(9) = 6.667
      STATES(2) = 2.24212874273495
      CONSTS(10) = 0.01
      CONSTS(11) = 70
      CONSTS(12) = 2
      STATES(3) = 1.50038322871398
      CONSTS(13) = 0.004
      CONSTS(14) = 10
      STATES(4) = 0.582873822068968
      CONSTS(15) = 0.004
      CONSTS(16) = 10
      CONSTS(17) = 1
      CONSTS(18) = 1
      CONSTS(19) = 1
      CONSTS(20) = 1
      STATES(5) = 0.274822042097049
      CONSTS(21) = 0.2
      CONSTS(22) = 200
      CONSTS(23) = 1
      STATES(6) = 0.274821056256821
      CONSTS(24) = 0.2
      CONSTS(25) = 200
      CONSTS(26) = 1
      STATES(7) = 1.43716360225001
      CONSTS(27) = 0.12
      CONSTS(28) = 80
      CONSTS(29) = 100
      STATES(8) = 0.658321340279823
      CONSTS(30) = 0.2
      CONSTS(31) = 200
      CONSTS(32) = 100
      CONSTS(33) = 100
      STATES(9) = 0.578637981993128
      CONSTS(34) = 0.2
      CONSTS(35) = 200
      CONSTS(36) = 100
      CONSTS(37) = 100
      CONSTS(38) = 100
      CONSTS(39) = 100
      STATES(10) = 0.609560299960373
      CONSTS(40) = 0.2
      CONSTS(41) = 200
      CONSTS(42) = 100
      CONSTS(43) = 100
      STATES(11) = 0.575700923676085
      CONSTS(44) = 0.2
      CONSTS(45) = 200
      CONSTS(46) = 100
      STATES(12) = 0.380415025779693
      CONSTS(47) = 0.06
      CONSTS(48) = 150
      CONSTS(49) = 100
      STATES(13) = 0.380415025779692
      CONSTS(50) = 0.06
      CONSTS(51) = 150
      CONSTS(52) = 100
      STATES(14) = 0.0726453866278764
      CONSTS(53) = 0.06
      STATES(15) = 0.0722595904606296
      CONSTS(54) = 0.06
      STATES(16) = 0.0717727262354559
      CONSTS(55) = 0.06
      STATES(17) = 0.0715328891139091
      CONSTS(56) = 0.06
      STATES(18) = 0.0714293557557069
      CONSTS(57) = 0.06
      STATES(19) = 0.0713497500336812
      CONSTS(58) = 0.06
      STATES(20) = 0.0712949341825069
      CONSTS(59) = 0.06
      STATES(21) = 0.0712201889071758
      CONSTS(60) = 0.06
      CONSTS(61) = 400
      CONSTS(62) = 400
      CONSTS(63) = 400
      CONSTS(64) = 400
      CONSTS(65) = 400
      CONSTS(66) = 400
      CONSTS(67) = 400
      CONSTS(68) = 400
      CONSTS(69) = 10000000
      CONSTS(70) = 100
      CONSTS(71) = 100
      CONSTS(72) = 100
      CONSTS(73) = 100
      CONSTS(74) = 100
      CONSTS(75) = 100
      CONSTS(76) = 100
      CONSTS(77) = 100
      STATES(22) = 0.226013520134949
      CONSTS(78) = 0.2
      CONSTS(79) = 200
      CONSTS(80) = 100
      STATES(23) = 0.226425525942081
      CONSTS(81) = 0.2
      CONSTS(82) = 200
      CONSTS(83) = 100
      CONSTS(84) = 100
      STATES(24) = 0.224729694488936
      CONSTS(85) = 0.2
      CONSTS(86) = 200
      CONSTS(87) = 100
      STATES(25) = 0.224899115493557
      CONSTS(88) = 0.2
      CONSTS(89) = 200
      CONSTS(90) = 100
      CONSTS(91) = 100
      STATES(26) = 0.0715428919998289
      CONSTS(92) = 0.06
      CONSTS(93) = 400
      CONSTS(94) = 0.2
      CONSTS(95) = 100
      STATES(27) = 0.071414457768405
      CONSTS(96) = 0.06
      CONSTS(97) = 400
      CONSTS(98) = 0.2
      CONSTS(99) = 100
      STATES(28) = 0.222843940557744
      CONSTS(100) = 0.2
      CONSTS(101) = 200
      CONSTS(102) = 100
      STATES(29) = 0.222892809988734
      CONSTS(103) = 0.2
      CONSTS(104) = 200
      CONSTS(105) = 100
      CONSTS(106) = 100
      STATES(30) = 0.222680450943523
      CONSTS(107) = 0.2
      CONSTS(108) = 200
      CONSTS(109) = 100
      STATES(31) = 0.222721287850087
      CONSTS(110) = 0.2
      CONSTS(111) = 200
      CONSTS(112) = 100
      CONSTS(113) = 100
      CONSTS(114) = 100
      STATES(32) = 0.105910755847991
      CONSTS(115) = 0.06
      CONSTS(116) = 100
      CONSTS(117) = 0.2
      CONSTS(118) = 2400
      CONSTS(119) = 100
      STATES(33) = 0.105910034847806
      CONSTS(120) = 0.06
      CONSTS(121) = 100
      CONSTS(122) = 0.2
      CONSTS(123) = 3000
      CONSTS(124) = 100
      CONSTS(125) = 100
      STATES(34) = 0.105787155209054
      CONSTS(126) = 0.06
      CONSTS(127) = 100
      CONSTS(128) = 0.2
      CONSTS(129) = 4000
      CONSTS(130) = 100
      STATES(35) = 0.105459895332545
      CONSTS(131) = 0.06
      CONSTS(132) = 100
      CONSTS(133) = 0.2
      CONSTS(134) = 1200
      CONSTS(135) = 100
      STATES(36) = 0.105317375726849
      CONSTS(136) = 0.06
      CONSTS(137) = 100
      CONSTS(138) = 0.2
      CONSTS(139) = 2000
      STATES(37) = 0.152013945600315
      CONSTS(140) = 0.06
      CONSTS(141) = 50
      CONSTS(142) = 0.2
      CONSTS(143) = 100
      CONSTS(144) = 800
      STATES(38) = 0.526188169704937
      CONSTS(145) = 0.06
      CONSTS(146) = 10
      CONSTS(147) = -0.5
      CONSTS(148) = 1
      STATES(39) = 0.0732527312969717
      CONSTS(149) = 0.06
      STATES(40) = 0.0724355014489091
      CONSTS(150) = 0.06
      STATES(41) = 0.0716607456049039
      CONSTS(151) = 0.06
      STATES(42) = 0.071501504166401
      CONSTS(152) = 0.06
      STATES(43) = 0.0714111536200381
      CONSTS(153) = 0.06
      STATES(44) = 0.0713861207574313
      CONSTS(154) = 0.06
      CONSTS(155) = 400
      CONSTS(156) = 400
      CONSTS(157) = 400
      CONSTS(158) = 400
      CONSTS(159) = 400
      CONSTS(160) = 400
      CONSTS(161) = 100
      CONSTS(162) = 100
      CONSTS(163) = 100
      CONSTS(164) = 100
      CONSTS(165) = 100
      CONSTS(166) = 100
      CONSTS(167) = 100
      STATES(45) = 0.0898390162613648
      CONSTS(168) = 0.06
      CONSTS(169) = 150
      CONSTS(170) = 4000
      CONSTS(171) = 100
      STATES(46) = 0.0898390162613648
      CONSTS(172) = 0.06
      CONSTS(173) = 150
      CONSTS(174) = 4000
      CONSTS(175) = ( CONSTS(3)* 3.14159265358979)/180.000
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) = VOI/CONSTS(4) - INT(VOI/CONSTS(4))
      ALGBRC(5) = TERNRY(ALGBRC(3).LT.0.400000,  (( (CONSTS(8) - CONSTS(9))*ALGBRC(3))/CONSTS(6))*EXP( 0.500000*(1.00000 - ALGBRC(3)/CONSTS(6) ** 2.00000))+CONSTS(9), TERNRY(ALGBRC(3).LT.0.600000,  (( (CONSTS(8) - CONSTS(9))*0.400000)/CONSTS(6))*EXP( 0.500000*(1.00000 - 0.400000/CONSTS(6) ** 2.00000))+CONSTS(9), CONSTS(9))
      ALGBRC(7) =  ALGBRC(5)*(STATES(1) - CONSTS(7))
      ALGBRC(22) =  CONSTS(28)*(STATES(7) - CONSTS(27))
      ALGBRC(25) = TERNRY(ALGBRC(7).GT.ALGBRC(22), (ALGBRC(7) - ALGBRC(22))/CONSTS(17), 0.00000)
      ALGBRC(13) =  CONSTS(14)*(STATES(3) - CONSTS(13))
      ALGBRC(19) = TERNRY(ALGBRC(13).GT.ALGBRC(7), (ALGBRC(13) - ALGBRC(7))/CONSTS(20), 0.00000)
      RATES(1) = ALGBRC(19) - ALGBRC(25)
      ALGBRC(9) = TERNRY(ALGBRC(3).LT.0.400000,  (( (CONSTS(11) - CONSTS(12))*ALGBRC(3))/CONSTS(6))*EXP( 0.500000*(1.00000 - ALGBRC(3)/CONSTS(6) ** 2.00000))+CONSTS(12), TERNRY(ALGBRC(3).LT.0.600000,  (( (CONSTS(11) - CONSTS(12))*0.400000)/CONSTS(6))*EXP( 0.500000*(1.00000 - 0.400000/CONSTS(6) ** 2.00000))+CONSTS(12), CONSTS(12))
      ALGBRC(11) =  ALGBRC(9)*(STATES(2) - CONSTS(10))
      ALGBRC(21) =  CONSTS(22)*(STATES(5) - CONSTS(21))
      ALGBRC(24) = TERNRY(ALGBRC(11).GT.ALGBRC(21), (ALGBRC(11) - ALGBRC(21))/CONSTS(18), 0.00000)
      ALGBRC(15) =  CONSTS(16)*(STATES(4) - CONSTS(15))
      ALGBRC(17) = TERNRY(ALGBRC(15).GT.ALGBRC(11), (ALGBRC(15) - ALGBRC(11))/CONSTS(19), 0.00000)
      RATES(2) = ALGBRC(17) - ALGBRC(24)
      ALGBRC(26) =  CONSTS(25)*(STATES(6) - CONSTS(24))
      ALGBRC(28) = (ALGBRC(21) - ALGBRC(26))/CONSTS(23)
      RATES(5) = ALGBRC(24) - ALGBRC(28)
      ALGBRC(30) = (ALGBRC(26) - ALGBRC(13))/CONSTS(26)
      RATES(3) = ALGBRC(30) - ALGBRC(19)
      RATES(6) = ALGBRC(28) - ALGBRC(30)
      ALGBRC(27) =  CONSTS(31)*(STATES(8) - CONSTS(30))
      ALGBRC(29) = (ALGBRC(22) - ALGBRC(27))/CONSTS(29)
      ALGBRC(31) =  CONSTS(35)*(STATES(9) - CONSTS(34))
      ALGBRC(32) = (ALGBRC(27) - ALGBRC(31))/CONSTS(32)
      ALGBRC(35) =  CONSTS(45)*(STATES(11) - CONSTS(44))
      ALGBRC(36) = (ALGBRC(27) - ALGBRC(35))/CONSTS(33)
      RATES(8) = (ALGBRC(29) - ALGBRC(32)) - ALGBRC(36)
      ALGBRC(37) = (ALGBRC(31) - ALGBRC(35))/CONSTS(37)
      ALGBRC(38) = (ALGBRC(22) - ALGBRC(35))/CONSTS(38)
      ALGBRC(33) =  CONSTS(41)*(STATES(10) - CONSTS(40))
      ALGBRC(39) = (ALGBRC(33) - ALGBRC(35))/CONSTS(43)
      ALGBRC(40) = (ALGBRC(35) - ALGBRC(15))/CONSTS(46)
      RATES(11) = (ALGBRC(37)+ALGBRC(36)+ALGBRC(38)+ALGBRC(39)) - ALGBRC(40)
      ALGBRC(34) = (ALGBRC(22) - ALGBRC(33))/CONSTS(39)
      ALGBRC(41) =  CONSTS(48)*(STATES(12) - CONSTS(47))
      ALGBRC(42) = (ALGBRC(33) - ALGBRC(41))/CONSTS(42)
      RATES(10) = (ALGBRC(34) - ALGBRC(42)) - ALGBRC(39)
      ALGBRC(43) = (ALGBRC(41) - ALGBRC(15))/CONSTS(49)
      RATES(12) = ALGBRC(42) - ALGBRC(43)
      ALGBRC(44) =  CONSTS(51)*(STATES(13) - CONSTS(50))
      ALGBRC(45) = (ALGBRC(31) - ALGBRC(44))/CONSTS(36)
      RATES(9) = (ALGBRC(32) - ALGBRC(45)) - ALGBRC(37)
      ALGBRC(46) = (ALGBRC(44) - ALGBRC(15))/CONSTS(52)
      RATES(13) = ALGBRC(45) - ALGBRC(46)
      ALGBRC(47) =  CONSTS(61)*(STATES(14) - CONSTS(53))
      ALGBRC(55) = (ALGBRC(22) - ALGBRC(47))/CONSTS(69)
      RATES(7) = (((ALGBRC(25) - ALGBRC(55)) - ALGBRC(29)) - ALGBRC(38)) - ALGBRC(34)
      ALGBRC(48) =  CONSTS(62)*(STATES(15) - CONSTS(54))
      ALGBRC(56) = (ALGBRC(47) - ALGBRC(48))/CONSTS(70)
      ALGBRC(63) =  CONSTS(79)*(STATES(22) - CONSTS(78))
      ALGBRC(64) = (ALGBRC(47) - ALGBRC(63))/CONSTS(77)
      RATES(14) = (ALGBRC(55) - ALGBRC(56)) - ALGBRC(64)
      ALGBRC(65) =  CONSTS(82)*(STATES(23) - CONSTS(81))
      ALGBRC(66) = (ALGBRC(63) - ALGBRC(65))/CONSTS(80)
      RATES(22) = ALGBRC(64) - ALGBRC(66)
      ALGBRC(49) =  CONSTS(63)*(STATES(16) - CONSTS(55))
      ALGBRC(57) = (ALGBRC(48) - ALGBRC(49))/CONSTS(71)
      ALGBRC(67) =  CONSTS(86)*(STATES(24) - CONSTS(85))
      ALGBRC(68) = (ALGBRC(48) - ALGBRC(67))/CONSTS(84)
      RATES(15) = (ALGBRC(56) - ALGBRC(57)) - ALGBRC(68)
      ALGBRC(69) =  CONSTS(89)*(STATES(25) - CONSTS(88))
      ALGBRC(70) = (ALGBRC(67) - ALGBRC(69))/CONSTS(87)
      RATES(24) = ALGBRC(68) - ALGBRC(70)
      ALGBRC(50) =  CONSTS(64)*(STATES(17) - CONSTS(56))
      ALGBRC(58) = (ALGBRC(49) - ALGBRC(50))/CONSTS(72)
      ALGBRC(71) =  CONSTS(93)*(STATES(26) - CONSTS(92))
      ALGBRC(72) = (ALGBRC(49) - ALGBRC(71))/CONSTS(91)
      RATES(16) = (ALGBRC(57) - ALGBRC(58)) - ALGBRC(72)
      ALGBRC(51) =  CONSTS(65)*(STATES(18) - CONSTS(57))
      ALGBRC(59) = (ALGBRC(50) - ALGBRC(51))/CONSTS(73)
      ALGBRC(73) =  CONSTS(97)*(STATES(27) - CONSTS(96))
      ALGBRC(74) = (ALGBRC(50) - ALGBRC(73))/CONSTS(95)
      RATES(17) = (ALGBRC(58) - ALGBRC(59)) - ALGBRC(74)
      ALGBRC(52) =  CONSTS(66)*(STATES(19) - CONSTS(58))
      ALGBRC(60) = (ALGBRC(51) - ALGBRC(52))/CONSTS(74)
      ALGBRC(75) =  CONSTS(101)*(STATES(28) - CONSTS(100))
      ALGBRC(76) = (ALGBRC(51) - ALGBRC(75))/CONSTS(99)
      RATES(18) = (ALGBRC(59) - ALGBRC(60)) - ALGBRC(76)
      ALGBRC(77) =  CONSTS(104)*(STATES(29) - CONSTS(103))
      ALGBRC(78) = (ALGBRC(75) - ALGBRC(77))/CONSTS(102)
      RATES(28) = ALGBRC(76) - ALGBRC(78)
      ALGBRC(53) =  CONSTS(67)*(STATES(20) - CONSTS(59))
      ALGBRC(61) = (ALGBRC(52) - ALGBRC(53))/CONSTS(75)
      ALGBRC(79) =  CONSTS(108)*(STATES(30) - CONSTS(107))
      ALGBRC(80) = (ALGBRC(52) - ALGBRC(79))/CONSTS(106)
      RATES(19) = (ALGBRC(60) - ALGBRC(61)) - ALGBRC(80)
      ALGBRC(81) =  CONSTS(111)*(STATES(31) - CONSTS(110))
      ALGBRC(82) = (ALGBRC(79) - ALGBRC(81))/CONSTS(109)
      RATES(30) = ALGBRC(80) - ALGBRC(82)
      ALGBRC(54) =  CONSTS(68)*(STATES(21) - CONSTS(60))
      ALGBRC(62) = (ALGBRC(53) - ALGBRC(54))/CONSTS(76)
      ALGBRC(92) =  CONSTS(137)*(STATES(36) - CONSTS(136))
      ALGBRC(93) = (ALGBRC(53) - ALGBRC(92))/CONSTS(113)
      RATES(20) = (ALGBRC(61) - ALGBRC(62)) - ALGBRC(93)
      ALGBRC(87) =  CONSTS(127)*(STATES(34) - CONSTS(126))
      ALGBRC(89) = (ALGBRC(73) - ALGBRC(87))/CONSTS(125)
      ALGBRC(90) =  CONSTS(132)*(STATES(35) - CONSTS(131))
      ALGBRC(91) = (ALGBRC(73) - ALGBRC(90))/CONSTS(130)
      ALGBRC(94) = (ALGBRC(73) - ALGBRC(92))/CONSTS(135)
      RATES(27) = ((ALGBRC(74) - ALGBRC(89)) - ALGBRC(91)) - ALGBRC(94)
      ALGBRC(83) =  CONSTS(116)*(STATES(32) - CONSTS(115))
      ALGBRC(84) = (ALGBRC(71) - ALGBRC(83))/CONSTS(114)
      ALGBRC(95) =  CONSTS(141)*(STATES(37) - CONSTS(140))
      ALGBRC(96) = (ALGBRC(83) - ALGBRC(95))/CONSTS(118)
      RATES(32) = ALGBRC(84) - ALGBRC(96)
      ALGBRC(85) =  CONSTS(121)*(STATES(33) - CONSTS(120))
      ALGBRC(86) = (ALGBRC(71) - ALGBRC(85))/CONSTS(119)
      ALGBRC(97) = (ALGBRC(85) - ALGBRC(95))/CONSTS(123)
      RATES(33) = ALGBRC(86) - ALGBRC(97)
      ALGBRC(88) = (ALGBRC(71) - ALGBRC(87))/CONSTS(124)
      ALGBRC(98) = (ALGBRC(87) - ALGBRC(95))/CONSTS(129)
      RATES(34) = (ALGBRC(88)+ALGBRC(89)) - ALGBRC(98)
      ALGBRC(99) = (ALGBRC(90) - ALGBRC(95))/CONSTS(134)
      RATES(35) = ALGBRC(91) - ALGBRC(99)
      ALGBRC(100) = (ALGBRC(92) - ALGBRC(95))/CONSTS(139)
      RATES(36) = (ALGBRC(94)+ALGBRC(93)) - ALGBRC(100)
      ALGBRC(101) =  CONSTS(146)*(STATES(38) - CONSTS(145))
      ALGBRC(102) = (ALGBRC(95) - ALGBRC(101))/CONSTS(143)
      RATES(37) = (ALGBRC(96)+ALGBRC(97)+ALGBRC(98)+ALGBRC(99)+ALGBRC(100)) - ALGBRC(102)
      ALGBRC(103) = (ALGBRC(71) - ALGBRC(101))/CONSTS(144)
      RATES(26) = (((ALGBRC(72) - ALGBRC(103)) - ALGBRC(84)) - ALGBRC(86)) - ALGBRC(88)
      ALGBRC(104) =  CONSTS(155)*(STATES(39) - CONSTS(149))
      ALGBRC(105) = (ALGBRC(65) - ALGBRC(104))/CONSTS(83)
      RATES(23) = ALGBRC(66) - ALGBRC(105)
      ALGBRC(106) =  CONSTS(156)*(STATES(40) - CONSTS(150))
      ALGBRC(107) = (ALGBRC(69) - ALGBRC(106))/CONSTS(90)
      RATES(25) = ALGBRC(70) - ALGBRC(107)
      ALGBRC(108) =  CONSTS(157)*(STATES(41) - CONSTS(151))
      ALGBRC(109) = (ALGBRC(101) - ALGBRC(108))/CONSTS(148)
      RATES(38) = (ALGBRC(103)+ALGBRC(102)) - ALGBRC(109)
      ALGBRC(110) =  CONSTS(158)*(STATES(42) - CONSTS(152))
      ALGBRC(111) = (ALGBRC(77) - ALGBRC(110))/CONSTS(105)
      RATES(29) = ALGBRC(78) - ALGBRC(111)
      ALGBRC(112) =  CONSTS(159)*(STATES(43) - CONSTS(153))
      ALGBRC(113) = (ALGBRC(81) - ALGBRC(112))/CONSTS(112)
      RATES(31) = ALGBRC(82) - ALGBRC(113)
      ALGBRC(115) = (ALGBRC(104) - ALGBRC(15))/CONSTS(161)
      RATES(4) = (ALGBRC(40)+ALGBRC(43)+ALGBRC(46)+ALGBRC(115)) - ALGBRC(17)
      ALGBRC(116) = (ALGBRC(106) - ALGBRC(104))/CONSTS(162)
      RATES(39) = (ALGBRC(116)+ALGBRC(105)) - ALGBRC(115)
      ALGBRC(117) = (ALGBRC(108) - ALGBRC(106))/CONSTS(163)
      RATES(40) = (ALGBRC(117)+ALGBRC(107)) - ALGBRC(116)
      ALGBRC(118) = (ALGBRC(110) - ALGBRC(108))/CONSTS(164)
      RATES(41) = (ALGBRC(118)+ALGBRC(109)) - ALGBRC(117)
      ALGBRC(119) = (ALGBRC(112) - ALGBRC(110))/CONSTS(165)
      RATES(42) = (ALGBRC(119)+ALGBRC(111)) - ALGBRC(118)
      ALGBRC(114) =  CONSTS(160)*(STATES(44) - CONSTS(154))
      ALGBRC(120) = (ALGBRC(114) - ALGBRC(112))/CONSTS(166)
      RATES(43) = (ALGBRC(120)+ALGBRC(113)) - ALGBRC(119)
      ALGBRC(121) =  CONSTS(169)*(STATES(45) - CONSTS(168))
      ALGBRC(122) = (ALGBRC(54) - ALGBRC(121))/CONSTS(167)
      ALGBRC(123) = (ALGBRC(121) - ALGBRC(114))/CONSTS(170)
      RATES(45) = ALGBRC(122) - ALGBRC(123)
      ALGBRC(124) =  CONSTS(173)*(STATES(46) - CONSTS(172))
      ALGBRC(125) = (ALGBRC(54) - ALGBRC(124))/CONSTS(171)
      RATES(21) = (ALGBRC(62) - ALGBRC(122)) - ALGBRC(125)
      ALGBRC(126) = (ALGBRC(124) - ALGBRC(114))/CONSTS(174)
      RATES(44) = (ALGBRC(123)+ALGBRC(126)) - ALGBRC(120)
      RATES(46) = ALGBRC(125) - ALGBRC(126)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) = VOI/CONSTS(4) - INT(VOI/CONSTS(4))
      ALGBRC(5) = TERNRY(ALGBRC(3).LT.0.400000,  (( (CONSTS(8) - CONSTS(9))*ALGBRC(3))/CONSTS(6))*EXP( 0.500000*(1.00000 - ALGBRC(3)/CONSTS(6) ** 2.00000))+CONSTS(9), TERNRY(ALGBRC(3).LT.0.600000,  (( (CONSTS(8) - CONSTS(9))*0.400000)/CONSTS(6))*EXP( 0.500000*(1.00000 - 0.400000/CONSTS(6) ** 2.00000))+CONSTS(9), CONSTS(9))
      ALGBRC(7) =  ALGBRC(5)*(STATES(1) - CONSTS(7))
      ALGBRC(22) =  CONSTS(28)*(STATES(7) - CONSTS(27))
      ALGBRC(25) = TERNRY(ALGBRC(7).GT.ALGBRC(22), (ALGBRC(7) - ALGBRC(22))/CONSTS(17), 0.00000)
      ALGBRC(13) =  CONSTS(14)*(STATES(3) - CONSTS(13))
      ALGBRC(19) = TERNRY(ALGBRC(13).GT.ALGBRC(7), (ALGBRC(13) - ALGBRC(7))/CONSTS(20), 0.00000)
      ALGBRC(9) = TERNRY(ALGBRC(3).LT.0.400000,  (( (CONSTS(11) - CONSTS(12))*ALGBRC(3))/CONSTS(6))*EXP( 0.500000*(1.00000 - ALGBRC(3)/CONSTS(6) ** 2.00000))+CONSTS(12), TERNRY(ALGBRC(3).LT.0.600000,  (( (CONSTS(11) - CONSTS(12))*0.400000)/CONSTS(6))*EXP( 0.500000*(1.00000 - 0.400000/CONSTS(6) ** 2.00000))+CONSTS(12), CONSTS(12))
      ALGBRC(11) =  ALGBRC(9)*(STATES(2) - CONSTS(10))
      ALGBRC(21) =  CONSTS(22)*(STATES(5) - CONSTS(21))
      ALGBRC(24) = TERNRY(ALGBRC(11).GT.ALGBRC(21), (ALGBRC(11) - ALGBRC(21))/CONSTS(18), 0.00000)
      ALGBRC(15) =  CONSTS(16)*(STATES(4) - CONSTS(15))
      ALGBRC(17) = TERNRY(ALGBRC(15).GT.ALGBRC(11), (ALGBRC(15) - ALGBRC(11))/CONSTS(19), 0.00000)
      ALGBRC(26) =  CONSTS(25)*(STATES(6) - CONSTS(24))
      ALGBRC(28) = (ALGBRC(21) - ALGBRC(26))/CONSTS(23)
      ALGBRC(30) = (ALGBRC(26) - ALGBRC(13))/CONSTS(26)
      ALGBRC(27) =  CONSTS(31)*(STATES(8) - CONSTS(30))
      ALGBRC(29) = (ALGBRC(22) - ALGBRC(27))/CONSTS(29)
      ALGBRC(31) =  CONSTS(35)*(STATES(9) - CONSTS(34))
      ALGBRC(32) = (ALGBRC(27) - ALGBRC(31))/CONSTS(32)
      ALGBRC(35) =  CONSTS(45)*(STATES(11) - CONSTS(44))
      ALGBRC(36) = (ALGBRC(27) - ALGBRC(35))/CONSTS(33)
      ALGBRC(37) = (ALGBRC(31) - ALGBRC(35))/CONSTS(37)
      ALGBRC(38) = (ALGBRC(22) - ALGBRC(35))/CONSTS(38)
      ALGBRC(33) =  CONSTS(41)*(STATES(10) - CONSTS(40))
      ALGBRC(39) = (ALGBRC(33) - ALGBRC(35))/CONSTS(43)
      ALGBRC(40) = (ALGBRC(35) - ALGBRC(15))/CONSTS(46)
      ALGBRC(34) = (ALGBRC(22) - ALGBRC(33))/CONSTS(39)
      ALGBRC(41) =  CONSTS(48)*(STATES(12) - CONSTS(47))
      ALGBRC(42) = (ALGBRC(33) - ALGBRC(41))/CONSTS(42)
      ALGBRC(43) = (ALGBRC(41) - ALGBRC(15))/CONSTS(49)
      ALGBRC(44) =  CONSTS(51)*(STATES(13) - CONSTS(50))
      ALGBRC(45) = (ALGBRC(31) - ALGBRC(44))/CONSTS(36)
      ALGBRC(46) = (ALGBRC(44) - ALGBRC(15))/CONSTS(52)
      ALGBRC(47) =  CONSTS(61)*(STATES(14) - CONSTS(53))
      ALGBRC(55) = (ALGBRC(22) - ALGBRC(47))/CONSTS(69)
      ALGBRC(48) =  CONSTS(62)*(STATES(15) - CONSTS(54))
      ALGBRC(56) = (ALGBRC(47) - ALGBRC(48))/CONSTS(70)
      ALGBRC(63) =  CONSTS(79)*(STATES(22) - CONSTS(78))
      ALGBRC(64) = (ALGBRC(47) - ALGBRC(63))/CONSTS(77)
      ALGBRC(65) =  CONSTS(82)*(STATES(23) - CONSTS(81))
      ALGBRC(66) = (ALGBRC(63) - ALGBRC(65))/CONSTS(80)
      ALGBRC(49) =  CONSTS(63)*(STATES(16) - CONSTS(55))
      ALGBRC(57) = (ALGBRC(48) - ALGBRC(49))/CONSTS(71)
      ALGBRC(67) =  CONSTS(86)*(STATES(24) - CONSTS(85))
      ALGBRC(68) = (ALGBRC(48) - ALGBRC(67))/CONSTS(84)
      ALGBRC(69) =  CONSTS(89)*(STATES(25) - CONSTS(88))
      ALGBRC(70) = (ALGBRC(67) - ALGBRC(69))/CONSTS(87)
      ALGBRC(50) =  CONSTS(64)*(STATES(17) - CONSTS(56))
      ALGBRC(58) = (ALGBRC(49) - ALGBRC(50))/CONSTS(72)
      ALGBRC(71) =  CONSTS(93)*(STATES(26) - CONSTS(92))
      ALGBRC(72) = (ALGBRC(49) - ALGBRC(71))/CONSTS(91)
      ALGBRC(51) =  CONSTS(65)*(STATES(18) - CONSTS(57))
      ALGBRC(59) = (ALGBRC(50) - ALGBRC(51))/CONSTS(73)
      ALGBRC(73) =  CONSTS(97)*(STATES(27) - CONSTS(96))
      ALGBRC(74) = (ALGBRC(50) - ALGBRC(73))/CONSTS(95)
      ALGBRC(52) =  CONSTS(66)*(STATES(19) - CONSTS(58))
      ALGBRC(60) = (ALGBRC(51) - ALGBRC(52))/CONSTS(74)
      ALGBRC(75) =  CONSTS(101)*(STATES(28) - CONSTS(100))
      ALGBRC(76) = (ALGBRC(51) - ALGBRC(75))/CONSTS(99)
      ALGBRC(77) =  CONSTS(104)*(STATES(29) - CONSTS(103))
      ALGBRC(78) = (ALGBRC(75) - ALGBRC(77))/CONSTS(102)
      ALGBRC(53) =  CONSTS(67)*(STATES(20) - CONSTS(59))
      ALGBRC(61) = (ALGBRC(52) - ALGBRC(53))/CONSTS(75)
      ALGBRC(79) =  CONSTS(108)*(STATES(30) - CONSTS(107))
      ALGBRC(80) = (ALGBRC(52) - ALGBRC(79))/CONSTS(106)
      ALGBRC(81) =  CONSTS(111)*(STATES(31) - CONSTS(110))
      ALGBRC(82) = (ALGBRC(79) - ALGBRC(81))/CONSTS(109)
      ALGBRC(54) =  CONSTS(68)*(STATES(21) - CONSTS(60))
      ALGBRC(62) = (ALGBRC(53) - ALGBRC(54))/CONSTS(76)
      ALGBRC(92) =  CONSTS(137)*(STATES(36) - CONSTS(136))
      ALGBRC(93) = (ALGBRC(53) - ALGBRC(92))/CONSTS(113)
      ALGBRC(87) =  CONSTS(127)*(STATES(34) - CONSTS(126))
      ALGBRC(89) = (ALGBRC(73) - ALGBRC(87))/CONSTS(125)
      ALGBRC(90) =  CONSTS(132)*(STATES(35) - CONSTS(131))
      ALGBRC(91) = (ALGBRC(73) - ALGBRC(90))/CONSTS(130)
      ALGBRC(94) = (ALGBRC(73) - ALGBRC(92))/CONSTS(135)
      ALGBRC(83) =  CONSTS(116)*(STATES(32) - CONSTS(115))
      ALGBRC(84) = (ALGBRC(71) - ALGBRC(83))/CONSTS(114)
      ALGBRC(95) =  CONSTS(141)*(STATES(37) - CONSTS(140))
      ALGBRC(96) = (ALGBRC(83) - ALGBRC(95))/CONSTS(118)
      ALGBRC(85) =  CONSTS(121)*(STATES(33) - CONSTS(120))
      ALGBRC(86) = (ALGBRC(71) - ALGBRC(85))/CONSTS(119)
      ALGBRC(97) = (ALGBRC(85) - ALGBRC(95))/CONSTS(123)
      ALGBRC(88) = (ALGBRC(71) - ALGBRC(87))/CONSTS(124)
      ALGBRC(98) = (ALGBRC(87) - ALGBRC(95))/CONSTS(129)
      ALGBRC(99) = (ALGBRC(90) - ALGBRC(95))/CONSTS(134)
      ALGBRC(100) = (ALGBRC(92) - ALGBRC(95))/CONSTS(139)
      ALGBRC(101) =  CONSTS(146)*(STATES(38) - CONSTS(145))
      ALGBRC(102) = (ALGBRC(95) - ALGBRC(101))/CONSTS(143)
      ALGBRC(103) = (ALGBRC(71) - ALGBRC(101))/CONSTS(144)
      ALGBRC(104) =  CONSTS(155)*(STATES(39) - CONSTS(149))
      ALGBRC(105) = (ALGBRC(65) - ALGBRC(104))/CONSTS(83)
      ALGBRC(106) =  CONSTS(156)*(STATES(40) - CONSTS(150))
      ALGBRC(107) = (ALGBRC(69) - ALGBRC(106))/CONSTS(90)
      ALGBRC(108) =  CONSTS(157)*(STATES(41) - CONSTS(151))
      ALGBRC(109) = (ALGBRC(101) - ALGBRC(108))/CONSTS(148)
      ALGBRC(110) =  CONSTS(158)*(STATES(42) - CONSTS(152))
      ALGBRC(111) = (ALGBRC(77) - ALGBRC(110))/CONSTS(105)
      ALGBRC(112) =  CONSTS(159)*(STATES(43) - CONSTS(153))
      ALGBRC(113) = (ALGBRC(81) - ALGBRC(112))/CONSTS(112)
      ALGBRC(115) = (ALGBRC(104) - ALGBRC(15))/CONSTS(161)
      ALGBRC(116) = (ALGBRC(106) - ALGBRC(104))/CONSTS(162)
      ALGBRC(117) = (ALGBRC(108) - ALGBRC(106))/CONSTS(163)
      ALGBRC(118) = (ALGBRC(110) - ALGBRC(108))/CONSTS(164)
      ALGBRC(119) = (ALGBRC(112) - ALGBRC(110))/CONSTS(165)
      ALGBRC(114) =  CONSTS(160)*(STATES(44) - CONSTS(154))
      ALGBRC(120) = (ALGBRC(114) - ALGBRC(112))/CONSTS(166)
      ALGBRC(121) =  CONSTS(169)*(STATES(45) - CONSTS(168))
      ALGBRC(122) = (ALGBRC(54) - ALGBRC(121))/CONSTS(167)
      ALGBRC(123) = (ALGBRC(121) - ALGBRC(114))/CONSTS(170)
      ALGBRC(124) =  CONSTS(173)*(STATES(46) - CONSTS(172))
      ALGBRC(125) = (ALGBRC(54) - ALGBRC(124))/CONSTS(171)
      ALGBRC(126) = (ALGBRC(124) - ALGBRC(114))/CONSTS(174)
      ALGBRC(1) = TERNRY(VOI.GT.100.000.AND.VOI.LT.101.000, 1.00000, 0.00000)
      ALGBRC(2) = (VOI/CONSTS(4) - CONSTS(5)) - INT((VOI -  CONSTS(5)*CONSTS(4))/CONSTS(4))
      ALGBRC(4) = STATES(4)+STATES(2)+STATES(3)+STATES(1)
      ALGBRC(6) = STATES(6)+STATES(5)
      ALGBRC(8) = STATES(8)+STATES(9)+STATES(10)+STATES(11)
      ALGBRC(10) = STATES(22)+STATES(23)+STATES(24)+STATES(25)
      ALGBRC(12) = STATES(7)+STATES(14)+STATES(15)+STATES(16)+STATES(17)+STATES(18)+STATES(19)+STATES(20)+STATES(21)
      ALGBRC(14) = STATES(28)+STATES(29)+STATES(30)+STATES(31)
      ALGBRC(16) = STATES(26)+STATES(27)+STATES(32)+STATES(33)+STATES(34)+STATES(35)+STATES(36)+STATES(37)+STATES(38)
      ALGBRC(18) = STATES(39)+STATES(40)+STATES(41)+STATES(42)+STATES(43)+STATES(44)
      ALGBRC(20) = STATES(12)+STATES(13)+STATES(45)+STATES(46)
      ALGBRC(23) = ALGBRC(4)+ALGBRC(6)+(ALGBRC(8))+ALGBRC(10)+ALGBRC(12)+ALGBRC(14)+ALGBRC(16)+ALGBRC(18)+ALGBRC(20)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END