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 43 entries in the algebraic variable array.
C There are a total of 7 entries in each of the rate and state variable arrays.
C There are a total of 41 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C CONSTS(1) is k0_12 in component parameters (per_M2_per_second).
C CONSTS(2) is k0_21 in component parameters (per_second).
C CONSTS(3) is k0_23 in component parameters (per_M_per_second).
C CONSTS(4) is k0_32 in component parameters (per_second).
C CONSTS(5) is k0_34 in component parameters (per_second).
C CONSTS(6) is k0_43 in component parameters (per_second).
C CONSTS(7) is k0_45 in component parameters (per_second).
C CONSTS(8) is k0_54 in component parameters (per_M_per_second).
C CONSTS(9) is k0_25 in component parameters (per_second).
C CONSTS(10) is k0_52 in component parameters (per_second).
C CONSTS(11) is k0_56 in component parameters (per_second).
C CONSTS(12) is k0_65_f in component parameters (per_M2_per_second).
C CONSTS(13) is k0_61_f in component parameters (per_second).
C CONSTS(14) is k0_16_f in component parameters (per_second).
C CONSTS(15) is k0_65_r in component parameters (per_M2_per_second).
C CONSTS(16) is k0_61_r in component parameters (per_second).
C CONSTS(17) is k0_16_r in component parameters (per_second).
C CONSTS(33) is k0_65 in component parameters (per_M2_per_second).
C CONSTS(34) is k0_61 in component parameters (per_second).
C CONSTS(35) is k0_16 in component parameters (per_second).
C CONSTS(18) is delta in component parameters (dimensionless).
C CONSTS(36) is alpha_p in component parameters (dimensionless).
C CONSTS(19) is alpha_pp in component parameters (dimensionless).
C CONSTS(20) is N_C in component parameters (dimensionless).
C CONSTS(21) is N_Avo in component parameters (per_mol).
C CONSTS(22) is area in component parameters (um2).
C CONSTS(37) is C_T in component parameters (umol).
C CONSTS(23) is n in component parameters (dimensionless).
C CONSTS(24) is z_c in component parameters (dimensionless).
C CONSTS(25) is z_Na in component parameters (dimensionless).
C CONSTS(26) is F in component parameters (C_per_mol).
C CONSTS(27) is R in component parameters (J_per_K_per_mol).
C CONSTS(28) is T in component parameters (kelvin).
C STATES(1) is V in component ion_concentrations (volt).
C ALGBRC(1) is mu in component parameters (dimensionless).
C CONSTS(29) is Na_o in component ion_concentrations (M).
C CONSTS(30) is Na_i in component ion_concentrations (M).
C CONSTS(31) is glucose_i in component ion_concentrations (M).
C CONSTS(32) is glucose_o in component ion_concentrations (M).
C ALGBRC(4) is k_12 in component rate_constants (per_second).
C ALGBRC(5) is k_21 in component rate_constants (per_second).
C CONSTS(38) is k_23 in component rate_constants (per_second).
C CONSTS(39) is k_32 in component rate_constants (per_second).
C ALGBRC(6) is k_34 in component rate_constants (per_second).
C ALGBRC(7) is k_43 in component rate_constants (per_second).
C CONSTS(40) is k_45 in component rate_constants (per_second).
C ALGBRC(15) is k_54 in component rate_constants (per_second).
C ALGBRC(8) is k_25 in component rate_constants (per_second).
C ALGBRC(13) is k_52 in component rate_constants (per_second).
C ALGBRC(9) is k_56 in component rate_constants (per_second).
C ALGBRC(10) is k_65 in component rate_constants (per_second).
C ALGBRC(11) is k_61 in component rate_constants (per_second).
C ALGBRC(12) is k_16 in component rate_constants (per_second).
C ALGBRC(2) is ks_12 in component rate_constants (per_M2_per_second).
C ALGBRC(14) is k0_54_temp in component rate_constants (per_M_per_second).
C ALGBRC(3) is k_52_temp in component rate_constants (per_second).
C STATES(2) is C_1 in component kinetic_equations (umol).
C STATES(3) is C_2 in component kinetic_equations (umol).
C STATES(4) is C_3 in component kinetic_equations (umol).
C STATES(5) is C_4 in component kinetic_equations (umol).
C STATES(6) is C_5 in component kinetic_equations (umol).
C ALGBRC(16) is C_6 in component kinetic_equations (umol).
C STATES(7) is C_6_temp in component kinetic_equations (umol).
C ALGBRC(17) is C1_sum in component king_altman_states (per_second5).
C ALGBRC(19) is C2_sum in component king_altman_states (per_second5).
C ALGBRC(23) is C3_sum in component king_altman_states (per_second5).
C ALGBRC(25) is C4_sum in component king_altman_states (per_second5).
C ALGBRC(28) is C5_sum in component king_altman_states (per_second5).
C ALGBRC(32) is C6_sum in component king_altman_states (per_second5).
C ALGBRC(36) is C_sum in component king_altman_states (per_second5).
C ALGBRC(37) is C1 in component king_altman_states (umol).
C ALGBRC(38) is C2 in component king_altman_states (umol).
C ALGBRC(39) is C3 in component king_altman_states (umol).
C ALGBRC(40) is C4 in component king_altman_states (umol).
C ALGBRC(41) is C5 in component king_altman_states (umol).
C ALGBRC(42) is C6 in component king_altman_states (umol).
C ALGBRC(20) is I_NaGl_pSS in component NBC_current (uA).
C ALGBRC(43) is I_NaGl_SS in component NBC_current (uA).
C ALGBRC(29) is epsilon in component phenomonological_constants (per_second).
C ALGBRC(18) is lambda in component phenomonological_constants (per_M3_per_second5).
C ALGBRC(21) is chi in component phenomonological_constants (M).
C ALGBRC(27) is alpha in component phenomonological_constants (M3).
C ALGBRC(24) is beta in component phenomonological_constants (M2).
C ALGBRC(22) is gamma in component phenomonological_constants (M3_per_second).
C ALGBRC(26) is phi in component phenomonological_constants (M_per_second).
C ALGBRC(33) is Imax_Na in component phenomonological_constants (uA).
C ALGBRC(34) is Imax_gluc in component phenomonological_constants (uA).
C ALGBRC(30) is Khalf_Na_sq in component phenomonological_constants (M2).
C ALGBRC(35) is Khalf_Na in component phenomonological_constants (M).
C ALGBRC(31) is Khalf_gluc in component phenomonological_constants (M).
C RATES(1) is d/dt V in component ion_concentrations (volt).
C RATES(2) is d/dt C_1 in component kinetic_equations (umol).
C RATES(3) is d/dt C_2 in component kinetic_equations (umol).
C RATES(4) is d/dt C_3 in component kinetic_equations (umol).
C RATES(5) is d/dt C_4 in component kinetic_equations (umol).
C RATES(6) is d/dt C_5 in component kinetic_equations (umol).
C RATES(7) is d/dt C_6_temp in component kinetic_equations (umol).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 80000
      CONSTS(2) = 500
      CONSTS(3) = 1e5
      CONSTS(4) = 20
      CONSTS(5) = 50
      CONSTS(6) = 50
      CONSTS(7) = 800
      CONSTS(8) = 1.219e4
      CONSTS(9) = 0.3
      CONSTS(10) = 9.1e-4
      CONSTS(11) = 10
      CONSTS(12) = 50
      CONSTS(13) = 5
      CONSTS(14) = 35
      CONSTS(15) = 4500
      CONSTS(16) = 3
      CONSTS(17) = 350
      CONSTS(18) = 0.7
      CONSTS(19) = 0
      CONSTS(20) = 3e6
      CONSTS(21) = 6.022e23
      CONSTS(22) = 1e6
      CONSTS(23) = 2
      CONSTS(24) = -2
      CONSTS(25) = 1
      CONSTS(26) = 96485.34
      CONSTS(27) = 8.314
      CONSTS(28) = 310
      STATES(1) = -150e-3
      CONSTS(29) = 10e-3
      CONSTS(30) = 500e-3
      CONSTS(31) = 100e-3
      CONSTS(32) = 0e-3
      STATES(2) = 0
      STATES(3) = 0
      STATES(4) = 0
      STATES(5) = 0
      STATES(6) = 0
      STATES(7) = 0
      CONSTS(33) = TERNRY(CONSTS(29).LT.0.100000, CONSTS(15), CONSTS(12))
      CONSTS(40) = 0.200000
      CONSTS(34) = TERNRY(CONSTS(29).LT.0.100000, CONSTS(16), CONSTS(13))
      CONSTS(35) = TERNRY(CONSTS(29).LT.0.100000, CONSTS(17), CONSTS(14))
      CONSTS(36) = (1.00000 - CONSTS(18)) - CONSTS(19)
      CONSTS(37) = ( 1.00000e+06*CONSTS(20))/CONSTS(21)
      CONSTS(38) =  CONSTS(3)*CONSTS(32)
      CONSTS(39) = CONSTS(4)
      CONSTS(40) = CONSTS(7)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(1) = CONSTS(40)
      ALGBRC(1) = ( CONSTS(26)*STATES(1))/( CONSTS(27)*CONSTS(28))
      ALGBRC(6) =  CONSTS(5)*EXP(( - (CONSTS(24)+CONSTS(23))*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(7) =  CONSTS(6)*EXP(( (CONSTS(24)+CONSTS(23))*CONSTS(18)*ALGBRC(1))/2.00000)
      RATES(4) = ( CONSTS(38)*STATES(3)+ ALGBRC(7)*STATES(5)) -  (CONSTS(39)+ALGBRC(6))*STATES(4)
      ALGBRC(9) =  CONSTS(11)*EXP(( - CONSTS(23)*CONSTS(25)*CONSTS(19)*ALGBRC(1))/2.00000)
      ALGBRC(10) =  CONSTS(33)*CONSTS(30) ** CONSTS(23)*EXP(( CONSTS(23)*CONSTS(25)*CONSTS(19)*ALGBRC(1))/2.00000)
      ALGBRC(11) =  CONSTS(34)*EXP(( CONSTS(24)*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(12) =  CONSTS(35)*EXP(( - CONSTS(24)*CONSTS(18)*ALGBRC(1))/2.00000)
      RATES(7) = ( ALGBRC(12)*STATES(2)+ ALGBRC(9)*STATES(6)) -  (ALGBRC(11)+ALGBRC(10))*STATES(7)
      ALGBRC(2) =  CONSTS(1)*EXP(( - CONSTS(23)*CONSTS(36)*ALGBRC(1))/2.00000)
      ALGBRC(4) =  ALGBRC(2)*CONSTS(29) ** CONSTS(23)
      ALGBRC(5) =  CONSTS(2)*EXP(( CONSTS(23)*CONSTS(25)*CONSTS(36)*ALGBRC(1))/2.00000)
      ALGBRC(8) =  CONSTS(9)*EXP(( - (CONSTS(24)+CONSTS(23))*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(13) = ( CONSTS(1)*ALGBRC(8)*CONSTS(11)*CONSTS(34))/( CONSTS(2)*CONSTS(35)*CONSTS(33))
      RATES(3) = ( ALGBRC(4)*STATES(2)+ CONSTS(39)*STATES(4)+ ALGBRC(13)*STATES(6)) -  (ALGBRC(5)+CONSTS(38)+ALGBRC(8))*STATES(3)
      ALGBRC(14) = ( CONSTS(3)*ALGBRC(6)*CONSTS(40)*ALGBRC(13))/( ALGBRC(7)*CONSTS(39)*ALGBRC(8))
      ALGBRC(15) =  ALGBRC(14)*CONSTS(31)
      RATES(5) = ( ALGBRC(6)*STATES(4)+ ALGBRC(15)*STATES(6)) -  (CONSTS(40)+ALGBRC(7))*STATES(5)
      ALGBRC(16) = CONSTS(37) - (STATES(2)+STATES(3)+STATES(4)+STATES(5)+STATES(6))
      RATES(2) = ( ALGBRC(5)*STATES(3)+ ALGBRC(11)*ALGBRC(16)) -  (ALGBRC(4)+ALGBRC(12))*STATES(2)
      RATES(6) = ( CONSTS(40)*STATES(5)+ ALGBRC(10)*ALGBRC(16)+ ALGBRC(8)*STATES(3)) -  (ALGBRC(15)+ALGBRC(13)+ALGBRC(9))*STATES(6)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(1) = ( CONSTS(26)*STATES(1))/( CONSTS(27)*CONSTS(28))
      ALGBRC(6) =  CONSTS(5)*EXP(( - (CONSTS(24)+CONSTS(23))*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(7) =  CONSTS(6)*EXP(( (CONSTS(24)+CONSTS(23))*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(9) =  CONSTS(11)*EXP(( - CONSTS(23)*CONSTS(25)*CONSTS(19)*ALGBRC(1))/2.00000)
      ALGBRC(10) =  CONSTS(33)*CONSTS(30) ** CONSTS(23)*EXP(( CONSTS(23)*CONSTS(25)*CONSTS(19)*ALGBRC(1))/2.00000)
      ALGBRC(11) =  CONSTS(34)*EXP(( CONSTS(24)*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(12) =  CONSTS(35)*EXP(( - CONSTS(24)*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(2) =  CONSTS(1)*EXP(( - CONSTS(23)*CONSTS(36)*ALGBRC(1))/2.00000)
      ALGBRC(4) =  ALGBRC(2)*CONSTS(29) ** CONSTS(23)
      ALGBRC(5) =  CONSTS(2)*EXP(( CONSTS(23)*CONSTS(25)*CONSTS(36)*ALGBRC(1))/2.00000)
      ALGBRC(8) =  CONSTS(9)*EXP(( - (CONSTS(24)+CONSTS(23))*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(13) = ( CONSTS(1)*ALGBRC(8)*CONSTS(11)*CONSTS(34))/( CONSTS(2)*CONSTS(35)*CONSTS(33))
      ALGBRC(14) = ( CONSTS(3)*ALGBRC(6)*CONSTS(40)*ALGBRC(13))/( ALGBRC(7)*CONSTS(39)*ALGBRC(8))
      ALGBRC(15) =  ALGBRC(14)*CONSTS(31)
      ALGBRC(16) = CONSTS(37) - (STATES(2)+STATES(3)+STATES(4)+STATES(5)+STATES(6))
      ALGBRC(3) =  CONSTS(10)*EXP(( (CONSTS(24)+CONSTS(23))*CONSTS(18)*ALGBRC(1))/2.00000)
      ALGBRC(17) =  ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ CONSTS(38)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(11)
      ALGBRC(18) =  ALGBRC(2)*CONSTS(3)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(2)*CONSTS(3)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(2)*CONSTS(3)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(2)*CONSTS(3)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(2)*CONSTS(3)*ALGBRC(6)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(2)*CONSTS(3)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ ALGBRC(2)*CONSTS(3)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*ALGBRC(6)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*ALGBRC(6)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*ALGBRC(7)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*ALGBRC(6)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*ALGBRC(6)*CONSTS(40)*ALGBRC(11)+ ALGBRC(2)*CONSTS(3)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)
      ALGBRC(19) =  ALGBRC(12)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(4)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*CONSTS(39)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(4)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(4)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(11)
      ALGBRC(20) =  - CONSTS(26)*( CONSTS(23)*CONSTS(25)*CONSTS(36)*( ALGBRC(4)*STATES(2) -  ALGBRC(5)*STATES(3))+ CONSTS(24)*CONSTS(18)*( ALGBRC(12)*STATES(2) -  ALGBRC(11)*ALGBRC(16))+ CONSTS(23)*CONSTS(25)*CONSTS(19)*( ALGBRC(9)*STATES(6) -  ALGBRC(10)*ALGBRC(16)))
      ALGBRC(21) =  (1.00000/ALGBRC(18))*( ALGBRC(2)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(8)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(8)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(8)*CONSTS(39)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(2)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(2)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(2)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(2)*CONSTS(39)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(2)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(2)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*ALGBRC(6)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*CONSTS(39)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(2)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(2)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(11)+ ALGBRC(2)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)+ ALGBRC(2)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(9)+ ALGBRC(2)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(9))
      ALGBRC(22) =  (1.00000/ALGBRC(18))*( ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10))
      ALGBRC(23) =  ALGBRC(12)*ALGBRC(8)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*CONSTS(38)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(4)*CONSTS(38)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*CONSTS(38)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(38)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*CONSTS(38)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(38)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(38)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*CONSTS(38)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*CONSTS(38)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*CONSTS(38)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(4)*CONSTS(38)*ALGBRC(7)*ALGBRC(13)*ALGBRC(11)
      ALGBRC(24) =  (1.00000/ALGBRC(18))*( CONSTS(3)*ALGBRC(12)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ CONSTS(3)*ALGBRC(12)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ CONSTS(3)*ALGBRC(12)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ CONSTS(3)*ALGBRC(12)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ CONSTS(3)*ALGBRC(12)*ALGBRC(6)*ALGBRC(13)*ALGBRC(10)+ CONSTS(3)*ALGBRC(12)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ CONSTS(3)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ CONSTS(3)*ALGBRC(12)*ALGBRC(6)*CONSTS(40)*ALGBRC(9))
      ALGBRC(25) =  ALGBRC(12)*ALGBRC(8)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*CONSTS(38)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(4)*CONSTS(38)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*CONSTS(39)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*CONSTS(38)*ALGBRC(6)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(38)*ALGBRC(6)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(4)*CONSTS(38)*ALGBRC(6)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(6)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*CONSTS(38)*ALGBRC(6)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(8)*CONSTS(39)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(4)*CONSTS(38)*ALGBRC(6)*ALGBRC(13)*ALGBRC(11)
      ALGBRC(26) =  (1.00000/ALGBRC(18))*(( - ALGBRC(2)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11) -  ALGBRC(2)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)) -  ALGBRC(2)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11))
      ALGBRC(27) =  (1.00000/ALGBRC(18))*( ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(7)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(15)*ALGBRC(10)+ ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)*ALGBRC(11)+ ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)*ALGBRC(11)+ ALGBRC(12)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(13)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(13))
      ALGBRC(28) =  ALGBRC(12)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*CONSTS(38)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ ALGBRC(4)*CONSTS(38)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(10)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(11)+ ALGBRC(4)*CONSTS(38)*ALGBRC(6)*CONSTS(40)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(11)+ ALGBRC(4)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(11)
      ALGBRC(29) =  (1.00000/ALGBRC(18))*- ALGBRC(2)*CONSTS(3)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)*ALGBRC(11)
      ALGBRC(30) = (ALGBRC(27)+ ALGBRC(24)*CONSTS(32))/(ALGBRC(21)+CONSTS(32))
      ALGBRC(31) = (ALGBRC(27)+ ALGBRC(21)*CONSTS(29) ** 2.00000)/(ALGBRC(24)+CONSTS(29) ** 2.00000)
      ALGBRC(32) =  ALGBRC(12)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)+ ALGBRC(4)*ALGBRC(8)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*CONSTS(38)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)+ ALGBRC(4)*CONSTS(38)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(9)+ ALGBRC(4)*ALGBRC(8)*CONSTS(39)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)+ ALGBRC(4)*ALGBRC(8)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(9)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(15)+ ALGBRC(12)*ALGBRC(5)*ALGBRC(6)*CONSTS(40)*ALGBRC(13)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*CONSTS(40)*ALGBRC(13)+ ALGBRC(12)*ALGBRC(5)*CONSTS(39)*ALGBRC(7)*ALGBRC(13)
      ALGBRC(33) = ( 2.00000*CONSTS(26)*CONSTS(37)*(ALGBRC(26)+ ALGBRC(29)*CONSTS(32)))/(ALGBRC(21)+CONSTS(32))
      ALGBRC(34) = ( 2.00000*CONSTS(26)*CONSTS(37)*ALGBRC(29)*CONSTS(29) ** 2.00000)/(ALGBRC(24)+CONSTS(29) ** 2.00000)
      ALGBRC(35) =  ALGBRC(30) ** (1.0 / 2)
      ALGBRC(36) = ALGBRC(17)+ALGBRC(19)+ALGBRC(23)+ALGBRC(25)+ALGBRC(28)+ALGBRC(32)
      ALGBRC(37) = ( CONSTS(37)*ALGBRC(17))/ALGBRC(36)
      ALGBRC(38) = ( CONSTS(37)*ALGBRC(19))/ALGBRC(36)
      ALGBRC(39) = ( CONSTS(37)*ALGBRC(23))/ALGBRC(36)
      ALGBRC(40) = ( CONSTS(37)*ALGBRC(25))/ALGBRC(36)
      ALGBRC(41) = ( CONSTS(37)*ALGBRC(28))/ALGBRC(36)
      ALGBRC(42) = ( CONSTS(37)*ALGBRC(32))/ALGBRC(36)
      ALGBRC(43) =  - CONSTS(26)*( CONSTS(24)*( ALGBRC(12)*ALGBRC(37) -  ALGBRC(11)*ALGBRC(42))+ (CONSTS(24)+ CONSTS(25)*CONSTS(23))*( ALGBRC(8)*ALGBRC(38) -  ALGBRC(13)*ALGBRC(41))+ (CONSTS(24)+ CONSTS(25)*CONSTS(23))*( ALGBRC(6)*ALGBRC(39) -  ALGBRC(7)*ALGBRC(40)))
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END
Source
Derived from workspace Eskandari, Wright, Loo 2005 at changeset 5d9d9d818d66.
Collaboration
To begin collaborating on this work, please use your git client and issue this command:
License
The terms of use/license for this work is unspecified.