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 56 entries in the algebraic variable array.
C There are a total of 4 entries in each of the rate and state variable arrays.
C There are a total of 578 entries in the constant variable array.
C
C
C CONSTS(1) is PA in component input_values (mmHg).
C CONSTS(2) is AUM in component input_values (dimensionless).
C CONSTS(3) is ANM in component input_values (dimensionless).
C CONSTS(4) is MYOGRS in component input_values (dimensionless).
C CONSTS(5) is PAMKRN in component input_values (dimensionless).
C CONSTS(6) is HM1 in component input_values (dimensionless).
C CONSTS(7) is PPC in component input_values (mmHg).
C CONSTS(8) is VTW in component input_values (litre).
C CONSTS(9) is CNA in component input_values (monovalent_mEq_per_litre).
C CONSTS(10) is ADHMK in component input_values (dimensionless).
C CONSTS(11) is AMNA in component input_values (dimensionless).
C CONSTS(12) is CKE in component input_values (monovalent_mEq_per_litre).
C CONSTS(13) is AMK in component input_values (dimensionless).
C CONSTS(14) is ANPX in component input_values (dimensionless).
C CONSTS(15) is A1K in component parameter_values (minute).
C CONSTS(16) is A2K in component parameter_values (minute).
C CONSTS(17) is A3K in component parameter_values (minute).
C CONSTS(18) is A4K in component parameter_values (minute).
C CONSTS(19) is A4K2 in component parameter_values (minute).
C CONSTS(20) is AARK in component parameter_values (mmHg_minute_per_L).
C CONSTS(21) is AARLL in component parameter_values (mmHg_minute_per_L).
C CONSTS(22) is ADHINF in component parameter_values (dimensionless).
C CONSTS(23) is ADHKLL in component parameter_values (dimensionless).
C CONSTS(24) is ADHKUL in component parameter_values (dimensionless).
C CONSTS(25) is ADHPAM in component parameter_values (per_mmHg2).
C CONSTS(26) is ADHPUL in component parameter_values (mmHg).
C CONSTS(27) is ADHTC in component parameter_values (minute).
C CONSTS(28) is ADHVLL in component parameter_values (dimensionless).
C CONSTS(29) is ADHVUL in component parameter_values (dimensionless).
C CONSTS(30) is AH10 in component parameter_values (dimensionless).
C CONSTS(31) is AH11 in component parameter_values (minute).
C CONSTS(32) is AH9 in component parameter_values (per_mmHg).
C CONSTS(33) is AHMNAR in component parameter_values (dimensionless).
C CONSTS(34) is AHTHM in component parameter_values (dimensionless).
C CONSTS(35) is ALCLK in component parameter_values (dimensionless).
C CONSTS(36) is ALDINF in component parameter_values (dimensionless).
C CONSTS(37) is ALDKNS in component parameter_values (dimensionless).
C CONSTS(38) is ALDMM in component parameter_values (dimensionless).
C CONSTS(39) is AM1LL in component parameter_values (dimensionless).
C CONSTS(40) is AM1UL in component parameter_values (dimensionless).
C CONSTS(41) is AMCSNS in component parameter_values (dimensionless).
C CONSTS(42) is AMKM in component parameter_values (dimensionless).
C CONSTS(43) is AMKMUL in component parameter_values (dimensionless).
C CONSTS(44) is AMM1LM in component parameter_values (dimensionless).
C CONSTS(45) is AMM4 in component parameter_values (dimensionless).
C CONSTS(46) is AMNALL in component parameter_values (dimensionless).
C CONSTS(47) is AMNAM in component parameter_values (dimensionless).
C CONSTS(48) is AMNAUL in component parameter_values (dimensionless).
C CONSTS(49) is AMT in component parameter_values (minute).
C CONSTS(50) is ANCSNS in component parameter_values (dimensionless).
C CONSTS(51) is ANGINF in component parameter_values (dimensionless).
C CONSTS(52) is ANGKNS in component parameter_values (dimensionless).
C CONSTS(53) is ANMALD in component parameter_values (dimensionless).
C CONSTS(54) is ANMAM in component parameter_values (dimensionless).
C CONSTS(55) is ANMARL in component parameter_values (dimensionless).
C CONSTS(56) is ANMEM in component parameter_values (dimensionless).
C CONSTS(57) is ANMKEL in component parameter_values (dimensionless).
C CONSTS(58) is ANMKEM in component parameter_values (dimensionless).
C CONSTS(59) is ANMLL in component parameter_values (dimensionless).
C CONSTS(60) is ANMNAM in component parameter_values (dimensionless).
C CONSTS(61) is ANMSLT in component parameter_values (dimensionless).
C CONSTS(62) is ANMTM in component parameter_values (dimensionless).
C CONSTS(63) is ANMUL in component parameter_values (dimensionless).
C CONSTS(64) is ANPTC in component parameter_values (minute).
C CONSTS(65) is ANPXAF in component parameter_values (mmHg_minute_per_L).
C CONSTS(66) is ANPXUL in component parameter_values (dimensionless).
C CONSTS(67) is ANT in component parameter_values (minute).
C CONSTS(68) is ANULL in component parameter_values (dimensionless).
C CONSTS(69) is ANUM in component parameter_values (dimensionless).
C CONSTS(70) is ANUVM in component parameter_values (dimensionless).
C CONSTS(71) is ANV in component parameter_values (minute).
C CONSTS(72) is ANXM in component parameter_values (dimensionless).
C CONSTS(73) is ANY in component parameter_values (litre).
C CONSTS(74) is AR1LM in component parameter_values (dimensionless).
C CONSTS(75) is AR2LM in component parameter_values (dimensionless).
C CONSTS(76) is ARF in component parameter_values (dimensionless).
C CONSTS(77) is AUC1 in component parameter_values (dimensionless).
C CONSTS(78) is AUDMP in component parameter_values (minute).
C CONSTS(79) is AUK in component parameter_values (per_minute).
C CONSTS(80) is AUL in component parameter_values (litre).
C CONSTS(81) is AULPM in component parameter_values (dimensionless).
C CONSTS(82) is AUM1 in component parameter_values (dimensionless).
C CONSTS(83) is AUM2 in component parameter_values (dimensionless).
C CONSTS(84) is AUMAX in component parameter_values (dimensionless).
C CONSTS(85) is AUMIN in component parameter_values (dimensionless).
C CONSTS(86) is AUMK1 in component parameter_values (dimensionless).
C CONSTS(87) is AUN1 in component parameter_values (dimensionless).
C CONSTS(88) is AURG in component parameter_values (dimensionless).
C CONSTS(89) is AUS in component parameter_values (dimensionless).
C CONSTS(90) is AUSLP in component parameter_values (dimensionless).
C CONSTS(91) is AUTO in component parameter_values (dimensionless).
C CONSTS(92) is AUTOGN in component parameter_values (dimensionless).
C CONSTS(93) is AUTOSN in component parameter_values (dimensionless).
C CONSTS(94) is AUX in component parameter_values (dimensionless).
C CONSTS(95) is BAROTC in component parameter_values (minute).
C CONSTS(96) is CFC in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(97) is CKEEX in component parameter_values (dimensionless).
C CONSTS(98) is CMPTSS in component parameter_values (dimensionless).
C CONSTS(99) is CN2 in component parameter_values (per_mmHg).
C CONSTS(100) is CN7 in component parameter_values (dimensionless).
C CONSTS(101) is CNR in component parameter_values (monovalent_mEq_per_litre).
C CONSTS(102) is CPF in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(103) is CPK in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(104) is CPR in component parameter_values (gram_per_L).
C CONSTS(105) is CV in component parameter_values (L_per_mmHg).
C CONSTS(106) is CRRFLX in component parameter_values (mmHg).
C CONSTS(107) is DIURET in component parameter_values (dimensionless).
C CONSTS(108) is DHDTR in component parameter_values (per_mmHg_per_minute).
C CONSTS(109) is DTNAR in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(110) is DTNARL in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(111) is EARK in component parameter_values (mmHg_minute_per_L).
C CONSTS(112) is EARLL in component parameter_values (mmHg_minute_per_L).
C CONSTS(113) is EFAFR in component parameter_values (dimensionless).
C CONSTS(114) is EXC in component parameter_values (dimensionless).
C CONSTS(115) is EXCML in component parameter_values (dimensionless).
C CONSTS(116) is EXCXP in component parameter_values (dimensionless).
C CONSTS(117) is EXCXP2 in component parameter_values (dimensionless).
C CONSTS(118) is EXE in component parameter_values (mmHg).
C CONSTS(119) is GCOPF in component parameter_values (per_mmHg).
C CONSTS(120) is GFLC in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(121) is GFNDMP in component parameter_values (dimensionless).
C CONSTS(122) is GFNLL in component parameter_values (L_per_minute).
C CONSTS(123) is GPPD in component parameter_values (dimensionless).
C CONSTS(124) is HM6 in component parameter_values (mmHg).
C CONSTS(125) is HM8 in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(126) is HTAUML in component parameter_values (dimensionless).
C CONSTS(127) is KORGN in component parameter_values (dimensionless).
C CONSTS(128) is KORNER in component parameter_values (dimensionless).
C CONSTS(129) is KORREN in component parameter_values (dimensionless).
C CONSTS(130) is KORTC in component parameter_values (minute).
C CONSTS(131) is LPDE in component parameter_values (dimensionless).
C CONSTS(132) is LPK in component parameter_values (L_per_minute).
C CONSTS(133) is LPPR in component parameter_values (gram_per_minute).
C CONSTS(134) is MDFL1 in component parameter_values (dimensionless).
C CONSTS(135) is MDFLKM in component parameter_values (monovalent_mEq_per_litre).
C CONSTS(136) is MDFLW2 in component parameter_values (dimensionless).
C CONSTS(137) is MDFLWX in component parameter_values (dimensionless).
C CONSTS(138) is MDMP in component parameter_values (dimensionless).
C CONSTS(139) is O2A in component parameter_values (dimensionless).
C CONSTS(140) is O2CHMO in component parameter_values (per_mmHg).
C CONSTS(141) is OMM in component parameter_values (mL_per_minute).
C CONSTS(142) is PA4 in component parameter_values (dimensionless).
C CONSTS(143) is PAEX in component parameter_values (dimensionless).
C CONSTS(144) is PCE in component parameter_values (dimensionless).
C CONSTS(145) is PCR in component parameter_values (mmHg).
C CONSTS(146) is PR1LL in component parameter_values (mmHg).
C CONSTS(147) is PGHF in component parameter_values (L_mmHg_per_gram).
C CONSTS(148) is PK2 in component parameter_values (mmHg_per_mL).
C CONSTS(149) is PL2 in component parameter_values (L_mL_per_minute_per_mmHg).
C CONSTS(150) is PLDF in component parameter_values (mmHg).
C CONSTS(151) is PM5 in component parameter_values (per_mmHg).
C CONSTS(152) is PO2ADV in component parameter_values (dimensionless).
C CONSTS(153) is POK in component parameter_values (dimensionless).
C CONSTS(154) is POM in component parameter_values (dimensionless).
C CONSTS(155) is POM2 in component parameter_values (dimensionless).
C CONSTS(156) is PON in component parameter_values (dimensionless).
C CONSTS(157) is POR in component parameter_values (mmHg).
C CONSTS(158) is POZ in component parameter_values (dimensionless).
C CONSTS(159) is PXTP in component parameter_values (mmHg).
C CONSTS(160) is QAOLM in component parameter_values (dimensionless).
C CONSTS(161) is QRF in component parameter_values (L_per_minute).
C CONSTS(162) is RABSC in component parameter_values (per_mmHg).
C CONSTS(163) is RAM in component parameter_values (mmHg_minute_per_L).
C CONSTS(164) is RAPRSP in component parameter_values (mmHg).
C CONSTS(165) is RAR in component parameter_values (mmHg_minute_per_L).
C CONSTS(166) is RCDFDP in component parameter_values (minute).
C CONSTS(167) is RCDFPC in component parameter_values (dimensionless).
C CONSTS(168) is REK in component parameter_values (dimensionless).
C CONSTS(169) is RFABKM in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(170) is RFABDM in component parameter_values (dimensionless).
C CONSTS(171) is RFABDP in component parameter_values (dimensionless).
C CONSTS(172) is RFCDFT in component parameter_values (dimensionless).
C CONSTS(173) is RKC in component parameter_values (per_minute).
C CONSTS(174) is RNAGTC in component parameter_values (minute).
C CONSTS(175) is RNAUAD in component parameter_values (per_minute).
C CONSTS(176) is RPS in component parameter_values (dimensionless).
C CONSTS(177) is RTPPR in component parameter_values (dimensionless).
C CONSTS(178) is RTPPRS in component parameter_values (mmHg).
C CONSTS(179) is RTSPRS in component parameter_values (mmHg).
C CONSTS(180) is RVRS in component parameter_values (mmHg_minute_per_L).
C CONSTS(181) is RVSM in component parameter_values (mmHg_minute_per_L).
C CONSTS(182) is SR in component parameter_values (dimensionless).
C CONSTS(183) is SR2 in component parameter_values (dimensionless).
C CONSTS(184) is SRK in component parameter_values (minute).
C CONSTS(185) is SRK2 in component parameter_values (minute).
C CONSTS(186) is STA in component parameter_values (dimensionless).
C CONSTS(187) is TENSGN in component parameter_values (dimensionless).
C CONSTS(188) is TENSTC in component parameter_values (minute).
C CONSTS(189) is TIMETR in component parameter_values (minute).
C CONSTS(190) is TRNSFS in component parameter_values (dimensionless).
C CONSTS(191) is TRPL in component parameter_values (L_per_minute).
C CONSTS(192) is TSSLML in component parameter_values (dimensionless).
C CONSTS(193) is TSSLTC in component parameter_values (per_minute).
C CONSTS(194) is TVDDL in component parameter_values (minute).
C CONSTS(195) is U in component parameter_values (dimensionless).
C CONSTS(196) is VID1 in component parameter_values (L_per_minute).
C CONSTS(197) is VIDML in component parameter_values (litre2_per_monovalent_mEq_per_minute).
C CONSTS(198) is VNTSTM in component parameter_values (dimensionless).
C CONSTS(199) is VP1 in component parameter_values (dimensionless).
C CONSTS(200) is VPTISS in component parameter_values (litre).
C CONSTS(201) is VTSF in component parameter_values (litre).
C CONSTS(202) is VV9 in component parameter_values (litre).
C CONSTS(203) is X in component parameter_values (dimensionless).
C CONSTS(204) is Z in component parameter_values (dimensionless).
C CONSTS(205) is Z4 in component parameter_values (dimensionless).
C CONSTS(206) is Z5 in component parameter_values (dimensionless).
C CONSTS(207) is Z6 in component parameter_values (dimensionless).
C CONSTS(208) is Z7 in component parameter_values (dimensionless).
C CONSTS(209) is Z10 in component parameter_values (mmHg).
C CONSTS(210) is Z11 in component parameter_values (per_mmHg2).
C CONSTS(211) is Z12 in component parameter_values (dimensionless).
C CONSTS(212) is Z13 in component parameter_values (dimensionless).
C CONSTS(213) is Z14 in component parameter_values (dimensionless).
C CONSTS(214) is Z16 in component parameter_values (dimensionless).
C CONSTS(215) is Z18 in component parameter_values (dimensionless).
C CONSTS(216) is Z19 in component parameter_values (dimensionless).
C CONSTS(217) is KID in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(218) is URFORM in component parameter_values (mOsm_per_minute).
C CONSTS(219) is AUV in component parameter_values (dimensionless).
C CONSTS(220) is AUY in component parameter_values (dimensionless).
C CONSTS(221) is O2M in component parameter_values (mL_per_minute).
C CONSTS(222) is NID in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(223) is ATRFBM in component parameter_values (dimensionless).
C CONSTS(224) is GBL in component parameter_values (mmHg).
C CONSTS(225) is RNAUGN in component parameter_values (minute_per_L).
C CONSTS(226) is RNAULL in component parameter_values (dimensionless).
C CONSTS(227) is RNAUUL in component parameter_values (dimensionless).
C CONSTS(228) is HSR in component parameter_values (dimensionless).
C CONSTS(229) is HSL in component parameter_values (dimensionless).
C CONSTS(230) is HMK in component parameter_values (dimensionless).
C CONSTS(231) is PO2AMB in component parameter_values (mmHg).
C CONSTS(232) is HKM in component parameter_values (dimensionless).
C CONSTS(233) is DR in component parameter_values (L_per_minute).
C CONSTS(234) is ANPKNS in component parameter_values (dimensionless).
C CONSTS(235) is RMULT1 in component parameter_values (dimensionless).
C CONSTS(236) is GLPCA in component parameter_values (mmHg).
C CONSTS(237) is RFN in component parameter_values (L_per_minute).
C CONSTS(238) is RFABX in component parameter_values (dimensionless).
C CONSTS(239) is ANPINF in component parameter_values (dimensionless).
C CONSTS(240) is FIS in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(241) is CPPD in component parameter_values (gram_per_L).
C CONSTS(242) is TRRBC in component parameter_values (L_per_minute).
C CONSTS(243) is HYL in component parameter_values (gram).
C CONSTS(244) is ATRVM in component parameter_values (litre).
C CONSTS(245) is ADHC_init in component state_variables (dimensionless).
C CONSTS(246) is AHY_init in component state_variables (dimensionless).
C CONSTS(247) is AMC_init in component state_variables (dimensionless).
C CONSTS(248) is AMM1T_init in component state_variables (dimensionless).
C CONSTS(249) is AMM2_init in component state_variables (dimensionless).
C CONSTS(250) is ANC_init in component state_variables (dimensionless).
C CONSTS(251) is ANPC_init in component state_variables (dimensionless).
C CONSTS(252) is ANX1_init in component state_variables (dimensionless).
C CONSTS(253) is AR1T_init in component state_variables (dimensionless).
C CONSTS(254) is AR2T_init in component state_variables (dimensionless).
C CONSTS(255) is AR3T_init in component state_variables (dimensionless).
C CONSTS(256) is AU1_init in component state_variables (dimensionless).
C CONSTS(257) is AU4_init in component state_variables (dimensionless).
C CONSTS(258) is AU6_init in component state_variables (dimensionless).
C CONSTS(259) is DTKA_init in component state_variables (monovalent_mEq_per_minute).
C CONSTS(260) is HMD1_init in component state_variables (dimensionless).
C CONSTS(261) is HPL_init in component state_variables (dimensionless).
C CONSTS(262) is HPR_init in component state_variables (dimensionless).
C CONSTS(263) is KTOT_init in component state_variables (monovalent_mEq).
C CONSTS(264) is NAE_init in component state_variables (monovalent_mEq).
C CONSTS(265) is O2VAD1_init in component state_variables (dimensionless).
C CONSTS(266) is OSV_init in component state_variables (dimensionless).
C CONSTS(267) is OVA_init in component state_variables (mL_per_L).
C CONSTS(268) is OVS_init in component state_variables (dimensionless).
C CONSTS(269) is PAR1_init in component state_variables (mmHg).
C CONSTS(270) is PC_init in component state_variables (mmHg).
C CONSTS(271) is PLUR_init in component state_variables (mOsm).
C CONSTS(272) is PPR1_init in component state_variables (gram).
C CONSTS(273) is PRP_init in component state_variables (gram).
C CONSTS(274) is QO2T_init in component state_variables (mL).
C CONSTS(275) is QOM1_init in component state_variables (mL).
C CONSTS(276) is RNAUG3_init in component state_variables (dimensionless).
C CONSTS(277) is TSP_init in component state_variables (gram).
C CONSTS(278) is TVD_init in component state_variables (L_per_minute).
C CONSTS(279) is VAS1_init in component state_variables (litre).
C CONSTS(280) is VIC_init in component state_variables (litre).
C CONSTS(281) is VLA1_init in component state_variables (litre).
C CONSTS(282) is VPA1_init in component state_variables (litre).
C CONSTS(283) is VPF1_init in component state_variables (litre).
C CONSTS(284) is VP_init in component state_variables (litre).
C CONSTS(285) is VRA1_init in component state_variables (litre).
C CONSTS(286) is VRC_init in component state_variables (litre).
C CONSTS(287) is VTS2_init in component state_variables (litre).
C CONSTS(288) is VTW_init in component state_variables (litre).
C CONSTS(289) is VV6_init in component state_variables (litre).
C CONSTS(290) is VV7_init in component state_variables (litre).
C CONSTS(291) is VVS1_init in component state_variables (litre).
C VOI is time in component environment (minute).
C ALGBRC(25) is RBF in component actual_renal_blood_flow (L_per_minute).
C ALGBRC(3) is MDFLW in component proximal_tubular_and_macula_densa_flow (L_per_minute).
C ALGBRC(45) is NOD in component actual_Na_excretion_rate (monovalent_mEq_per_minute).
C ALGBRC(48) is KOD in component actual_K_excretion_rate (monovalent_mEq_per_minute).
C ALGBRC(56) is VUD in component actual_urine_volume (L_per_minute).
C CONSTS(292) is A1K in component parameter_values (minute).
C CONSTS(293) is A2K in component parameter_values (minute).
C CONSTS(294) is A3K in component parameter_values (minute).
C CONSTS(295) is A4K in component parameter_values (minute).
C CONSTS(296) is A4K2 in component parameter_values (minute).
C CONSTS(297) is AARK in component parameter_values (mmHg_minute_per_L).
C CONSTS(298) is AARLL in component parameter_values (mmHg_minute_per_L).
C CONSTS(299) is ADHINF in component parameter_values (dimensionless).
C CONSTS(300) is ADHKLL in component parameter_values (dimensionless).
C CONSTS(301) is ADHKUL in component parameter_values (dimensionless).
C CONSTS(302) is ADHPAM in component parameter_values (per_mmHg2).
C CONSTS(303) is ADHPUL in component parameter_values (mmHg).
C CONSTS(304) is ADHTC in component parameter_values (minute).
C CONSTS(305) is ADHVLL in component parameter_values (dimensionless).
C CONSTS(306) is ADHVUL in component parameter_values (dimensionless).
C CONSTS(307) is AH10 in component parameter_values (dimensionless).
C CONSTS(308) is AH11 in component parameter_values (minute).
C CONSTS(309) is AH9 in component parameter_values (per_mmHg).
C CONSTS(310) is AHMNAR in component parameter_values (dimensionless).
C CONSTS(311) is AHTHM in component parameter_values (dimensionless).
C CONSTS(312) is ALCLK in component parameter_values (dimensionless).
C CONSTS(313) is ALDINF in component parameter_values (dimensionless).
C CONSTS(314) is ALDKNS in component parameter_values (dimensionless).
C CONSTS(315) is ALDMM in component parameter_values (dimensionless).
C CONSTS(316) is AM1LL in component parameter_values (dimensionless).
C CONSTS(317) is AM1UL in component parameter_values (dimensionless).
C CONSTS(318) is AMCSNS in component parameter_values (dimensionless).
C CONSTS(319) is AMKM in component parameter_values (dimensionless).
C CONSTS(320) is AMKMUL in component parameter_values (dimensionless).
C CONSTS(321) is AMM1LM in component parameter_values (dimensionless).
C CONSTS(322) is AMM4 in component parameter_values (dimensionless).
C CONSTS(323) is AMNALL in component parameter_values (dimensionless).
C CONSTS(324) is AMNAM in component parameter_values (dimensionless).
C CONSTS(325) is AMNAUL in component parameter_values (dimensionless).
C CONSTS(326) is AMT in component parameter_values (minute).
C CONSTS(327) is ANCSNS in component parameter_values (dimensionless).
C CONSTS(328) is ANGINF in component parameter_values (dimensionless).
C CONSTS(329) is ANGKNS in component parameter_values (dimensionless).
C CONSTS(330) is ANMALD in component parameter_values (dimensionless).
C CONSTS(331) is ANMAM in component parameter_values (dimensionless).
C CONSTS(332) is ANMARL in component parameter_values (dimensionless).
C CONSTS(333) is ANMEM in component parameter_values (dimensionless).
C CONSTS(334) is ANMKEL in component parameter_values (dimensionless).
C CONSTS(335) is ANMKEM in component parameter_values (dimensionless).
C CONSTS(336) is ANMLL in component parameter_values (dimensionless).
C CONSTS(337) is ANMNAM in component parameter_values (dimensionless).
C CONSTS(338) is ANMSLT in component parameter_values (dimensionless).
C CONSTS(339) is ANMTM in component parameter_values (dimensionless).
C CONSTS(340) is ANMUL in component parameter_values (dimensionless).
C CONSTS(341) is ANPTC in component parameter_values (minute).
C CONSTS(342) is ANPXAF in component parameter_values (mmHg_minute_per_L).
C CONSTS(343) is ANPXUL in component parameter_values (dimensionless).
C CONSTS(344) is ANT in component parameter_values (minute).
C CONSTS(345) is ANULL in component parameter_values (dimensionless).
C CONSTS(346) is ANUM in component parameter_values (dimensionless).
C CONSTS(347) is ANUVM in component parameter_values (dimensionless).
C CONSTS(348) is ANV in component parameter_values (minute).
C CONSTS(349) is ANXM in component parameter_values (dimensionless).
C CONSTS(350) is ANY in component parameter_values (litre).
C CONSTS(351) is AR1LM in component parameter_values (dimensionless).
C CONSTS(352) is AR2LM in component parameter_values (dimensionless).
C CONSTS(353) is ARF in component parameter_values (dimensionless).
C CONSTS(354) is AUC1 in component parameter_values (dimensionless).
C CONSTS(355) is AUDMP in component parameter_values (minute).
C CONSTS(356) is AUK in component parameter_values (per_minute).
C CONSTS(357) is AUL in component parameter_values (litre).
C CONSTS(358) is AULPM in component parameter_values (dimensionless).
C CONSTS(359) is AUM1 in component parameter_values (dimensionless).
C CONSTS(360) is AUM2 in component parameter_values (dimensionless).
C CONSTS(361) is AUMAX in component parameter_values (dimensionless).
C CONSTS(362) is AUMIN in component parameter_values (dimensionless).
C CONSTS(363) is AUMK1 in component parameter_values (dimensionless).
C CONSTS(364) is AUN1 in component parameter_values (dimensionless).
C CONSTS(365) is AURG in component parameter_values (dimensionless).
C CONSTS(366) is AUS in component parameter_values (dimensionless).
C CONSTS(367) is AUSLP in component parameter_values (dimensionless).
C CONSTS(368) is AUTO in component parameter_values (dimensionless).
C CONSTS(369) is AUTOGN in component parameter_values (dimensionless).
C CONSTS(370) is AUTOSN in component parameter_values (dimensionless).
C CONSTS(371) is AUX in component parameter_values (dimensionless).
C CONSTS(372) is BAROTC in component parameter_values (minute).
C CONSTS(373) is CFC in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(374) is CKEEX in component parameter_values (dimensionless).
C CONSTS(375) is CMPTSS in component parameter_values (dimensionless).
C CONSTS(376) is CN2 in component parameter_values (per_mmHg).
C CONSTS(377) is CN7 in component parameter_values (dimensionless).
C CONSTS(378) is CNR in component parameter_values (monovalent_mEq_per_litre).
C CONSTS(379) is CPF in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(380) is CPK in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(381) is CPR in component parameter_values (gram_per_L).
C CONSTS(382) is CV in component parameter_values (L_per_mmHg).
C CONSTS(383) is CRRFLX in component parameter_values (mmHg).
C CONSTS(384) is DIURET in component parameter_values (dimensionless).
C CONSTS(385) is DHDTR in component parameter_values (per_mmHg_per_minute).
C CONSTS(386) is DTNAR in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(387) is DTNARL in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(388) is EARK in component parameter_values (mmHg_minute_per_L).
C CONSTS(389) is EARLL in component parameter_values (mmHg_minute_per_L).
C CONSTS(390) is EFAFR in component parameter_values (dimensionless).
C CONSTS(391) is EXC in component parameter_values (dimensionless).
C CONSTS(392) is EXCML in component parameter_values (dimensionless).
C CONSTS(393) is EXCXP in component parameter_values (dimensionless).
C CONSTS(394) is EXCXP2 in component parameter_values (dimensionless).
C CONSTS(395) is EXE in component parameter_values (mmHg).
C CONSTS(396) is GCOPF in component parameter_values (per_mmHg).
C CONSTS(397) is GFLC in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(398) is GFNDMP in component parameter_values (dimensionless).
C CONSTS(399) is GFNLL in component parameter_values (L_per_minute).
C CONSTS(400) is GPPD in component parameter_values (dimensionless).
C CONSTS(401) is HM6 in component parameter_values (mmHg).
C CONSTS(402) is HM8 in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(403) is HTAUML in component parameter_values (dimensionless).
C CONSTS(404) is KORGN in component parameter_values (dimensionless).
C CONSTS(405) is KORNER in component parameter_values (dimensionless).
C CONSTS(406) is KORREN in component parameter_values (dimensionless).
C CONSTS(407) is KORTC in component parameter_values (minute).
C CONSTS(408) is LPDE in component parameter_values (dimensionless).
C CONSTS(409) is LPK in component parameter_values (L_per_minute).
C CONSTS(410) is LPPR in component parameter_values (gram_per_minute).
C CONSTS(411) is MDFL1 in component parameter_values (dimensionless).
C CONSTS(412) is MDFLKM in component parameter_values (monovalent_mEq_per_litre).
C CONSTS(413) is MDFLW2 in component parameter_values (dimensionless).
C CONSTS(414) is MDFLWX in component parameter_values (dimensionless).
C CONSTS(415) is MDMP in component parameter_values (dimensionless).
C CONSTS(416) is O2A in component parameter_values (dimensionless).
C CONSTS(417) is O2CHMO in component parameter_values (per_mmHg).
C CONSTS(418) is OMM in component parameter_values (mL_per_minute).
C CONSTS(419) is PA4 in component parameter_values (dimensionless).
C CONSTS(420) is PAEX in component parameter_values (dimensionless).
C CONSTS(421) is PCE in component parameter_values (dimensionless).
C CONSTS(422) is PCR in component parameter_values (mmHg).
C CONSTS(423) is PR1LL in component parameter_values (mmHg).
C CONSTS(424) is PGHF in component parameter_values (L_mmHg_per_gram).
C CONSTS(425) is PK2 in component parameter_values (mmHg_per_mL).
C CONSTS(426) is PL2 in component parameter_values (L_mL_per_minute_per_mmHg).
C CONSTS(427) is PLDF in component parameter_values (mmHg).
C CONSTS(428) is PM5 in component parameter_values (per_mmHg).
C CONSTS(429) is PO2ADV in component parameter_values (dimensionless).
C CONSTS(430) is POK in component parameter_values (dimensionless).
C CONSTS(431) is POM in component parameter_values (dimensionless).
C CONSTS(432) is POM2 in component parameter_values (dimensionless).
C CONSTS(433) is PON in component parameter_values (dimensionless).
C CONSTS(434) is POR in component parameter_values (mmHg).
C CONSTS(435) is POZ in component parameter_values (dimensionless).
C CONSTS(436) is PXTP in component parameter_values (mmHg).
C CONSTS(437) is QAOLM in component parameter_values (dimensionless).
C CONSTS(438) is QRF in component parameter_values (L_per_minute).
C CONSTS(439) is RABSC in component parameter_values (per_mmHg).
C CONSTS(440) is RAM in component parameter_values (mmHg_minute_per_L).
C CONSTS(441) is RAPRSP in component parameter_values (mmHg).
C CONSTS(442) is RAR in component parameter_values (mmHg_minute_per_L).
C CONSTS(443) is RCDFDP in component parameter_values (minute).
C CONSTS(444) is RCDFPC in component parameter_values (dimensionless).
C CONSTS(445) is REK in component parameter_values (dimensionless).
C CONSTS(446) is RFABKM in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(447) is RFABDM in component parameter_values (dimensionless).
C CONSTS(448) is RFABDP in component parameter_values (dimensionless).
C CONSTS(449) is RFCDFT in component parameter_values (dimensionless).
C CONSTS(450) is RKC in component parameter_values (per_minute).
C CONSTS(451) is RNAGTC in component parameter_values (minute).
C CONSTS(452) is RNAUAD in component parameter_values (per_minute).
C CONSTS(453) is RPS in component parameter_values (dimensionless).
C CONSTS(454) is RTPPR in component parameter_values (dimensionless).
C CONSTS(455) is RTPPRS in component parameter_values (mmHg).
C CONSTS(456) is RTSPRS in component parameter_values (mmHg).
C CONSTS(457) is RVRS in component parameter_values (mmHg_minute_per_L).
C CONSTS(458) is RVSM in component parameter_values (mmHg_minute_per_L).
C CONSTS(459) is SR in component parameter_values (dimensionless).
C CONSTS(460) is SR2 in component parameter_values (dimensionless).
C CONSTS(461) is SRK in component parameter_values (minute).
C CONSTS(462) is SRK2 in component parameter_values (minute).
C CONSTS(463) is STA in component parameter_values (dimensionless).
C CONSTS(464) is TENSGN in component parameter_values (dimensionless).
C CONSTS(465) is TENSTC in component parameter_values (minute).
C CONSTS(466) is TIMETR in component parameter_values (minute).
C CONSTS(467) is TRNSFS in component parameter_values (dimensionless).
C CONSTS(468) is TRPL in component parameter_values (L_per_minute).
C CONSTS(469) is TSSLML in component parameter_values (dimensionless).
C CONSTS(470) is TSSLTC in component parameter_values (per_minute).
C CONSTS(471) is TVDDL in component parameter_values (minute).
C CONSTS(472) is U in component parameter_values (dimensionless).
C CONSTS(473) is VID1 in component parameter_values (L_per_minute).
C CONSTS(474) is VIDML in component parameter_values (litre2_per_monovalent_mEq_per_minute).
C CONSTS(475) is VNTSTM in component parameter_values (dimensionless).
C CONSTS(476) is VP1 in component parameter_values (dimensionless).
C CONSTS(477) is VPTISS in component parameter_values (litre).
C CONSTS(478) is VTSF in component parameter_values (litre).
C CONSTS(479) is VV9 in component parameter_values (litre).
C CONSTS(480) is X in component parameter_values (dimensionless).
C CONSTS(481) is Z in component parameter_values (dimensionless).
C CONSTS(482) is Z4 in component parameter_values (dimensionless).
C CONSTS(483) is Z5 in component parameter_values (dimensionless).
C CONSTS(484) is Z6 in component parameter_values (dimensionless).
C CONSTS(485) is Z7 in component parameter_values (dimensionless).
C CONSTS(486) is Z10 in component parameter_values (mmHg).
C CONSTS(487) is Z11 in component parameter_values (per_mmHg2).
C CONSTS(488) is Z12 in component parameter_values (dimensionless).
C CONSTS(489) is Z13 in component parameter_values (dimensionless).
C CONSTS(490) is Z14 in component parameter_values (dimensionless).
C CONSTS(491) is Z16 in component parameter_values (dimensionless).
C CONSTS(492) is Z18 in component parameter_values (dimensionless).
C CONSTS(493) is Z19 in component parameter_values (dimensionless).
C CONSTS(494) is KID in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(495) is URFORM in component parameter_values (mOsm_per_minute).
C CONSTS(496) is AUV in component parameter_values (dimensionless).
C CONSTS(497) is AUY in component parameter_values (dimensionless).
C CONSTS(498) is O2M in component parameter_values (mL_per_minute).
C CONSTS(499) is NID in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(500) is ATRFBM in component parameter_values (dimensionless).
C CONSTS(501) is GBL in component parameter_values (mmHg).
C CONSTS(502) is RNAUGN in component parameter_values (minute_per_L).
C CONSTS(503) is RNAULL in component parameter_values (dimensionless).
C CONSTS(504) is RNAUUL in component parameter_values (dimensionless).
C CONSTS(505) is HSR in component parameter_values (dimensionless).
C CONSTS(506) is HSL in component parameter_values (dimensionless).
C CONSTS(507) is HMK in component parameter_values (dimensionless).
C CONSTS(508) is PO2AMB in component parameter_values (mmHg).
C CONSTS(509) is HKM in component parameter_values (dimensionless).
C CONSTS(510) is DR in component parameter_values (L_per_minute).
C CONSTS(511) is ANPKNS in component parameter_values (dimensionless).
C CONSTS(512) is RMULT1 in component parameter_values (dimensionless).
C CONSTS(513) is GLPCA in component parameter_values (mmHg).
C CONSTS(514) is RFN in component parameter_values (L_per_minute).
C CONSTS(515) is RFABX in component parameter_values (dimensionless).
C CONSTS(516) is ANPINF in component parameter_values (dimensionless).
C CONSTS(517) is FIS in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(518) is CPPD in component parameter_values (gram_per_L).
C CONSTS(519) is TRRBC in component parameter_values (L_per_minute).
C CONSTS(520) is HYL in component parameter_values (gram).
C CONSTS(521) is ATRVM in component parameter_values (litre).
C CONSTS(522) is ADHC_init in component state_variables (dimensionless).
C CONSTS(523) is AHY_init in component state_variables (dimensionless).
C CONSTS(524) is AMC_init in component state_variables (dimensionless).
C CONSTS(525) is AMM1T_init in component state_variables (dimensionless).
C CONSTS(526) is AMM2_init in component state_variables (dimensionless).
C CONSTS(527) is ANC_init in component state_variables (dimensionless).
C CONSTS(528) is ANPC_init in component state_variables (dimensionless).
C CONSTS(529) is ANX1_init in component state_variables (dimensionless).
C CONSTS(530) is AR1T_init in component state_variables (dimensionless).
C CONSTS(531) is AR2T_init in component state_variables (dimensionless).
C CONSTS(532) is AR3T_init in component state_variables (dimensionless).
C CONSTS(533) is AU1_init in component state_variables (dimensionless).
C CONSTS(534) is AU4_init in component state_variables (dimensionless).
C CONSTS(535) is AU6_init in component state_variables (dimensionless).
C CONSTS(536) is DTKA_init in component state_variables (monovalent_mEq_per_minute).
C CONSTS(537) is HMD1_init in component state_variables (dimensionless).
C CONSTS(538) is HPL_init in component state_variables (dimensionless).
C CONSTS(539) is HPR_init in component state_variables (dimensionless).
C CONSTS(540) is KTOT_init in component state_variables (monovalent_mEq).
C CONSTS(541) is NAE_init in component state_variables (monovalent_mEq).
C CONSTS(542) is O2VAD1_init in component state_variables (dimensionless).
C CONSTS(543) is OSV_init in component state_variables (dimensionless).
C CONSTS(544) is OVA_init in component state_variables (mL_per_L).
C CONSTS(545) is OVS_init in component state_variables (dimensionless).
C CONSTS(546) is PAR1_init in component state_variables (mmHg).
C CONSTS(547) is PC_init in component state_variables (mmHg).
C CONSTS(548) is PLUR_init in component state_variables (mOsm).
C CONSTS(549) is PPR1_init in component state_variables (gram).
C CONSTS(550) is PRP_init in component state_variables (gram).
C CONSTS(551) is QO2T_init in component state_variables (mL).
C CONSTS(552) is QOM1_init in component state_variables (mL).
C CONSTS(553) is RNAUG3_init in component state_variables (dimensionless).
C CONSTS(554) is TSP_init in component state_variables (gram).
C CONSTS(555) is TVD_init in component state_variables (L_per_minute).
C CONSTS(556) is VAS1_init in component state_variables (litre).
C CONSTS(557) is VIC_init in component state_variables (litre).
C CONSTS(558) is VLA1_init in component state_variables (litre).
C CONSTS(559) is VPA1_init in component state_variables (litre).
C CONSTS(560) is VPF1_init in component state_variables (litre).
C CONSTS(561) is VP_init in component state_variables (litre).
C CONSTS(562) is VRA1_init in component state_variables (litre).
C CONSTS(563) is VRC_init in component state_variables (litre).
C CONSTS(564) is VTS2_init in component state_variables (litre).
C CONSTS(565) is VTW_init in component state_variables (litre).
C CONSTS(566) is VV6_init in component state_variables (litre).
C CONSTS(567) is VV7_init in component state_variables (litre).
C CONSTS(568) is VVS1_init in component state_variables (litre).
C ALGBRC(1) is PAR in component perfusion_pressure (mmHg).
C STATES(1) is PAR1 in component perfusion_pressure (mmHg).
C ALGBRC(4) is RNAUG2 in component renal_autoregulatory_feedback_factor (dimensionless).
C ALGBRC(5) is RNAUG1 in component renal_autoregulatory_feedback_factor (dimensionless).
C ALGBRC(6) is RNAUG1T in component renal_autoregulatory_feedback_factor (dimensionless).
C STATES(2) is RNAUG3 in component renal_autoregulatory_feedback_factor (dimensionless).
C CONSTS(570) is AUMK in component autonomic_effect_on_AAR (dimensionless).
C CONSTS(569) is AUMKT in component autonomic_effect_on_AAR (dimensionless).
C CONSTS(572) is ANMAR in component angiotensin_effect_on_AAR (dimensionless).
C CONSTS(571) is ANMAR1 in component angiotensin_effect_on_AAR (dimensionless).
C ALGBRC(7) is AAR1 in component AAR_calculation (mmHg_minute_per_L).
C ALGBRC(8) is AAR in component atrial_natriuretic_peptide_effect_on_AAR (mmHg_minute_per_L).
C ALGBRC(9) is AART in component atrial_natriuretic_peptide_effect_on_AAR (mmHg_minute_per_L).
C CONSTS(573) is AUMK2 in component autonomic_effect_on_EAR (dimensionless).
C CONSTS(574) is ANMER in component angiotensin_effect_on_EAR (dimensionless).
C ALGBRC(10) is RNAUG4 in component effect_of_renal_autoregulatory_feedback_on_EAR (dimensionless).
C ALGBRC(11) is EAR in component EAR_calculation (mmHg_minute_per_L).
C ALGBRC(12) is EAR1 in component EAR_calculation (mmHg_minute_per_L).
C ALGBRC(13) is RR in component total_renal_resistance (mmHg_minute_per_L).
C ALGBRC(14) is RFN in component normal_renal_blood_flow (L_per_minute).
C ALGBRC(15) is GFN in component glomerular_filtration_rate (L_per_minute).
C ALGBRC(16) is GLPC in component glomerular_colloid_osmotic_pressure (mmHg).
C ALGBRC(17) is EFAFPR in component glomerular_colloid_osmotic_pressure (dimensionless).
C ALGBRC(18) is EFAFPR1 in component glomerular_colloid_osmotic_pressure (dimensionless).
C ALGBRC(19) is GLP in component glomerular_pressure (mmHg).
C ALGBRC(20) is APD in component glomerular_pressure (mmHg).
C ALGBRC(26) is GFR in component glomerular_filtration_rate (L_per_minute).
C ALGBRC(21) is PFL in component glomerular_filtration_rate (mmHg).
C ALGBRC(22) is GFN1 in component glomerular_filtration_rate (L_per_minute).
C ALGBRC(50) is UROD in component actual_urea_excretion_rate (mOsm_per_minute).
C STATES(3) is PLUR in component glomerular_urea_concentration (mOsm).
C ALGBRC(23) is PTFL in component proximal_tubular_and_macula_densa_flow (L_per_minute).
C ALGBRC(24) is MDFLWT in component proximal_tubular_and_macula_densa_flow (L_per_minute).
C ALGBRC(28) is RTSPPC in component renal_tissue_osmotic_pressure (mmHg).
C ALGBRC(27) is RTSPPC1 in component renal_tissue_osmotic_pressure (mmHg).
C ALGBRC(2) is PLURC in component plasma_urea_concentration (mOsm_per_litre).
C ALGBRC(29) is RCPRS in component peritubular_capillary_pressure (mmHg).
C ALGBRC(34) is RFABD in component peritubular_capillary_reabsorption_factor (dimensionless).
C ALGBRC(30) is RABSPR in component peritubular_capillary_reabsorption_factor (mmHg).
C ALGBRC(31) is RFAB1 in component peritubular_capillary_reabsorption_factor (dimensionless).
C ALGBRC(32) is RFAB in component peritubular_capillary_reabsorption_factor (dimensionless).
C ALGBRC(33) is RFABD1 in component peritubular_capillary_reabsorption_factor (dimensionless).
C ALGBRC(35) is DTNAI in component distal_tubular_Na_delivery (monovalent_mEq_per_minute).
C ALGBRC(37) is DTNARA in component Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C ALGBRC(36) is DTNARA1 in component Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C CONSTS(576) is DTNANG in component angiotensin_induced_Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C CONSTS(575) is DTNANG1 in component angiotensin_induced_Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C ALGBRC(38) is DTKI in component distal_tubular_K_delivery (monovalent_mEq_per_minute).
C ALGBRC(39) is RFABK in component effect_of_physical_forces_on_distal_K_reabsorption (monovalent_mEq_per_minute).
C ALGBRC(41) is MDFLK in component effect_of_fluid_flow_on_distal_K_reabsorption (monovalent_mEq_per_minute).
C ALGBRC(40) is MDFLK1 in component effect_of_fluid_flow_on_distal_K_reabsorption (monovalent_mEq_per_minute).
C ALGBRC(47) is KODN in component normal_K_excretion (monovalent_mEq_per_minute).
C ALGBRC(55) is VUDN in component normal_urine_volume (L_per_minute).
C STATES(4) is DTKA in component K_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C ALGBRC(42) is DTKSC in component K_secretion_from_distal_tubules (monovalent_mEq_per_minute).
C CONSTS(577) is ANMKE1 in component K_secretion_from_distal_tubules (dimensionless).
C CONSTS(578) is ANMKE in component K_secretion_from_distal_tubules (dimensionless).
C ALGBRC(44) is NODN in component normal_Na_excretion (monovalent_mEq_per_minute).
C ALGBRC(43) is NODN1 in component normal_Na_excretion (monovalent_mEq_per_minute).
C ALGBRC(46) is KODN1 in component normal_K_excretion (monovalent_mEq_per_minute).
C ALGBRC(49) is DTURI in component normal_urea_excretion (mOsm_per_minute).
C ALGBRC(51) is OSMOPN1 in component normal_osmolar_and_water_excretion (mOsm_per_minute).
C ALGBRC(52) is OSMOPN in component normal_osmolar_and_water_excretion (mOsm_per_minute).
C ALGBRC(53) is OSMOP1T in component normal_urine_volume (mOsm_per_minute).
C ALGBRC(54) is OSMOP1 in component normal_urine_volume (mOsm_per_minute).
C RATES(1) is d/dt PAR1 in component perfusion_pressure (mmHg).
C RATES(2) is d/dt RNAUG3 in component renal_autoregulatory_feedback_factor (dimensionless).
C RATES(3) is d/dt PLUR in component glomerular_urea_concentration (mOsm).
C RATES(4) is d/dt DTKA in component K_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 103.525
      CONSTS(2) = 1.00066
      CONSTS(3) = 0.987545
      CONSTS(4) = 1.0
      CONSTS(5) = 1.0
      CONSTS(6) = 0.39984739
      CONSTS(7) = 29.9941
      CONSTS(8) = 39.8952
      CONSTS(9) = 142.035
      CONSTS(10) = 1.0
      CONSTS(11) = 1.0
      CONSTS(12) = 4.44092
      CONSTS(13) = 1.037
      CONSTS(14) = 1.0
      CONSTS(15) = 0.5
      CONSTS(16) = 60
      CONSTS(17) = 40000
      CONSTS(18) = 0.1
      CONSTS(19) = 40000
      CONSTS(20) = 1
      CONSTS(21) = 4
      CONSTS(22) = 0
      CONSTS(23) = 0.2
      CONSTS(24) = 5
      CONSTS(25) = 0.3
      CONSTS(26) = 85
      CONSTS(27) = 15
      CONSTS(28) = 0.93617
      CONSTS(29) = 2.5
      CONSTS(30) = 0.333
      CONSTS(31) = 1000
      CONSTS(32) = 1
      CONSTS(33) = 0.3
      CONSTS(34) = 2
      CONSTS(35) = 0.3
      CONSTS(36) = 0
      CONSTS(37) = 0
      CONSTS(38) = 2.5
      CONSTS(39) = 0
      CONSTS(40) = 5
      CONSTS(41) = 0.65
      CONSTS(42) = 0.5
      CONSTS(43) = 12
      CONSTS(44) = 0.005
      CONSTS(45) = 0.005
      CONSTS(46) = 0.04
      CONSTS(47) = 0.8
      CONSTS(48) = 15
      CONSTS(49) = 60
      CONSTS(50) = 0.4
      CONSTS(51) = 0
      CONSTS(52) = 0
      CONSTS(53) = 2.5
      CONSTS(54) = 1.4
      CONSTS(55) = 0.86
      CONSTS(56) = 1.6
      CONSTS(57) = 0.3
      CONSTS(58) = 2
      CONSTS(59) = 0.7
      CONSTS(60) = 1
      CONSTS(61) = 2
      CONSTS(62) = 1.5
      CONSTS(63) = 1.8
      CONSTS(64) = 4
      CONSTS(65) = 1.5
      CONSTS(66) = 10
      CONSTS(67) = 12
      CONSTS(68) = 0.8
      CONSTS(69) = 6
      CONSTS(70) = 0
      CONSTS(71) = 5000
      CONSTS(72) = 0
      CONSTS(73) = -0.2
      CONSTS(74) = 0.005
      CONSTS(75) = 0.005
      CONSTS(76) = 0.5
      CONSTS(77) = 0.3
      CONSTS(78) = 0.3
      CONSTS(79) = 0.004
      CONSTS(80) = 1.5
      CONSTS(81) = 0
      CONSTS(82) = 3
      CONSTS(83) = 1
      CONSTS(84) = 5.0
      CONSTS(85) = 0.4
      CONSTS(86) = 0.3
      CONSTS(87) = 0.5
      CONSTS(88) = 1.0
      CONSTS(89) = 1
      CONSTS(90) = 1.5
      CONSTS(91) = 0
      CONSTS(92) = 7
      CONSTS(93) = 0.9
      CONSTS(94) = 1
      CONSTS(95) = 0.16
      CONSTS(96) = 0.01167
      CONSTS(97) = 4
      CONSTS(98) = 2
      CONSTS(99) = 0.0212
      CONSTS(100) = 0.2
      CONSTS(101) = 139
      CONSTS(102) = 0.0003
      CONSTS(103) = 0.000253
      CONSTS(104) = 40
      CONSTS(105) = 0.1
      CONSTS(106) = 0
      CONSTS(107) = 1
      CONSTS(108) = 0.05
      CONSTS(109) = 0.675
      CONSTS(110) = 1e-06
      CONSTS(111) = 1
      CONSTS(112) = 24
      CONSTS(113) = 0
      CONSTS(114) = 1
      CONSTS(115) = 0.01
      CONSTS(116) = 1.0
      CONSTS(117) = 0.17
      CONSTS(118) = 0
      CONSTS(119) = 0.8092
      CONSTS(120) = 0.0208333
      CONSTS(121) = 3
      CONSTS(122) = 0.001
      CONSTS(123) = 1.0
      CONSTS(124) = 1850
      CONSTS(125) = 4.714e-08
      CONSTS(126) = 0.4
      CONSTS(127) = 1
      CONSTS(128) = 0
      CONSTS(129) = 0
      CONSTS(130) = 10000
      CONSTS(131) = 8
      CONSTS(132) = 2.728e-14
      CONSTS(133) = 0.03
      CONSTS(134) = 10
      CONSTS(135) = 0.667
      CONSTS(136) = 0.005
      CONSTS(137) = 0.01
      CONSTS(138) = 0
      CONSTS(139) = 0.1
      CONSTS(140) = 0.01
      CONSTS(141) = 57.1
      CONSTS(142) = 1
      CONSTS(143) = 2
      CONSTS(144) = 1
      CONSTS(145) = 15
      CONSTS(146) = 0
      CONSTS(147) = -2
      CONSTS(148) = 0.79167
      CONSTS(149) = 1.8
      CONSTS(150) = 4.2
      CONSTS(151) = 30
      CONSTS(152) = 0.001
      CONSTS(153) = 0.1
      CONSTS(154) = 0.04
      CONSTS(155) = 2
      CONSTS(156) = 0.1
      CONSTS(157) = 35
      CONSTS(158) = 2
      CONSTS(159) = 8
      CONSTS(160) = 0.2
      CONSTS(161) = 0.15
      CONSTS(162) = 0.5
      CONSTS(163) = 96.3
      CONSTS(164) = 0
      CONSTS(165) = 30.52
      CONSTS(166) = 2000
      CONSTS(167) = 0
      CONSTS(168) = 1
      CONSTS(169) = 0.03
      CONSTS(170) = 0.3
      CONSTS(171) = 1
      CONSTS(172) = 0
      CONSTS(173) = 5.8e-06
      CONSTS(174) = 15
      CONSTS(175) = 0
      CONSTS(176) = 0
      CONSTS(177) = 0.9
      CONSTS(178) = 15.2
      CONSTS(179) = 6
      CONSTS(180) = 19.167
      CONSTS(181) = 1
      CONSTS(182) = 1
      CONSTS(183) = 1
      CONSTS(184) = 5
      CONSTS(185) = 10000
      CONSTS(186) = 0
      CONSTS(187) = 0
      CONSTS(188) = 0.05
      CONSTS(189) = 0
      CONSTS(190) = 0
      CONSTS(191) = 0
      CONSTS(192) = 0.15
      CONSTS(193) = 0.005
      CONSTS(194) = 30
      CONSTS(195) = 4
      CONSTS(196) = 0.01
      CONSTS(197) = 0.01
      CONSTS(198) = 1
      CONSTS(199) = 0.01
      CONSTS(200) = 0.0175
      CONSTS(201) = 6
      CONSTS(202) = 2.51
      CONSTS(203) = 4
      CONSTS(204) = 1
      CONSTS(205) = 0.5
      CONSTS(206) = 0.3
      CONSTS(207) = 20
      CONSTS(208) = 5
      CONSTS(209) = 45
      CONSTS(210) = 0.01
      CONSTS(211) = 5
      CONSTS(212) = 0.625
      CONSTS(213) = 5
      CONSTS(214) = 5
      CONSTS(215) = 5
      CONSTS(216) = 5
      CONSTS(217) = 0.08
      CONSTS(218) = 0.24
      CONSTS(219) = 0.55
      CONSTS(220) = 0
      CONSTS(221) = 164
      CONSTS(222) = 0.1
      CONSTS(223) = 0
      CONSTS(224) = 0
      CONSTS(225) = 0.6
      CONSTS(226) = 0.3
      CONSTS(227) = 10
      CONSTS(228) = 1
      CONSTS(229) = 1
      CONSTS(230) = 90
      CONSTS(231) = 150
      CONSTS(232) = 0.53333
      CONSTS(233) = 0
      CONSTS(234) = 0
      CONSTS(235) = 1
      CONSTS(236) = 1.0
      CONSTS(237) = 1.22057
      CONSTS(238) = 0.8
      CONSTS(239) = 0
      CONSTS(240) = 0
      CONSTS(241) = 31.9719
      CONSTS(242) = 0
      CONSTS(243) = 60
      CONSTS(244) = 0
      CONSTS(245) = 1.0
      CONSTS(246) = 0.301963
      CONSTS(247) = 1.0
      CONSTS(248) = 1.00269
      CONSTS(249) = 1.09071
      CONSTS(250) = 0.859476
      CONSTS(251) = 1.0
      CONSTS(252) = 0
      CONSTS(253) = 1.02127
      CONSTS(254) = 1.01179
      CONSTS(255) = 1.1448
      CONSTS(256) = 1.00007
      CONSTS(257) = -0.060024
      CONSTS(258) = 1.00132
      CONSTS(259) = 0.0367573
      CONSTS(260) = 1
      CONSTS(261) = 1.00163
      CONSTS(262) = 1.00237
      CONSTS(263) = 3622.54
      CONSTS(264) = 2109.91
      CONSTS(265) = 2.368e-07
      CONSTS(266) = 0.694163
      CONSTS(267) = 204.497
      CONSTS(268) = 0.699673
      CONSTS(269) = 103.525
      CONSTS(270) = 16.9144
      CONSTS(271) = 159.549
      CONSTS(272) = 0.419998
      CONSTS(273) = 216.243
      CONSTS(274) = 72.2362
      CONSTS(275) = 48.0839
      CONSTS(276) = 0
      CONSTS(277) = 279.945
      CONSTS(278) = 0.000980838
      CONSTS(279) = 0.862514
      CONSTS(280) = 25.0404
      CONSTS(281) = 0.379883
      CONSTS(282) = 0.38131
      CONSTS(283) = 0.0123238
      CONSTS(284) = 3.00449
      CONSTS(285) = 0.100043
      CONSTS(286) = 2.00439
      CONSTS(287) = 0
      CONSTS(288) = 39.8952
      CONSTS(289) = 0.0101913
      CONSTS(290) = 0.00366525
      CONSTS(291) = 3.28246
      CONSTS(292) = 0.5
      CONSTS(293) = 60
      CONSTS(294) = 40000
      CONSTS(295) = 0.1
      CONSTS(296) = 40000
      CONSTS(297) = 1
      CONSTS(298) = 4
      CONSTS(299) = 0
      CONSTS(300) = 0.2
      CONSTS(301) = 5
      CONSTS(302) = 0.3
      CONSTS(303) = 85
      CONSTS(304) = 15
      CONSTS(305) = 0.93617
      CONSTS(306) = 2.5
      CONSTS(307) = 0.333
      CONSTS(308) = 1000
      CONSTS(309) = 1
      CONSTS(310) = 0.3
      CONSTS(311) = 2
      CONSTS(312) = 0.3
      CONSTS(313) = 0
      CONSTS(314) = 0
      CONSTS(315) = 2.5
      CONSTS(316) = 0
      CONSTS(317) = 5
      CONSTS(318) = 0.65
      CONSTS(319) = 0.5
      CONSTS(320) = 12
      CONSTS(321) = 0.005
      CONSTS(322) = 0.005
      CONSTS(323) = 0.04
      CONSTS(324) = 0.8
      CONSTS(325) = 15
      CONSTS(326) = 60
      CONSTS(327) = 0.4
      CONSTS(328) = 0
      CONSTS(329) = 0
      CONSTS(330) = 2.5
      CONSTS(331) = 1.4
      CONSTS(332) = 0.86
      CONSTS(333) = 1.6
      CONSTS(334) = 0.3
      CONSTS(335) = 2
      CONSTS(336) = 0.7
      CONSTS(337) = 1
      CONSTS(338) = 2
      CONSTS(339) = 1.5
      CONSTS(340) = 1.8
      CONSTS(341) = 4
      CONSTS(342) = 1.5
      CONSTS(343) = 10
      CONSTS(344) = 12
      CONSTS(345) = 0.8
      CONSTS(346) = 6
      CONSTS(347) = 0
      CONSTS(348) = 5000
      CONSTS(349) = 0
      CONSTS(350) = -0.2
      CONSTS(351) = 0.005
      CONSTS(352) = 0.005
      CONSTS(353) = 0.5
      CONSTS(354) = 0.3
      CONSTS(355) = 0.3
      CONSTS(356) = 0.004
      CONSTS(357) = 1.5
      CONSTS(358) = 0
      CONSTS(359) = 3
      CONSTS(360) = 1
      CONSTS(361) = 5.0
      CONSTS(362) = 0.4
      CONSTS(363) = 0.3
      CONSTS(364) = 0.5
      CONSTS(365) = 1.0
      CONSTS(366) = 1
      CONSTS(367) = 1.5
      CONSTS(368) = 0
      CONSTS(369) = 7
      CONSTS(370) = 0.9
      CONSTS(371) = 1
      CONSTS(372) = 0.16
      CONSTS(373) = 0.01167
      CONSTS(374) = 4
      CONSTS(375) = 2
      CONSTS(376) = 0.0212
      CONSTS(377) = 0.2
      CONSTS(378) = 139
      CONSTS(379) = 0.0003
      CONSTS(380) = 0.000253
      CONSTS(381) = 40
      CONSTS(382) = 0.1
      CONSTS(383) = 0
      CONSTS(384) = 1
      CONSTS(385) = 0.05
      CONSTS(386) = 0.675
      CONSTS(387) = 1e-06
      CONSTS(388) = 1
      CONSTS(389) = 24
      CONSTS(390) = 0
      CONSTS(391) = 1
      CONSTS(392) = 0.01
      CONSTS(393) = 1.0
      CONSTS(394) = 0.17
      CONSTS(395) = 0
      CONSTS(396) = 0.8092
      CONSTS(397) = 0.0208333
      CONSTS(398) = 3
      CONSTS(399) = 0.001
      CONSTS(400) = 1.0
      CONSTS(401) = 1850
      CONSTS(402) = 4.714e-08
      CONSTS(403) = 0.4
      CONSTS(404) = 1
      CONSTS(405) = 0
      CONSTS(406) = 0
      CONSTS(407) = 10000
      CONSTS(408) = 8
      CONSTS(409) = 2.728e-14
      CONSTS(410) = 0.03
      CONSTS(411) = 10
      CONSTS(412) = 0.667
      CONSTS(413) = 0.005
      CONSTS(414) = 0.01
      CONSTS(415) = 0
      CONSTS(416) = 0.1
      CONSTS(417) = 0.01
      CONSTS(418) = 57.1
      CONSTS(419) = 1
      CONSTS(420) = 2
      CONSTS(421) = 1
      CONSTS(422) = 15
      CONSTS(423) = 0
      CONSTS(424) = -2
      CONSTS(425) = 0.79167
      CONSTS(426) = 1.8
      CONSTS(427) = 4.2
      CONSTS(428) = 30
      CONSTS(429) = 0.001
      CONSTS(430) = 0.1
      CONSTS(431) = 0.04
      CONSTS(432) = 2
      CONSTS(433) = 0.1
      CONSTS(434) = 35
      CONSTS(435) = 2
      CONSTS(436) = 8
      CONSTS(437) = 0.2
      CONSTS(438) = 0.15
      CONSTS(439) = 0.5
      CONSTS(440) = 96.3
      CONSTS(441) = 0
      CONSTS(442) = 30.52
      CONSTS(443) = 2000
      CONSTS(444) = 0
      CONSTS(445) = 1
      CONSTS(446) = 0.03
      CONSTS(447) = 0.3
      CONSTS(448) = 1
      CONSTS(449) = 0
      CONSTS(450) = 5.8e-06
      CONSTS(451) = 15
      CONSTS(452) = 0
      CONSTS(453) = 0
      CONSTS(454) = 0.9
      CONSTS(455) = 15.2
      CONSTS(456) = 6
      CONSTS(457) = 19.167
      CONSTS(458) = 1
      CONSTS(459) = 1
      CONSTS(460) = 1
      CONSTS(461) = 5
      CONSTS(462) = 10000
      CONSTS(463) = 0
      CONSTS(464) = 0
      CONSTS(465) = 0.05
      CONSTS(466) = 0
      CONSTS(467) = 0
      CONSTS(468) = 0
      CONSTS(469) = 0.15
      CONSTS(470) = 0.005
      CONSTS(471) = 30
      CONSTS(472) = 4
      CONSTS(473) = 0.01
      CONSTS(474) = 0.01
      CONSTS(475) = 1
      CONSTS(476) = 0.01
      CONSTS(477) = 0.0175
      CONSTS(478) = 6
      CONSTS(479) = 2.51
      CONSTS(480) = 4
      CONSTS(481) = 1
      CONSTS(482) = 0.5
      CONSTS(483) = 0.3
      CONSTS(484) = 20
      CONSTS(485) = 5
      CONSTS(486) = 45
      CONSTS(487) = 0.01
      CONSTS(488) = 5
      CONSTS(489) = 0.625
      CONSTS(490) = 5
      CONSTS(491) = 5
      CONSTS(492) = 5
      CONSTS(493) = 5
      CONSTS(494) = 0.08
      CONSTS(495) = 0.24
      CONSTS(496) = 0.55
      CONSTS(497) = 0
      CONSTS(498) = 164
      CONSTS(499) = 0.1
      CONSTS(500) = 0
      CONSTS(501) = 0
      CONSTS(502) = 0.6
      CONSTS(503) = 0.3
      CONSTS(504) = 10
      CONSTS(505) = 1
      CONSTS(506) = 1
      CONSTS(507) = 90
      CONSTS(508) = 150
      CONSTS(509) = 0.53333
      CONSTS(510) = 0
      CONSTS(511) = 0
      CONSTS(512) = 1
      CONSTS(513) = 1.0
      CONSTS(514) = 1.22057
      CONSTS(515) = 0.8
      CONSTS(516) = 0
      CONSTS(517) = 0
      CONSTS(518) = 31.9719
      CONSTS(519) = 0
      CONSTS(520) = 60
      CONSTS(521) = 0
      CONSTS(522) = 1.0
      CONSTS(523) = 0.301963
      CONSTS(524) = 1.0
      CONSTS(525) = 1.00269
      CONSTS(526) = 1.09071
      CONSTS(527) = 0.859476
      CONSTS(528) = 1.0
      CONSTS(529) = 0
      CONSTS(530) = 1.02127
      CONSTS(531) = 1.01179
      CONSTS(532) = 1.1448
      CONSTS(533) = 1.00007
      CONSTS(534) = -0.060024
      CONSTS(535) = 1.00132
      CONSTS(536) = 0.0367573
      CONSTS(537) = 1
      CONSTS(538) = 1.00163
      CONSTS(539) = 1.00237
      CONSTS(540) = 3622.54
      CONSTS(541) = 2109.91
      CONSTS(542) = 2.368e-07
      CONSTS(543) = 0.694163
      CONSTS(544) = 204.497
      CONSTS(545) = 0.699673
      CONSTS(546) = 103.525
      CONSTS(547) = 16.9144
      CONSTS(548) = 159.549
      CONSTS(549) = 0.419998
      CONSTS(550) = 216.243
      CONSTS(551) = 72.2362
      CONSTS(552) = 48.0839
      CONSTS(553) = 0
      CONSTS(554) = 279.945
      CONSTS(555) = 0.000980838
      CONSTS(556) = 0.862514
      CONSTS(557) = 25.0404
      CONSTS(558) = 0.379883
      CONSTS(559) = 0.38131
      CONSTS(560) = 0.0123238
      CONSTS(561) = 3.00449
      CONSTS(562) = 0.100043
      CONSTS(563) = 2.00439
      CONSTS(564) = 0
      CONSTS(565) = 39.8952
      CONSTS(566) = 0.0101913
      CONSTS(567) = 0.00366525
      CONSTS(568) = 3.28246
      CONSTS(569) =  (CONSTS(2) - 1.00000)*CONSTS(353)+1.00000
      CONSTS(570) = TERNRY(CONSTS(569).LT.0.800000, 0.800000, CONSTS(569))
      CONSTS(571) =  (CONSTS(3) - 1.00000)*CONSTS(331)+1.00000
      CONSTS(572) = TERNRY(CONSTS(571).LT.CONSTS(332), CONSTS(332), CONSTS(571))
      CONSTS(573) =  (CONSTS(570) - 1.00000)*CONSTS(363)+1.00000
      CONSTS(574) =  (CONSTS(3) - 1.00000)*CONSTS(333)+1.00000
      CONSTS(575) =  ( (CONSTS(3) - 1.00000)*CONSTS(337)+1.00000)*0.100000
      CONSTS(576) = TERNRY(CONSTS(575).LT.0.00000, 0.00000, CONSTS(575))
      CONSTS(577) =  (CONSTS(3) - 1.00000)*CONSTS(335)+1.00000
      CONSTS(578) = TERNRY(CONSTS(577).LT.CONSTS(334), CONSTS(334), CONSTS(577))
      STATES(1) = CONSTS(546)
      STATES(2) = CONSTS(553)
      STATES(3) = CONSTS(548)
      STATES(4) = CONSTS(536)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = ((100.000+ (CONSTS(1) - 100.000)*CONSTS(444)) - STATES(1))/CONSTS(443)
      ALGBRC(1) = TERNRY(CONSTS(441).GT.0.00000.AND.CONSTS(449).LE.0.00000, CONSTS(441), TERNRY(CONSTS(449).GT.0.00000, STATES(1), CONSTS(1) - CONSTS(501))
CALL minimize(minfunc_0, CONSTS, VARIABLES, <VAR>)
      RATES(2) =  (ALGBRC(4) - 1.00000)*CONSTS(452)
      ALGBRC(2) = STATES(3)/CONSTS(8)
      ALGBRC(49) =  ALGBRC(15) ** 2.00000*ALGBRC(2)*3.84000
      ALGBRC(50) =  ALGBRC(49)*CONSTS(445)
      RATES(3) = CONSTS(495) - ALGBRC(50)
      ALGBRC(35) =  ALGBRC(3)*CONSTS(9)*0.00616190
      ALGBRC(38) = ( ALGBRC(35)*CONSTS(12))/CONSTS(9)
      ALGBRC(27) =  ALGBRC(16)*CONSTS(454) - CONSTS(455)
      ALGBRC(28) = TERNRY(ALGBRC(27).LT.1.00000, 1.00000, ALGBRC(27))
      ALGBRC(29) =  ( (ALGBRC(14) - 1.20000)*CONSTS(515)+1.20000)*CONSTS(457)
      ALGBRC(30) = ((ALGBRC(16)+CONSTS(456)) - ALGBRC(29)) - ALGBRC(28)
      ALGBRC(31) =  ALGBRC(30)*CONSTS(439)
      ALGBRC(32) = ALGBRC(31)
      ALGBRC(33) =  (ALGBRC(32) - 1.00000)*CONSTS(447)+1.00000
      ALGBRC(34) = TERNRY(ALGBRC(33).LT.0.000100000, 0.000100000, ALGBRC(33))
      ALGBRC(39) =  (ALGBRC(34) - 1.00000)*CONSTS(446)
      ALGBRC(40) =  (ALGBRC(3) - 1.00000)*CONSTS(412)+1.00000
      ALGBRC(41) = TERNRY(ALGBRC(40).LT.0.100000, 0.100000, ALGBRC(40))
      ALGBRC(42) = ( CONSTS(12)/4.40000 ** CONSTS(374)*CONSTS(13)*0.0800000*ALGBRC(41))/CONSTS(578)
      ALGBRC(46) = ((ALGBRC(38)+ALGBRC(42)) - STATES(4)) - ALGBRC(39)
      ALGBRC(47) = TERNRY(ALGBRC(46).LT.0.00000, 0.00000, ALGBRC(46))
      ALGBRC(36) =  (( CONSTS(11)*ALGBRC(34)*CONSTS(386))/CONSTS(384))*( (CONSTS(10) - 1.00000)*CONSTS(310)+1.00000)
      ALGBRC(37) = TERNRY(ALGBRC(36).LT.CONSTS(387), CONSTS(387), ALGBRC(36))
      ALGBRC(43) = (ALGBRC(35) - ALGBRC(37)) - CONSTS(576)
      ALGBRC(44) = TERNRY(ALGBRC(43).LT.1.00000e-08, 1.00000e-08, ALGBRC(43))
      ALGBRC(51) = ALGBRC(49)+ 2.00000*(ALGBRC(44)+ALGBRC(47))
      ALGBRC(52) = TERNRY(ALGBRC(51).GT.0.600000, 0.600000, ALGBRC(51))
      ALGBRC(53) = ALGBRC(51) - 0.600000
      ALGBRC(54) = TERNRY(ALGBRC(53).LT.0.00000, 0.00000, ALGBRC(53))
      ALGBRC(55) = ALGBRC(52)/( 600.000*CONSTS(10))+ALGBRC(54)/360.000
      RATES(4) =  ( (ALGBRC(47)/ALGBRC(55))*0.000451800 - STATES(4))*1.00000
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = TERNRY(CONSTS(441).GT.0.00000.AND.CONSTS(449).LE.0.00000, CONSTS(441), TERNRY(CONSTS(449).GT.0.00000, STATES(1), CONSTS(1) - CONSTS(501))
      ALGBRC(2) = STATES(3)/CONSTS(8)
      ALGBRC(49) =  ALGBRC(15) ** 2.00000*ALGBRC(2)*3.84000
      ALGBRC(50) =  ALGBRC(49)*CONSTS(445)
      ALGBRC(35) =  ALGBRC(3)*CONSTS(9)*0.00616190
      ALGBRC(38) = ( ALGBRC(35)*CONSTS(12))/CONSTS(9)
      ALGBRC(27) =  ALGBRC(16)*CONSTS(454) - CONSTS(455)
      ALGBRC(28) = TERNRY(ALGBRC(27).LT.1.00000, 1.00000, ALGBRC(27))
      ALGBRC(29) =  ( (ALGBRC(14) - 1.20000)*CONSTS(515)+1.20000)*CONSTS(457)
      ALGBRC(30) = ((ALGBRC(16)+CONSTS(456)) - ALGBRC(29)) - ALGBRC(28)
      ALGBRC(31) =  ALGBRC(30)*CONSTS(439)
      ALGBRC(32) = ALGBRC(31)
      ALGBRC(33) =  (ALGBRC(32) - 1.00000)*CONSTS(447)+1.00000
      ALGBRC(34) = TERNRY(ALGBRC(33).LT.0.000100000, 0.000100000, ALGBRC(33))
      ALGBRC(39) =  (ALGBRC(34) - 1.00000)*CONSTS(446)
      ALGBRC(40) =  (ALGBRC(3) - 1.00000)*CONSTS(412)+1.00000
      ALGBRC(41) = TERNRY(ALGBRC(40).LT.0.100000, 0.100000, ALGBRC(40))
      ALGBRC(42) = ( CONSTS(12)/4.40000 ** CONSTS(374)*CONSTS(13)*0.0800000*ALGBRC(41))/CONSTS(578)
      ALGBRC(46) = ((ALGBRC(38)+ALGBRC(42)) - STATES(4)) - ALGBRC(39)
      ALGBRC(47) = TERNRY(ALGBRC(46).LT.0.00000, 0.00000, ALGBRC(46))
      ALGBRC(36) =  (( CONSTS(11)*ALGBRC(34)*CONSTS(386))/CONSTS(384))*( (CONSTS(10) - 1.00000)*CONSTS(310)+1.00000)
      ALGBRC(37) = TERNRY(ALGBRC(36).LT.CONSTS(387), CONSTS(387), ALGBRC(36))
      ALGBRC(43) = (ALGBRC(35) - ALGBRC(37)) - CONSTS(576)
      ALGBRC(44) = TERNRY(ALGBRC(43).LT.1.00000e-08, 1.00000e-08, ALGBRC(43))
      ALGBRC(51) = ALGBRC(49)+ 2.00000*(ALGBRC(44)+ALGBRC(47))
      ALGBRC(52) = TERNRY(ALGBRC(51).GT.0.600000, 0.600000, ALGBRC(51))
      ALGBRC(53) = ALGBRC(51) - 0.600000
      ALGBRC(54) = TERNRY(ALGBRC(53).LT.0.00000, 0.00000, ALGBRC(53))
      ALGBRC(55) = ALGBRC(52)/( 600.000*CONSTS(10))+ALGBRC(54)/360.000
      ALGBRC(25) =  CONSTS(445)*ALGBRC(14)
      ALGBRC(26) =  ALGBRC(15)*CONSTS(445)
      ALGBRC(45) =  ALGBRC(44)*CONSTS(445)
      ALGBRC(48) =  ALGBRC(47)*CONSTS(445)
      ALGBRC(56) =  ALGBRC(55)*CONSTS(445)
      RETURN
      END
REAL FUNCTION minfunc_0(CONSTS, VARIABLES)
            REAL CONSTS(*), VARIABLES(*) 
      minfunc_0[1] = abs(ALGBRC(6) - ( (ALGBRC(3) - 1.00000)*CONSTS(502)+1.00000))
      minfunc_0[2] = abs(ALGBRC(5) - (TERNRY(ALGBRC(6).LT.CONSTS(503), CONSTS(503), TERNRY(ALGBRC(6).GT.CONSTS(504), CONSTS(504), ALGBRC(6))))
      minfunc_0[3] = abs(ALGBRC(4) - (ALGBRC(5) - STATES(2)))
      minfunc_0[4] = abs(ALGBRC(7) -  CONSTS(297)*CONSTS(5)*CONSTS(570)*ALGBRC(4)*CONSTS(572)*40.0000*CONSTS(4))
      minfunc_0[5] = abs(ALGBRC(9) - ((ALGBRC(7) -  CONSTS(14)*CONSTS(342))+CONSTS(342)))
      minfunc_0[6] = abs(ALGBRC(8) - (TERNRY(ALGBRC(9).LT.CONSTS(298), CONSTS(298), ALGBRC(9))))
      minfunc_0[7] = abs(ALGBRC(10) - ( (ALGBRC(4) - 1.00000)*CONSTS(390)+1.00000))
      minfunc_0[8] = abs(ALGBRC(12) -  43.3330*CONSTS(388)*CONSTS(574)*ALGBRC(10)*CONSTS(4)*CONSTS(573))
      minfunc_0[9] = abs(ALGBRC(11) - (TERNRY(ALGBRC(12).LT.CONSTS(389), CONSTS(389), ALGBRC(12))))
      minfunc_0[10] = abs(ALGBRC(13) - (ALGBRC(8)+ALGBRC(11)))
      minfunc_0[11] = abs(ALGBRC(14) - ALGBRC(1)/ALGBRC(13))
      minfunc_0[12] = abs(ALGBRC(18) - ( ALGBRC(14)*(1.00000 - CONSTS(6)))/( ALGBRC(14)*(1.00000 - CONSTS(6)) - ALGBRC(15)))
      minfunc_0[13] = abs(ALGBRC(17) - (TERNRY(ALGBRC(18).LT.1.00000, 1.00000, ALGBRC(18))))
      minfunc_0[14] = abs(ALGBRC(16) - (TERNRY(CONSTS(513).GT.0.00000,  ALGBRC(17) ** 1.35000*CONSTS(7)*0.980000, CONSTS(7)+4.00000)))
      minfunc_0[15] = abs(ALGBRC(20) -  ALGBRC(8)*ALGBRC(14))
      minfunc_0[16] = abs(ALGBRC(19) - (ALGBRC(1) - ALGBRC(20)))
      minfunc_0[17] = abs(ALGBRC(21) - ((ALGBRC(19) - ALGBRC(16)) - CONSTS(436)))
      minfunc_0[18] = abs(ALGBRC(22) -  ALGBRC(21)*CONSTS(397))
      minfunc_0[19] = abs(ALGBRC(15) - (TERNRY(ALGBRC(22).LT.CONSTS(399), CONSTS(399), ALGBRC(22))))
      minfunc_0[20] = abs(ALGBRC(23) -  ALGBRC(15)*8.00000)
      minfunc_0[21] = abs(ALGBRC(24) - ( (ALGBRC(23) - 1.00000)*CONSTS(411)+1.00000))
      minfunc_0[22] = abs(ALGBRC(3) - (TERNRY(ALGBRC(24).LT.0.00000, 0.00000, ALGBRC(24))))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END