Generated Code
The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
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