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 93 entries in the algebraic variable array.
C There are a total of 33 entries in each of the rate and state variable arrays.
C There are a total of 139 entries in the constant variable array.
C
C
C VOI is t in component main (second).
C CONSTS(1) is rho in component main (Js2_per_m5).
C CONSTS(2) is g in component main (m_per_s2).
C CONSTS(3) is theta_deg in component main (dim).
C CONSTS(139) is theta_rad in component main (radian).
C ALGBRC(1) is z in component main (dim).
C CONSTS(4) is T in component main (second).
C ALGBRC(3) is mt in component main (dim).
C ALGBRC(2) is mta in component main (dim).
C CONSTS(5) is delay in component main (dim).
C CONSTS(6) is To_pk in component main (dim).
C ALGBRC(5) is E_B_LV in component main (kPa_per_L).
C CONSTS(7) is E_B_LVmax in component main (kPa_per_L).
C CONSTS(8) is E_B_LVmin in component main (kPa_per_L).
C STATES(1) is q_B_LV in component main (litre).
C CONSTS(9) is q_B_LV_us in component main (litre).
C ALGBRC(7) is u_B_LV in component main (kPa).
C ALGBRC(9) is E_B_RV in component main (kPa_per_L).
C CONSTS(10) is E_B_RVmax in component main (kPa_per_L).
C CONSTS(11) is E_B_RVmin in component main (kPa_per_L).
C STATES(2) is q_B_RV in component main (litre).
C CONSTS(12) is q_B_RV_us in component main (litre).
C ALGBRC(11) is u_B_RV in component main (kPa).
C CONSTS(13) is E_B_LA in component main (kPa_per_L).
C STATES(3) is q_B_LA in component main (litre).
C CONSTS(14) is q_B_LA_us in component main (litre).
C ALGBRC(13) is u_B_LA in component main (kPa).
C CONSTS(15) is E_B_RA in component main (kPa_per_L).
C STATES(4) is q_B_RA in component main (litre).
C CONSTS(16) is q_B_RA_us in component main (litre).
C ALGBRC(15) is u_B_RA in component main (kPa).
C CONSTS(17) is R_B_AV in component main (kPa_s_per_L).
C ALGBRC(30) is v_B_AV in component main (L_per_s).
C CONSTS(18) is R_B_PV in component main (kPa_s_per_L).
C ALGBRC(21) is v_B_PV in component main (L_per_s).
C CONSTS(19) is R_B_TV in component main (kPa_s_per_L).
C ALGBRC(17) is v_B_TV in component main (L_per_s).
C CONSTS(20) is R_B_MV in component main (kPa_s_per_L).
C ALGBRC(18) is v_B_MV in component main (L_per_s).
C STATES(5) is q_B_PA in component main (litre).
C CONSTS(21) is q_B_PA_us in component main (litre).
C ALGBRC(19) is u_B_PA in component main (kPa).
C CONSTS(22) is E_B_PA in component main (kPa_per_L).
C STATES(6) is q_B_lung in component main (litre).
C CONSTS(23) is q_B_lung_us in component main (litre).
C ALGBRC(23) is u_B_lung in component main (kPa).
C CONSTS(24) is E_B_lung in component main (kPa_per_L).
C CONSTS(25) is z_lung in component main (meter).
C ALGBRC(26) is v_B_lung1 in component main (L_per_s).
C ALGBRC(29) is v_B_lung2 in component main (L_per_s).
C CONSTS(26) is R_B_lung1 in component main (kPa_s_per_L).
C CONSTS(27) is R_B_lung2 in component main (kPa_s_per_L).
C CONSTS(28) is E_B_pulmVein in component main (kPa_per_L).
C STATES(7) is q_B_pulmVein in component main (litre).
C CONSTS(29) is q_B_pulmVein_us in component main (litre).
C ALGBRC(28) is u_B_pulmVein in component main (kPa).
C ALGBRC(31) is v_B_pulmVein in component main (L_per_s).
C CONSTS(30) is R_B_pulmVein in component main (kPa_s_per_L).
C STATES(8) is q_B_brain in component main (litre).
C CONSTS(31) is q_B_brain_us in component main (litre).
C STATES(9) is q_B_brainVein in component main (litre).
C CONSTS(32) is q_B_brainVein_us in component main (litre).
C ALGBRC(20) is u_B_brain in component main (kPa).
C ALGBRC(22) is u_B_brainVein in component main (kPa).
C CONSTS(33) is E_B_brain in component main (kPa_per_L).
C CONSTS(34) is E_B_brainVein in component main (kPa_per_L).
C ALGBRC(32) is v_B_brain1 in component main (L_per_s).
C ALGBRC(24) is v_B_brain2 in component main (L_per_s).
C ALGBRC(25) is v_B_brain3 in component main (L_per_s).
C CONSTS(35) is R_B_brain1 in component main (kPa_s_per_L).
C CONSTS(36) is R_B_brain2 in component main (kPa_s_per_L).
C CONSTS(37) is R_B_brain3 in component main (kPa_s_per_L).
C STATES(10) is q_B_AA in component main (litre).
C CONSTS(38) is q_B_AA_us in component main (litre).
C ALGBRC(27) is u_B_AA in component main (kPa).
C CONSTS(39) is E_B_AA in component main (kPa_per_L).
C STATES(11) is q_B_celiac in component main (litre).
C CONSTS(40) is q_B_celiac_us in component main (litre).
C ALGBRC(33) is u_B_celiac in component main (kPa).
C CONSTS(41) is E_B_celiac in component main (kPa_per_L).
C CONSTS(42) is z_celiac in component main (meter).
C ALGBRC(34) is v_B_celiac in component main (L_per_s).
C CONSTS(43) is R_B_celiac in component main (kPa_s_per_L).
C STATES(12) is q_B_supMes in component main (litre).
C CONSTS(44) is q_B_supMes_us in component main (litre).
C ALGBRC(35) is u_B_supMes in component main (kPa).
C CONSTS(45) is E_B_supMes in component main (kPa_per_L).
C CONSTS(46) is z_supMes in component main (meter).
C ALGBRC(36) is v_B_supMes in component main (L_per_s).
C CONSTS(47) is R_B_supMes in component main (kPa_s_per_L).
C ALGBRC(47) is v_B_infMes in component main (L_per_s).
C CONSTS(48) is R_B_infMes in component main (kPa_s_per_L).
C STATES(13) is q_B_stomach in component main (litre).
C CONSTS(49) is q_B_stomach_us in component main (litre).
C ALGBRC(37) is u_B_stomach in component main (kPa).
C CONSTS(50) is E_B_stomach in component main (kPa_per_L).
C CONSTS(51) is z_stomach in component main (meter).
C ALGBRC(38) is v_B_stomach1 in component main (L_per_s).
C ALGBRC(50) is v_B_stomach2 in component main (L_per_s).
C CONSTS(52) is R_B_stomach1 in component main (kPa_s_per_L).
C CONSTS(53) is R_B_stomach2 in component main (kPa_s_per_L).
C STATES(14) is q_B_spleen in component main (litre).
C CONSTS(54) is q_B_spleen_us in component main (litre).
C ALGBRC(39) is u_B_spleen in component main (kPa).
C CONSTS(55) is E_B_spleen in component main (kPa_per_L).
C CONSTS(56) is z_spleen in component main (meter).
C ALGBRC(40) is v_B_spleen1 in component main (L_per_s).
C ALGBRC(51) is v_B_spleen2 in component main (L_per_s).
C CONSTS(57) is R_B_spleen1 in component main (kPa_s_per_L).
C CONSTS(58) is R_B_spleen2 in component main (kPa_s_per_L).
C STATES(15) is q_B_pancreas in component main (litre).
C CONSTS(59) is q_B_pancreas_us in component main (litre).
C ALGBRC(41) is u_B_pancreas in component main (kPa).
C CONSTS(60) is E_B_pancreas in component main (kPa_per_L).
C CONSTS(61) is z_pancreas in component main (meter).
C ALGBRC(42) is v_B_pancreas1 in component main (L_per_s).
C ALGBRC(43) is v_B_pancreas2 in component main (L_per_s).
C ALGBRC(52) is v_B_pancreas3 in component main (L_per_s).
C CONSTS(62) is R_B_pancreas1 in component main (kPa_s_per_L).
C CONSTS(63) is R_B_pancreas2 in component main (kPa_s_per_L).
C CONSTS(64) is R_B_pancreas3 in component main (kPa_s_per_L).
C STATES(16) is q_B_intestine in component main (litre).
C CONSTS(65) is q_B_intestine_us in component main (litre).
C ALGBRC(44) is u_B_intestine in component main (kPa).
C CONSTS(66) is E_B_intestine in component main (kPa_per_L).
C CONSTS(67) is z_ in component main (meter).
C ALGBRC(45) is v_B_intestine1 in component main (L_per_s).
C ALGBRC(53) is v_B_intestine2 in component main (L_per_s).
C CONSTS(68) is R_B_intestine1 in component main (kPa_s_per_L).
C CONSTS(69) is R_B_intestine2 in component main (kPa_s_per_L).
C STATES(17) is q_B_colon in component main (litre).
C CONSTS(70) is q_B_colon_us in component main (litre).
C ALGBRC(46) is u_B_colon in component main (kPa).
C CONSTS(71) is E_B_colon in component main (kPa_per_L).
C CONSTS(72) is z_colon in component main (meter).
C ALGBRC(48) is v_B_colon1 in component main (L_per_s).
C ALGBRC(54) is v_B_colon2 in component main (L_per_s).
C CONSTS(73) is R_B_colon1 in component main (kPa_s_per_L).
C CONSTS(74) is R_B_colon2 in component main (kPa_s_per_L).
C STATES(18) is q_B_portalVein in component main (litre).
C CONSTS(75) is q_B_portalVein_us in component main (litre).
C ALGBRC(49) is u_B_portalVein in component main (kPa).
C CONSTS(76) is E_B_portalVein in component main (kPa_per_L).
C CONSTS(77) is z_portalVein in component main (meter).
C ALGBRC(56) is v_B_portalVein in component main (L_per_s).
C CONSTS(78) is R_B_portalVein in component main (kPa_s_per_L).
C STATES(19) is q_B_liver in component main (litre).
C CONSTS(79) is q_B_liver_us in component main (litre).
C ALGBRC(55) is u_B_liver in component main (kPa).
C CONSTS(80) is E_B_liver in component main (kPa_per_L).
C CONSTS(81) is z_liver in component main (meter).
C ALGBRC(57) is v_B_liver1 in component main (L_per_s).
C ALGBRC(58) is v_B_liver2 in component main (L_per_s).
C CONSTS(82) is R_B_liver1 in component main (kPa_s_per_L).
C CONSTS(83) is R_B_liver2 in component main (kPa_s_per_L).
C STATES(20) is q_B_lGlomerulus in component main (litre).
C STATES(21) is q_B_rGlomerulus in component main (litre).
C STATES(22) is q_B_lPeritubular in component main (litre).
C STATES(23) is q_B_rPeritubular in component main (litre).
C CONSTS(84) is q_B_lGlomerulus_us in component main (litre).
C CONSTS(85) is q_B_rGlomerulus_us in component main (litre).
C CONSTS(86) is q_B_lPeritubular_us in component main (litre).
C CONSTS(87) is q_B_rPeritubular_us in component main (litre).
C ALGBRC(59) is u_B_lGlomerulus in component main (kPa).
C ALGBRC(60) is u_B_rGlomerulus in component main (kPa).
C ALGBRC(61) is u_B_lPeritubular in component main (kPa).
C ALGBRC(62) is u_B_rPeritubular in component main (kPa).
C CONSTS(88) is E_B_lGlomerulus in component main (kPa_per_L).
C CONSTS(89) is E_B_rGlomerulus in component main (kPa_per_L).
C CONSTS(90) is E_B_lPeritubular in component main (kPa_per_L).
C CONSTS(91) is E_B_rPeritubular in component main (kPa_per_L).
C CONSTS(92) is z_kidney in component main (meter).
C ALGBRC(64) is v_B_lRenal1 in component main (L_per_s).
C ALGBRC(65) is v_B_rRenal1 in component main (L_per_s).
C ALGBRC(66) is v_B_lRenal2 in component main (L_per_s).
C ALGBRC(63) is v_B_rRenal2 in component main (L_per_s).
C ALGBRC(67) is v_B_lRenal3 in component main (L_per_s).
C ALGBRC(68) is v_B_rRenal3 in component main (L_per_s).
C CONSTS(93) is R_B_lRenal1 in component main (kPa_s_per_L).
C CONSTS(94) is R_B_rRenal1 in component main (kPa_s_per_L).
C CONSTS(95) is R_B_lRenal2 in component main (kPa_s_per_L).
C CONSTS(96) is R_B_rRenal2 in component main (kPa_s_per_L).
C CONSTS(97) is R_B_lRenal3 in component main (kPa_s_per_L).
C CONSTS(98) is R_B_rRenal3 in component main (kPa_s_per_L).
C STATES(24) is q_B_musTorso in component main (litre).
C STATES(25) is q_B_musTorsoVein in component main (litre).
C CONSTS(99) is q_B_musTorso_us in component main (litre).
C CONSTS(100) is q_B_musTorsoVein_us in component main (litre).
C ALGBRC(69) is u_B_musTorso in component main (kPa).
C ALGBRC(70) is u_B_musTorsoVein in component main (kPa).
C CONSTS(101) is E_B_musTorso in component main (kPa_per_L).
C CONSTS(102) is E_B_musTorsoVein in component main (kPa_per_L).
C CONSTS(103) is z_musTorso in component main (meter).
C ALGBRC(71) is v_B_musTorso1 in component main (L_per_s).
C ALGBRC(72) is v_B_musTorso2 in component main (L_per_s).
C ALGBRC(73) is v_B_musTorso3 in component main (L_per_s).
C CONSTS(104) is R_B_musTorso1 in component main (kPa_s_per_L).
C CONSTS(105) is R_B_musTorso2 in component main (kPa_s_per_L).
C CONSTS(106) is R_B_musTorso3 in component main (kPa_s_per_L).
C STATES(26) is q_B_lArm in component main (litre).
C STATES(27) is q_B_lArmVein in component main (litre).
C CONSTS(107) is q_B_lArm_us in component main (litre).
C CONSTS(108) is q_B_lArmVein_us in component main (litre).
C ALGBRC(74) is u_B_lArm in component main (kPa).
C ALGBRC(75) is u_B_lArmVein in component main (kPa).
C CONSTS(109) is E_B_lArm in component main (kPa_per_L).
C CONSTS(110) is E_B_lArmVein in component main (kPa_per_L).
C CONSTS(111) is z_lArm in component main (meter).
C ALGBRC(76) is v_B_lArm1 in component main (L_per_s).
C ALGBRC(77) is v_B_lArm2 in component main (L_per_s).
C ALGBRC(78) is v_B_lArm3 in component main (L_per_s).
C CONSTS(112) is R_B_lArm1 in component main (kPa_s_per_L).
C CONSTS(113) is R_B_lArm2 in component main (kPa_s_per_L).
C CONSTS(114) is R_B_lArm3 in component main (kPa_s_per_L).
C STATES(28) is q_B_rArm in component main (litre).
C STATES(29) is q_B_rArmVein in component main (litre).
C CONSTS(115) is q_B_rArm_us in component main (litre).
C CONSTS(116) is q_B_rArmVein_us in component main (litre).
C ALGBRC(79) is u_B_rArm in component main (kPa).
C ALGBRC(80) is u_B_rArmVein in component main (kPa).
C CONSTS(117) is E_B_rArm in component main (kPa_per_L).
C CONSTS(118) is E_B_rArmVein in component main (kPa_per_L).
C CONSTS(119) is z_rArm in component main (meter).
C ALGBRC(81) is v_B_rArm1 in component main (L_per_s).
C ALGBRC(82) is v_B_rArm2 in component main (L_per_s).
C ALGBRC(83) is v_B_rArm3 in component main (L_per_s).
C CONSTS(120) is R_B_rArm1 in component main (kPa_s_per_L).
C CONSTS(121) is R_B_rArm2 in component main (kPa_s_per_L).
C CONSTS(122) is R_B_rArm3 in component main (kPa_s_per_L).
C STATES(30) is q_B_lLeg in component main (litre).
C STATES(31) is q_B_lLegVein in component main (litre).
C CONSTS(123) is q_B_lLeg_us in component main (litre).
C CONSTS(124) is q_B_lLegVein_us in component main (litre).
C ALGBRC(84) is u_B_lLeg in component main (kPa).
C ALGBRC(85) is u_B_lLegVein in component main (kPa).
C CONSTS(125) is E_B_lLeg in component main (kPa_per_L).
C CONSTS(126) is E_B_lLegVein in component main (kPa_per_L).
C CONSTS(127) is z_lLeg in component main (meter).
C ALGBRC(86) is v_B_lLeg1 in component main (L_per_s).
C ALGBRC(87) is v_B_lLeg2 in component main (L_per_s).
C ALGBRC(88) is v_B_lLeg3 in component main (L_per_s).
C CONSTS(128) is R_B_lLeg1 in component main (kPa_s_per_L).
C CONSTS(129) is R_B_lLeg2 in component main (kPa_s_per_L).
C CONSTS(130) is R_B_lLeg3 in component main (kPa_s_per_L).
C STATES(32) is q_B_rLeg in component main (litre).
C STATES(33) is q_B_rLegVein in component main (litre).
C CONSTS(131) is q_B_rLeg_us in component main (litre).
C CONSTS(132) is q_B_rLegVein_us in component main (litre).
C ALGBRC(89) is u_B_rLeg in component main (kPa).
C ALGBRC(90) is u_B_rLegVein in component main (kPa).
C CONSTS(133) is E_B_rLeg in component main (kPa_per_L).
C CONSTS(134) is E_B_rLegVein in component main (kPa_per_L).
C CONSTS(135) is z_rLeg in component main (meter).
C ALGBRC(91) is v_B_rLeg1 in component main (L_per_s).
C ALGBRC(92) is v_B_rLeg2 in component main (L_per_s).
C ALGBRC(93) is v_B_rLeg3 in component main (L_per_s).
C CONSTS(136) is R_B_rLeg1 in component main (kPa_s_per_L).
C CONSTS(137) is R_B_rLeg2 in component main (kPa_s_per_L).
C CONSTS(138) is R_B_rLeg3 in component main (kPa_s_per_L).
C ALGBRC(4) is q_B_heartTot in component main (litre).
C ALGBRC(6) is q_B_lungTot in component main (litre).
C ALGBRC(8) is q_B_renalTot in component main (litre).
C ALGBRC(10) is q_B_gutTot in component main (litre).
C ALGBRC(12) is q_B_brainTot in component main (litre).
C ALGBRC(14) is q_B_limbTot in component main (litre).
C ALGBRC(16) is q_Blood in component main (litre).
C RATES(1) is d/dt q_B_LV in component main (litre).
C RATES(2) is d/dt q_B_RV in component main (litre).
C RATES(3) is d/dt q_B_LA in component main (litre).
C RATES(4) is d/dt q_B_RA in component main (litre).
C RATES(5) is d/dt q_B_PA in component main (litre).
C RATES(6) is d/dt q_B_lung in component main (litre).
C RATES(7) is d/dt q_B_pulmVein in component main (litre).
C RATES(8) is d/dt q_B_brain in component main (litre).
C RATES(9) is d/dt q_B_brainVein in component main (litre).
C RATES(10) is d/dt q_B_AA in component main (litre).
C RATES(11) is d/dt q_B_celiac in component main (litre).
C RATES(12) is d/dt q_B_supMes in component main (litre).
C RATES(13) is d/dt q_B_stomach in component main (litre).
C RATES(14) is d/dt q_B_spleen in component main (litre).
C RATES(15) is d/dt q_B_pancreas in component main (litre).
C RATES(16) is d/dt q_B_intestine in component main (litre).
C RATES(17) is d/dt q_B_colon in component main (litre).
C RATES(18) is d/dt q_B_portalVein in component main (litre).
C RATES(19) is d/dt q_B_liver in component main (litre).
C RATES(20) is d/dt q_B_lGlomerulus in component main (litre).
C RATES(21) is d/dt q_B_rGlomerulus in component main (litre).
C RATES(22) is d/dt q_B_lPeritubular in component main (litre).
C RATES(23) is d/dt q_B_rPeritubular in component main (litre).
C RATES(24) is d/dt q_B_musTorso in component main (litre).
C RATES(25) is d/dt q_B_musTorsoVein in component main (litre).
C RATES(26) is d/dt q_B_lArm in component main (litre).
C RATES(27) is d/dt q_B_lArmVein in component main (litre).
C RATES(28) is d/dt q_B_rArm in component main (litre).
C RATES(29) is d/dt q_B_rArmVein in component main (litre).
C RATES(30) is d/dt q_B_lLeg in component main (litre).
C RATES(31) is d/dt q_B_lLegVein in component main (litre).
C RATES(32) is d/dt q_B_rLeg in component main (litre).
C RATES(33) is d/dt q_B_rLegVein in component main (litre).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 1e3
CONSTS(2) = 9.81
CONSTS(3) = 90
CONSTS(4) = 1
CONSTS(5) = 0.1
CONSTS(6) = 0.1
CONSTS(7) = 320
CONSTS(8) = 6.667
STATES(1) = 1
CONSTS(9) = 0.005
CONSTS(10) = 70
CONSTS(11) = 2
STATES(2) = 0.2
CONSTS(12) = 0.01
CONSTS(13) = 10
STATES(3) = 0.02
CONSTS(14) = 0.004
CONSTS(15) = 10
STATES(4) = 0.02
CONSTS(16) = 0.004
CONSTS(17) = 1
CONSTS(18) = 1
CONSTS(19) = 1
CONSTS(20) = 1
STATES(5) = 0.3
CONSTS(21) = 0.2
CONSTS(22) = 200
STATES(6) = 0.2
CONSTS(23) = 0.1
CONSTS(24) = 100
CONSTS(25) = 0.20
CONSTS(26) = 3
CONSTS(27) = 10
CONSTS(28) = 20
STATES(7) = 0.02
CONSTS(29) = 0.012
CONSTS(30) = 1
STATES(8) = 0.1
CONSTS(31) = 0.06
STATES(9) = 0.1
CONSTS(32) = 0.06
CONSTS(33) = 20
CONSTS(34) = 20
CONSTS(35) = 100
CONSTS(36) = 600
CONSTS(37) = 1
STATES(10) = 0.2
CONSTS(38) = 0.12
CONSTS(39) = 80
STATES(11) = 0.1
CONSTS(40) = 0.06
CONSTS(41) = 400
CONSTS(42) = 0.20
CONSTS(43) = 100
STATES(12) = 0.1
CONSTS(44) = 0.06
CONSTS(45) = 400
CONSTS(46) = 0.20
CONSTS(47) = 100
CONSTS(48) = 100
STATES(13) = 0.1
CONSTS(49) = 0.06
CONSTS(50) = 100
CONSTS(51) = 0.20
CONSTS(52) = 100
CONSTS(53) = 2400
STATES(14) = 0.1
CONSTS(54) = 0.06
CONSTS(55) = 100
CONSTS(56) = 0.20
CONSTS(57) = 100
CONSTS(58) = 3000
STATES(15) = 0.1
CONSTS(59) = 0.06
CONSTS(60) = 100
CONSTS(61) = 0.20
CONSTS(62) = 100
CONSTS(63) = 100
CONSTS(64) = 4000
STATES(16) = 0.1
CONSTS(65) = 0.06
CONSTS(66) = 100
CONSTS(67) = 0.20
CONSTS(68) = 100
CONSTS(69) = 1200
STATES(17) = 0.1
CONSTS(70) = 0.06
CONSTS(71) = 100
CONSTS(72) = 0.20
CONSTS(73) = 100
CONSTS(74) = 2000
STATES(18) = 0.1
CONSTS(75) = 0.06
CONSTS(76) = 50
CONSTS(77) = 0.20
CONSTS(78) = 100
STATES(19) = 0.2
CONSTS(79) = 0.06
CONSTS(80) = 10
CONSTS(81) = -0.50
CONSTS(82) = 800
CONSTS(83) = 1
STATES(20) = 0.1
STATES(21) = 0.1
STATES(22) = 0.1
STATES(23) = 0.1
CONSTS(84) = 0.06
CONSTS(85) = 0.06
CONSTS(86) = 0.06
CONSTS(87) = 0.06
CONSTS(88) = 150
CONSTS(89) = 150
CONSTS(90) = 50
CONSTS(91) = 50
CONSTS(92) = -0.25
CONSTS(93) = 100
CONSTS(94) = 100
CONSTS(95) = 2000
CONSTS(96) = 2000
CONSTS(97) = 1
CONSTS(98) = 1
STATES(24) = 0.2
STATES(25) = 0.2
CONSTS(99) = 0.12
CONSTS(100) = 0.12
CONSTS(101) = 60
CONSTS(102) = 10
CONSTS(103) = -0.50
CONSTS(104) = 100
CONSTS(105) = 1000
CONSTS(106) = 1
STATES(26) = 0.1
STATES(27) = 0.1
CONSTS(107) = 0.06
CONSTS(108) = 0.06
CONSTS(109) = 150
CONSTS(110) = 60
CONSTS(111) = -0.50
CONSTS(112) = 100
CONSTS(113) = 6000
CONSTS(114) = 1
STATES(28) = 0.1
STATES(29) = 0.1
CONSTS(115) = 0.06
CONSTS(116) = 0.06
CONSTS(117) = 150
CONSTS(118) = 60
CONSTS(119) = -0.50
CONSTS(120) = 100
CONSTS(121) = 6000
CONSTS(122) = 1
STATES(30) = 0.1
STATES(31) = 0.1
CONSTS(123) = 0.06
CONSTS(124) = 0.06
CONSTS(125) = 150
CONSTS(126) = 60
CONSTS(127) = -0.50
CONSTS(128) = 100
CONSTS(129) = 4000
CONSTS(130) = 1
STATES(32) = 0.1
STATES(33) = 0.1
CONSTS(131) = 0.06
CONSTS(132) = 0.06
CONSTS(133) = 150
CONSTS(134) = 60
CONSTS(135) = -0.50
CONSTS(136) = 100
CONSTS(137) = 4000
CONSTS(138) = 1
CONSTS(139) = ( CONSTS(3)* 3.14159265358979)/180.000
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(3) = VOI/CONSTS(4) - INT(VOI/CONSTS(4))
ALGBRC(9) = TERNRY(ALGBRC(3).LT.0.400000, (( (CONSTS(10) - CONSTS(11))*ALGBRC(3))/CONSTS(6))*EXP( 0.500000*(1.00000 - ALGBRC(3)/CONSTS(6) ** 2.00000))+CONSTS(11), TERNRY(ALGBRC(3).LT.0.600000, (( (CONSTS(10) - CONSTS(11))*0.400000)/CONSTS(6))*EXP( 0.500000*(1.00000 - 0.400000/CONSTS(6) ** 2.00000))+CONSTS(11), CONSTS(11))
ALGBRC(11) = ALGBRC(9)*(STATES(2) - CONSTS(12))
ALGBRC(19) = CONSTS(22)*(STATES(5) - CONSTS(21))
ALGBRC(21) = TERNRY(ALGBRC(11).GT.ALGBRC(19), (ALGBRC(11) - ALGBRC(19))/CONSTS(18), 0.00000)
ALGBRC(15) = CONSTS(15)*(STATES(4) - CONSTS(16))
ALGBRC(17) = TERNRY(ALGBRC(15).GT.ALGBRC(11), (ALGBRC(15) - ALGBRC(11))/CONSTS(19), 0.00000)
RATES(2) = ALGBRC(17) - ALGBRC(21)
ALGBRC(23) = CONSTS(24)*(STATES(6) - CONSTS(23)) - CONSTS(1)*CONSTS(2)*CONSTS(25)*cos(CONSTS(139))
ALGBRC(26) = (ALGBRC(19) - ALGBRC(23))/CONSTS(26)
RATES(5) = ALGBRC(21) - ALGBRC(26)
ALGBRC(20) = CONSTS(33)*(STATES(8) - CONSTS(31))
ALGBRC(22) = CONSTS(34)*(STATES(9) - CONSTS(32))
ALGBRC(24) = (ALGBRC(20) - ALGBRC(22))/CONSTS(36)
ALGBRC(25) = (ALGBRC(22) - ALGBRC(15))/CONSTS(37)
RATES(9) = ALGBRC(24) - ALGBRC(25)
ALGBRC(5) = TERNRY(ALGBRC(3).LT.0.400000, (( (CONSTS(7) - CONSTS(8))*ALGBRC(3))/CONSTS(6))*EXP( 0.500000*(1.00000 - ALGBRC(3)/CONSTS(6) ** 2.00000))+CONSTS(8), TERNRY(ALGBRC(3).LT.0.600000, (( (CONSTS(7) - CONSTS(8))*0.400000)/CONSTS(6))*EXP( 0.500000*(1.00000 - 0.400000/CONSTS(6) ** 2.00000))+CONSTS(8), CONSTS(8))
ALGBRC(7) = ALGBRC(5)*(STATES(1) - CONSTS(9))
ALGBRC(27) = CONSTS(39)*(STATES(10) - CONSTS(38))
ALGBRC(30) = TERNRY(ALGBRC(7).GT.ALGBRC(27), (ALGBRC(7) - ALGBRC(27))/CONSTS(17), 0.00000)
ALGBRC(13) = CONSTS(13)*(STATES(3) - CONSTS(14))
ALGBRC(18) = TERNRY(ALGBRC(13).GT.ALGBRC(7), (ALGBRC(13) - ALGBRC(7))/CONSTS(20), 0.00000)
RATES(1) = ALGBRC(18) - ALGBRC(30)
ALGBRC(28) = CONSTS(28)*(STATES(7) - CONSTS(29))
ALGBRC(29) = (ALGBRC(23) - ALGBRC(28))/CONSTS(27)
RATES(6) = ALGBRC(26) - ALGBRC(29)
ALGBRC(31) = (ALGBRC(28) - ALGBRC(13))/CONSTS(30)
RATES(3) = ALGBRC(31) - ALGBRC(18)
RATES(7) = ALGBRC(29) - ALGBRC(31)
ALGBRC(32) = (ALGBRC(27) - ALGBRC(20))/CONSTS(35)
RATES(8) = ALGBRC(32) - ALGBRC(24)
ALGBRC(35) = CONSTS(45)*(STATES(12) - CONSTS(44)) - CONSTS(1)*CONSTS(2)*CONSTS(46)*cos(CONSTS(139))
ALGBRC(36) = (ALGBRC(27) - ALGBRC(35))/CONSTS(47)
ALGBRC(41) = CONSTS(60)*(STATES(15) - CONSTS(59)) - CONSTS(1)*CONSTS(2)*CONSTS(61)*cos(CONSTS(139))
ALGBRC(43) = (ALGBRC(35) - ALGBRC(41))/CONSTS(63)
ALGBRC(44) = CONSTS(66)*(STATES(16) - CONSTS(65))
ALGBRC(45) = (ALGBRC(35) - ALGBRC(44))/CONSTS(68)
ALGBRC(46) = CONSTS(71)*(STATES(17) - CONSTS(70)) - CONSTS(1)*CONSTS(2)*CONSTS(72)*cos(CONSTS(139))
ALGBRC(48) = (ALGBRC(35) - ALGBRC(46))/CONSTS(73)
RATES(12) = ((ALGBRC(36) - ALGBRC(43)) - ALGBRC(45)) - ALGBRC(48)
ALGBRC(33) = CONSTS(41)*(STATES(11) - CONSTS(40)) - CONSTS(1)*CONSTS(2)*CONSTS(42)*cos(CONSTS(139))
ALGBRC(37) = CONSTS(50)*(STATES(13) - CONSTS(49)) - CONSTS(1)*CONSTS(2)*CONSTS(51)*cos(CONSTS(139))
ALGBRC(38) = (ALGBRC(33) - ALGBRC(37))/CONSTS(52)
ALGBRC(49) = CONSTS(76)*(STATES(18) - CONSTS(75)) - CONSTS(1)*CONSTS(2)*CONSTS(77)*cos(CONSTS(139))
ALGBRC(50) = (ALGBRC(37) - ALGBRC(49))/CONSTS(53)
RATES(13) = ALGBRC(38) - ALGBRC(50)
ALGBRC(39) = CONSTS(71)*(STATES(14) - CONSTS(54)) - CONSTS(1)*CONSTS(2)*CONSTS(56)*cos(CONSTS(139))
ALGBRC(40) = (ALGBRC(33) - ALGBRC(39))/CONSTS(57)
ALGBRC(51) = (ALGBRC(39) - ALGBRC(49))/CONSTS(58)
RATES(14) = ALGBRC(40) - ALGBRC(51)
ALGBRC(42) = (ALGBRC(33) - ALGBRC(41))/CONSTS(62)
ALGBRC(52) = (ALGBRC(41) - ALGBRC(49))/CONSTS(64)
RATES(15) = (ALGBRC(42)+ALGBRC(43)) - ALGBRC(52)
ALGBRC(53) = (ALGBRC(44) - ALGBRC(49))/CONSTS(69)
RATES(16) = ALGBRC(45) - ALGBRC(53)
ALGBRC(47) = (ALGBRC(27) - ALGBRC(46))/CONSTS(48)
ALGBRC(54) = (ALGBRC(46) - ALGBRC(49))/CONSTS(74)
RATES(17) = (ALGBRC(48)+ALGBRC(47)) - ALGBRC(54)
ALGBRC(55) = CONSTS(80)*(STATES(19) - CONSTS(79)) - CONSTS(1)*CONSTS(2)*CONSTS(81)*cos(CONSTS(139))
ALGBRC(56) = (ALGBRC(49) - ALGBRC(55))/CONSTS(78)
RATES(18) = (ALGBRC(50)+ALGBRC(51)+ALGBRC(52)+ALGBRC(53)+ALGBRC(54)) - ALGBRC(56)
ALGBRC(34) = (ALGBRC(27) - ALGBRC(33))/CONSTS(43)
ALGBRC(57) = (ALGBRC(33) - ALGBRC(55))/CONSTS(82)
RATES(11) = (((ALGBRC(34) - ALGBRC(57)) - ALGBRC(38)) - ALGBRC(40)) - ALGBRC(42)
ALGBRC(58) = (ALGBRC(55) - ALGBRC(15))/CONSTS(83)
RATES(19) = (ALGBRC(57)+ALGBRC(56)) - ALGBRC(58)
ALGBRC(59) = CONSTS(88)*(STATES(20) - CONSTS(84)) - CONSTS(1)*CONSTS(2)*CONSTS(92)*cos(CONSTS(139))
ALGBRC(64) = (ALGBRC(27) - ALGBRC(59))/CONSTS(93)
ALGBRC(61) = CONSTS(90)*(STATES(22) - CONSTS(86)) - CONSTS(1)*CONSTS(2)*CONSTS(92)*cos(CONSTS(139))
ALGBRC(66) = (ALGBRC(59) - ALGBRC(61))/CONSTS(95)
RATES(20) = ALGBRC(64) - ALGBRC(66)
RATES(21) = ALGBRC(64) - ALGBRC(66)
ALGBRC(67) = (ALGBRC(61) - ALGBRC(15))/CONSTS(97)
RATES(22) = ALGBRC(66) - ALGBRC(67)
RATES(23) = ALGBRC(66) - ALGBRC(67)
ALGBRC(69) = CONSTS(101)*(STATES(24) - CONSTS(99)) - CONSTS(1)*CONSTS(2)*CONSTS(103)*cos(CONSTS(139))
ALGBRC(71) = (ALGBRC(27) - ALGBRC(69))/CONSTS(104)
ALGBRC(70) = CONSTS(102)*(STATES(25) - CONSTS(100)) - CONSTS(1)*CONSTS(2)*CONSTS(103)*cos(CONSTS(139))
ALGBRC(72) = (ALGBRC(69) - ALGBRC(70))/CONSTS(105)
RATES(24) = ALGBRC(71) - ALGBRC(72)
ALGBRC(73) = (ALGBRC(70) - ALGBRC(15))/CONSTS(106)
RATES(25) = ALGBRC(72) - ALGBRC(73)
ALGBRC(74) = CONSTS(109)*(STATES(26) - CONSTS(107)) - CONSTS(1)*CONSTS(2)*CONSTS(111)*cos(CONSTS(139))
ALGBRC(76) = (ALGBRC(27) - ALGBRC(74))/CONSTS(112)
ALGBRC(75) = CONSTS(110)*(STATES(27) - CONSTS(108)) - CONSTS(1)*CONSTS(2)*CONSTS(111)*cos(CONSTS(139))
ALGBRC(77) = (ALGBRC(74) - ALGBRC(75))/CONSTS(113)
RATES(26) = ALGBRC(76) - ALGBRC(77)
ALGBRC(78) = (ALGBRC(75) - ALGBRC(15))/CONSTS(114)
RATES(27) = ALGBRC(77) - ALGBRC(78)
ALGBRC(79) = CONSTS(117)*(STATES(28) - CONSTS(115)) - CONSTS(1)*CONSTS(2)*CONSTS(119)*cos(CONSTS(139))
ALGBRC(81) = (ALGBRC(27) - ALGBRC(79))/CONSTS(120)
ALGBRC(80) = CONSTS(118)*(STATES(29) - CONSTS(115)) - CONSTS(1)*CONSTS(2)*CONSTS(119)*cos(CONSTS(139))
ALGBRC(82) = (ALGBRC(79) - ALGBRC(80))/CONSTS(121)
RATES(28) = ALGBRC(81) - ALGBRC(82)
ALGBRC(83) = (ALGBRC(80) - ALGBRC(15))/CONSTS(122)
RATES(29) = ALGBRC(82) - ALGBRC(83)
ALGBRC(84) = CONSTS(125)*(STATES(30) - CONSTS(123)) - CONSTS(1)*CONSTS(2)*CONSTS(127)*cos(CONSTS(139))
ALGBRC(86) = (ALGBRC(27) - ALGBRC(84))/CONSTS(128)
ALGBRC(85) = CONSTS(126)*(STATES(31) - CONSTS(123)) - CONSTS(1)*CONSTS(2)*CONSTS(127)*cos(CONSTS(139))
ALGBRC(87) = (ALGBRC(84) - ALGBRC(85))/CONSTS(129)
RATES(30) = ALGBRC(86) - ALGBRC(87)
ALGBRC(88) = (ALGBRC(85) - ALGBRC(15))/CONSTS(130)
RATES(31) = ALGBRC(87) - ALGBRC(88)
ALGBRC(60) = CONSTS(89)*(STATES(21) - CONSTS(85)) - CONSTS(1)*CONSTS(2)*CONSTS(92)*cos(CONSTS(139))
ALGBRC(65) = (ALGBRC(27) - ALGBRC(60))/CONSTS(94)
ALGBRC(89) = CONSTS(133)*(STATES(32) - CONSTS(131)) - CONSTS(1)*CONSTS(2)*CONSTS(135)*cos(CONSTS(139))
ALGBRC(91) = (ALGBRC(27) - ALGBRC(89))/CONSTS(136)
RATES(10) = ((((((((((ALGBRC(30) - ALGBRC(32)) - ALGBRC(34)) - ALGBRC(36)) - ALGBRC(47)) - ALGBRC(64)) - ALGBRC(65)) - ALGBRC(71)) - ALGBRC(76)) - ALGBRC(81)) - ALGBRC(86)) - ALGBRC(91)
ALGBRC(90) = CONSTS(134)*(STATES(33) - CONSTS(131)) - CONSTS(1)*CONSTS(2)*CONSTS(135)*cos(CONSTS(139))
ALGBRC(92) = (ALGBRC(89) - ALGBRC(90))/CONSTS(137)
RATES(32) = ALGBRC(91) - ALGBRC(92)
ALGBRC(62) = CONSTS(91)*(STATES(23) - CONSTS(87)) - CONSTS(1)*CONSTS(2)*CONSTS(92)*cos(CONSTS(139))
ALGBRC(68) = (ALGBRC(62) - ALGBRC(15))/CONSTS(98)
ALGBRC(93) = (ALGBRC(90) - ALGBRC(15))/CONSTS(138)
RATES(4) = (ALGBRC(25)+ALGBRC(58)+ALGBRC(67)+ALGBRC(68)+ALGBRC(73)+ALGBRC(78)+ALGBRC(83)+ALGBRC(88)+ALGBRC(93)) - ALGBRC(17)
RATES(33) = ALGBRC(92) - ALGBRC(93)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(3) = VOI/CONSTS(4) - INT(VOI/CONSTS(4))
ALGBRC(9) = TERNRY(ALGBRC(3).LT.0.400000, (( (CONSTS(10) - CONSTS(11))*ALGBRC(3))/CONSTS(6))*EXP( 0.500000*(1.00000 - ALGBRC(3)/CONSTS(6) ** 2.00000))+CONSTS(11), TERNRY(ALGBRC(3).LT.0.600000, (( (CONSTS(10) - CONSTS(11))*0.400000)/CONSTS(6))*EXP( 0.500000*(1.00000 - 0.400000/CONSTS(6) ** 2.00000))+CONSTS(11), CONSTS(11))
ALGBRC(11) = ALGBRC(9)*(STATES(2) - CONSTS(12))
ALGBRC(19) = CONSTS(22)*(STATES(5) - CONSTS(21))
ALGBRC(21) = TERNRY(ALGBRC(11).GT.ALGBRC(19), (ALGBRC(11) - ALGBRC(19))/CONSTS(18), 0.00000)
ALGBRC(15) = CONSTS(15)*(STATES(4) - CONSTS(16))
ALGBRC(17) = TERNRY(ALGBRC(15).GT.ALGBRC(11), (ALGBRC(15) - ALGBRC(11))/CONSTS(19), 0.00000)
ALGBRC(23) = CONSTS(24)*(STATES(6) - CONSTS(23)) - CONSTS(1)*CONSTS(2)*CONSTS(25)*cos(CONSTS(139))
ALGBRC(26) = (ALGBRC(19) - ALGBRC(23))/CONSTS(26)
ALGBRC(20) = CONSTS(33)*(STATES(8) - CONSTS(31))
ALGBRC(22) = CONSTS(34)*(STATES(9) - CONSTS(32))
ALGBRC(24) = (ALGBRC(20) - ALGBRC(22))/CONSTS(36)
ALGBRC(25) = (ALGBRC(22) - ALGBRC(15))/CONSTS(37)
ALGBRC(5) = TERNRY(ALGBRC(3).LT.0.400000, (( (CONSTS(7) - CONSTS(8))*ALGBRC(3))/CONSTS(6))*EXP( 0.500000*(1.00000 - ALGBRC(3)/CONSTS(6) ** 2.00000))+CONSTS(8), TERNRY(ALGBRC(3).LT.0.600000, (( (CONSTS(7) - CONSTS(8))*0.400000)/CONSTS(6))*EXP( 0.500000*(1.00000 - 0.400000/CONSTS(6) ** 2.00000))+CONSTS(8), CONSTS(8))
ALGBRC(7) = ALGBRC(5)*(STATES(1) - CONSTS(9))
ALGBRC(27) = CONSTS(39)*(STATES(10) - CONSTS(38))
ALGBRC(30) = TERNRY(ALGBRC(7).GT.ALGBRC(27), (ALGBRC(7) - ALGBRC(27))/CONSTS(17), 0.00000)
ALGBRC(13) = CONSTS(13)*(STATES(3) - CONSTS(14))
ALGBRC(18) = TERNRY(ALGBRC(13).GT.ALGBRC(7), (ALGBRC(13) - ALGBRC(7))/CONSTS(20), 0.00000)
ALGBRC(28) = CONSTS(28)*(STATES(7) - CONSTS(29))
ALGBRC(29) = (ALGBRC(23) - ALGBRC(28))/CONSTS(27)
ALGBRC(31) = (ALGBRC(28) - ALGBRC(13))/CONSTS(30)
ALGBRC(32) = (ALGBRC(27) - ALGBRC(20))/CONSTS(35)
ALGBRC(35) = CONSTS(45)*(STATES(12) - CONSTS(44)) - CONSTS(1)*CONSTS(2)*CONSTS(46)*cos(CONSTS(139))
ALGBRC(36) = (ALGBRC(27) - ALGBRC(35))/CONSTS(47)
ALGBRC(41) = CONSTS(60)*(STATES(15) - CONSTS(59)) - CONSTS(1)*CONSTS(2)*CONSTS(61)*cos(CONSTS(139))
ALGBRC(43) = (ALGBRC(35) - ALGBRC(41))/CONSTS(63)
ALGBRC(44) = CONSTS(66)*(STATES(16) - CONSTS(65))
ALGBRC(45) = (ALGBRC(35) - ALGBRC(44))/CONSTS(68)
ALGBRC(46) = CONSTS(71)*(STATES(17) - CONSTS(70)) - CONSTS(1)*CONSTS(2)*CONSTS(72)*cos(CONSTS(139))
ALGBRC(48) = (ALGBRC(35) - ALGBRC(46))/CONSTS(73)
ALGBRC(33) = CONSTS(41)*(STATES(11) - CONSTS(40)) - CONSTS(1)*CONSTS(2)*CONSTS(42)*cos(CONSTS(139))
ALGBRC(37) = CONSTS(50)*(STATES(13) - CONSTS(49)) - CONSTS(1)*CONSTS(2)*CONSTS(51)*cos(CONSTS(139))
ALGBRC(38) = (ALGBRC(33) - ALGBRC(37))/CONSTS(52)
ALGBRC(49) = CONSTS(76)*(STATES(18) - CONSTS(75)) - CONSTS(1)*CONSTS(2)*CONSTS(77)*cos(CONSTS(139))
ALGBRC(50) = (ALGBRC(37) - ALGBRC(49))/CONSTS(53)
ALGBRC(39) = CONSTS(71)*(STATES(14) - CONSTS(54)) - CONSTS(1)*CONSTS(2)*CONSTS(56)*cos(CONSTS(139))
ALGBRC(40) = (ALGBRC(33) - ALGBRC(39))/CONSTS(57)
ALGBRC(51) = (ALGBRC(39) - ALGBRC(49))/CONSTS(58)
ALGBRC(42) = (ALGBRC(33) - ALGBRC(41))/CONSTS(62)
ALGBRC(52) = (ALGBRC(41) - ALGBRC(49))/CONSTS(64)
ALGBRC(53) = (ALGBRC(44) - ALGBRC(49))/CONSTS(69)
ALGBRC(47) = (ALGBRC(27) - ALGBRC(46))/CONSTS(48)
ALGBRC(54) = (ALGBRC(46) - ALGBRC(49))/CONSTS(74)
ALGBRC(55) = CONSTS(80)*(STATES(19) - CONSTS(79)) - CONSTS(1)*CONSTS(2)*CONSTS(81)*cos(CONSTS(139))
ALGBRC(56) = (ALGBRC(49) - ALGBRC(55))/CONSTS(78)
ALGBRC(34) = (ALGBRC(27) - ALGBRC(33))/CONSTS(43)
ALGBRC(57) = (ALGBRC(33) - ALGBRC(55))/CONSTS(82)
ALGBRC(58) = (ALGBRC(55) - ALGBRC(15))/CONSTS(83)
ALGBRC(59) = CONSTS(88)*(STATES(20) - CONSTS(84)) - CONSTS(1)*CONSTS(2)*CONSTS(92)*cos(CONSTS(139))
ALGBRC(64) = (ALGBRC(27) - ALGBRC(59))/CONSTS(93)
ALGBRC(61) = CONSTS(90)*(STATES(22) - CONSTS(86)) - CONSTS(1)*CONSTS(2)*CONSTS(92)*cos(CONSTS(139))
ALGBRC(66) = (ALGBRC(59) - ALGBRC(61))/CONSTS(95)
ALGBRC(67) = (ALGBRC(61) - ALGBRC(15))/CONSTS(97)
ALGBRC(69) = CONSTS(101)*(STATES(24) - CONSTS(99)) - CONSTS(1)*CONSTS(2)*CONSTS(103)*cos(CONSTS(139))
ALGBRC(71) = (ALGBRC(27) - ALGBRC(69))/CONSTS(104)
ALGBRC(70) = CONSTS(102)*(STATES(25) - CONSTS(100)) - CONSTS(1)*CONSTS(2)*CONSTS(103)*cos(CONSTS(139))
ALGBRC(72) = (ALGBRC(69) - ALGBRC(70))/CONSTS(105)
ALGBRC(73) = (ALGBRC(70) - ALGBRC(15))/CONSTS(106)
ALGBRC(74) = CONSTS(109)*(STATES(26) - CONSTS(107)) - CONSTS(1)*CONSTS(2)*CONSTS(111)*cos(CONSTS(139))
ALGBRC(76) = (ALGBRC(27) - ALGBRC(74))/CONSTS(112)
ALGBRC(75) = CONSTS(110)*(STATES(27) - CONSTS(108)) - CONSTS(1)*CONSTS(2)*CONSTS(111)*cos(CONSTS(139))
ALGBRC(77) = (ALGBRC(74) - ALGBRC(75))/CONSTS(113)
ALGBRC(78) = (ALGBRC(75) - ALGBRC(15))/CONSTS(114)
ALGBRC(79) = CONSTS(117)*(STATES(28) - CONSTS(115)) - CONSTS(1)*CONSTS(2)*CONSTS(119)*cos(CONSTS(139))
ALGBRC(81) = (ALGBRC(27) - ALGBRC(79))/CONSTS(120)
ALGBRC(80) = CONSTS(118)*(STATES(29) - CONSTS(115)) - CONSTS(1)*CONSTS(2)*CONSTS(119)*cos(CONSTS(139))
ALGBRC(82) = (ALGBRC(79) - ALGBRC(80))/CONSTS(121)
ALGBRC(83) = (ALGBRC(80) - ALGBRC(15))/CONSTS(122)
ALGBRC(84) = CONSTS(125)*(STATES(30) - CONSTS(123)) - CONSTS(1)*CONSTS(2)*CONSTS(127)*cos(CONSTS(139))
ALGBRC(86) = (ALGBRC(27) - ALGBRC(84))/CONSTS(128)
ALGBRC(85) = CONSTS(126)*(STATES(31) - CONSTS(123)) - CONSTS(1)*CONSTS(2)*CONSTS(127)*cos(CONSTS(139))
ALGBRC(87) = (ALGBRC(84) - ALGBRC(85))/CONSTS(129)
ALGBRC(88) = (ALGBRC(85) - ALGBRC(15))/CONSTS(130)
ALGBRC(60) = CONSTS(89)*(STATES(21) - CONSTS(85)) - CONSTS(1)*CONSTS(2)*CONSTS(92)*cos(CONSTS(139))
ALGBRC(65) = (ALGBRC(27) - ALGBRC(60))/CONSTS(94)
ALGBRC(89) = CONSTS(133)*(STATES(32) - CONSTS(131)) - CONSTS(1)*CONSTS(2)*CONSTS(135)*cos(CONSTS(139))
ALGBRC(91) = (ALGBRC(27) - ALGBRC(89))/CONSTS(136)
ALGBRC(90) = CONSTS(134)*(STATES(33) - CONSTS(131)) - CONSTS(1)*CONSTS(2)*CONSTS(135)*cos(CONSTS(139))
ALGBRC(92) = (ALGBRC(89) - ALGBRC(90))/CONSTS(137)
ALGBRC(62) = CONSTS(91)*(STATES(23) - CONSTS(87)) - CONSTS(1)*CONSTS(2)*CONSTS(92)*cos(CONSTS(139))
ALGBRC(68) = (ALGBRC(62) - ALGBRC(15))/CONSTS(98)
ALGBRC(93) = (ALGBRC(90) - ALGBRC(15))/CONSTS(138)
ALGBRC(1) = TERNRY(VOI.GT.100.000.AND.VOI.LT.101.000, 1.00000, 0.00000)
ALGBRC(2) = (VOI/CONSTS(4) - CONSTS(5)) - INT((VOI - CONSTS(5)*CONSTS(4))/CONSTS(4))
ALGBRC(4) = STATES(1)+STATES(2)+STATES(3)+STATES(4)+STATES(10)+STATES(5)
ALGBRC(6) = STATES(6)+STATES(7)
ALGBRC(8) = STATES(20)+STATES(21)+STATES(22)+STATES(23)
ALGBRC(10) = STATES(11)+STATES(12)+STATES(13)+STATES(14)+STATES(15)+STATES(16)+STATES(17)+STATES(18)+STATES(19)+STATES(24)+STATES(25)
ALGBRC(12) = STATES(8)+STATES(9)
ALGBRC(14) = STATES(26)+STATES(27)+STATES(28)+STATES(29)+STATES(30)+STATES(31)+STATES(32)+STATES(33)
ALGBRC(16) = ALGBRC(4)+ALGBRC(6)+ALGBRC(8)+ALGBRC(10)+ALGBRC(12)+ALGBRC(14)
ALGBRC(63) = (ALGBRC(60) - ALGBRC(62))/CONSTS(96)
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END