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 72 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is ADHMK in component kidney (dimensionless).
C CONSTS(2) is AMK in component kidney (dimensionless).
C CONSTS(3) is AMNA in component kidney (dimensionless).
C CONSTS(4) is ANM in component kidney (dimensionless).
C CONSTS(5) is ANPX in component kidney (dimensionless).
C CONSTS(6) is AUM in component kidney (dimensionless).
C CONSTS(7) is CKE in component kidney (monovalent_mEq_per_litre).
C CONSTS(8) is CNA in component kidney (monovalent_mEq_per_litre).
C CONSTS(9) is HM1 in component kidney (dimensionless).
C CONSTS(10) is MYOGRS in component kidney (dimensionless).
C CONSTS(11) is PA in component kidney (mmHg).
C CONSTS(12) is PAMKRN in component kidney (dimensionless).
C CONSTS(13) is PPC in component kidney (mmHg).
C CONSTS(14) is VTW in component kidney (litre).
C ALGBRC(1) is PAR in component perfusion_pressure (mmHg).
C CONSTS(15) is GBL in component parameter_values (mmHg).
C CONSTS(16) is RAPRSP in component parameter_values (mmHg).
C CONSTS(17) is RFCDFT in component parameter_values (dimensionless).
C CONSTS(18) is RCDFPC in component parameter_values (dimensionless).
C CONSTS(19) is RCDFDP in component parameter_values (minute).
C STATES(1) is PAR1 in component perfusion_pressure (mmHg).
C ALGBRC(3) is MDFLW in component proximal_tubular_and_macula_densa_flow (L_per_minute).
C ALGBRC(4) is RNAUG2 in component renal_autoregulatory_feedback_factor (dimensionless).
C CONSTS(20) is RNAUGN in component parameter_values (minute_per_L).
C CONSTS(21) is RNAULL in component parameter_values (dimensionless).
C CONSTS(22) is RNAUUL in component parameter_values (dimensionless).
C CONSTS(23) is RNAUAD in component parameter_values (per_minute).
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(64) is AUMK in component autonomic_effect_on_AAR (dimensionless).
C CONSTS(24) is ARF in component parameter_values (dimensionless).
C CONSTS(63) is AUMKT in component autonomic_effect_on_AAR (dimensionless).
C CONSTS(66) is ANMAR in component angiotensin_effect_on_AAR (dimensionless).
C CONSTS(25) is ANMAM in component parameter_values (dimensionless).
C CONSTS(26) is ANMARL in component parameter_values (dimensionless).
C CONSTS(65) is ANMAR1 in component angiotensin_effect_on_AAR (dimensionless).
C ALGBRC(7) is AAR1 in component AAR_calculation (mmHg_minute_per_L).
C CONSTS(27) is AARK in component parameter_values (mmHg_minute_per_L).
C ALGBRC(8) is AAR in component atrial_natriuretic_peptide_effect_on_AAR (mmHg_minute_per_L).
C CONSTS(28) is ANPXAF in component parameter_values (mmHg_minute_per_L).
C CONSTS(29) is AARLL in component parameter_values (mmHg_minute_per_L).
C ALGBRC(9) is AART in component atrial_natriuretic_peptide_effect_on_AAR (mmHg_minute_per_L).
C CONSTS(67) is AUMK2 in component autonomic_effect_on_EAR (dimensionless).
C CONSTS(30) is AUMK1 in component parameter_values (dimensionless).
C CONSTS(68) is ANMER in component angiotensin_effect_on_EAR (dimensionless).
C CONSTS(31) is ANMEM in component parameter_values (dimensionless).
C ALGBRC(10) is RNAUG4 in component effect_of_renal_autoregulatory_feedback_on_EAR (dimensionless).
C CONSTS(32) is EFAFR in component parameter_values (dimensionless).
C ALGBRC(11) is EAR in component EAR_calculation (mmHg_minute_per_L).
C CONSTS(33) is EARK in component parameter_values (mmHg_minute_per_L).
C CONSTS(34) is EARLL in component parameter_values (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(25) is RBF in component actual_renal_blood_flow (L_per_minute).
C CONSTS(35) is REK in component parameter_values (dimensionless).
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 CONSTS(36) is GPPD in component parameter_values (dimensionless).
C CONSTS(37) is GLPCA in component parameter_values (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 CONSTS(38) is PXTP in component parameter_values (mmHg).
C CONSTS(39) is GFLC in component parameter_values (L_per_minute_per_mmHg).
C CONSTS(40) is GFNLL in component parameter_values (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 CONSTS(41) is MDFL1 in component parameter_values (dimensionless).
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 CONSTS(42) is RTPPR in component parameter_values (dimensionless).
C CONSTS(43) is RTPPRS in component parameter_values (mmHg).
C ALGBRC(27) is RTSPPC1 in component renal_tissue_osmotic_pressure (mmHg).
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 CONSTS(44) is URFORM in component parameter_values (mOsm_per_minute).
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 CONSTS(45) is RFABX in component parameter_values (dimensionless).
C CONSTS(46) is RVRS in component parameter_values (mmHg_minute_per_L).
C ALGBRC(34) is RFABD in component peritubular_capillary_reabsorption_factor (dimensionless).
C CONSTS(47) is RTSPRS in component parameter_values (mmHg).
C CONSTS(48) is RABSC in component parameter_values (per_mmHg).
C CONSTS(49) is RFABDP in component parameter_values (dimensionless).
C CONSTS(50) is RFABDM in component parameter_values (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 CONSTS(51) is DTNAR in component parameter_values (monovalent_mEq_per_minute).
C CONSTS(52) is DIURET in component parameter_values (dimensionless).
C CONSTS(53) is AHMNAR in component parameter_values (dimensionless).
C CONSTS(54) is DTNARL in component parameter_values (monovalent_mEq_per_minute).
C ALGBRC(36) is DTNARA1 in component Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C CONSTS(70) is DTNANG in component angiotensin_induced_Na_reabsorption_into_distal_tubules (monovalent_mEq_per_minute).
C CONSTS(55) is ANMNAM in component parameter_values (dimensionless).
C CONSTS(69) 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 CONSTS(56) is RFABKM in component parameter_values (monovalent_mEq_per_minute).
C ALGBRC(41) is MDFLK in component effect_of_fluid_flow_on_distal_K_reabsorption (monovalent_mEq_per_minute).
C CONSTS(57) is MDFLKM in component parameter_values (monovalent_mEq_per_litre).
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(58) is ANMKEM in component parameter_values (dimensionless).
C CONSTS(59) is ANMKEL in component parameter_values (dimensionless).
C CONSTS(60) is CKEEX in component parameter_values (dimensionless).
C CONSTS(71) is ANMKE1 in component K_secretion_from_distal_tubules (dimensionless).
C CONSTS(72) 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(45) is KODN1 in component normal_K_excretion (monovalent_mEq_per_minute).
C ALGBRC(48) 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 ALGBRC(46) is NOD in component actual_Na_excretion_rate (monovalent_mEq_per_minute).
C ALGBRC(49) 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(61) is RNAGTC in component parameter_values (minute).
C CONSTS(62) is GFNDMP in component parameter_values (dimensionless).
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) = 1.0
CONSTS(2) = 1.037
CONSTS(3) = 1.0
CONSTS(4) = 0.987545
CONSTS(5) = 1.0
CONSTS(6) = 1.00066
CONSTS(7) = 4.44092
CONSTS(8) = 142.035
CONSTS(9) = 0.39984739
CONSTS(10) = 1.0
CONSTS(11) = 103.525
CONSTS(12) = 1.0
CONSTS(13) = 29.9941
CONSTS(14) = 39.8952
CONSTS(15) = 0
CONSTS(16) = 0
CONSTS(17) = 0
CONSTS(18) = 0
CONSTS(19) = 2000
STATES(1) = 103.525
CONSTS(20) = 0.6
CONSTS(21) = 0.3
CONSTS(22) = 10
CONSTS(23) = 0
STATES(2) = 0.0
CONSTS(24) = 0.5
CONSTS(25) = 1.4
CONSTS(26) = 0.86
CONSTS(27) = 1
CONSTS(28) = 1.5
CONSTS(29) = 4
CONSTS(30) = 0.3
CONSTS(31) = 1.6
CONSTS(32) = 0
CONSTS(33) = 1
CONSTS(34) = 24
CONSTS(35) = 1
CONSTS(36) = 1.0
CONSTS(37) = 1.0
CONSTS(38) = 8
CONSTS(39) = 0.0208333
CONSTS(40) = 0.001
CONSTS(41) = 10
CONSTS(42) = 0.9
CONSTS(43) = 15.2
STATES(3) = 159.549
CONSTS(44) = 0.24
CONSTS(45) = 0.8
CONSTS(46) = 19.167
CONSTS(47) = 6
CONSTS(48) = 0.5
CONSTS(49) = 1
CONSTS(50) = 0.3
CONSTS(51) = 0.675
CONSTS(52) = 1
CONSTS(53) = 0.3
CONSTS(54) = 1e-06
CONSTS(55) = 1
CONSTS(56) = 0.03
CONSTS(57) = 0.667
STATES(4) = 0.0367573
CONSTS(58) = 2
CONSTS(59) = 0.3
CONSTS(60) = 4
CONSTS(61) = 15
CONSTS(62) = 3
CONSTS(63) = (CONSTS(6) - 1.00000)*CONSTS(24)+1.00000
CONSTS(64) = TERNRY(CONSTS(63).LT.0.800000, 0.800000, CONSTS(63))
CONSTS(65) = (CONSTS(4) - 1.00000)*CONSTS(25)+1.00000
CONSTS(66) = TERNRY(CONSTS(65).LT.CONSTS(26), CONSTS(26), CONSTS(65))
CONSTS(67) = (CONSTS(64) - 1.00000)*CONSTS(30)+1.00000
CONSTS(68) = (CONSTS(4) - 1.00000)*CONSTS(31)+1.00000
CONSTS(69) = ( (CONSTS(4) - 1.00000)*CONSTS(55)+1.00000)*0.100000
CONSTS(70) = TERNRY(CONSTS(69).LT.0.00000, 0.00000, CONSTS(69))
CONSTS(71) = (CONSTS(4) - 1.00000)*CONSTS(58)+1.00000
CONSTS(72) = TERNRY(CONSTS(71).LT.CONSTS(59), CONSTS(59), CONSTS(71))
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(1) = ((100.000+ (CONSTS(11) - 100.000)*CONSTS(18)) - STATES(1))/CONSTS(19)
ALGBRC(1) = TERNRY(CONSTS(16).GT.0.00000.AND.CONSTS(17).LE.0.00000, CONSTS(16), TERNRY(CONSTS(17).GT.0.00000, STATES(1), CONSTS(11) - CONSTS(15))
CALL minimize(minfunc_0, CONSTS, VARIABLES, <VAR>)
RATES(2) = (ALGBRC(4) - 1.00000)*CONSTS(23)
ALGBRC(2) = STATES(3)/CONSTS(14)
ALGBRC(48) = ALGBRC(15) ** 2.00000*ALGBRC(2)*3.84000
ALGBRC(50) = ALGBRC(48)*CONSTS(35)
RATES(3) = CONSTS(44) - ALGBRC(50)
ALGBRC(35) = ALGBRC(3)*CONSTS(8)*0.00616190
ALGBRC(38) = ( ALGBRC(35)*CONSTS(7))/CONSTS(8)
ALGBRC(27) = ALGBRC(16)*CONSTS(42) - CONSTS(43)
ALGBRC(28) = TERNRY(ALGBRC(27).LT.1.00000, 1.00000, ALGBRC(27))
ALGBRC(29) = ( (ALGBRC(14) - 1.20000)*CONSTS(45)+1.20000)*CONSTS(46)
ALGBRC(30) = ((ALGBRC(16)+CONSTS(47)) - ALGBRC(29)) - ALGBRC(28)
ALGBRC(31) = ALGBRC(30)*CONSTS(48)
ALGBRC(32) = ALGBRC(31)
ALGBRC(33) = (ALGBRC(32) - 1.00000)*CONSTS(50)+1.00000
ALGBRC(34) = TERNRY(ALGBRC(33).LT.0.000100000, 0.000100000, ALGBRC(33))
ALGBRC(39) = (ALGBRC(34) - 1.00000)*CONSTS(56)
ALGBRC(40) = (ALGBRC(3) - 1.00000)*CONSTS(57)+1.00000
ALGBRC(41) = TERNRY(ALGBRC(40).LT.0.100000, 0.100000, ALGBRC(40))
ALGBRC(42) = ( CONSTS(7)/4.40000 ** CONSTS(60)*CONSTS(2)*0.0800000*ALGBRC(41))/CONSTS(72)
ALGBRC(45) = ((ALGBRC(38)+ALGBRC(42)) - STATES(4)) - ALGBRC(39)
ALGBRC(47) = TERNRY(ALGBRC(45).LT.0.00000, 0.00000, ALGBRC(45))
ALGBRC(36) = (( CONSTS(3)*ALGBRC(34)*CONSTS(51))/CONSTS(52))*( (CONSTS(1) - 1.00000)*CONSTS(53)+1.00000)
ALGBRC(37) = TERNRY(ALGBRC(36).LT.CONSTS(54), CONSTS(54), ALGBRC(36))
ALGBRC(43) = (ALGBRC(35) - ALGBRC(37)) - CONSTS(70)
ALGBRC(44) = TERNRY(ALGBRC(43).LT.1.00000e-08, 1.00000e-08, ALGBRC(43))
ALGBRC(51) = ALGBRC(48)+ 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(1))+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(16).GT.0.00000.AND.CONSTS(17).LE.0.00000, CONSTS(16), TERNRY(CONSTS(17).GT.0.00000, STATES(1), CONSTS(11) - CONSTS(15))
ALGBRC(2) = STATES(3)/CONSTS(14)
ALGBRC(48) = ALGBRC(15) ** 2.00000*ALGBRC(2)*3.84000
ALGBRC(50) = ALGBRC(48)*CONSTS(35)
ALGBRC(35) = ALGBRC(3)*CONSTS(8)*0.00616190
ALGBRC(38) = ( ALGBRC(35)*CONSTS(7))/CONSTS(8)
ALGBRC(27) = ALGBRC(16)*CONSTS(42) - CONSTS(43)
ALGBRC(28) = TERNRY(ALGBRC(27).LT.1.00000, 1.00000, ALGBRC(27))
ALGBRC(29) = ( (ALGBRC(14) - 1.20000)*CONSTS(45)+1.20000)*CONSTS(46)
ALGBRC(30) = ((ALGBRC(16)+CONSTS(47)) - ALGBRC(29)) - ALGBRC(28)
ALGBRC(31) = ALGBRC(30)*CONSTS(48)
ALGBRC(32) = ALGBRC(31)
ALGBRC(33) = (ALGBRC(32) - 1.00000)*CONSTS(50)+1.00000
ALGBRC(34) = TERNRY(ALGBRC(33).LT.0.000100000, 0.000100000, ALGBRC(33))
ALGBRC(39) = (ALGBRC(34) - 1.00000)*CONSTS(56)
ALGBRC(40) = (ALGBRC(3) - 1.00000)*CONSTS(57)+1.00000
ALGBRC(41) = TERNRY(ALGBRC(40).LT.0.100000, 0.100000, ALGBRC(40))
ALGBRC(42) = ( CONSTS(7)/4.40000 ** CONSTS(60)*CONSTS(2)*0.0800000*ALGBRC(41))/CONSTS(72)
ALGBRC(45) = ((ALGBRC(38)+ALGBRC(42)) - STATES(4)) - ALGBRC(39)
ALGBRC(47) = TERNRY(ALGBRC(45).LT.0.00000, 0.00000, ALGBRC(45))
ALGBRC(36) = (( CONSTS(3)*ALGBRC(34)*CONSTS(51))/CONSTS(52))*( (CONSTS(1) - 1.00000)*CONSTS(53)+1.00000)
ALGBRC(37) = TERNRY(ALGBRC(36).LT.CONSTS(54), CONSTS(54), ALGBRC(36))
ALGBRC(43) = (ALGBRC(35) - ALGBRC(37)) - CONSTS(70)
ALGBRC(44) = TERNRY(ALGBRC(43).LT.1.00000e-08, 1.00000e-08, ALGBRC(43))
ALGBRC(51) = ALGBRC(48)+ 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(1))+ALGBRC(54)/360.000
ALGBRC(25) = CONSTS(35)*ALGBRC(14)
ALGBRC(26) = ALGBRC(15)*CONSTS(35)
ALGBRC(46) = ALGBRC(44)*CONSTS(35)
ALGBRC(49) = ALGBRC(47)*CONSTS(35)
ALGBRC(56) = ALGBRC(55)*CONSTS(35)
RETURN
END
REAL FUNCTION minfunc_0(CONSTS, VARIABLES)
REAL CONSTS(*), VARIABLES(*)
minfunc_0[1] = abs(ALGBRC(6) - ( (ALGBRC(3) - 1.00000)*CONSTS(20)+1.00000))
minfunc_0[2] = abs(ALGBRC(5) - (TERNRY(ALGBRC(6).LT.CONSTS(21), CONSTS(21), TERNRY(ALGBRC(6).GT.CONSTS(22), CONSTS(22), ALGBRC(6))))
minfunc_0[3] = abs(ALGBRC(4) - (ALGBRC(5) - STATES(2)))
minfunc_0[4] = abs(ALGBRC(7) - CONSTS(27)*CONSTS(12)*CONSTS(64)*ALGBRC(4)*CONSTS(66)*40.0000*CONSTS(10))
minfunc_0[5] = abs(ALGBRC(9) - ((ALGBRC(7) - CONSTS(5)*CONSTS(28))+CONSTS(28)))
minfunc_0[6] = abs(ALGBRC(8) - (TERNRY(ALGBRC(9).LT.CONSTS(29), CONSTS(29), ALGBRC(9))))
minfunc_0[7] = abs(ALGBRC(10) - ( (ALGBRC(4) - 1.00000)*CONSTS(32)+1.00000))
minfunc_0[8] = abs(ALGBRC(12) - 43.3330*CONSTS(33)*CONSTS(68)*ALGBRC(10)*CONSTS(10)*CONSTS(67))
minfunc_0[9] = abs(ALGBRC(11) - (TERNRY(ALGBRC(12).LT.CONSTS(34), CONSTS(34), 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(9)))/( ALGBRC(14)*(1.00000 - CONSTS(9)) - 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(37).GT.0.00000, ALGBRC(17) ** 1.35000*CONSTS(13)*0.980000, CONSTS(13)+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(38)))
minfunc_0[18] = abs(ALGBRC(22) - ALGBRC(21)*CONSTS(39))
minfunc_0[19] = abs(ALGBRC(15) - (TERNRY(ALGBRC(22).LT.CONSTS(40), CONSTS(40), ALGBRC(22))))
minfunc_0[20] = abs(ALGBRC(23) - ALGBRC(15)*8.00000)
minfunc_0[21] = abs(ALGBRC(24) - ( (ALGBRC(23) - 1.00000)*CONSTS(41)+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