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