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 45 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C CONSTS(1) is ADHMV in component circulatory_dynamics (dimensionless).
C CONSTS(2) is AMM in component circulatory_dynamics (dimensionless).
C CONSTS(3) is ANU in component circulatory_dynamics (dimensionless).
C CONSTS(4) is ANUVN in component circulatory_dynamics (dimensionless).
C CONSTS(5) is ARM in component circulatory_dynamics (dimensionless).
C CONSTS(6) is ATRRFB in component circulatory_dynamics (dimensionless).
C CONSTS(7) is ATRVFB in component circulatory_dynamics (litre).
C CONSTS(8) is AU in component circulatory_dynamics (dimensionless).
C CONSTS(9) is AUH in component circulatory_dynamics (dimensionless).
C CONSTS(10) is AUM in component circulatory_dynamics (dimensionless).
C CONSTS(11) is AVE in component circulatory_dynamics (dimensionless).
C CONSTS(12) is HMD in component circulatory_dynamics (dimensionless).
C CONSTS(13) is HPL in component circulatory_dynamics (dimensionless).
C CONSTS(14) is HPR in component circulatory_dynamics (dimensionless).
C CONSTS(15) is MYOGRS in component circulatory_dynamics (dimensionless).
C CONSTS(16) is OSA in component circulatory_dynamics (dimensionless).
C CONSTS(17) is PAMK in component circulatory_dynamics (dimensionless).
C CONSTS(18) is PC in component circulatory_dynamics (mmHg).
C CONSTS(19) is RBF in component circulatory_dynamics (L_per_minute).
C CONSTS(20) is VIM in component circulatory_dynamics (dimensionless).
C CONSTS(21) is VP in component circulatory_dynamics (litre).
C CONSTS(22) is VRC in component circulatory_dynamics (litre).
C CONSTS(23) is VV6 in component circulatory_dynamics (litre).
C CONSTS(24) is VV7 in component circulatory_dynamics (litre).
C CONSTS(25) is VVR in component circulatory_dynamics (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(46) is QRO in component right_ventricular_output (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(4) is PRA in component right_atrial_pressure (mmHg).
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 CONSTS(26) is HTAUML in component parameter_values (dimensionless).
C ALGBRC(9) is PPA in component pulmonary_vasculature_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(42) is QLO in component left_ventricular_output (L_per_minute).
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 CONSTS(27) is QRF in component parameter_values (L_per_minute).
C CONSTS(28) is HSR in component parameter_values (dimensionless).
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(15) is RPA in component pulmonary_arterial_resistance (mmHg_minute_per_L).
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(18) is PLA in component left_atrial_pressure (mmHg).
C ALGBRC(20) is RPV in component pulmonary_venous_resistance (mmHg_minute_per_L).
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(37) is PA in component arterial_pressure_and_pressure_gradient (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 CONSTS(29) is HSL in component parameter_values (dimensionless).
C ALGBRC(41) is QLO1 in component left_ventricular_output (L_per_minute).
C ALGBRC(59) is QAO in component systemic_blood_flow (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(41) is VVA in component angiotensin_induced_venous_constriction (litre).
C CONSTS(30) is ANY in component parameter_values (litre).
C ALGBRC(28) is VVE in component venous_excess_volume (litre).
C ALGBRC(27) is VVE1 in component venous_excess_volume (litre).
C ALGBRC(30) is PVS in component venous_average_pressure (mmHg).
C CONSTS(31) is CV in component parameter_values (L_per_mmHg).
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 CONSTS(32) is PR1LL in component parameter_values (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(44) is RVS in component venous_resistance (mmHg_minute_per_L).
C CONSTS(33) is CN7 in component parameter_values (dimensionless).
C CONSTS(34) is CN2 in component parameter_values (per_mmHg).
C CONSTS(35) is RVSM in component parameter_values (mmHg_minute_per_L).
C CONSTS(42) is CN3 in component venous_resistance (dimensionless).
C CONSTS(43) is RV1 in component venous_resistance (mmHg_minute_per_L).
C CONSTS(45) 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 CONSTS(36) is PAEX in component parameter_values (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 CONSTS(37) is RAR in component parameter_values (mmHg_minute_per_L).
C CONSTS(38) is RMULT1 in component parameter_values (dimensionless).
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 CONSTS(39) is RAM in component parameter_values (mmHg_minute_per_L).
C ALGBRC(54) is RSN in component total_NM_NR_systemic_resistance (mmHg_minute_per_L).
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(57) is FISFLO in component blood_flow_through_AV_fistulas (L_per_minute).
C CONSTS(40) is FIS in component parameter_values (L_per_minute_per_mmHg).
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) = 1.0
      CONSTS(2) = 1.0
      CONSTS(3) = 0.925271
      CONSTS(4) = 1.0
      CONSTS(5) = 1.16463
      CONSTS(6) = 1.0
      CONSTS(7) = 0.0
      CONSTS(8) = 1.00022
      CONSTS(9) = 1.00012
      CONSTS(10) = 1.00066
      CONSTS(11) = 1.0
      CONSTS(12) = 1.0
      CONSTS(13) = 1.00163
      CONSTS(14) = 1.00237
      CONSTS(15) = 1.0
      CONSTS(16) = 0.97287
      CONSTS(17) = 1.0
      CONSTS(18) = 16.9144
      CONSTS(19) = 1.22057
      CONSTS(20) = 1.00076
      CONSTS(21) = 3.00449
      CONSTS(22) = 2.00439
      CONSTS(23) = 0.0101913
      CONSTS(24) = 0.00366525
      CONSTS(25) = 2.50967
      STATES(1) = 3.28246
      STATES(2) = 0.862514
      STATES(3) = 0.379883
      STATES(4) = 0.38131
      STATES(5) = 0.100043
      CONSTS(26) = 0.4
      CONSTS(27) = 0.15
      CONSTS(28) = 1
      CONSTS(29) = 1
      CONSTS(30) = -0.2
      CONSTS(31) = 0.1
      CONSTS(32) = 0
      CONSTS(33) = 0.2
      CONSTS(34) = 0.0212
      CONSTS(35) = 1
      CONSTS(36) = 2
      CONSTS(37) = 30.52
      CONSTS(38) = 1
      CONSTS(39) = 96.3
      CONSTS(40) = 0
      CONSTS(41) =  (CONSTS(3) - 1.00000)*CONSTS(30)
      CONSTS(42) =  ( (CONSTS(18) - 17.0000)*CONSTS(33)+17.0000)*CONSTS(34)
      CONSTS(43) = CONSTS(35)/CONSTS(42)
      CONSTS(44) =  CONSTS(11)*CONSTS(43)*CONSTS(20)*CONSTS(4)
      CONSTS(45) =  CONSTS(44)*1.79000
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = ((((((CONSTS(21)+CONSTS(22)) - 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(26)*(CONSTS(8) - 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(9)*CONSTS(16))
      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(9)*CONSTS(29)*CONSTS(12)*CONSTS(13)
      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(26) = STATES(1)+ ALGBRC(1)*0.398600
      ALGBRC(27) = ((((ALGBRC(26) - CONSTS(25)) - CONSTS(41)) - CONSTS(24)) - CONSTS(23)) - CONSTS(7)
      ALGBRC(28) = TERNRY(ALGBRC(27).LT.0.000100000, 0.000100000, ALGBRC(27))
      ALGBRC(29) = 3.70000+(ALGBRC(28) - 0.740000)/CONSTS(31)
      ALGBRC(30) = TERNRY(ALGBRC(29).LT.0.000100000, 0.000100000, ALGBRC(29))
      ALGBRC(32) = 0.740000/ALGBRC(30)/( CONSTS(20)*3.70000) ** 0.500000
      ALGBRC(2) = STATES(5)+ ALGBRC(1)*0.0574000
      ALGBRC(3) = ALGBRC(2) - 0.100000
      ALGBRC(4) = ALGBRC(3)/0.00500000
      ALGBRC(31) = TERNRY(ALGBRC(4).LT.CONSTS(32), CONSTS(32), ALGBRC(4))
      ALGBRC(33) = ALGBRC(30) - ALGBRC(31)
      ALGBRC(34) = ALGBRC(33)/ALGBRC(32)
      ALGBRC(10) = (ALGBRC(9)/CONSTS(9))/CONSTS(16)
      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(27))*CONSTS(9)*ALGBRC(11)*CONSTS(28)*CONSTS(12)*CONSTS(14)+( CONSTS(27)*ALGBRC(42))/ALGBRC(25)
      ALGBRC(5) =  (ALGBRC(4)+8.00000)*( CONSTS(26)*(CONSTS(8) - 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(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(40)
      ALGBRC(52) = ALGBRC(37) - ALGBRC(30)
      ALGBRC(47) = ALGBRC(37)/100.000 ** CONSTS(36)
      ALGBRC(50) = (( CONSTS(3)*CONSTS(1)*CONSTS(10)*CONSTS(20)*CONSTS(17))/ALGBRC(47))/CONSTS(6)
      ALGBRC(53) =  CONSTS(39)*CONSTS(2)*ALGBRC(50)*CONSTS(15)*CONSTS(38)
      ALGBRC(55) = ALGBRC(52)/ALGBRC(53)
      ALGBRC(51) =  CONSTS(37)*CONSTS(5)*ALGBRC(50)*CONSTS(15)*CONSTS(38)
      ALGBRC(54) = ALGBRC(51)+CONSTS(45)
      ALGBRC(56) = ALGBRC(52)/ALGBRC(54)
      ALGBRC(58) = ALGBRC(55)+ALGBRC(56)+CONSTS(19)
      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(21)+CONSTS(22)) - 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(26)*(CONSTS(8) - 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(9)*CONSTS(16))
      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(9)*CONSTS(29)*CONSTS(12)*CONSTS(13)
      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(26) = STATES(1)+ ALGBRC(1)*0.398600
      ALGBRC(27) = ((((ALGBRC(26) - CONSTS(25)) - CONSTS(41)) - CONSTS(24)) - CONSTS(23)) - CONSTS(7)
      ALGBRC(28) = TERNRY(ALGBRC(27).LT.0.000100000, 0.000100000, ALGBRC(27))
      ALGBRC(29) = 3.70000+(ALGBRC(28) - 0.740000)/CONSTS(31)
      ALGBRC(30) = TERNRY(ALGBRC(29).LT.0.000100000, 0.000100000, ALGBRC(29))
      ALGBRC(32) = 0.740000/ALGBRC(30)/( CONSTS(20)*3.70000) ** 0.500000
      ALGBRC(2) = STATES(5)+ ALGBRC(1)*0.0574000
      ALGBRC(3) = ALGBRC(2) - 0.100000
      ALGBRC(4) = ALGBRC(3)/0.00500000
      ALGBRC(31) = TERNRY(ALGBRC(4).LT.CONSTS(32), CONSTS(32), ALGBRC(4))
      ALGBRC(33) = ALGBRC(30) - ALGBRC(31)
      ALGBRC(34) = ALGBRC(33)/ALGBRC(32)
      ALGBRC(10) = (ALGBRC(9)/CONSTS(9))/CONSTS(16)
      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(27))*CONSTS(9)*ALGBRC(11)*CONSTS(28)*CONSTS(12)*CONSTS(14)+( CONSTS(27)*ALGBRC(42))/ALGBRC(25)
      ALGBRC(5) =  (ALGBRC(4)+8.00000)*( CONSTS(26)*(CONSTS(8) - 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(48) = ALGBRC(34) - ALGBRC(46)
      ALGBRC(49) = ALGBRC(46) - ALGBRC(23)
      ALGBRC(45) = ALGBRC(37) - ALGBRC(4)
      ALGBRC(57) =  ALGBRC(45)*CONSTS(40)
      ALGBRC(52) = ALGBRC(37) - ALGBRC(30)
      ALGBRC(47) = ALGBRC(37)/100.000 ** CONSTS(36)
      ALGBRC(50) = (( CONSTS(3)*CONSTS(1)*CONSTS(10)*CONSTS(20)*CONSTS(17))/ALGBRC(47))/CONSTS(6)
      ALGBRC(53) =  CONSTS(39)*CONSTS(2)*ALGBRC(50)*CONSTS(15)*CONSTS(38)
      ALGBRC(55) = ALGBRC(52)/ALGBRC(53)
      ALGBRC(51) =  CONSTS(37)*CONSTS(5)*ALGBRC(50)*CONSTS(15)*CONSTS(38)
      ALGBRC(54) = ALGBRC(51)+CONSTS(45)
      ALGBRC(56) = ALGBRC(52)/ALGBRC(54)
      ALGBRC(58) = ALGBRC(55)+ALGBRC(56)+CONSTS(19)
      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