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 62 entries in the algebraic variable array.
C There are a total of 5 entries in each of the rate and state variable arrays.
C There are a total of 584 entries in the constant variable array.
C
C
C CONSTS(1) is ANU in component input_values (dimensionless).
C CONSTS(2) is ADHMV in component input_values (dimensionless).
C CONSTS(3) is PAMK in component input_values (dimensionless).
C CONSTS(4) is ATRRFB in component input_values (dimensionless).
C CONSTS(5) is MYOGRS in component input_values (dimensionless).
C CONSTS(6) is ARM in component input_values (dimensionless).
C CONSTS(7) is VIM in component input_values (dimensionless).
C CONSTS(8) is AVE in component input_values (dimensionless).
C CONSTS(9) is ANUVN in component input_values (dimensionless).
C CONSTS(10) is PC in component input_values (mmHg).
C CONSTS(11) is RBF in component input_values (L_per_minute).
C CONSTS(12) is VVR in component input_values (litre).
C CONSTS(13) is VV6 in component input_values (litre).
C CONSTS(14) is VV7 in component input_values (litre).
C CONSTS(15) is ATRVFB in component input_values (litre).
C CONSTS(16) is VP in component input_values (litre).
C CONSTS(17) is VRC in component input_values (litre).
C CONSTS(18) is HPR in component input_values (dimensionless).
C CONSTS(19) is HMD in component input_values (dimensionless).
C CONSTS(20) is AUH in component input_values (dimensionless).
C CONSTS(21) is OSA in component input_values (dimensionless).
C CONSTS(22) is HPL in component input_values (dimensionless).
C CONSTS(23) is AMM in component input_values (dimensionless).
C CONSTS(24) is AUM in component input_values (dimensionless).
C CONSTS(25) is AU in component input_values (dimensionless).
C VOI is time in component environment (minute).
C CONSTS(26) is A1K in component parameter_values (minute).
C CONSTS(27) is A2K in component parameter_values (minute).
C CONSTS(28) is A3K in component parameter_values (minute).
C CONSTS(29) is A4K in component parameter_values (minute).
C CONSTS(30) is A4K2 in component parameter_values (minute).
C CONSTS(31) is AARK in component parameter_values (mmHg_minute_per_L).
C CONSTS(32) is AARLL in component parameter_values (mmHg_minute_per_L).
C CONSTS(33) is ADHINF in component parameter_values (dimensionless).
C CONSTS(34) is ADHKLL in component parameter_values (dimensionless).
C CONSTS(35) is ADHKUL in component parameter_values (dimensionless).
C CONSTS(36) is ADHPAM in component parameter_values (per_mmHg2).
C CONSTS(37) is ADHPUL in component parameter_values (mmHg).
C CONSTS(38) is ADHTC in component parameter_values (minute).
C CONSTS(39) is ADHVLL in component parameter_values (dimensionless).
C CONSTS(40) is ADHVUL in component parameter_values (dimensionless).
C CONSTS(41) is AH10 in component parameter_values (dimensionless).
C CONSTS(42) is AH11 in component parameter_values (minute).
C CONSTS(43) is AH9 in component parameter_values (per_mmHg).
C CONSTS(44) is AHMNAR in component parameter_values (dimensionless).
C CONSTS(45) is AHTHM in component parameter_values (dimensionless).
C CONSTS(46) is ALCLK in component parameter_values (dimensionless).
C CONSTS(47) is ALDINF in component parameter_values (dimensionless).
C CONSTS(48) is ALDKNS in component parameter_values (dimensionless).
C CONSTS(49) is ALDMM in component parameter_values (dimensionless).
C CONSTS(50) is AM1LL in component parameter_values (dimensionless).
C CONSTS(51) is AM1UL in component parameter_values (dimensionless).
C CONSTS(52) is AMCSNS in component parameter_values (dimensionless).
C CONSTS(53) is AMKM in component parameter_values (dimensionless).
C CONSTS(54) is AMKMUL in component parameter_values (dimensionless).
C CONSTS(55) is AMM1LM in component parameter_values (dimensionless).
C CONSTS(56) is AMM4 in component parameter_values (dimensionless).
C CONSTS(57) is AMNALL in component parameter_values (dimensionless).
C CONSTS(58) is AMNAM in component parameter_values (dimensionless).
C CONSTS(59) is AMNAUL in component parameter_values (dimensionless).
C CONSTS(60) is AMT in component parameter_values (minute).
C CONSTS(61) is ANCSNS in component parameter_values (dimensionless).
C CONSTS(62) is ANGINF in component parameter_values (dimensionless).
C CONSTS(63) is ANGKNS in component parameter_values (dimensionless).
C CONSTS(64) is ANMALD in component parameter_values (dimensionless).
C CONSTS(65) is ANMAM in component parameter_values (dimensionless).
C CONSTS(66) is ANMARL in component parameter_values (dimensionless).
C CONSTS(67) is ANMEM in component parameter_values (dimensionless).
C CONSTS(68) is ANMKEL in component parameter_values (dimensionless).
C CONSTS(69) is ANMKEM in component parameter_values (dimensionless).
C CONSTS(70) is ANMLL in component parameter_values (dimensionless).
C CONSTS(71) is ANMNAM in component parameter_values (dimensionless).
C CONSTS(72) is ANMSLT in component parameter_values (dimensionless).
C CONSTS(73) is ANMTM in component parameter_values (dimensionless).
C CONSTS(74) is ANMUL in component parameter_values (dimensionless).
C CONSTS(75) is ANPTC in component parameter_values (minute).
C CONSTS(76) is ANPXAF in component parameter_values (mmHg_minute_per_L).
C CONSTS(77) is ANPXUL in component parameter_values (dimensionless).
C CONSTS(78) is ANT in component parameter_values (minute).
C CONSTS(79) is ANULL in component parameter_values (dimensionless).
C CONSTS(80) is ANUM in component parameter_values (dimensionless).
C CONSTS(81) is ANUVM in component parameter_values (dimensionless).
C CONSTS(82) is ANV in component parameter_values (minute).
C CONSTS(83) is ANXM in component parameter_values (dimensionless).
C CONSTS(84) is ANY in component parameter_values (litre).
C CONSTS(85) is AR1LM in component parameter_values (dimensionless).
C CONSTS(86) is AR2LM in component parameter_values (dimensionless).
C CONSTS(87) is ARF in component parameter_values (dimensionless).
C CONSTS(88) is AUC1 in component parameter_values (dimensionless).
C CONSTS(89) is AUDMP in component parameter_values (minute).
C CONSTS(90) is AUK in component parameter_values (per_minute).
C CONSTS(91) is AUL in component parameter_values (litre).
C CONSTS(92) is AULPM in component parameter_values (dimensionless).
C CONSTS(93) is AUM1 in component parameter_values (dimensionless).
C CONSTS(94) is AUM2 in component parameter_values (dimensionless).
C CONSTS(95) is AUMAX in component parameter_values (dimensionless).
C CONSTS(96) is AUMIN in component parameter_values (dimensionless).
C CONSTS(97) is AUMK1 in component parameter_values (dimensionless).
C CONSTS(98) is AUN1 in component parameter_values (dimensionless).
C CONSTS(99) is AURG in component parameter_values (dimensionless).
C CONSTS(100) is AUS in component parameter_values (dimensionless).
C CONSTS(101) is AUSLP in component parameter_values (dimensionless).
C CONSTS(102) is AUTO in component parameter_values (dimensionless).
C CONSTS(103) is AUTOGN in component parameter_values (dimensionless).
C CONSTS(104) is AUTOSN in component parameter_values (dimensionless).
C CONSTS(105) is AUX in component parameter_values (dimensionless).
C CONSTS(106) is BAROTC in component parameter_values (minute).
C CONSTS(107) is CFC in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(108) is CKEEX in component parameter_values (dimensionless).
C CONSTS(109) is CMPTSS in component parameter_values (dimensionless).
C CONSTS(110) is CN2 in component parameter_values (per_mmHg).
C CONSTS(111) is CN7 in component parameter_values (dimensionless).
C CONSTS(112) is CNR in component parameter_values (monovalent_mEq_per_litre).
C CONSTS(113) is CPF in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(114) is CPK in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(115) is CPR in component parameter_values (gram_per_L).
C CONSTS(116) is CV in component parameter_values (L_per_mmHg).
C CONSTS(117) is CRRFLX in component parameter_values (mmHg).
C CONSTS(118) is DIURET in component parameter_values (dimensionless).
C CONSTS(119) is DHDTR in component parameter_values (per_mmHg_per_minute).
C CONSTS(120) is DTNAR in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(121) is DTNARL in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(122) is EARK in component parameter_values (mmHg_minute_per_L).
C CONSTS(123) is EARLL in component parameter_values (mmHg_minute_per_L).
C CONSTS(124) is EFAFR in component parameter_values (dimensionless).
C CONSTS(125) is EXC in component parameter_values (dimensionless).
C CONSTS(126) is EXCML in component parameter_values (dimensionless).
C CONSTS(127) is EXCXP in component parameter_values (dimensionless).
C CONSTS(128) is EXCXP2 in component parameter_values (dimensionless).
C CONSTS(129) is EXE in component parameter_values (mmHg).
C CONSTS(130) is GCOPF in component parameter_values (per_mmHg).
C CONSTS(131) is GFLC in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(132) is GFNDMP in component parameter_values (dimensionless).
C CONSTS(133) is GFNLL in component parameter_values (L_per_minute).
C CONSTS(134) is GPPD in component parameter_values (dimensionless).
C CONSTS(135) is HM6 in component parameter_values (mmHg).
C CONSTS(136) is HM8 in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(137) is HTAUML in component parameter_values (dimensionless).
C CONSTS(138) is KORGN in component parameter_values (dimensionless).
C CONSTS(139) is KORNER in component parameter_values (dimensionless).
C CONSTS(140) is KORREN in component parameter_values (dimensionless).
C CONSTS(141) is KORTC in component parameter_values (minute).
C CONSTS(142) is LPDE in component parameter_values (dimensionless).
C CONSTS(143) is LPK in component parameter_values (L_per_minute).
C CONSTS(144) is LPPR in component parameter_values (gram_per_minute).
C CONSTS(145) is MDFL1 in component parameter_values (dimensionless).
C CONSTS(146) is MDFLKM in component parameter_values (monovalent_mEq_per_litre).
C CONSTS(147) is MDFLW2 in component parameter_values (dimensionless).
C CONSTS(148) is MDFLWX in component parameter_values (dimensionless).
C CONSTS(149) is MDMP in component parameter_values (dimensionless).
C CONSTS(150) is O2A in component parameter_values (dimensionless).
C CONSTS(151) is O2CHMO in component parameter_values (per_mmHg).
C CONSTS(152) is OMM in component parameter_values (mL_per_minute).
C CONSTS(153) is PA4 in component parameter_values (dimensionless).
C CONSTS(154) is PAEX in component parameter_values (dimensionless).
C CONSTS(155) is PCE in component parameter_values (dimensionless).
C CONSTS(156) is PCR in component parameter_values (mmHg).
C CONSTS(157) is PR1LL in component parameter_values (mmHg).
C CONSTS(158) is PGHF in component parameter_values (L_mmHg_per_gram).
C CONSTS(159) is PK2 in component parameter_values (mmHg_per_mL).
C CONSTS(160) is PL2 in component parameter_values (L_mL_per_minute_per_mmHg).
C CONSTS(161) is PLDF in component parameter_values (mmHg).
C CONSTS(162) is PM5 in component parameter_values (per_mmHg).
C CONSTS(163) is PO2ADV in component parameter_values (dimensionless).
C CONSTS(164) is POK in component parameter_values (dimensionless).
C CONSTS(165) is POM in component parameter_values (dimensionless).
C CONSTS(166) is POM2 in component parameter_values (dimensionless).
C CONSTS(167) is PON in component parameter_values (dimensionless).
C CONSTS(168) is POR in component parameter_values (mmHg).
C CONSTS(169) is POZ in component parameter_values (dimensionless).
C CONSTS(170) is PXTP in component parameter_values (mmHg).
C CONSTS(171) is QAOLM in component parameter_values (dimensionless).
C CONSTS(172) is QRF in component parameter_values (L_per_minute).
C CONSTS(173) is RABSC in component parameter_values (per_mmHg).
C CONSTS(174) is RAM in component parameter_values (mmHg_minute_per_L).
C CONSTS(175) is RAPRSP in component parameter_values (mmHg).
C CONSTS(176) is RAR in component parameter_values (mmHg_minute_per_L).
C CONSTS(177) is RCDFDP in component parameter_values (minute).
C CONSTS(178) is RCDFPC in component parameter_values (dimensionless).
C CONSTS(179) is REK in component parameter_values (dimensionless).
C CONSTS(180) is RFABKM in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(181) is RFABDM in component parameter_values (dimensionless).
C CONSTS(182) is RFABDP in component parameter_values (dimensionless).
C CONSTS(183) is RFCDFT in component parameter_values (dimensionless).
C CONSTS(184) is RKC in component parameter_values (per_minute).
C CONSTS(185) is RNAGTC in component parameter_values (minute).
C CONSTS(186) is RNAUAD in component parameter_values (per_minute).
C CONSTS(187) is RPS in component parameter_values (dimensionless).
C CONSTS(188) is RTPPR in component parameter_values (dimensionless).
C CONSTS(189) is RTPPRS in component parameter_values (mmHg).
C CONSTS(190) is RTSPRS in component parameter_values (mmHg).
C CONSTS(191) is RVRS in component parameter_values (mmHg_minute_per_L).
C CONSTS(192) is RVSM in component parameter_values (mmHg_minute_per_L).
C CONSTS(193) is SR in component parameter_values (dimensionless).
C CONSTS(194) is SR2 in component parameter_values (dimensionless).
C CONSTS(195) is SRK in component parameter_values (minute).
C CONSTS(196) is SRK2 in component parameter_values (minute).
C CONSTS(197) is STA in component parameter_values (dimensionless).
C CONSTS(198) is TENSGN in component parameter_values (dimensionless).
C CONSTS(199) is TENSTC in component parameter_values (minute).
C CONSTS(200) is TIMETR in component parameter_values (minute).
C CONSTS(201) is TRNSFS in component parameter_values (dimensionless).
C CONSTS(202) is TRPL in component parameter_values (L_per_minute).
C CONSTS(203) is TSSLML in component parameter_values (dimensionless).
C CONSTS(204) is TSSLTC in component parameter_values (per_minute).
C CONSTS(205) is TVDDL in component parameter_values (minute).
C CONSTS(206) is U in component parameter_values (dimensionless).
C CONSTS(207) is VID1 in component parameter_values (L_per_minute).
C CONSTS(208) is VIDML in component parameter_values (litre2_per_monovalent_mEq_per_minute).
C CONSTS(209) is VNTSTM in component parameter_values (dimensionless).
C CONSTS(210) is VP1 in component parameter_values (dimensionless).
C CONSTS(211) is VPTISS in component parameter_values (litre).
C CONSTS(212) is VTSF in component parameter_values (litre).
C CONSTS(213) is VV9 in component parameter_values (litre).
C CONSTS(214) is X in component parameter_values (dimensionless).
C CONSTS(215) is Z in component parameter_values (dimensionless).
C CONSTS(216) is Z4 in component parameter_values (dimensionless).
C CONSTS(217) is Z5 in component parameter_values (dimensionless).
C CONSTS(218) is Z6 in component parameter_values (dimensionless).
C CONSTS(219) is Z7 in component parameter_values (dimensionless).
C CONSTS(220) is Z10 in component parameter_values (mmHg).
C CONSTS(221) is Z11 in component parameter_values (per_mmHg2).
C CONSTS(222) is Z12 in component parameter_values (dimensionless).
C CONSTS(223) is Z13 in component parameter_values (dimensionless).
C CONSTS(224) is Z14 in component parameter_values (dimensionless).
C CONSTS(225) is Z16 in component parameter_values (dimensionless).
C CONSTS(226) is Z18 in component parameter_values (dimensionless).
C CONSTS(227) is Z19 in component parameter_values (dimensionless).
C CONSTS(228) is KID in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(229) is URFORM in component parameter_values (mOsm_per_minute).
C CONSTS(230) is AUV in component parameter_values (dimensionless).
C CONSTS(231) is AUY in component parameter_values (dimensionless).
C CONSTS(232) is O2M in component parameter_values (mL_per_minute).
C CONSTS(233) is NID in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(234) is ATRFBM in component parameter_values (dimensionless).
C CONSTS(235) is GBL in component parameter_values (mmHg).
C CONSTS(236) is RNAUGN in component parameter_values (minute_per_L).
C CONSTS(237) is RNAULL in component parameter_values (dimensionless).
C CONSTS(238) is RNAUUL in component parameter_values (dimensionless).
C CONSTS(239) is HSR in component parameter_values (dimensionless).
C CONSTS(240) is HSL in component parameter_values (dimensionless).
C CONSTS(241) is HMK in component parameter_values (dimensionless).
C CONSTS(242) is PO2AMB in component parameter_values (mmHg).
C CONSTS(243) is HKM in component parameter_values (dimensionless).
C CONSTS(244) is DR in component parameter_values (L_per_minute).
C CONSTS(245) is ANPKNS in component parameter_values (dimensionless).
C CONSTS(246) is RMULT1 in component parameter_values (dimensionless).
C CONSTS(247) is GLPCA in component parameter_values (mmHg).
C CONSTS(248) is RFN in component parameter_values (L_per_minute).
C CONSTS(249) is RFABX in component parameter_values (dimensionless).
C CONSTS(250) is ANPINF in component parameter_values (dimensionless).
C CONSTS(251) is FIS in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(252) is CPPD in component parameter_values (gram_per_L).
C CONSTS(253) is TRRBC in component parameter_values (L_per_minute).
C CONSTS(254) is HYL in component parameter_values (gram).
C CONSTS(255) is ATRVM in component parameter_values (litre).
C CONSTS(256) is ADHC_init in component state_variables (dimensionless).
C CONSTS(257) is AHY_init in component state_variables (dimensionless).
C CONSTS(258) is AMC_init in component state_variables (dimensionless).
C CONSTS(259) is AMM1T_init in component state_variables (dimensionless).
C CONSTS(260) is AMM2_init in component state_variables (dimensionless).
C CONSTS(261) is ANC_init in component state_variables (dimensionless).
C CONSTS(262) is ANPC_init in component state_variables (dimensionless).
C CONSTS(263) is ANX1_init in component state_variables (dimensionless).
C CONSTS(264) is AR1T_init in component state_variables (dimensionless).
C CONSTS(265) is AR2T_init in component state_variables (dimensionless).
C CONSTS(266) is AR3T_init in component state_variables (dimensionless).
C CONSTS(267) is AU1_init in component state_variables (dimensionless).
C CONSTS(268) is AU4_init in component state_variables (dimensionless).
C CONSTS(269) is AU6_init in component state_variables (dimensionless).
C CONSTS(270) is DTKA_init in component state_variables (monovalent_mEq_per_minute).
C CONSTS(271) is HMD1_init in component state_variables (dimensionless).
C CONSTS(272) is HPL_init in component state_variables (dimensionless).
C CONSTS(273) is HPR_init in component state_variables (dimensionless).
C CONSTS(274) is KTOT_init in component state_variables (monovalent_mEq).
C CONSTS(275) is NAE_init in component state_variables (monovalent_mEq).
C CONSTS(276) is O2VAD1_init in component state_variables (dimensionless).
C CONSTS(277) is OSV_init in component state_variables (dimensionless).
C CONSTS(278) is OVA_init in component state_variables (mL_per_L).
C CONSTS(279) is OVS_init in component state_variables (dimensionless).
C CONSTS(280) is PAR1_init in component state_variables (mmHg).
C CONSTS(281) is PC_init in component state_variables (mmHg).
C CONSTS(282) is PLUR_init in component state_variables (mOsm).
C CONSTS(283) is PPR1_init in component state_variables (gram).
C CONSTS(284) is PRP_init in component state_variables (gram).
C CONSTS(285) is QO2T_init in component state_variables (mL).
C CONSTS(286) is QOM1_init in component state_variables (mL).
C CONSTS(287) is RNAUG3_init in component state_variables (dimensionless).
C CONSTS(288) is TSP_init in component state_variables (gram).
C CONSTS(289) is TVD_init in component state_variables (L_per_minute).
C CONSTS(290) is VAS1_init in component state_variables (litre).
C CONSTS(291) is VIC_init in component state_variables (litre).
C CONSTS(292) is VLA1_init in component state_variables (litre).
C CONSTS(293) is VPA1_init in component state_variables (litre).
C CONSTS(294) is VPF1_init in component state_variables (litre).
C CONSTS(295) is VP_init in component state_variables (litre).
C CONSTS(296) is VRA1_init in component state_variables (litre).
C CONSTS(297) is VRC_init in component state_variables (litre).
C CONSTS(298) is VTS2_init in component state_variables (litre).
C CONSTS(299) is VTW_init in component state_variables (litre).
C CONSTS(300) is VV6_init in component state_variables (litre).
C CONSTS(301) is VV7_init in component state_variables (litre).
C CONSTS(302) is VVS1_init in component state_variables (litre).
C ALGBRC(55) is BFM in component blood_flow_through_M_tissues (L_per_minute).
C ALGBRC(56) is BFN in component blood_flow_through_NM_NR_tissues (L_per_minute).
C ALGBRC(37) is PA in component arterial_pressure_and_pressure_gradient (mmHg).
C ALGBRC(18) is PLA in component left_atrial_pressure (mmHg).
C ALGBRC(9) is PPA in component pulmonary_vasculature_pressure (mmHg).
C ALGBRC(4) is PRA in component right_atrial_pressure (mmHg).
C ALGBRC(30) is PVS in component venous_average_pressure (mmHg).
C ALGBRC(59) is QAO in component systemic_blood_flow (L_per_minute).
C ALGBRC(42) is QLO in component left_ventricular_output (L_per_minute).
C ALGBRC(46) is QRO in component right_ventricular_output (L_per_minute).
C ALGBRC(15) is RPA in component pulmonary_arterial_resistance (mmHg_minute_per_L).
C ALGBRC(20) is RPV in component pulmonary_venous_resistance (mmHg_minute_per_L).
C CONSTS(583) is RVS in component venous_resistance (mmHg_minute_per_L).
C ALGBRC(28) is VVE in component venous_excess_volume (litre).
C CONSTS(303) is A1K in component parameter_values (minute).
C CONSTS(304) is A2K in component parameter_values (minute).
C CONSTS(305) is A3K in component parameter_values (minute).
C CONSTS(306) is A4K in component parameter_values (minute).
C CONSTS(307) is A4K2 in component parameter_values (minute).
C CONSTS(308) is AARK in component parameter_values (mmHg_minute_per_L).
C CONSTS(309) is AARLL in component parameter_values (mmHg_minute_per_L).
C CONSTS(310) is ADHINF in component parameter_values (dimensionless).
C CONSTS(311) is ADHKLL in component parameter_values (dimensionless).
C CONSTS(312) is ADHKUL in component parameter_values (dimensionless).
C CONSTS(313) is ADHPAM in component parameter_values (per_mmHg2).
C CONSTS(314) is ADHPUL in component parameter_values (mmHg).
C CONSTS(315) is ADHTC in component parameter_values (minute).
C CONSTS(316) is ADHVLL in component parameter_values (dimensionless).
C CONSTS(317) is ADHVUL in component parameter_values (dimensionless).
C CONSTS(318) is AH10 in component parameter_values (dimensionless).
C CONSTS(319) is AH11 in component parameter_values (minute).
C CONSTS(320) is AH9 in component parameter_values (per_mmHg).
C CONSTS(321) is AHMNAR in component parameter_values (dimensionless).
C CONSTS(322) is AHTHM in component parameter_values (dimensionless).
C CONSTS(323) is ALCLK in component parameter_values (dimensionless).
C CONSTS(324) is ALDINF in component parameter_values (dimensionless).
C CONSTS(325) is ALDKNS in component parameter_values (dimensionless).
C CONSTS(326) is ALDMM in component parameter_values (dimensionless).
C CONSTS(327) is AM1LL in component parameter_values (dimensionless).
C CONSTS(328) is AM1UL in component parameter_values (dimensionless).
C CONSTS(329) is AMCSNS in component parameter_values (dimensionless).
C CONSTS(330) is AMKM in component parameter_values (dimensionless).
C CONSTS(331) is AMKMUL in component parameter_values (dimensionless).
C CONSTS(332) is AMM1LM in component parameter_values (dimensionless).
C CONSTS(333) is AMM4 in component parameter_values (dimensionless).
C CONSTS(334) is AMNALL in component parameter_values (dimensionless).
C CONSTS(335) is AMNAM in component parameter_values (dimensionless).
C CONSTS(336) is AMNAUL in component parameter_values (dimensionless).
C CONSTS(337) is AMT in component parameter_values (minute).
C CONSTS(338) is ANCSNS in component parameter_values (dimensionless).
C CONSTS(339) is ANGINF in component parameter_values (dimensionless).
C CONSTS(340) is ANGKNS in component parameter_values (dimensionless).
C CONSTS(341) is ANMALD in component parameter_values (dimensionless).
C CONSTS(342) is ANMAM in component parameter_values (dimensionless).
C CONSTS(343) is ANMARL in component parameter_values (dimensionless).
C CONSTS(344) is ANMEM in component parameter_values (dimensionless).
C CONSTS(345) is ANMKEL in component parameter_values (dimensionless).
C CONSTS(346) is ANMKEM in component parameter_values (dimensionless).
C CONSTS(347) is ANMLL in component parameter_values (dimensionless).
C CONSTS(348) is ANMNAM in component parameter_values (dimensionless).
C CONSTS(349) is ANMSLT in component parameter_values (dimensionless).
C CONSTS(350) is ANMTM in component parameter_values (dimensionless).
C CONSTS(351) is ANMUL in component parameter_values (dimensionless).
C CONSTS(352) is ANPTC in component parameter_values (minute).
C CONSTS(353) is ANPXAF in component parameter_values (mmHg_minute_per_L).
C CONSTS(354) is ANPXUL in component parameter_values (dimensionless).
C CONSTS(355) is ANT in component parameter_values (minute).
C CONSTS(356) is ANULL in component parameter_values (dimensionless).
C CONSTS(357) is ANUM in component parameter_values (dimensionless).
C CONSTS(358) is ANUVM in component parameter_values (dimensionless).
C CONSTS(359) is ANV in component parameter_values (minute).
C CONSTS(360) is ANXM in component parameter_values (dimensionless).
C CONSTS(361) is ANY in component parameter_values (litre).
C CONSTS(362) is AR1LM in component parameter_values (dimensionless).
C CONSTS(363) is AR2LM in component parameter_values (dimensionless).
C CONSTS(364) is ARF in component parameter_values (dimensionless).
C CONSTS(365) is AUC1 in component parameter_values (dimensionless).
C CONSTS(366) is AUDMP in component parameter_values (minute).
C CONSTS(367) is AUK in component parameter_values (per_minute).
C CONSTS(368) is AUL in component parameter_values (litre).
C CONSTS(369) is AULPM in component parameter_values (dimensionless).
C CONSTS(370) is AUM1 in component parameter_values (dimensionless).
C CONSTS(371) is AUM2 in component parameter_values (dimensionless).
C CONSTS(372) is AUMAX in component parameter_values (dimensionless).
C CONSTS(373) is AUMIN in component parameter_values (dimensionless).
C CONSTS(374) is AUMK1 in component parameter_values (dimensionless).
C CONSTS(375) is AUN1 in component parameter_values (dimensionless).
C CONSTS(376) is AURG in component parameter_values (dimensionless).
C CONSTS(377) is AUS in component parameter_values (dimensionless).
C CONSTS(378) is AUSLP in component parameter_values (dimensionless).
C CONSTS(379) is AUTO in component parameter_values (dimensionless).
C CONSTS(380) is AUTOGN in component parameter_values (dimensionless).
C CONSTS(381) is AUTOSN in component parameter_values (dimensionless).
C CONSTS(382) is AUX in component parameter_values (dimensionless).
C CONSTS(383) is BAROTC in component parameter_values (minute).
C CONSTS(384) is CFC in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(385) is CKEEX in component parameter_values (dimensionless).
C CONSTS(386) is CMPTSS in component parameter_values (dimensionless).
C CONSTS(387) is CN2 in component parameter_values (per_mmHg).
C CONSTS(388) is CN7 in component parameter_values (dimensionless).
C CONSTS(389) is CNR in component parameter_values (monovalent_mEq_per_litre).
C CONSTS(390) is CPF in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(391) is CPK in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(392) is CPR in component parameter_values (gram_per_L).
C CONSTS(393) is CV in component parameter_values (L_per_mmHg).
C CONSTS(394) is CRRFLX in component parameter_values (mmHg).
C CONSTS(395) is DIURET in component parameter_values (dimensionless).
C CONSTS(396) is DHDTR in component parameter_values (per_mmHg_per_minute).
C CONSTS(397) is DTNAR in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(398) is DTNARL in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(399) is EARK in component parameter_values (mmHg_minute_per_L).
C CONSTS(400) is EARLL in component parameter_values (mmHg_minute_per_L).
C CONSTS(401) is EFAFR in component parameter_values (dimensionless).
C CONSTS(402) is EXC in component parameter_values (dimensionless).
C CONSTS(403) is EXCML in component parameter_values (dimensionless).
C CONSTS(404) is EXCXP in component parameter_values (dimensionless).
C CONSTS(405) is EXCXP2 in component parameter_values (dimensionless).
C CONSTS(406) is EXE in component parameter_values (mmHg).
C CONSTS(407) is GCOPF in component parameter_values (per_mmHg).
C CONSTS(408) is GFLC in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(409) is GFNDMP in component parameter_values (dimensionless).
C CONSTS(410) is GFNLL in component parameter_values (L_per_minute).
C CONSTS(411) is GPPD in component parameter_values (dimensionless).
C CONSTS(412) is HM6 in component parameter_values (mmHg).
C CONSTS(413) is HM8 in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(414) is HTAUML in component parameter_values (dimensionless).
C CONSTS(415) is KORGN in component parameter_values (dimensionless).
C CONSTS(416) is KORNER in component parameter_values (dimensionless).
C CONSTS(417) is KORREN in component parameter_values (dimensionless).
C CONSTS(418) is KORTC in component parameter_values (minute).
C CONSTS(419) is LPDE in component parameter_values (dimensionless).
C CONSTS(420) is LPK in component parameter_values (L_per_minute).
C CONSTS(421) is LPPR in component parameter_values (gram_per_minute).
C CONSTS(422) is MDFL1 in component parameter_values (dimensionless).
C CONSTS(423) is MDFLKM in component parameter_values (monovalent_mEq_per_litre).
C CONSTS(424) is MDFLW2 in component parameter_values (dimensionless).
C CONSTS(425) is MDFLWX in component parameter_values (dimensionless).
C CONSTS(426) is MDMP in component parameter_values (dimensionless).
C CONSTS(427) is O2A in component parameter_values (dimensionless).
C CONSTS(428) is O2CHMO in component parameter_values (per_mmHg).
C CONSTS(429) is OMM in component parameter_values (mL_per_minute).
C CONSTS(430) is PA4 in component parameter_values (dimensionless).
C CONSTS(431) is PAEX in component parameter_values (dimensionless).
C CONSTS(432) is PCE in component parameter_values (dimensionless).
C CONSTS(433) is PCR in component parameter_values (mmHg).
C CONSTS(434) is PR1LL in component parameter_values (mmHg).
C CONSTS(435) is PGHF in component parameter_values (L_mmHg_per_gram).
C CONSTS(436) is PK2 in component parameter_values (mmHg_per_mL).
C CONSTS(437) is PL2 in component parameter_values (L_mL_per_minute_per_mmHg).
C CONSTS(438) is PLDF in component parameter_values (mmHg).
C CONSTS(439) is PM5 in component parameter_values (per_mmHg).
C CONSTS(440) is PO2ADV in component parameter_values (dimensionless).
C CONSTS(441) is POK in component parameter_values (dimensionless).
C CONSTS(442) is POM in component parameter_values (dimensionless).
C CONSTS(443) is POM2 in component parameter_values (dimensionless).
C CONSTS(444) is PON in component parameter_values (dimensionless).
C CONSTS(445) is POR in component parameter_values (mmHg).
C CONSTS(446) is POZ in component parameter_values (dimensionless).
C CONSTS(447) is PXTP in component parameter_values (mmHg).
C CONSTS(448) is QAOLM in component parameter_values (dimensionless).
C CONSTS(449) is QRF in component parameter_values (L_per_minute).
C CONSTS(450) is RABSC in component parameter_values (per_mmHg).
C CONSTS(451) is RAM in component parameter_values (mmHg_minute_per_L).
C CONSTS(452) is RAPRSP in component parameter_values (mmHg).
C CONSTS(453) is RAR in component parameter_values (mmHg_minute_per_L).
C CONSTS(454) is RCDFDP in component parameter_values (minute).
C CONSTS(455) is RCDFPC in component parameter_values (dimensionless).
C CONSTS(456) is REK in component parameter_values (dimensionless).
C CONSTS(457) is RFABKM in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(458) is RFABDM in component parameter_values (dimensionless).
C CONSTS(459) is RFABDP in component parameter_values (dimensionless).
C CONSTS(460) is RFCDFT in component parameter_values (dimensionless).
C CONSTS(461) is RKC in component parameter_values (per_minute).
C CONSTS(462) is RNAGTC in component parameter_values (minute).
C CONSTS(463) is RNAUAD in component parameter_values (per_minute).
C CONSTS(464) is RPS in component parameter_values (dimensionless).
C CONSTS(465) is RTPPR in component parameter_values (dimensionless).
C CONSTS(466) is RTPPRS in component parameter_values (mmHg).
C CONSTS(467) is RTSPRS in component parameter_values (mmHg).
C CONSTS(468) is RVRS in component parameter_values (mmHg_minute_per_L).
C CONSTS(469) is RVSM in component parameter_values (mmHg_minute_per_L).
C CONSTS(470) is SR in component parameter_values (dimensionless).
C CONSTS(471) is SR2 in component parameter_values (dimensionless).
C CONSTS(472) is SRK in component parameter_values (minute).
C CONSTS(473) is SRK2 in component parameter_values (minute).
C CONSTS(474) is STA in component parameter_values (dimensionless).
C CONSTS(475) is TENSGN in component parameter_values (dimensionless).
C CONSTS(476) is TENSTC in component parameter_values (minute).
C CONSTS(477) is TIMETR in component parameter_values (minute).
C CONSTS(478) is TRNSFS in component parameter_values (dimensionless).
C CONSTS(479) is TRPL in component parameter_values (L_per_minute).
C CONSTS(480) is TSSLML in component parameter_values (dimensionless).
C CONSTS(481) is TSSLTC in component parameter_values (per_minute).
C CONSTS(482) is TVDDL in component parameter_values (minute).
C CONSTS(483) is U in component parameter_values (dimensionless).
C CONSTS(484) is VID1 in component parameter_values (L_per_minute).
C CONSTS(485) is VIDML in component parameter_values (litre2_per_monovalent_mEq_per_minute).
C CONSTS(486) is VNTSTM in component parameter_values (dimensionless).
C CONSTS(487) is VP1 in component parameter_values (dimensionless).
C CONSTS(488) is VPTISS in component parameter_values (litre).
C CONSTS(489) is VTSF in component parameter_values (litre).
C CONSTS(490) is VV9 in component parameter_values (litre).
C CONSTS(491) is X in component parameter_values (dimensionless).
C CONSTS(492) is Z in component parameter_values (dimensionless).
C CONSTS(493) is Z4 in component parameter_values (dimensionless).
C CONSTS(494) is Z5 in component parameter_values (dimensionless).
C CONSTS(495) is Z6 in component parameter_values (dimensionless).
C CONSTS(496) is Z7 in component parameter_values (dimensionless).
C CONSTS(497) is Z10 in component parameter_values (mmHg).
C CONSTS(498) is Z11 in component parameter_values (per_mmHg2).
C CONSTS(499) is Z12 in component parameter_values (dimensionless).
C CONSTS(500) is Z13 in component parameter_values (dimensionless).
C CONSTS(501) is Z14 in component parameter_values (dimensionless).
C CONSTS(502) is Z16 in component parameter_values (dimensionless).
C CONSTS(503) is Z18 in component parameter_values (dimensionless).
C CONSTS(504) is Z19 in component parameter_values (dimensionless).
C CONSTS(505) is KID in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(506) is URFORM in component parameter_values (mOsm_per_minute).
C CONSTS(507) is AUV in component parameter_values (dimensionless).
C CONSTS(508) is AUY in component parameter_values (dimensionless).
C CONSTS(509) is O2M in component parameter_values (mL_per_minute).
C CONSTS(510) is NID in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(511) is ATRFBM in component parameter_values (dimensionless).
C CONSTS(512) is GBL in component parameter_values (mmHg).
C CONSTS(513) is RNAUGN in component parameter_values (minute_per_L).
C CONSTS(514) is RNAULL in component parameter_values (dimensionless).
C CONSTS(515) is RNAUUL in component parameter_values (dimensionless).
C CONSTS(516) is HSR in component parameter_values (dimensionless).
C CONSTS(517) is HSL in component parameter_values (dimensionless).
C CONSTS(518) is HMK in component parameter_values (dimensionless).
C CONSTS(519) is PO2AMB in component parameter_values (mmHg).
C CONSTS(520) is HKM in component parameter_values (dimensionless).
C CONSTS(521) is DR in component parameter_values (L_per_minute).
C CONSTS(522) is ANPKNS in component parameter_values (dimensionless).
C CONSTS(523) is RMULT1 in component parameter_values (dimensionless).
C CONSTS(524) is GLPCA in component parameter_values (mmHg).
C CONSTS(525) is RFN in component parameter_values (L_per_minute).
C CONSTS(526) is RFABX in component parameter_values (dimensionless).
C CONSTS(527) is ANPINF in component parameter_values (dimensionless).
C CONSTS(528) is FIS in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(529) is CPPD in component parameter_values (gram_per_L).
C CONSTS(530) is TRRBC in component parameter_values (L_per_minute).
C CONSTS(531) is HYL in component parameter_values (gram).
C CONSTS(532) is ATRVM in component parameter_values (litre).
C CONSTS(533) is ADHC_init in component state_variables (dimensionless).
C CONSTS(534) is AHY_init in component state_variables (dimensionless).
C CONSTS(535) is AMC_init in component state_variables (dimensionless).
C CONSTS(536) is AMM1T_init in component state_variables (dimensionless).
C CONSTS(537) is AMM2_init in component state_variables (dimensionless).
C CONSTS(538) is ANC_init in component state_variables (dimensionless).
C CONSTS(539) is ANPC_init in component state_variables (dimensionless).
C CONSTS(540) is ANX1_init in component state_variables (dimensionless).
C CONSTS(541) is AR1T_init in component state_variables (dimensionless).
C CONSTS(542) is AR2T_init in component state_variables (dimensionless).
C CONSTS(543) is AR3T_init in component state_variables (dimensionless).
C CONSTS(544) is AU1_init in component state_variables (dimensionless).
C CONSTS(545) is AU4_init in component state_variables (dimensionless).
C CONSTS(546) is AU6_init in component state_variables (dimensionless).
C CONSTS(547) is DTKA_init in component state_variables (monovalent_mEq_per_minute).
C CONSTS(548) is HMD1_init in component state_variables (dimensionless).
C CONSTS(549) is HPL_init in component state_variables (dimensionless).
C CONSTS(550) is HPR_init in component state_variables (dimensionless).
C CONSTS(551) is KTOT_init in component state_variables (monovalent_mEq).
C CONSTS(552) is NAE_init in component state_variables (monovalent_mEq).
C CONSTS(553) is O2VAD1_init in component state_variables (dimensionless).
C CONSTS(554) is OSV_init in component state_variables (dimensionless).
C CONSTS(555) is OVA_init in component state_variables (mL_per_L).
C CONSTS(556) is OVS_init in component state_variables (dimensionless).
C CONSTS(557) is PAR1_init in component state_variables (mmHg).
C CONSTS(558) is PC_init in component state_variables (mmHg).
C CONSTS(559) is PLUR_init in component state_variables (mOsm).
C CONSTS(560) is PPR1_init in component state_variables (gram).
C CONSTS(561) is PRP_init in component state_variables (gram).
C CONSTS(562) is QO2T_init in component state_variables (mL).
C CONSTS(563) is QOM1_init in component state_variables (mL).
C CONSTS(564) is RNAUG3_init in component state_variables (dimensionless).
C CONSTS(565) is TSP_init in component state_variables (gram).
C CONSTS(566) is TVD_init in component state_variables (L_per_minute).
C CONSTS(567) is VAS1_init in component state_variables (litre).
C CONSTS(568) is VIC_init in component state_variables (litre).
C CONSTS(569) is VLA1_init in component state_variables (litre).
C CONSTS(570) is VPA1_init in component state_variables (litre).
C CONSTS(571) is VPF1_init in component state_variables (litre).
C CONSTS(572) is VP_init in component state_variables (litre).
C CONSTS(573) is VRA1_init in component state_variables (litre).
C CONSTS(574) is VRC_init in component state_variables (litre).
C CONSTS(575) is VTS2_init in component state_variables (litre).
C CONSTS(576) is VTW_init in component state_variables (litre).
C CONSTS(577) is VV6_init in component state_variables (litre).
C CONSTS(578) is VV7_init in component state_variables (litre).
C CONSTS(579) is VVS1_init in component state_variables (litre).
C STATES(1) is VVS1 in component venous_blood_volume (litre).
C STATES(2) is VAS1 in component arterial_blood_volume (litre).
C STATES(3) is VLA1 in component left_atrial_blood_volume (litre).
C STATES(4) is VPA1 in component pulmonary_vasculature_blood_volume (litre).
C STATES(5) is VRA1 in component right_atrial_blood_volume (litre).
C ALGBRC(1) is VBD in component total_blood_volume_change (litre).
C ALGBRC(34) is QVO in component rate_of_blood_flow_from_veins_to_right_atrium (L_per_minute).
C ALGBRC(2) is VRA in component right_atrial_blood_volume (litre).
C ALGBRC(48) is DRA in component right_atrial_blood_volume (L_per_minute).
C ALGBRC(3) is VRE in component right_atrial_pressure (litre).
C ALGBRC(5) is PRA1 in component autonomic_stimulation_effect_on_right_atrial_pressure (mmHg).
C ALGBRC(11) is RVM in component pressure_effect_on_right_ventricular_pumping (dimensionless).
C ALGBRC(10) is PP2 in component pressure_effect_on_right_ventricular_pumping (mmHg).
C ALGBRC(25) is QLN in component left_ventricular_output (L_per_minute).
C ALGBRC(43) is HPEF in component pumping_effectiveness_of_right_ventricle (L_per_minute).
C ALGBRC(6) is QRN in component right_ventricular_output (dimensionless).
C ALGBRC(23) is QPO in component rate_of_blood_flow_from_pulmonary_veins_to_left_atrium (L_per_minute).
C ALGBRC(7) is VPA in component pulmonary_vasculature_blood_volume (litre).
C ALGBRC(49) is DPA in component pulmonary_vasculature_blood_volume (L_per_minute).
C ALGBRC(8) is VPE in component pulmonary_vasculature_pressure (litre).
C ALGBRC(12) is PP1T in component pulmonary_arterial_resistance (L_per_minute_per_mmHg).
C ALGBRC(13) is PP1 in component pulmonary_arterial_resistance (L_per_minute_per_mmHg).
C ALGBRC(14) is CPA in component pulmonary_arterial_resistance (L_per_minute_per_mmHg).
C ALGBRC(19) is PL1 in component pulmonary_venous_resistance (mmHg).
C ALGBRC(21) is RPT in component total_pulmonary_vascular_resistance (mmHg_minute_per_L).
C ALGBRC(22) is PGL in component pressure_gradient_through_the_lungs (mmHg).
C ALGBRC(16) is VLA in component left_atrial_blood_volume (litre).
C ALGBRC(44) is DLA in component left_atrial_blood_volume (L_per_minute).
C ALGBRC(17) is VLE in component left_atrial_pressure (litre).
C ALGBRC(24) is PLA1 in component autonomic_stimulation_effect_on_left_atrial_pressure (mmHg).
C ALGBRC(39) is LVM in component pumping_effectiveness_of_left_ventricle (dimensionless).
C ALGBRC(38) is PA2 in component pumping_effectiveness_of_left_ventricle (mmHg).
C ALGBRC(40) is QLOT in component left_ventricular_output (L_per_minute).
C ALGBRC(41) is QLO1 in component left_ventricular_output (L_per_minute).
C ALGBRC(26) is VVS in component venous_blood_volume (litre).
C ALGBRC(60) is DVS in component venous_blood_volume (L_per_minute).
C CONSTS(580) is VVA in component angiotensin_induced_venous_constriction (litre).
C ALGBRC(27) is VVE1 in component venous_excess_volume (litre).
C ALGBRC(29) is PVS1 in component venous_average_pressure (mmHg).
C ALGBRC(31) is PR1 in component venous_outflow_pressure_into_heart (mmHg).
C ALGBRC(32) is RVG in component resistance_from_veins_to_right_atrium (mmHg_minute_per_L).
C ALGBRC(33) is PGV in component rate_of_blood_flow_from_veins_to_right_atrium (mmHg).
C CONSTS(581) is CN3 in component venous_resistance (dimensionless).
C CONSTS(582) is RV1 in component venous_resistance (mmHg_minute_per_L).
C CONSTS(584) is NNRVR in component NM_NR_venous_resistance (mmHg_minute_per_L).
C ALGBRC(35) is VAS in component arterial_blood_volume (litre).
C ALGBRC(61) is DAS in component arterial_blood_volume (L_per_minute).
C ALGBRC(45) is PAG in component arterial_pressure_and_pressure_gradient (mmHg).
C ALGBRC(36) is VAE in component arterial_pressure_and_pressure_gradient (litre).
C ALGBRC(47) is PAM in component pressure_effect_on_arterial_distention (dimensionless).
C ALGBRC(50) is R1 in component non_renal_systemic_arterial_resistance_multiplier (dimensionless).
C ALGBRC(51) is NNRAR in component NM_NR_arterial_resistance (mmHg_minute_per_L).
C ALGBRC(52) is PGS in component pressure_gradient_from_arteries_to_veins (mmHg).
C ALGBRC(53) is RSM in component M_systemic_resistance (mmHg_minute_per_L).
C ALGBRC(54) is RSN in component total_NM_NR_systemic_resistance (mmHg_minute_per_L).
C ALGBRC(57) is FISFLO in component blood_flow_through_AV_fistulas (L_per_minute).
C ALGBRC(58) is SYSFLO in component systemic_blood_flow (L_per_minute).
C ALGBRC(62) is RTP in component total_peripheral_resistance (mmHg_minute_per_L).
C RATES(5) is d/dt VRA1 in component right_atrial_blood_volume (litre).
C RATES(4) is d/dt VPA1 in component pulmonary_vasculature_blood_volume (litre).
C RATES(3) is d/dt VLA1 in component left_atrial_blood_volume (litre).
C RATES(1) is d/dt VVS1 in component venous_blood_volume (litre).
C RATES(2) is d/dt VAS1 in component arterial_blood_volume (litre).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 0.925271
CONSTS(2) = 1.0
CONSTS(3) = 1.0
CONSTS(4) = 1.0
CONSTS(5) = 1.0
CONSTS(6) = 1.16463
CONSTS(7) = 1.00076
CONSTS(8) = 1.0
CONSTS(9) = 1.0
CONSTS(10) = 16.9144
CONSTS(11) = 1.22057
CONSTS(12) = 2.50967
CONSTS(13) = 0.0101913
CONSTS(14) = 0.00366525
CONSTS(15) = 0.0
CONSTS(16) = 3.00449
CONSTS(17) = 2.00439
CONSTS(18) = 1.00237
CONSTS(19) = 1.0
CONSTS(20) = 1.00012
CONSTS(21) = 0.97287
CONSTS(22) = 1.00163
CONSTS(23) = 1.0
CONSTS(24) = 1.00066
CONSTS(25) = 1.00022
CONSTS(26) = 0.5
CONSTS(27) = 60
CONSTS(28) = 40000
CONSTS(29) = 0.1
CONSTS(30) = 40000
CONSTS(31) = 1
CONSTS(32) = 4
CONSTS(33) = 0
CONSTS(34) = 0.2
CONSTS(35) = 5
CONSTS(36) = 0.3
CONSTS(37) = 85
CONSTS(38) = 15
CONSTS(39) = 0.93617
CONSTS(40) = 2.5
CONSTS(41) = 0.333
CONSTS(42) = 1000
CONSTS(43) = 1
CONSTS(44) = 0.3
CONSTS(45) = 2
CONSTS(46) = 0.3
CONSTS(47) = 0
CONSTS(48) = 0
CONSTS(49) = 2.5
CONSTS(50) = 0
CONSTS(51) = 5
CONSTS(52) = 0.65
CONSTS(53) = 0.5
CONSTS(54) = 12
CONSTS(55) = 0.005
CONSTS(56) = 0.005
CONSTS(57) = 0.04
CONSTS(58) = 0.8
CONSTS(59) = 15
CONSTS(60) = 60
CONSTS(61) = 0.4
CONSTS(62) = 0
CONSTS(63) = 0
CONSTS(64) = 2.5
CONSTS(65) = 1.4
CONSTS(66) = 0.86
CONSTS(67) = 1.6
CONSTS(68) = 0.3
CONSTS(69) = 2
CONSTS(70) = 0.7
CONSTS(71) = 1
CONSTS(72) = 2
CONSTS(73) = 1.5
CONSTS(74) = 1.8
CONSTS(75) = 4
CONSTS(76) = 1.5
CONSTS(77) = 10
CONSTS(78) = 12
CONSTS(79) = 0.8
CONSTS(80) = 6
CONSTS(81) = 0
CONSTS(82) = 5000
CONSTS(83) = 0
CONSTS(84) = -0.2
CONSTS(85) = 0.005
CONSTS(86) = 0.005
CONSTS(87) = 0.5
CONSTS(88) = 0.3
CONSTS(89) = 0.3
CONSTS(90) = 0.004
CONSTS(91) = 1.5
CONSTS(92) = 0
CONSTS(93) = 3
CONSTS(94) = 1
CONSTS(95) = 5.0
CONSTS(96) = 0.4
CONSTS(97) = 0.3
CONSTS(98) = 0.5
CONSTS(99) = 1.0
CONSTS(100) = 1
CONSTS(101) = 1.5
CONSTS(102) = 0
CONSTS(103) = 7
CONSTS(104) = 0.9
CONSTS(105) = 1
CONSTS(106) = 0.16
CONSTS(107) = 0.01167
CONSTS(108) = 4
CONSTS(109) = 2
CONSTS(110) = 0.0212
CONSTS(111) = 0.2
CONSTS(112) = 139
CONSTS(113) = 0.0003
CONSTS(114) = 0.000253
CONSTS(115) = 40
CONSTS(116) = 0.1
CONSTS(117) = 0
CONSTS(118) = 1
CONSTS(119) = 0.05
CONSTS(120) = 0.675
CONSTS(121) = 1e-06
CONSTS(122) = 1
CONSTS(123) = 24
CONSTS(124) = 0
CONSTS(125) = 1
CONSTS(126) = 0.01
CONSTS(127) = 1.0
CONSTS(128) = 0.17
CONSTS(129) = 0
CONSTS(130) = 0.8092
CONSTS(131) = 0.0208333
CONSTS(132) = 3
CONSTS(133) = 0.001
CONSTS(134) = 1.0
CONSTS(135) = 1850
CONSTS(136) = 4.714e-08
CONSTS(137) = 0.4
CONSTS(138) = 1
CONSTS(139) = 0
CONSTS(140) = 0
CONSTS(141) = 10000
CONSTS(142) = 8
CONSTS(143) = 2.728e-14
CONSTS(144) = 0.03
CONSTS(145) = 10
CONSTS(146) = 0.667
CONSTS(147) = 0.005
CONSTS(148) = 0.01
CONSTS(149) = 0
CONSTS(150) = 0.1
CONSTS(151) = 0.01
CONSTS(152) = 57.1
CONSTS(153) = 1
CONSTS(154) = 2
CONSTS(155) = 1
CONSTS(156) = 15
CONSTS(157) = 0
CONSTS(158) = -2
CONSTS(159) = 0.79167
CONSTS(160) = 1.8
CONSTS(161) = 4.2
CONSTS(162) = 30
CONSTS(163) = 0.001
CONSTS(164) = 0.1
CONSTS(165) = 0.04
CONSTS(166) = 2
CONSTS(167) = 0.1
CONSTS(168) = 35
CONSTS(169) = 2
CONSTS(170) = 8
CONSTS(171) = 0.2
CONSTS(172) = 0.15
CONSTS(173) = 0.5
CONSTS(174) = 96.3
CONSTS(175) = 0
CONSTS(176) = 30.52
CONSTS(177) = 2000
CONSTS(178) = 0
CONSTS(179) = 1
CONSTS(180) = 0.03
CONSTS(181) = 0.3
CONSTS(182) = 1
CONSTS(183) = 0
CONSTS(184) = 5.8e-06
CONSTS(185) = 15
CONSTS(186) = 0
CONSTS(187) = 0
CONSTS(188) = 0.9
CONSTS(189) = 15.2
CONSTS(190) = 6
CONSTS(191) = 19.167
CONSTS(192) = 1
CONSTS(193) = 1
CONSTS(194) = 1
CONSTS(195) = 5
CONSTS(196) = 10000
CONSTS(197) = 0
CONSTS(198) = 0
CONSTS(199) = 0.05
CONSTS(200) = 0
CONSTS(201) = 0
CONSTS(202) = 0
CONSTS(203) = 0.15
CONSTS(204) = 0.005
CONSTS(205) = 30
CONSTS(206) = 4
CONSTS(207) = 0.01
CONSTS(208) = 0.01
CONSTS(209) = 1
CONSTS(210) = 0.01
CONSTS(211) = 0.0175
CONSTS(212) = 6
CONSTS(213) = 2.51
CONSTS(214) = 4
CONSTS(215) = 1
CONSTS(216) = 0.5
CONSTS(217) = 0.3
CONSTS(218) = 20
CONSTS(219) = 5
CONSTS(220) = 45
CONSTS(221) = 0.01
CONSTS(222) = 5
CONSTS(223) = 0.625
CONSTS(224) = 5
CONSTS(225) = 5
CONSTS(226) = 5
CONSTS(227) = 5
CONSTS(228) = 0.08
CONSTS(229) = 0.24
CONSTS(230) = 0.55
CONSTS(231) = 0
CONSTS(232) = 164
CONSTS(233) = 0.1
CONSTS(234) = 0
CONSTS(235) = 0
CONSTS(236) = 0.6
CONSTS(237) = 0.3
CONSTS(238) = 10
CONSTS(239) = 1
CONSTS(240) = 1
CONSTS(241) = 90
CONSTS(242) = 150
CONSTS(243) = 0.53333
CONSTS(244) = 0
CONSTS(245) = 0
CONSTS(246) = 1
CONSTS(247) = 1.0
CONSTS(248) = 1.22057
CONSTS(249) = 0.8
CONSTS(250) = 0
CONSTS(251) = 0
CONSTS(252) = 31.9719
CONSTS(253) = 0
CONSTS(254) = 60
CONSTS(255) = 0
CONSTS(256) = 1.0
CONSTS(257) = 0.301963
CONSTS(258) = 1.0
CONSTS(259) = 1.00269
CONSTS(260) = 1.09071
CONSTS(261) = 0.859476
CONSTS(262) = 1.0
CONSTS(263) = 0
CONSTS(264) = 1.02127
CONSTS(265) = 1.01179
CONSTS(266) = 1.1448
CONSTS(267) = 1.00007
CONSTS(268) = -0.060024
CONSTS(269) = 1.00132
CONSTS(270) = 0.0367573
CONSTS(271) = 1
CONSTS(272) = 1.00163
CONSTS(273) = 1.00237
CONSTS(274) = 3622.54
CONSTS(275) = 2109.91
CONSTS(276) = 2.368e-07
CONSTS(277) = 0.694163
CONSTS(278) = 204.497
CONSTS(279) = 0.699673
CONSTS(280) = 103.525
CONSTS(281) = 16.9144
CONSTS(282) = 159.549
CONSTS(283) = 0.419998
CONSTS(284) = 216.243
CONSTS(285) = 72.2362
CONSTS(286) = 48.0839
CONSTS(287) = 0
CONSTS(288) = 279.945
CONSTS(289) = 0.000980838
CONSTS(290) = 0.862514
CONSTS(291) = 25.0404
CONSTS(292) = 0.379883
CONSTS(293) = 0.38131
CONSTS(294) = 0.0123238
CONSTS(295) = 3.00449
CONSTS(296) = 0.100043
CONSTS(297) = 2.00439
CONSTS(298) = 0
CONSTS(299) = 39.8952
CONSTS(300) = 0.0101913
CONSTS(301) = 0.00366525
CONSTS(302) = 3.28246
CONSTS(303) = 0.5
CONSTS(304) = 60
CONSTS(305) = 40000
CONSTS(306) = 0.1
CONSTS(307) = 40000
CONSTS(308) = 1
CONSTS(309) = 4
CONSTS(310) = 0
CONSTS(311) = 0.2
CONSTS(312) = 5
CONSTS(313) = 0.3
CONSTS(314) = 85
CONSTS(315) = 15
CONSTS(316) = 0.93617
CONSTS(317) = 2.5
CONSTS(318) = 0.333
CONSTS(319) = 1000
CONSTS(320) = 1
CONSTS(321) = 0.3
CONSTS(322) = 2
CONSTS(323) = 0.3
CONSTS(324) = 0
CONSTS(325) = 0
CONSTS(326) = 2.5
CONSTS(327) = 0
CONSTS(328) = 5
CONSTS(329) = 0.65
CONSTS(330) = 0.5
CONSTS(331) = 12
CONSTS(332) = 0.005
CONSTS(333) = 0.005
CONSTS(334) = 0.04
CONSTS(335) = 0.8
CONSTS(336) = 15
CONSTS(337) = 60
CONSTS(338) = 0.4
CONSTS(339) = 0
CONSTS(340) = 0
CONSTS(341) = 2.5
CONSTS(342) = 1.4
CONSTS(343) = 0.86
CONSTS(344) = 1.6
CONSTS(345) = 0.3
CONSTS(346) = 2
CONSTS(347) = 0.7
CONSTS(348) = 1
CONSTS(349) = 2
CONSTS(350) = 1.5
CONSTS(351) = 1.8
CONSTS(352) = 4
CONSTS(353) = 1.5
CONSTS(354) = 10
CONSTS(355) = 12
CONSTS(356) = 0.8
CONSTS(357) = 6
CONSTS(358) = 0
CONSTS(359) = 5000
CONSTS(360) = 0
CONSTS(361) = -0.2
CONSTS(362) = 0.005
CONSTS(363) = 0.005
CONSTS(364) = 0.5
CONSTS(365) = 0.3
CONSTS(366) = 0.3
CONSTS(367) = 0.004
CONSTS(368) = 1.5
CONSTS(369) = 0
CONSTS(370) = 3
CONSTS(371) = 1
CONSTS(372) = 5.0
CONSTS(373) = 0.4
CONSTS(374) = 0.3
CONSTS(375) = 0.5
CONSTS(376) = 1.0
CONSTS(377) = 1
CONSTS(378) = 1.5
CONSTS(379) = 0
CONSTS(380) = 7
CONSTS(381) = 0.9
CONSTS(382) = 1
CONSTS(383) = 0.16
CONSTS(384) = 0.01167
CONSTS(385) = 4
CONSTS(386) = 2
CONSTS(387) = 0.0212
CONSTS(388) = 0.2
CONSTS(389) = 139
CONSTS(390) = 0.0003
CONSTS(391) = 0.000253
CONSTS(392) = 40
CONSTS(393) = 0.1
CONSTS(394) = 0
CONSTS(395) = 1
CONSTS(396) = 0.05
CONSTS(397) = 0.675
CONSTS(398) = 1e-06
CONSTS(399) = 1
CONSTS(400) = 24
CONSTS(401) = 0
CONSTS(402) = 1
CONSTS(403) = 0.01
CONSTS(404) = 1.0
CONSTS(405) = 0.17
CONSTS(406) = 0
CONSTS(407) = 0.8092
CONSTS(408) = 0.0208333
CONSTS(409) = 3
CONSTS(410) = 0.001
CONSTS(411) = 1.0
CONSTS(412) = 1850
CONSTS(413) = 4.714e-08
CONSTS(414) = 0.4
CONSTS(415) = 1
CONSTS(416) = 0
CONSTS(417) = 0
CONSTS(418) = 10000
CONSTS(419) = 8
CONSTS(420) = 2.728e-14
CONSTS(421) = 0.03
CONSTS(422) = 10
CONSTS(423) = 0.667
CONSTS(424) = 0.005
CONSTS(425) = 0.01
CONSTS(426) = 0
CONSTS(427) = 0.1
CONSTS(428) = 0.01
CONSTS(429) = 57.1
CONSTS(430) = 1
CONSTS(431) = 2
CONSTS(432) = 1
CONSTS(433) = 15
CONSTS(434) = 0
CONSTS(435) = -2
CONSTS(436) = 0.79167
CONSTS(437) = 1.8
CONSTS(438) = 4.2
CONSTS(439) = 30
CONSTS(440) = 0.001
CONSTS(441) = 0.1
CONSTS(442) = 0.04
CONSTS(443) = 2
CONSTS(444) = 0.1
CONSTS(445) = 35
CONSTS(446) = 2
CONSTS(447) = 8
CONSTS(448) = 0.2
CONSTS(449) = 0.15
CONSTS(450) = 0.5
CONSTS(451) = 96.3
CONSTS(452) = 0
CONSTS(453) = 30.52
CONSTS(454) = 2000
CONSTS(455) = 0
CONSTS(456) = 1
CONSTS(457) = 0.03
CONSTS(458) = 0.3
CONSTS(459) = 1
CONSTS(460) = 0
CONSTS(461) = 5.8e-06
CONSTS(462) = 15
CONSTS(463) = 0
CONSTS(464) = 0
CONSTS(465) = 0.9
CONSTS(466) = 15.2
CONSTS(467) = 6
CONSTS(468) = 19.167
CONSTS(469) = 1
CONSTS(470) = 1
CONSTS(471) = 1
CONSTS(472) = 5
CONSTS(473) = 10000
CONSTS(474) = 0
CONSTS(475) = 0
CONSTS(476) = 0.05
CONSTS(477) = 0
CONSTS(478) = 0
CONSTS(479) = 0
CONSTS(480) = 0.15
CONSTS(481) = 0.005
CONSTS(482) = 30
CONSTS(483) = 4
CONSTS(484) = 0.01
CONSTS(485) = 0.01
CONSTS(486) = 1
CONSTS(487) = 0.01
CONSTS(488) = 0.0175
CONSTS(489) = 6
CONSTS(490) = 2.51
CONSTS(491) = 4
CONSTS(492) = 1
CONSTS(493) = 0.5
CONSTS(494) = 0.3
CONSTS(495) = 20
CONSTS(496) = 5
CONSTS(497) = 45
CONSTS(498) = 0.01
CONSTS(499) = 5
CONSTS(500) = 0.625
CONSTS(501) = 5
CONSTS(502) = 5
CONSTS(503) = 5
CONSTS(504) = 5
CONSTS(505) = 0.08
CONSTS(506) = 0.24
CONSTS(507) = 0.55
CONSTS(508) = 0
CONSTS(509) = 164
CONSTS(510) = 0.1
CONSTS(511) = 0
CONSTS(512) = 0
CONSTS(513) = 0.6
CONSTS(514) = 0.3
CONSTS(515) = 10
CONSTS(516) = 1
CONSTS(517) = 1
CONSTS(518) = 90
CONSTS(519) = 150
CONSTS(520) = 0.53333
CONSTS(521) = 0
CONSTS(522) = 0
CONSTS(523) = 1
CONSTS(524) = 1.0
CONSTS(525) = 1.22057
CONSTS(526) = 0.8
CONSTS(527) = 0
CONSTS(528) = 0
CONSTS(529) = 31.9719
CONSTS(530) = 0
CONSTS(531) = 60
CONSTS(532) = 0
CONSTS(533) = 1.0
CONSTS(534) = 0.301963
CONSTS(535) = 1.0
CONSTS(536) = 1.00269
CONSTS(537) = 1.09071
CONSTS(538) = 0.859476
CONSTS(539) = 1.0
CONSTS(540) = 0
CONSTS(541) = 1.02127
CONSTS(542) = 1.01179
CONSTS(543) = 1.1448
CONSTS(544) = 1.00007
CONSTS(545) = -0.060024
CONSTS(546) = 1.00132
CONSTS(547) = 0.0367573
CONSTS(548) = 1
CONSTS(549) = 1.00163
CONSTS(550) = 1.00237
CONSTS(551) = 3622.54
CONSTS(552) = 2109.91
CONSTS(553) = 2.368e-07
CONSTS(554) = 0.694163
CONSTS(555) = 204.497
CONSTS(556) = 0.699673
CONSTS(557) = 103.525
CONSTS(558) = 16.9144
CONSTS(559) = 159.549
CONSTS(560) = 0.419998
CONSTS(561) = 216.243
CONSTS(562) = 72.2362
CONSTS(563) = 48.0839
CONSTS(564) = 0
CONSTS(565) = 279.945
CONSTS(566) = 0.000980838
CONSTS(567) = 0.862514
CONSTS(568) = 25.0404
CONSTS(569) = 0.379883
CONSTS(570) = 0.38131
CONSTS(571) = 0.0123238
CONSTS(572) = 3.00449
CONSTS(573) = 0.100043
CONSTS(574) = 2.00439
CONSTS(575) = 0
CONSTS(576) = 39.8952
CONSTS(577) = 0.0101913
CONSTS(578) = 0.00366525
CONSTS(579) = 3.28246
CONSTS(580) = (CONSTS(1) - 1.00000)*CONSTS(361)
CONSTS(581) = ( (CONSTS(10) - 17.0000)*CONSTS(388)+17.0000)*CONSTS(387)
CONSTS(582) = CONSTS(469)/CONSTS(581)
CONSTS(583) = CONSTS(8)*CONSTS(582)*CONSTS(7)*CONSTS(9)
CONSTS(584) = CONSTS(583)*1.79000
STATES(1) = CONSTS(579)
STATES(2) = CONSTS(567)
STATES(3) = CONSTS(569)
STATES(4) = CONSTS(570)
STATES(5) = CONSTS(573)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = ((((((CONSTS(16)+CONSTS(17)) - STATES(1)) - STATES(2)) - STATES(3)) - STATES(4)) - STATES(5))/2.00000
ALGBRC(16) = STATES(3)+ ALGBRC(1)*0.128000
ALGBRC(17) = ALGBRC(16) - 0.380000
ALGBRC(18) = ALGBRC(17)/0.0100000
ALGBRC(24) = (ALGBRC(18)+4.00000)*( CONSTS(414)*(CONSTS(25) - 1.00000)+1.00000) - 4.00000
ALGBRC(25) = TERNRY(ALGBRC(24).LE.- 2.00000, 0.0100000, TERNRY(ALGBRC(24).GT.- 2.00000.AND.ALGBRC(24).LE.1.00000, 0.0100000+( (3.60000 - 0.0100000)*(ALGBRC(24) - - 2.00000))/(1.00000 - - 2.00000), TERNRY(ALGBRC(24).GT.1.00000.AND.ALGBRC(24).LE.5.00000, 3.60000+( (9.40000 - 3.60000)*(ALGBRC(24) - 1.00000))/(5.00000 - 1.00000), TERNRY(ALGBRC(24).GT.5.00000.AND.ALGBRC(24).LE.8.00000, 9.40000+( (11.6000 - 9.40000)*(ALGBRC(24) - 5.00000))/(8.00000 - 5.00000), TERNRY(ALGBRC(24).GT.8.00000.AND.ALGBRC(24).LE.12.0000, 11.6000+( (13.5000 - 11.6000)*(ALGBRC(24) - 8.00000))/(12.0000 - 8.00000), 13.5000)
ALGBRC(35) = STATES(2)+ ALGBRC(1)*0.261000
ALGBRC(36) = ALGBRC(35) - 0.495000
ALGBRC(37) = ALGBRC(36)/0.00355000
ALGBRC(38) = ALGBRC(37)/( CONSTS(20)*CONSTS(21))
ALGBRC(39) = TERNRY(ALGBRC(38).LE.0.00000, 1.04000, TERNRY(ALGBRC(38).GT.0.00000.AND.ALGBRC(38).LE.60.0000, 1.04000+( (1.02500 - 1.04000)*(ALGBRC(38) - 0.00000))/(60.0000 - 0.00000), TERNRY(ALGBRC(38).GT.60.0000.AND.ALGBRC(38).LE.125.000, 1.02500+( (0.970000 - 1.02500)*(ALGBRC(38) - 60.0000))/(125.000 - 60.0000), TERNRY(ALGBRC(38).GT.125.000.AND.ALGBRC(38).LE.160.000, 0.970000+( (0.880000 - 0.970000)*(ALGBRC(38) - 125.000))/(160.000 - 125.000), TERNRY(ALGBRC(38).GT.160.000.AND.ALGBRC(38).LE.200.000, 0.880000+( (0.590000 - 0.880000)*(ALGBRC(38) - 160.000))/(200.000 - 160.000), TERNRY(ALGBRC(38).GT.200.000.AND.ALGBRC(38).LE.240.000, 0.590000+( (0.00000 - 0.590000)*(ALGBRC(38) - 200.000))/(240.000 - 200.000), 0.00000)
ALGBRC(40) = ALGBRC(39)*ALGBRC(25)*CONSTS(20)*CONSTS(517)*CONSTS(19)*CONSTS(22)
ALGBRC(41) = (ALGBRC(18) - ALGBRC(37))/3.00000
ALGBRC(42) = TERNRY(ALGBRC(41).GT.0.00000, ALGBRC(40)+ALGBRC(41), ALGBRC(40))
ALGBRC(7) = STATES(4)+ ALGBRC(1)*0.155000
ALGBRC(8) = ALGBRC(7) - 0.306250
ALGBRC(9) = ALGBRC(8)/0.00480000
ALGBRC(12) = 0.0260000*ALGBRC(9)
ALGBRC(13) = TERNRY(ALGBRC(12).LT.1.00000e-05, 1.00000e-05, ALGBRC(12))
ALGBRC(14) = ALGBRC(13) ** 0.500000
ALGBRC(15) = 1.00000/ALGBRC(14)
ALGBRC(19) = ALGBRC(18)+18.0000
ALGBRC(20) = 1.00000/( ALGBRC(19)*0.0357000)
ALGBRC(21) = ALGBRC(20)+ALGBRC(15)
ALGBRC(22) = ALGBRC(9) - ALGBRC(18)
ALGBRC(23) = ALGBRC(22)/ALGBRC(21)
ALGBRC(44) = ALGBRC(23) - ALGBRC(42)
RATES(3) = ALGBRC(44)
ALGBRC(10) = (ALGBRC(9)/CONSTS(20))/CONSTS(21)
ALGBRC(11) = TERNRY(ALGBRC(10).LE.0.00000, 1.06000, TERNRY(ALGBRC(10).GT.0.00000.AND.ALGBRC(10).LE.32.0000, 1.06000+( (0.970000 - 1.06000)*(ALGBRC(10) - 0.00000))/(32.0000 - 0.00000), TERNRY(ALGBRC(10).GT.32.0000.AND.ALGBRC(10).LE.38.4000, 0.970000+( (0.930000 - 0.970000)*(ALGBRC(10) - 32.0000))/(38.4000 - 32.0000), TERNRY(ALGBRC(10).GT.38.4000.AND.ALGBRC(10).LE.48.0000, 0.930000+( (0.800000 - 0.930000)*(ALGBRC(10) - 38.4000))/(48.0000 - 38.4000), TERNRY(ALGBRC(10).GT.48.0000.AND.ALGBRC(10).LE.60.8000, 0.800000+( (0.460000 - 0.800000)*(ALGBRC(10) - 48.0000))/(60.8000 - 48.0000), TERNRY(ALGBRC(10).GT.60.8000.AND.ALGBRC(10).LE.72.0000, 0.460000+( (0.00000 - 0.460000)*(ALGBRC(10) - 60.8000))/(72.0000 - 60.8000), 0.00000)
ALGBRC(43) = (1.00000 - CONSTS(449))*CONSTS(20)*ALGBRC(11)*CONSTS(516)*CONSTS(19)*CONSTS(18)+( CONSTS(449)*ALGBRC(42))/ALGBRC(25)
ALGBRC(2) = STATES(5)+ ALGBRC(1)*0.0574000
ALGBRC(3) = ALGBRC(2) - 0.100000
ALGBRC(4) = ALGBRC(3)/0.00500000
ALGBRC(5) = (ALGBRC(4)+8.00000)*( CONSTS(414)*(CONSTS(25) - 1.00000)+1.00000) - 8.00000
ALGBRC(6) = TERNRY(ALGBRC(5).LE.- 8.00000, 0.00000, TERNRY(ALGBRC(5).GT.- 8.00000.AND.ALGBRC(5).LE.- 6.00000, 0.00000+( (0.750000 - 0.00000)*(ALGBRC(5) - - 8.00000))/(- 6.00000 - - 8.00000), TERNRY(ALGBRC(5).GT.- 6.00000.AND.ALGBRC(5).LE.- 2.00000, 0.750000+( (2.60000 - 0.750000)*(ALGBRC(5) - - 6.00000))/(- 2.00000 - - 6.00000), TERNRY(ALGBRC(5).GT.- 2.00000.AND.ALGBRC(5).LE.4.00000, 2.60000+( (9.80000 - 2.60000)*(ALGBRC(5) - - 2.00000))/(4.00000 - - 2.00000), TERNRY(ALGBRC(5).GT.4.00000.AND.ALGBRC(5).LE.12.0000, 9.80000+( (13.5000 - 9.80000)*(ALGBRC(5) - 4.00000))/(12.0000 - 4.00000), 13.5000)
ALGBRC(46) = ALGBRC(6)*ALGBRC(43)
ALGBRC(26) = STATES(1)+ ALGBRC(1)*0.398600
ALGBRC(27) = ((((ALGBRC(26) - CONSTS(12)) - CONSTS(580)) - CONSTS(14)) - CONSTS(13)) - CONSTS(15)
ALGBRC(28) = TERNRY(ALGBRC(27).LT.0.000100000, 0.000100000, ALGBRC(27))
ALGBRC(29) = 3.70000+(ALGBRC(28) - 0.740000)/CONSTS(393)
ALGBRC(30) = TERNRY(ALGBRC(29).LT.0.000100000, 0.000100000, ALGBRC(29))
ALGBRC(32) = 0.740000/ALGBRC(30)/( CONSTS(7)*3.70000) ** 0.500000
ALGBRC(31) = TERNRY(ALGBRC(4).LT.CONSTS(434), CONSTS(434), ALGBRC(4))
ALGBRC(33) = ALGBRC(30) - ALGBRC(31)
ALGBRC(34) = ALGBRC(33)/ALGBRC(32)
ALGBRC(48) = ALGBRC(34) - ALGBRC(46)
RATES(5) = ALGBRC(48)
ALGBRC(49) = ALGBRC(46) - ALGBRC(23)
RATES(4) = ALGBRC(49)
ALGBRC(45) = ALGBRC(37) - ALGBRC(4)
ALGBRC(57) = ALGBRC(45)*CONSTS(528)
ALGBRC(52) = ALGBRC(37) - ALGBRC(30)
ALGBRC(47) = ALGBRC(37)/100.000 ** CONSTS(431)
ALGBRC(50) = (( CONSTS(1)*CONSTS(2)*CONSTS(24)*CONSTS(7)*CONSTS(3))/ALGBRC(47))/CONSTS(4)
ALGBRC(53) = CONSTS(451)*CONSTS(23)*ALGBRC(50)*CONSTS(5)*CONSTS(523)
ALGBRC(55) = ALGBRC(52)/ALGBRC(53)
ALGBRC(51) = CONSTS(453)*CONSTS(6)*ALGBRC(50)*CONSTS(5)*CONSTS(523)
ALGBRC(54) = ALGBRC(51)+CONSTS(584)
ALGBRC(56) = ALGBRC(52)/ALGBRC(54)
ALGBRC(58) = ALGBRC(55)+ALGBRC(56)+CONSTS(11)
ALGBRC(59) = ALGBRC(58)+ALGBRC(57)
ALGBRC(60) = ALGBRC(59) - ALGBRC(34)
RATES(1) = ALGBRC(60)
ALGBRC(61) = ALGBRC(42) - ALGBRC(59)
RATES(2) = ALGBRC(61)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = ((((((CONSTS(16)+CONSTS(17)) - STATES(1)) - STATES(2)) - STATES(3)) - STATES(4)) - STATES(5))/2.00000
ALGBRC(16) = STATES(3)+ ALGBRC(1)*0.128000
ALGBRC(17) = ALGBRC(16) - 0.380000
ALGBRC(18) = ALGBRC(17)/0.0100000
ALGBRC(24) = (ALGBRC(18)+4.00000)*( CONSTS(414)*(CONSTS(25) - 1.00000)+1.00000) - 4.00000
ALGBRC(25) = TERNRY(ALGBRC(24).LE.- 2.00000, 0.0100000, TERNRY(ALGBRC(24).GT.- 2.00000.AND.ALGBRC(24).LE.1.00000, 0.0100000+( (3.60000 - 0.0100000)*(ALGBRC(24) - - 2.00000))/(1.00000 - - 2.00000), TERNRY(ALGBRC(24).GT.1.00000.AND.ALGBRC(24).LE.5.00000, 3.60000+( (9.40000 - 3.60000)*(ALGBRC(24) - 1.00000))/(5.00000 - 1.00000), TERNRY(ALGBRC(24).GT.5.00000.AND.ALGBRC(24).LE.8.00000, 9.40000+( (11.6000 - 9.40000)*(ALGBRC(24) - 5.00000))/(8.00000 - 5.00000), TERNRY(ALGBRC(24).GT.8.00000.AND.ALGBRC(24).LE.12.0000, 11.6000+( (13.5000 - 11.6000)*(ALGBRC(24) - 8.00000))/(12.0000 - 8.00000), 13.5000)
ALGBRC(35) = STATES(2)+ ALGBRC(1)*0.261000
ALGBRC(36) = ALGBRC(35) - 0.495000
ALGBRC(37) = ALGBRC(36)/0.00355000
ALGBRC(38) = ALGBRC(37)/( CONSTS(20)*CONSTS(21))
ALGBRC(39) = TERNRY(ALGBRC(38).LE.0.00000, 1.04000, TERNRY(ALGBRC(38).GT.0.00000.AND.ALGBRC(38).LE.60.0000, 1.04000+( (1.02500 - 1.04000)*(ALGBRC(38) - 0.00000))/(60.0000 - 0.00000), TERNRY(ALGBRC(38).GT.60.0000.AND.ALGBRC(38).LE.125.000, 1.02500+( (0.970000 - 1.02500)*(ALGBRC(38) - 60.0000))/(125.000 - 60.0000), TERNRY(ALGBRC(38).GT.125.000.AND.ALGBRC(38).LE.160.000, 0.970000+( (0.880000 - 0.970000)*(ALGBRC(38) - 125.000))/(160.000 - 125.000), TERNRY(ALGBRC(38).GT.160.000.AND.ALGBRC(38).LE.200.000, 0.880000+( (0.590000 - 0.880000)*(ALGBRC(38) - 160.000))/(200.000 - 160.000), TERNRY(ALGBRC(38).GT.200.000.AND.ALGBRC(38).LE.240.000, 0.590000+( (0.00000 - 0.590000)*(ALGBRC(38) - 200.000))/(240.000 - 200.000), 0.00000)
ALGBRC(40) = ALGBRC(39)*ALGBRC(25)*CONSTS(20)*CONSTS(517)*CONSTS(19)*CONSTS(22)
ALGBRC(41) = (ALGBRC(18) - ALGBRC(37))/3.00000
ALGBRC(42) = TERNRY(ALGBRC(41).GT.0.00000, ALGBRC(40)+ALGBRC(41), ALGBRC(40))
ALGBRC(7) = STATES(4)+ ALGBRC(1)*0.155000
ALGBRC(8) = ALGBRC(7) - 0.306250
ALGBRC(9) = ALGBRC(8)/0.00480000
ALGBRC(12) = 0.0260000*ALGBRC(9)
ALGBRC(13) = TERNRY(ALGBRC(12).LT.1.00000e-05, 1.00000e-05, ALGBRC(12))
ALGBRC(14) = ALGBRC(13) ** 0.500000
ALGBRC(15) = 1.00000/ALGBRC(14)
ALGBRC(19) = ALGBRC(18)+18.0000
ALGBRC(20) = 1.00000/( ALGBRC(19)*0.0357000)
ALGBRC(21) = ALGBRC(20)+ALGBRC(15)
ALGBRC(22) = ALGBRC(9) - ALGBRC(18)
ALGBRC(23) = ALGBRC(22)/ALGBRC(21)
ALGBRC(44) = ALGBRC(23) - ALGBRC(42)
ALGBRC(10) = (ALGBRC(9)/CONSTS(20))/CONSTS(21)
ALGBRC(11) = TERNRY(ALGBRC(10).LE.0.00000, 1.06000, TERNRY(ALGBRC(10).GT.0.00000.AND.ALGBRC(10).LE.32.0000, 1.06000+( (0.970000 - 1.06000)*(ALGBRC(10) - 0.00000))/(32.0000 - 0.00000), TERNRY(ALGBRC(10).GT.32.0000.AND.ALGBRC(10).LE.38.4000, 0.970000+( (0.930000 - 0.970000)*(ALGBRC(10) - 32.0000))/(38.4000 - 32.0000), TERNRY(ALGBRC(10).GT.38.4000.AND.ALGBRC(10).LE.48.0000, 0.930000+( (0.800000 - 0.930000)*(ALGBRC(10) - 38.4000))/(48.0000 - 38.4000), TERNRY(ALGBRC(10).GT.48.0000.AND.ALGBRC(10).LE.60.8000, 0.800000+( (0.460000 - 0.800000)*(ALGBRC(10) - 48.0000))/(60.8000 - 48.0000), TERNRY(ALGBRC(10).GT.60.8000.AND.ALGBRC(10).LE.72.0000, 0.460000+( (0.00000 - 0.460000)*(ALGBRC(10) - 60.8000))/(72.0000 - 60.8000), 0.00000)
ALGBRC(43) = (1.00000 - CONSTS(449))*CONSTS(20)*ALGBRC(11)*CONSTS(516)*CONSTS(19)*CONSTS(18)+( CONSTS(449)*ALGBRC(42))/ALGBRC(25)
ALGBRC(2) = STATES(5)+ ALGBRC(1)*0.0574000
ALGBRC(3) = ALGBRC(2) - 0.100000
ALGBRC(4) = ALGBRC(3)/0.00500000
ALGBRC(5) = (ALGBRC(4)+8.00000)*( CONSTS(414)*(CONSTS(25) - 1.00000)+1.00000) - 8.00000
ALGBRC(6) = TERNRY(ALGBRC(5).LE.- 8.00000, 0.00000, TERNRY(ALGBRC(5).GT.- 8.00000.AND.ALGBRC(5).LE.- 6.00000, 0.00000+( (0.750000 - 0.00000)*(ALGBRC(5) - - 8.00000))/(- 6.00000 - - 8.00000), TERNRY(ALGBRC(5).GT.- 6.00000.AND.ALGBRC(5).LE.- 2.00000, 0.750000+( (2.60000 - 0.750000)*(ALGBRC(5) - - 6.00000))/(- 2.00000 - - 6.00000), TERNRY(ALGBRC(5).GT.- 2.00000.AND.ALGBRC(5).LE.4.00000, 2.60000+( (9.80000 - 2.60000)*(ALGBRC(5) - - 2.00000))/(4.00000 - - 2.00000), TERNRY(ALGBRC(5).GT.4.00000.AND.ALGBRC(5).LE.12.0000, 9.80000+( (13.5000 - 9.80000)*(ALGBRC(5) - 4.00000))/(12.0000 - 4.00000), 13.5000)
ALGBRC(46) = ALGBRC(6)*ALGBRC(43)
ALGBRC(26) = STATES(1)+ ALGBRC(1)*0.398600
ALGBRC(27) = ((((ALGBRC(26) - CONSTS(12)) - CONSTS(580)) - CONSTS(14)) - CONSTS(13)) - CONSTS(15)
ALGBRC(28) = TERNRY(ALGBRC(27).LT.0.000100000, 0.000100000, ALGBRC(27))
ALGBRC(29) = 3.70000+(ALGBRC(28) - 0.740000)/CONSTS(393)
ALGBRC(30) = TERNRY(ALGBRC(29).LT.0.000100000, 0.000100000, ALGBRC(29))
ALGBRC(32) = 0.740000/ALGBRC(30)/( CONSTS(7)*3.70000) ** 0.500000
ALGBRC(31) = TERNRY(ALGBRC(4).LT.CONSTS(434), CONSTS(434), ALGBRC(4))
ALGBRC(33) = ALGBRC(30) - ALGBRC(31)
ALGBRC(34) = ALGBRC(33)/ALGBRC(32)
ALGBRC(48) = ALGBRC(34) - ALGBRC(46)
ALGBRC(49) = ALGBRC(46) - ALGBRC(23)
ALGBRC(45) = ALGBRC(37) - ALGBRC(4)
ALGBRC(57) = ALGBRC(45)*CONSTS(528)
ALGBRC(52) = ALGBRC(37) - ALGBRC(30)
ALGBRC(47) = ALGBRC(37)/100.000 ** CONSTS(431)
ALGBRC(50) = (( CONSTS(1)*CONSTS(2)*CONSTS(24)*CONSTS(7)*CONSTS(3))/ALGBRC(47))/CONSTS(4)
ALGBRC(53) = CONSTS(451)*CONSTS(23)*ALGBRC(50)*CONSTS(5)*CONSTS(523)
ALGBRC(55) = ALGBRC(52)/ALGBRC(53)
ALGBRC(51) = CONSTS(453)*CONSTS(6)*ALGBRC(50)*CONSTS(5)*CONSTS(523)
ALGBRC(54) = ALGBRC(51)+CONSTS(584)
ALGBRC(56) = ALGBRC(52)/ALGBRC(54)
ALGBRC(58) = ALGBRC(55)+ALGBRC(56)+CONSTS(11)
ALGBRC(59) = ALGBRC(58)+ALGBRC(57)
ALGBRC(60) = ALGBRC(59) - ALGBRC(34)
ALGBRC(61) = ALGBRC(42) - ALGBRC(59)
ALGBRC(62) = ALGBRC(45)/ALGBRC(59)
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END