From 01cc637f10e8b2079059514f609d244152b18eca Mon Sep 17 00:00:00 2001 From: Piyush Sharda <34922596+psharda@users.noreply.github.com> Date: Wed, 18 Sep 2024 00:14:38 +0200 Subject: [PATCH] Correctly handle ices in metal chemistry (#1650) Minor modifications to ensure we handle ices correctly in the metal chemistry network (following the methodology used in KROME - Grassi et al. 2014). Also updates the reference solution outputs (minor differences from what we already have). --- .github/workflows/burn_cell_metal_chem.yml | 8 +- EOS/metal_chem/_parameters | 4 +- networks/metal_chem/actual_rhs.H | 978 +++++++++--------- .../reference_solution_1.out | 74 +- .../reference_solution_1_z10.out | 78 +- .../reference_solution_1e-1.out | 74 +- .../reference_solution_1e-2.out | 74 +- .../reference_solution_1e-3.out | 74 +- .../reference_solution_1e-4.out | 74 +- .../reference_solution_1e-5.out | 76 +- .../reference_solution_1e-6.out | 84 +- 11 files changed, 800 insertions(+), 798 deletions(-) diff --git a/.github/workflows/burn_cell_metal_chem.yml b/.github/workflows/burn_cell_metal_chem.yml index 03386c464..0631ca629 100644 --- a/.github/workflows/burn_cell_metal_chem.yml +++ b/.github/workflows/burn_cell_metal_chem.yml @@ -39,14 +39,14 @@ jobs: cd unit_test/burn_cell_metal_chem declare -A line_numbers_map=( - ["Z=1"]="4461 4463 4465 4467 4468 4471 4472 4475 4481 4486" + ["Z=1"]="4450 4452 4454 4456 4457 4460 4461 4464 4470 4475" ["Z=1_z10"]="1613 1615 4456 4458 4460 4463 4470 4476 4481" ["Z=1e-1"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463" ["Z=1e-2"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463" ["Z=1e-3"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463" ["Z=1e-4"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463" ["Z=1e-5"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463" - ["Z=1e-6"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463" + ["Z=1e-6"]="4451 4453 4455 4457 4458 4461 4462 4465 4471 4476" ) declare -A ref_map=( @@ -92,9 +92,11 @@ jobs: # Adjust the line number for the reference file if [[ "$Z" == "Z=1" ]]; then - reference_line_number=$((line_number - 4460)) + reference_line_number=$((line_number - 4449)) elif [[ "$Z" == "Z=1_z10" ]]; then reference_line_number=${ref_line_number_z10[$index]} + elif [[ "$Z" == "Z=1e-6" ]]; then + reference_line_number=$((line_number - 4450)) else reference_line_number=$((line_number - 4437)) fi diff --git a/EOS/metal_chem/_parameters b/EOS/metal_chem/_parameters index 8fa847fd8..238fad957 100644 --- a/EOS/metal_chem/_parameters +++ b/EOS/metal_chem/_parameters @@ -3,11 +3,11 @@ eos_gamma_default real 1.4 # define the specie names, and their masses and gammas -species_1_name string "co_ice" +species_1_name string "co_total" species_1_gamma real 5./3. species_1_mass real 0.0 -species_2_name string "h2o_ice" +species_2_name string "h2o_total" species_2_gamma real 5./3. species_2_mass real 0.0 diff --git a/networks/metal_chem/actual_rhs.H b/networks/metal_chem/actual_rhs.H index 3e96e81cd..bf64a1e17 100644 --- a/networks/metal_chem/actual_rhs.H +++ b/networks/metal_chem/actual_rhs.H @@ -32,968 +32,968 @@ void rhs_specie(const burn_t& state, Real T = state.T; - Real x0 = 2.0e-99*X(0); + Real x0 = std::max(T, 2.73*z + 2.73); - Real x1 = std::max(T, 2.73*z + 2.73); + Real x1 = ((x0)*(x0)); - Real x2 = ((x1)*(x1)); + Real x2 = 0.002*x0 + 1.0; - Real x3 = 0.002*x1 + 1.0; + Real x3 = std::sqrt(x0); - Real x4 = std::sqrt(x1); + Real x4 = user_dust2gas_ratio*x3; - Real x5 = user_dust2gas_ratio*x4; + Real x5 = x4*(3.3470650363800003e-24*X(10) + 5.0186540981899997e-24*X(11) + 5.01956503638e-24*X(12) + 6.6902431600000005e-24*X(13) + 6.6911540981899994e-24*X(14) + 6.6920650363799998e-24*X(15) + 2.007528417094e-23*X(16) + 2.0076195109127999e-23*X(17) + 2.1749727627315999e-23*X(18) + 2.3423260145503998e-23*X(19) + 9.1093818800000008e-28*X(2) + 2.5096792663692001e-23*X(20) + 2.6767349207315999e-23*X(21) + 2.6768260145504001e-23*X(22) + 2.6770325181879997e-23*X(23) + 2.8440881725504001e-23*X(24) + 2.8441792663692003e-23*X(25) + 3.0114414243691998e-23*X(26) + 3.011532518188e-23*X(27) + 3.178794676188e-23*X(28) + 4.6843544316443998e-23*X(29) + 1.6726215800000001e-24*X(3) + 4.6844455254632e-23*X(30) + 5.3535609352819999e-23*X(31) + 5.3536520291008001e-23*X(32) + 7.3612715400136001e-23*X(33) + 1.6735325181900001e-24*X(4) + 1.6744434563800001e-24*X(5) + 3.3451215800000003e-24*X(6) + 3.3460325181899999e-24*X(7) + 3.3461540981899999e-24*X(8) + 3.3469434563800003e-24*X(9))/(8.0000000000000013e-6*x1 + x2 + 0.040000000000000008*std::sqrt(x0 + 1)); - Real x6 = x5*(3.3470650363800003e-24*X(10) + 5.0186540981899997e-24*X(11) + 5.01956503638e-24*X(12) + 6.6902431600000005e-24*X(13) + 6.6911540981899994e-24*X(14) + 6.6920650363799998e-24*X(15) + 2.007528417094e-23*X(16) + 2.0076195109127999e-23*X(17) + 2.1749727627315999e-23*X(18) + 2.3423260145503998e-23*X(19) + 9.1093818800000008e-28*X(2) + 2.5096792663692001e-23*X(20) + 2.6767349207315999e-23*X(21) + 2.6768260145504001e-23*X(22) + 2.6770325181879997e-23*X(23) + 2.8440881725504001e-23*X(24) + 2.8441792663692003e-23*X(25) + 3.0114414243691998e-23*X(26) + 3.011532518188e-23*X(27) + 3.178794676188e-23*X(28) + 4.6843544316443998e-23*X(29) + 1.6726215800000001e-24*X(3) + 4.6844455254632e-23*X(30) + 5.3535609352819999e-23*X(31) + 5.3536520291008001e-23*X(32) + 7.3612715400136001e-23*X(33) + 1.6735325181900001e-24*X(4) + 1.6744434563800001e-24*X(5) + 3.3451215800000003e-24*X(6) + 3.3460325181899999e-24*X(7) + 3.3461540981899999e-24*X(8) + 3.3469434563800003e-24*X(9))/(8.0000000000000013e-6*x2 + x3 + 0.040000000000000008*std::sqrt(x1 + 1)); + Real x6 = 1195783.1580871048*x5 + 2.0e-99; - Real x7 = 1195783.1580871048*x6 + 2.0e-99; + Real x7 = X(16)*X(32); - Real x8 = 2.0e-99*X(1); + Real x8 = 4.5299999999999999e-10*x7; - Real x9 = 1491378.0618531078*x6 + 2.0e-99; + Real x9 = std::sqrt(x0); - Real x10 = 1.0/x4; + Real x10 = 1.0/x0; - Real x11 = 5.7884371785482823e-10*x10*std::pow(0.00060040841663220993*x4 + 1.0, -1.7524)*std::pow(0.32668576019240059*x4 + 1.0, -0.24759999999999999) + 1.0e-99; + Real x11 = 1.6974097914175e-15*x9*std::exp(-7550.0*x10) + 1.0e-99; - Real x12 = X(13)*X(2); + Real x12 = X(10)*X(33); - Real x13 = -x11*x12; + Real x13 = x11*x12; - Real x14 = std::pow(x1, -0.69999999999999996); + Real x14 = 7.5e-10*X(29)*X(4); - Real x15 = X(2)*X(31); + Real x15 = 1.0e-99 + 2.8e-11*std::exp(-26500.0*x10); - Real x16 = x15*(-1.0568649366703787e-5*x14 - 1.0e-99); + Real x16 = X(22)*X(33); - Real x17 = 7.1999999999999996e-8*x10 + 1.0e-99; + Real x17 = 1.0e-99 + 4.1999999999999999e-12*std::exp(-24000.0*x10); - Real x18 = X(11)*X(2); + Real x18 = X(30)*X(32); - Real x19 = -x17*x18; + Real x19 = x15*x16 - x17*x18; - Real x20 = 9.7355109978586414e-8*std::pow(x1, -0.47999999999999998) + 1.0e-99; + Real x20 = 1.0e-13*X(25)*X(30); - Real x21 = X(2)*X(24); + Real x21 = 1.0e-99 + 2.54e-10*std::exp(-13300.0*x10); - Real x22 = -x20*x21; + Real x22 = X(33)*X(4); - Real x23 = std::log(x1); + Real x23 = -x20 + x21*x22; - Real x24 = 0.10684732509875319*x23 - 1; + Real x24 = std::pow(x0, -0.34000000000000002); - Real x25 = ((x24)*(x24)); + Real x25 = ((x0 <= 1052.0) ? ( + 3.2682787072746432e-10*x24 +) +: ( + 0 +)) + 1.0e-99; - Real x26 = ((x24)*(x24)*(x24)); + Real x26 = X(17)*X(32); - Real x27 = ((((x24)*(x24)))*(((x24)*(x24)))); + Real x27 = ((x0 > 1052.0) ? ( + 3.7991369449339675e-16*std::pow(x0, 1.54)*std::exp(613.0*x10) +) +: ( + 0 +)) + 1.0e-99; - Real x28 = ((x24)*(x24)*(x24)*(x24)*(x24)); + Real x28 = x25*x26 + x26*x27; - Real x29 = std::exp((6)*std::log(std::abs(x24))); + Real x29 = x0 <= 10.0; - Real x30 = ((x24)*(x24)*(x24)*(x24)*(x24)*(x24)*(x24)); + Real x30 = ((x29) ? ( + 7.0509999999999997e-11 +) +: ( + 0 +)) + 1.0e-99; - Real x31 = std::exp((8)*std::log(std::abs(x24))); + Real x31 = X(17)*X(25); - Real x32 = std::pow(x1, 23.915965629999999)*std::exp(-941.91483008144996*x25 + 2506.9866529060901*x26 - 4361.9927099007555*x27 + 4879.7345146260486*x28 - 3366.4639698826941*x29 + 1300.3028484326148*x30 - 214.82451513312137*x31); + Real x32 = x0 > 10.0; - Real x33 = 4.3524079114767552e-117*x32 + 1.0e-99; + Real x33 = ((x32) ? ( + 1.5557027645121538e-10*std::pow(x0, -0.33900000000000002)*std::exp(-0.108*x10) +) +: ( + 0 +)) + 1.0e-99; - Real x34 = X(15)*X(2); + Real x34 = 1.0e-99 + 1.02e-10*std::exp(-914.0*x10); - Real x35 = -x33*x34; + Real x35 = X(18)*X(22); - Real x36 = std::pow(x1, -0.75); + Real x36 = x30*x31 + x31*x33 + x34*x35; - Real x37 = 2.5950363272655348e-10*x36 + 1.0e-99; + Real x37 = -X(30)*x6 + x13 + x14 + x19 + x23 + x28 + x36 + x8; - Real x38 = 2*X(2)*X(6); + Real x38 = 1491378.0618531078*x5 + 2.0e-99; - Real x39 = -x37*x38; + Real x39 = ((X(25))*(X(25))); - Real x40 = 6.3350403285625223e-6*std::pow(x1, -0.55000000000000004) + 1.0e-99; + Real x40 = std::pow(x0, 1.1399999999999999)*std::exp(-50.0*x10); - Real x41 = X(2)*X(29); + Real x41 = x39*(2.4749459212614312e-15*x40 + 1.0e-99); - Real x42 = -x40*x41; + Real x42 = ((x29) ? ( + 4.2020000000000002e-8 +) +: ( + 0 +)) + 1.0e-99; - Real x43 = std::pow(x1, 2.360852208681)*std::exp(-24.766609674457612*x25 + 13.307984239358756*x26 - 258.18559308467116*x27 + 846.15238706523724*x28 - 1113.0879095147111*x29 + 671.95094388835207*x30 - 154.90262957142161*x31); + Real x43 = X(27)*X(3); - Real x44 = 3.7903999274394518e-18*x43 + 1.0e-99; + Real x44 = 1.0/x9; - Real x45 = X(2)*X(5); + Real x45 = ((x32) ? ( + 1.2730571616977901e-7*x44 + 1.767e-9 +) +: ( + 0 +)) + 1.0e-99; - Real x46 = -x44*x45; + Real x46 = -x42*x43 - x43*x45; - Real x47 = 1.3300135414628029e-18*std::pow(x1, 0.94999999999999996)*std::exp(-0.00010729613733905579*x1) + 1.0e-99; + Real x47 = 3.5999999999999998e-16*std::pow(T, 1.52)*std::exp(-1740.0*x10) + 1.0e-99; - Real x48 = X(2)*X(7); + Real x48 = X(10)*X(25); - Real x49 = -x47*x48; + Real x49 = 1.8706148721743872e-6*x44 + 1.0e-99; - Real x50 = std::pow(x1, 0.25); + Real x50 = X(2)*X(28); - Real x51 = 1.0/x50; + Real x51 = 6.3600000000000002e-14*x0*std::exp(-9610.0*x10) + 1.0e-99; - Real x52 = 1.0/x1; + Real x52 = X(27)*X(4); - Real x53 = 130666.72639118577*x52; + Real x53 = x47*x48 + x49*x50 - x51*x52; - Real x54 = x51*std::exp(-x53)/(x53 + 0.193); + Real x54 = -X(27)*x38 + x13 + x41 + x46 + x53; - Real x55 = 1.3023623244766063e-6*x54 + 1.0e-99; + Real x55 = 1.0/x3; - Real x56 = X(17)*X(2); + Real x56 = 5.7884371785482823e-10*x55*std::pow(0.00060040841663220993*x3 + 1.0, -1.7524)*std::pow(0.32668576019240059*x3 + 1.0, -0.24759999999999999) + 1.0e-99; - Real x57 = -x55*x56; + Real x57 = X(13)*X(2); - Real x58 = std::pow(x1, -0.34000000000000002); + Real x58 = -x56*x57; - Real x59 = 157821.26811013557*x52; + Real x59 = std::pow(x0, -0.69999999999999996); - Real x60 = x58*std::exp(-x59)/(x59 + 0.072999999999999995); + Real x60 = X(2)*X(31); - Real x61 = 2.1012161411986507e-6*x60 + 1.0e-99; + Real x61 = x60*(-1.0568649366703787e-5*x59 - 1.0e-99); - Real x62 = X(2)*X(22); + Real x62 = 7.1999999999999996e-8*x55 + 1.0e-99; - Real x63 = -x61*x62; + Real x63 = X(11)*X(2); - Real x64 = std::pow(x1, 0.34999999999999998)*std::exp(-102000.0*x52); + Real x64 = -x62*x63; - Real x65 = 4.3799999999999999e-10*x64 + 1.0e-99; + Real x65 = 9.7355109978586414e-8*std::pow(x0, -0.47999999999999998) + 1.0e-99; - Real x66 = X(10)*X(2); + Real x66 = X(2)*X(24); - Real x67 = std::pow(x1, 43.933476326349997)*std::exp(-1618.789587733125*x25 + 3854.4033653120221*x26 - 5902.1601240760483*x27 + 5825.9326359379538*x28 - 3578.1439181805955*x29 + 1242.7294446825149*x30 - 186.35635455381879*x31); + Real x67 = -x65*x66; - Real x68 = X(14)*X(2); + Real x68 = std::log(x0); - Real x69 = std::pow(x1, 13.536555999999999)*std::exp(-502.72883252679094*x25 + 1281.477767828706*x26 - 2207.4643501257692*x27 + 2500.8077583366976*x28 - 1768.8867461266502*x29 + 704.19926629500367*x30 - 120.0438480494693*x31); + Real x69 = 0.10684732509875319*x68 - 1; - Real x70 = X(2)*X(4); + Real x70 = ((x69)*(x69)); - Real x71 = ((x1 <= 1160.0) ? ( - 1.4643482606109061e-16*std::pow(x1, 1.78186) -) -: ( - 0 -)); + Real x71 = ((x69)*(x69)*(x69)); + + Real x72 = ((((x69)*(x69)))*(((x69)*(x69)))); + + Real x73 = ((x69)*(x69)*(x69)*(x69)*(x69)); + + Real x74 = std::exp((6)*std::log(std::abs(x69))); + + Real x75 = ((x69)*(x69)*(x69)*(x69)*(x69)*(x69)*(x69)); + + Real x76 = std::exp((8)*std::log(std::abs(x69))); + + Real x77 = std::pow(x0, 23.915965629999999)*std::exp(-941.91483008144996*x70 + 2506.9866529060901*x71 - 4361.9927099007555*x72 + 4879.7345146260486*x73 - 3366.4639698826941*x74 + 1300.3028484326148*x75 - 214.82451513312137*x76); + + Real x78 = 4.3524079114767552e-117*x77 + 1.0e-99; + + Real x79 = X(15)*X(2); + + Real x80 = -x78*x79; + + Real x81 = std::pow(x0, -0.75); + + Real x82 = 2.5950363272655348e-10*x81 + 1.0e-99; + + Real x83 = 2*X(2)*X(6); + + Real x84 = -x82*x83; + + Real x85 = 6.3350403285625223e-6*std::pow(x0, -0.55000000000000004) + 1.0e-99; + + Real x86 = X(2)*X(29); + + Real x87 = -x85*x86; + + Real x88 = std::pow(x0, 2.360852208681)*std::exp(-24.766609674457612*x70 + 13.307984239358756*x71 - 258.18559308467116*x72 + 846.15238706523724*x73 - 1113.0879095147111*x74 + 671.95094388835207*x75 - 154.90262957142161*x76); + + Real x89 = 3.7903999274394518e-18*x88 + 1.0e-99; + + Real x90 = X(2)*X(5); + + Real x91 = -x89*x90; + + Real x92 = 1.3300135414628029e-18*std::pow(x0, 0.94999999999999996)*std::exp(-0.00010729613733905579*x0) + 1.0e-99; + + Real x93 = X(2)*X(7); + + Real x94 = -x92*x93; + + Real x95 = std::pow(x0, 0.25); + + Real x96 = 1.0/x95; + + Real x97 = 130666.72639118577*x10; + + Real x98 = x96*std::exp(-x97)/(x97 + 0.193); + + Real x99 = 1.3023623244766063e-6*x98 + 1.0e-99; + + Real x100 = X(17)*X(2); + + Real x101 = -x100*x99; + + Real x102 = 157821.26811013557*x10; + + Real x103 = x24*std::exp(-x102)/(x102 + 0.072999999999999995); + + Real x104 = 2.1012161411986507e-6*x103 + 1.0e-99; + + Real x105 = X(2)*X(22); + + Real x106 = -x104*x105; + + Real x107 = std::pow(x0, 0.34999999999999998)*std::exp(-102000.0*x10); + + Real x108 = 4.3799999999999999e-10*x107 + 1.0e-99; + + Real x109 = X(10)*X(2); + + Real x110 = std::pow(x0, 43.933476326349997)*std::exp(-1618.789587733125*x70 + 3854.4033653120221*x71 - 5902.1601240760483*x72 + 5825.9326359379538*x73 - 3578.1439181805955*x74 + 1242.7294446825149*x75 - 186.35635455381879*x76); - Real x72 = x71 + 1.0e-99; + Real x111 = X(14)*X(2); - Real x73 = X(4)*X(5); + Real x112 = std::pow(x0, 13.536555999999999)*std::exp(-502.72883252679094*x70 + 1281.477767828706*x71 - 2207.4643501257692*x72 + 2500.8077583366976*x73 - 1768.8867461266502*x74 + 704.19926629500367*x75 - 120.0438480494693*x76); - Real x74 = ((x24)*(x24)*(x24)*(x24)*(x24)*(x24)*(x24)*(x24)*(x24)); + Real x113 = X(2)*X(4); - Real x75 = ((x1 > 1160.0) ? ( - 3.3178155742407601e-14*std::pow(x1, 1.139449335841631)*std::exp(-12.447178055372776*x25 + 6.9391784778399117*x26 - 10.993097527150175*x27 + 14.449862906216714*x28 + 58.228375789703179*x29 - 162.59852239006702*x30 + 144.55426734953477*x31 - 44.454280878123605*x74) + Real x114 = ((x0 <= 1160.0) ? ( + 1.4643482606109061e-16*std::pow(x0, 1.78186) ) : ( 0 )); - Real x76 = x75 + 1.0e-99; + Real x115 = x114 + 1.0e-99; - Real x77 = std::sqrt(x1); + Real x116 = X(4)*X(5); - Real x78 = 1.0/x77; + Real x117 = ((x69)*(x69)*(x69)*(x69)*(x69)*(x69)*(x69)*(x69)*(x69)); - Real x79 = 4.4686910835277033e-6*x78 + 1.0e-99; + Real x118 = ((x0 > 1160.0) ? ( + 3.3178155742407601e-14*std::pow(x0, 1.139449335841631)*std::exp(-44.454280878123605*x117 - 12.447178055372776*x70 + 6.9391784778399117*x71 - 10.993097527150175*x72 + 14.449862906216714*x73 + 58.228375789703179*x74 - 162.59852239006702*x75 + 144.55426734953477*x76) +) +: ( + 0 +)); - Real x80 = X(2)*X(28); + Real x119 = x118 + 1.0e-99; - Real x81 = 1.8706148721743872e-6*x78 + 1.0e-99; + Real x120 = 4.4686910835277033e-6*x44 + 1.0e-99; - Real x82 = -x79*x80 - x80*x81; + Real x121 = -x120*x50 - x49*x50; - Real x83 = 1.4895636945092344e-6*x78 + 1.0e-99; + Real x122 = 1.4895636945092344e-6*x44 + 1.0e-99; - Real x84 = X(2)*X(26); + Real x123 = X(2)*X(26); - Real x85 = 6.7549981495186206e-7*x78 + 1.0e-99; + Real x124 = 6.7549981495186206e-7*x44 + 1.0e-99; - Real x86 = -x83*x84 - x84*x85; + Real x125 = -x122*x123 - x123*x124; - Real x87 = x1 <= 10000.0; + Real x126 = x0 <= 10000.0; - Real x88 = ((x87) ? ( - -5.5279999999999998e-28*((x1)*(x1)*(x1)*(x1)*(x1)) + 3.3467999999999999e-23*((((x1)*(x1)))*(((x1)*(x1)))) - 7.5474000000000004e-19*((x1)*(x1)*(x1)) - 2.3088e-11*x1 + 7.3427999999999993e-15*x2 + 4.2277999999999996e-8 + Real x127 = ((x126) ? ( + -5.5279999999999998e-28*((x0)*(x0)*(x0)*(x0)*(x0)) + 3.3467999999999999e-23*((((x0)*(x0)))*(((x0)*(x0)))) - 7.5474000000000004e-19*((x0)*(x0)*(x0)) - 2.3088e-11*x0 + 7.3427999999999993e-15*x1 + 4.2277999999999996e-8 ) : ( 0 )); - Real x89 = X(2)*X(8); + Real x128 = X(2)*X(8); - Real x90 = 1.0e-8*std::pow(x1, -0.40000000000000002) + 1.0e-99; + Real x129 = 1.0e-8*std::pow(x0, -0.40000000000000002) + 1.0e-99; - Real x91 = X(3)*X(5); + Real x130 = X(3)*X(5); - Real x92 = x89*(-x88 - 1.0e-99) + x90*x91; + Real x131 = x128*(-x127 - 1.0e-99) + x129*x130; - Real x93 = ((x1 <= 400.0) ? ( - 1.2999999999999999e-10*std::pow(x1, -0.64000000000000001) + Real x132 = ((x0 <= 400.0) ? ( + 1.2999999999999999e-10*std::pow(x0, -0.64000000000000001) ) : ( 0 )) + 1.0e-99; - Real x94 = X(2)*X(21); + Real x133 = X(2)*X(21); - Real x95 = std::pow(x1, -1.5); + Real x134 = std::pow(x0, -1.5); - Real x96 = ((x1 > 400.0) ? ( - 1.41e-10*std::pow(x1, -0.66000000000000003) + 0.00073999999999999999*x95*(1.0 + 0.062*std::exp(-145000.0*x52))*std::exp(-175000.0*x52) + Real x135 = ((x0 > 400.0) ? ( + 1.41e-10*std::pow(x0, -0.66000000000000003) + 0.00073999999999999999*x134*(1.0 + 0.062*std::exp(-145000.0*x10))*std::exp(-175000.0*x10) ) : ( 0 )) + 1.0e-99; - Real x97 = -x93*x94 - x94*x96; + Real x136 = -x132*x133 - x133*x135; - Real x98 = 1.4981088130721367e-10*std::pow(x1, -0.63529999999999998); + Real x137 = 1.4981088130721367e-10*std::pow(x0, -0.63529999999999998); - Real x99 = x1 > 2.73; + Real x138 = x0 > 2.73; - Real x100 = ((x99 && x1 <= 9280.0) ? ( - x98 + Real x139 = ((x138 && x0 <= 9280.0) ? ( + x137 ) : ( 0 )) + 1.0e-99; - Real x101 = x1 <= 100000000.0; + Real x140 = x0 <= 100000000.0; - Real x102 = ((x101 && x1 > 9280.0) ? ( - 1250086.112245841*x95*(1.5400000000000001e-9 + 4.6200000000000001e-10*std::exp(-93988.701501924661*x52))*std::exp(-469943.50750964211*x52) + x98 + Real x141 = ((x140 && x0 > 9280.0) ? ( + 1250086.112245841*x134*(1.5400000000000001e-9 + 4.6200000000000001e-10*std::exp(-93988.701501924661*x10))*std::exp(-469943.50750964211*x10) + x137 ) : ( 0 )) + 1.0e-99; - Real x103 = 3.8571873359681448e-209*x67 + 1.0e-99; + Real x142 = 3.8571873359681448e-209*x110 + 1.0e-99; - Real x104 = -x100*x68 - x102*x68 - x103*x68; + Real x143 = -x111*x139 - x111*x141 - x111*x142; - Real x105 = 2.6534040307116387e-9*std::pow(x1, -0.10000000000000001) + 1.0e-99; + Real x144 = 2.6534040307116387e-9*std::pow(x0, -0.10000000000000001) + 1.0e-99; - Real x106 = X(4)*X(9); + Real x145 = X(4)*X(9); - Real x107 = X(5)*X(7); + Real x146 = X(5)*X(7); - Real x108 = x105*x106 + x105*x107; + Real x147 = x144*x145 + x144*x146; - Real x109 = ((x1 <= 7950.0) ? ( - 1.4308352583277889e-10*std::pow(x1, -0.59999999999999998) + Real x148 = ((x0 <= 7950.0) ? ( + 1.4308352583277889e-10*std::pow(x0, -0.59999999999999998) ) : ( 0 )) + 1.0e-99; - Real x110 = X(16)*X(2); + Real x149 = X(16)*X(2); - Real x111 = ((x1 > 21140.0) ? ( - 0.00023813936486848612*std::pow(x1, -1.3700000000000001)*std::exp(-115786.2*x52) + Real x150 = ((x0 > 21140.0) ? ( + 0.00023813936486848612*std::pow(x0, -1.3700000000000001)*std::exp(-115786.2*x10) ) : ( 0 )) + 1.0e-99; - Real x112 = ((x1 <= 21140.0 && x1 > 7950.0) ? ( - 8.3535906003219973e-24*std::pow(x1, 2.4900000000000002)*std::exp(21845.599999999999*x52) + Real x151 = ((x0 <= 21140.0 && x0 > 7950.0) ? ( + 8.3535906003219973e-24*std::pow(x0, 2.4900000000000002)*std::exp(21845.599999999999*x10) ) : ( 0 )) + 1.0e-99; - Real x113 = -x109*x110 - x110*x111 - x110*x112; + Real x152 = -x148*x149 - x149*x150 - x149*x151; - Real x114 = ((x99 && x1 <= 5500.0) ? ( - x98 + Real x153 = ((x138 && x0 <= 5500.0) ? ( + x137 ) : ( 0 )) + 1.0e-99; - Real x115 = X(2)*X(3); + Real x154 = X(2)*X(3); - Real x116 = ((x101 && x1 > 5500.0) ? ( - 3.2867337024382733e-10*std::pow(x1, -0.7241125657826851)*std::exp(-1.7746868094247411*x25 - 1.9518356165136789*x26 - 2.4649195146505537*x27 - 1.0207737270119371*x28 + 3.3530579587656565*x29 + 3.6203127646377788*x30 - 1.0930705283186734*x31 - 1.6921001126637108*x74) + Real x155 = ((x140 && x0 > 5500.0) ? ( + 3.2867337024382733e-10*std::pow(x0, -0.7241125657826851)*std::exp(-1.6921001126637108*x117 - 1.7746868094247411*x70 - 1.9518356165136789*x71 - 2.4649195146505537*x72 - 1.0207737270119371*x73 + 3.3530579587656565*x74 + 3.6203127646377788*x75 - 1.0930705283186734*x76) ) : ( 0 )) + 1.0e-99; - Real x117 = -x114*x115 - x115*x116; - - Real x118 = 35.5*std::pow(x1, -2.2799999999999998)*std::exp(-46707.0*x52) + 1.0e-99; + Real x156 = -x153*x154 - x154*x155; - Real x119 = (1.3500000000000001e-9*std::pow(x1, 0.098492999999999997) + 4.4350199999999998e-10*std::pow(x1, 0.55610000000000004) + 3.7408500000000004e-16*std::pow(x1, 2.1825999999999999))/(0.0061910000000000003*std::pow(x1, 1.0461) + 8.9711999999999997e-11*std::pow(x1, 3.0424000000000002) + 3.2575999999999999e-14*std::pow(x1, 3.7740999999999998) + 1.0) + 1.0e-99; + Real x157 = 35.5*std::pow(x0, -2.2799999999999998)*std::exp(-46707.0*x10) + 1.0e-99; - Real x120 = -x118*x66 + x119*x73 - x65*x66; + Real x158 = (1.3500000000000001e-9*std::pow(x0, 0.098492999999999997) + 4.4350199999999998e-10*std::pow(x0, 0.55610000000000004) + 3.7408500000000004e-16*std::pow(x0, 2.1825999999999999))/(0.0061910000000000003*std::pow(x0, 1.0461) + 8.9711999999999997e-11*std::pow(x0, 3.0424000000000002) + 3.2575999999999999e-14*std::pow(x0, 3.7740999999999998) + 1.0) + 1.0e-99; - Real x121 = 5.9082438637265071e-70*x69 + 1.0e-99; + Real x159 = -x108*x109 - x109*x157 + x116*x158; - Real x122 = 1.4000000000000001e-18*std::pow(x1, 0.92800000000000005)*std::exp(-6.1728395061728397e-5*x1) + 1.0e-99; + Real x160 = 5.9082438637265071e-70*x112 + 1.0e-99; - Real x123 = -x121*x70 - x122*x70; + Real x161 = 1.4000000000000001e-18*std::pow(x0, 0.92800000000000005)*std::exp(-6.1728395061728397e-5*x0) + 1.0e-99; - Real x124 = 7.9674337148168363e-7*x78 + 1.0e-99; + Real x162 = -x113*x160 - x113*x161; - Real x125 = X(3)*X(9); + Real x163 = 7.9674337148168363e-7*x44 + 1.0e-99; - Real x126 = -x124*x125; + Real x164 = X(3)*X(9); - Real x127 = 2.0000000000000001e-9*X(3)*X(32); + Real x165 = -x163*x164; - Real x128 = -x127; + Real x166 = 2.0000000000000001e-9*X(3)*X(32); - Real x129 = 7.5e-10*X(29)*X(4); + Real x167 = -x166; - Real x130 = 2.8833736969617052e-16*x50 + 1.0e-99; + Real x168 = 2.8833736969617052e-16*x95 + 1.0e-99; - Real x131 = X(14)*X(4); + Real x169 = X(14)*X(4); - Real x132 = ((x1 > 10000.0) ? ( - 4.0000000000000003e-37*std::pow(x1, 4.7400000000000002) + Real x170 = ((x0 > 10000.0) ? ( + 4.0000000000000003e-37*std::pow(x0, 4.7400000000000002) ) : ( 0 )) + 1.0e-99; - Real x133 = X(15)*X(3); + Real x171 = X(15)*X(3); - Real x134 = ((x87) ? ( - 1.26e-9*x36*std::exp(-127500.0*x52) + Real x172 = ((x126) ? ( + 1.26e-9*x81*std::exp(-127500.0*x10) ) : ( 0 )) + 1.0e-99; - Real x135 = x130*x131 - x132*x133 - x133*x134; - - Real x136 = 6.3999999999999996e-10*X(11)*X(4); - - Real x137 = M_LN10; - - Real x138 = 1.0/x137; - - Real x139 = x138*x23; + Real x173 = x168*x169 - x170*x171 - x171*x172; - Real x140 = ((x23)*(x23)); + Real x174 = 6.3999999999999996e-10*X(11)*X(4); - Real x141 = std::exp((-2)*std::log(std::abs(x137))); + Real x175 = M_LN10; - Real x142 = x140*x141; + Real x176 = 1.0/x175; - Real x143 = 8.4600000000000008e-10*x139 - 1.3700000000000002e-10*x142 + 4.1700000000000001e-10; + Real x177 = x176*x68; - Real x144 = X(10)*X(6); + Real x178 = ((x68)*(x68)); - Real x145 = 1.0e-99 + 1.0000000000000001e-9*std::exp(-457.0*x52); + Real x179 = std::exp((-2)*std::log(std::abs(x175))); - Real x146 = X(12)*X(3); + Real x180 = x178*x179; - Real x147 = x136 + x143*x144 - x145*x146; + Real x181 = 8.4600000000000008e-10*x177 - 1.3700000000000002e-10*x180 + 4.1700000000000001e-10; - Real x148 = x1 <= 10.0; + Real x182 = X(10)*X(6); - Real x149 = ((x148) ? ( - 4.2020000000000002e-8 -) -: ( - 0 -)) + 1.0e-99; + Real x183 = 1.0e-99 + 1.0000000000000001e-9*std::exp(-457.0*x10); - Real x150 = X(27)*X(3); + Real x184 = X(12)*X(3); - Real x151 = x1 > 10.0; + Real x185 = x174 + x181*x182 - x183*x184; - Real x152 = ((x151) ? ( - 1.2730571616977901e-7*x78 + 1.767e-9 -) -: ( - 0 -)) + 1.0e-99; - - Real x153 = -x149*x150 - x150*x152; - - Real x154 = ((x151 && x1 <= 100000.0) ? ( - -7.7700000000000002e-13*x1 + 2.96e-6*x10 + 2.5000000000000002e-10*x4 - 1.73e-9 + Real x186 = ((x32 && x0 <= 100000.0) ? ( + -7.7700000000000002e-13*x0 + 2.5000000000000002e-10*x3 + 2.96e-6*x55 - 1.73e-9 ) : ( 0 )); - Real x155 = -x90*x91 + x91*(-x154 - 1.0e-99); + Real x187 = -x129*x130 + x130*(-x186 - 1.0e-99); - Real x156 = std::pow(10.0, -0.12690000000000001*amrex::Math::powi<-3>(x137)*((x23)*(x23)*(x23)) + 1.1180000000000001*amrex::Math::powi<-2>(x137)*((x23)*(x23)) - 1.5229999999999999*x139 - 19.379999999999999) + 1.0e-99; + Real x188 = std::pow(10.0, -0.12690000000000001*amrex::Math::powi<-3>(x175)*((x68)*(x68)*(x68)) + 1.1180000000000001*amrex::Math::powi<-2>(x175)*((x68)*(x68)) - 1.5229999999999999*x177 - 19.379999999999999) + 1.0e-99; - Real x157 = -x156; + Real x189 = -x188; - Real x158 = X(3)*X(7); + Real x190 = X(3)*X(7); - Real x159 = x1 > 50.0; + Real x191 = x0 > 50.0; - Real x160 = ((x159) ? ( - 2.0000000000000001e-10*std::pow(x1, 0.40200000000000002)*std::exp(-37.100000000000001*x52) - 3.3099999999999998e-17*std::pow(x1, 1.48) + Real x192 = ((x191) ? ( + 2.0000000000000001e-10*std::pow(x0, 0.40200000000000002)*std::exp(-37.100000000000001*x10) - 3.3099999999999998e-17*std::pow(x0, 1.48) ) : ( 0 )) + 1.0e-99; - Real x161 = ((x159) ? ( - 2.0299999999999998e-9*std::pow(x1, -0.33200000000000002) + 2.0600000000000001e-10*std::pow(x1, 0.39600000000000002)*std::exp(-33.0*x52) + Real x193 = ((x191) ? ( + 2.0299999999999998e-9*std::pow(x0, -0.33200000000000002) + 2.0600000000000001e-10*std::pow(x0, 0.39600000000000002)*std::exp(-33.0*x10) ) : ( 0 )) + 1.0e-99; - Real x162 = X(4)*X(6); + Real x194 = X(4)*X(6); - Real x163 = x157*x158 - x158*x160 + x161*x162; + Real x195 = x189*x190 - x190*x192 + x193*x194; - Real x164 = 8.7882738061352424e-22*std::pow(x1, 1.96)*std::exp(-170000.0*x52) + 1.0e-99; + Real x196 = 8.7882738061352424e-22*std::pow(x0, 1.96)*std::exp(-170000.0*x10) + 1.0e-99; - Real x165 = X(16)*X(4); + Real x197 = X(16)*X(4); - Real x166 = 3.8999999999999998e-16*std::pow(x1, 0.21299999999999999) + 1.0e-99; + Real x198 = 3.8999999999999998e-16*std::pow(x0, 0.21299999999999999) + 1.0e-99; - Real x167 = X(17)*X(3); + Real x199 = X(17)*X(3); - Real x168 = x164*x165 - x166*x167; + Real x200 = x196*x197 - x198*x199; - Real x169 = ((x148) ? ( + Real x201 = ((x29) ? ( 3.7450000000000001e-8 ) : ( 0 )) + 1.0e-99; - Real x170 = X(25)*X(3); + Real x202 = X(25)*X(3); - Real x171 = ((x151) ? ( - 1.1352943524561152e-7*x78 + 1.55e-9 + Real x203 = ((x32) ? ( + 1.1352943524561152e-7*x44 + 1.55e-9 ) : ( 0 )) + 1.0e-99; - Real x172 = -x169*x170 - x170*x171; + Real x204 = -x201*x202 - x202*x203; - Real x173 = 7.5399999999999998e-10*std::pow(x1, -0.45800000000000002) + 4.9899999999999997e-11*std::pow(x1, 0.40500000000000003) + 1.0e-99; + Real x205 = 7.5399999999999998e-10*std::pow(x0, -0.45800000000000002) + 4.9899999999999997e-11*std::pow(x0, 0.40500000000000003) + 1.0e-99; - Real x174 = X(21)*X(4); + Real x206 = X(21)*X(4); - Real x175 = (4.0000000000000001e-10*std::pow(x1, 0.0066899999999999998) + 1.0799999999999999e-11*std::pow(x1, 0.51700000000000002))*std::exp(-227.0*x52) + 1.0e-99; + Real x207 = (4.0000000000000001e-10*std::pow(x0, 0.0066899999999999998) + 1.0799999999999999e-11*std::pow(x0, 0.51700000000000002))*std::exp(-227.0*x10) + 1.0e-99; - Real x176 = X(22)*X(3); + Real x208 = X(22)*X(3); - Real x177 = x173*x174 - x175*x176; + Real x209 = x205*x206 - x207*x208; - Real x178 = 6.0e-10*X(4)*X(8); + Real x210 = 6.0e-10*X(4)*X(8); - Real x179 = ((x23)*(x23)*(x23)); + Real x211 = ((x68)*(x68)*(x68)); - Real x180 = ((((x23)*(x23)))*(((x23)*(x23)))); + Real x212 = ((((x68)*(x68)))*(((x68)*(x68)))); - Real x181 = ((x23)*(x23)*(x23)*(x23)*(x23)); + Real x213 = ((x68)*(x68)*(x68)*(x68)*(x68)); - Real x182 = ((x1 <= 30000.0 && x1 > 100.0) ? ( - (-1.4491368e-7*x140 + 3.4172804999999998e-8*x179 - 4.7813727999999997e-9*x180 + 3.9731542e-10*x181 + 3.5311931999999998e-13*((x23)*(x23)*(x23)*(x23)*(x23)*(x23)*(x23)) - 1.8171411000000001e-11*std::exp((6)*std::log(std::abs(x23))) + 3.3735381999999997e-7*x23 - 3.3232183000000002e-7)*std::exp(-21237.150000000001*x52) + Real x214 = ((x0 <= 30000.0 && x0 > 100.0) ? ( + (-1.4491368e-7*x178 + 3.4172804999999998e-8*x211 - 4.7813727999999997e-9*x212 + 3.9731542e-10*x213 + 3.5311931999999998e-13*((x68)*(x68)*(x68)*(x68)*(x68)*(x68)*(x68)) - 1.8171411000000001e-11*std::exp((6)*std::log(std::abs(x68))) + 3.3735381999999997e-7*x68 - 3.3232183000000002e-7)*std::exp(-21237.150000000001*x10) ) : ( 0 )) + 1.0e-99; - Real x183 = X(10)*X(3); + Real x215 = X(10)*X(3); - Real x184 = x178 - x182*x183; + Real x216 = x210 - x214*x215; - Real x185 = x179/((x137)*(x137)*(x137)); + Real x217 = x211/((x175)*(x175)*(x175)); - Real x186 = ((x1 > 30.0) ? ( - std::pow(10, -3.194*x139 + 1.786*x140*x141 - 0.2072*x185 - 18.199999999999999) + Real x218 = ((x0 > 30.0) ? ( + std::pow(10, -3.194*x177 + 1.786*x178*x179 - 0.2072*x217 - 18.199999999999999) ) : ( 0 )) + 1.0e-99; - Real x187 = X(3)*X(4); + Real x219 = X(3)*X(4); - Real x188 = std::pow(x1, -0.14999999999999999); + Real x220 = std::pow(x0, -0.14999999999999999); - Real x189 = ((x1 <= 30.0) ? ( - 3.4977396723747635e-20*x188 + Real x221 = ((x0 <= 30.0) ? ( + 3.4977396723747635e-20*x220 ) : ( 0 )) + 1.0e-99; - Real x190 = -x186*x187 - x187*x189; - - Real x191 = 6.3999999999999996e-10*X(10)*X(26); - - Real x192 = -x136; - - Real x193 = 1.0000000000000001e-17*X(17); + Real x222 = -x218*x219 - x219*x221; - Real x194 = X(4)*x193; + Real x223 = 6.3999999999999996e-10*X(10)*X(26); - Real x195 = -x194; + Real x224 = -x174; - Real x196 = ((X(10))*(X(10))); + Real x225 = 1.0000000000000001e-17*X(17); - Real x197 = std::log(T); + Real x226 = X(4)*x225; - Real x198 = 2.0*X(10) + X(11) + X(12) + X(18) + 2.0*X(19) + 3.0*X(20) + 4.0*X(23) + X(24) + X(25) + 2.0*X(26) + 2.0*X(27) + 3.0*X(28) + X(3) + X(4) + X(5) + 2.0*X(8); + Real x227 = -x226; - Real x199 = 1.0/(std::pow(10.0, 1.3*x138*(x197 - 9.2103403719761836) - 137.42519902360013*x141*((0.10857362047581294*x197 - 1)*(0.10857362047581294*x197 - 1)) - 4.8449999999999998)*x198 + 1.0); + Real x228 = ((X(10))*(X(10))); - Real x200 = std::pow(1.1800000000000001e-10*std::exp(-69500.0*x52), x199)*std::pow(8.1250000000000003e-8*std::pow(T, -0.5)*(1.0 - std::exp(-6000.0*x52))*std::exp(-52000.0*x52), 1.0 - x199); + Real x229 = std::log(T); - Real x201 = 2*x200 + 2.0e-99; + Real x230 = 2.0*X(1) + 2.0*X(10) + X(11) + X(12) + X(18) + 2.0*X(19) + 3.0*X(20) + 4.0*X(23) + X(24) + X(25) + 2.0*X(26) + 3.0*X(28) + X(3) + X(4) + X(5) + 2.0*X(8); - Real x202 = ((X(4))*(X(4))*(X(4))); + Real x231 = 1.0/(std::pow(10.0, 1.3*x176*(x229 - 9.2103403719761836) - 137.42519902360013*x179*((0.10857362047581294*x229 - 1)*(0.10857362047581294*x229 - 1)) - 4.8449999999999998)*x230 + 1.0); - Real x203 = 0.25820969999999999*x185; + Real x232 = std::pow(1.1800000000000001e-10*std::exp(-69500.0*x10), x231)*std::pow(8.1250000000000003e-8*std::pow(T, -0.5)*(1.0 - std::exp(-6000.0*x10))*std::exp(-52000.0*x10), 1.0 - x231); - Real x204 = std::log(27535.310000000001*x52 + 1.0); + Real x233 = 2*x232 + 2.0e-99; - Real x205 = -2.4640089999999999*x138*x23 + 0.19859550000000001*x142 + 743.05999999999995*x52; + Real x234 = ((X(4))*(X(4))*(X(4))); - Real x206 = std::exp(-0.0022727272727272726*x1); + Real x235 = 0.25820969999999999*x217; - Real x207 = std::exp(-0.00054054054054054055*x1); + Real x236 = std::log(27535.310000000001*x10 + 1.0); - Real x208 = 2.9375070000000001*x206 + 0.23588480000000001*x207 + 0.75022860000000002; + Real x237 = 743.05999999999995*x10 - 2.4640089999999999*x176*x68 + 0.19859550000000001*x180; - Real x209 = std::pow(10.0, 21.360939999999999*x138*x204 + 42.707410000000003*x138*x23 - 2.0273650000000001*x142 - x203 - 21467.790000000001*x52 - 1657.4099999999999*x52/(std::pow(std::pow(10.0, -x205 - 8.1313220000000008)*x198, x208) + 1.0) - 142.7664 - (21.360939999999999*x138*x204 + 11.28215*x138*std::log(14254.549999999999*x52 + 1.0) + 70.138370000000009*x139 - 4.7035149999999994*x142 - x203 - 203.11568)/(std::pow(std::pow(10.0, -x205 - 9.3055640000000004)*x198, x208) + 1.0)); + Real x238 = std::exp(-0.0022727272727272726*x0); - Real x210 = 4.6331670000000003*x185; + Real x239 = std::exp(-0.00054054054054054055*x0); - Real x211 = std::log(40870.379999999997*x52 + 1.0); + Real x240 = 2.9375070000000001*x238 + 0.23588480000000001*x239 + 0.75022860000000002; - Real x212 = -4.8909149999999997*x138*x23 + 0.47490300000000002*x142 - 133.82830000000001*x52; + Real x241 = std::pow(10.0, -21467.790000000001*x10 - 1657.4099999999999*x10/(std::pow(std::pow(10.0, -x237 - 8.1313220000000008)*x230, x240) + 1.0) + 21.360939999999999*x176*x236 + 42.707410000000003*x176*x68 - 2.0273650000000001*x180 - x235 - 142.7664 - (21.360939999999999*x176*x236 + 11.28215*x176*std::log(14254.549999999999*x10 + 1.0) + 70.138370000000009*x177 - 4.7035149999999994*x180 - x235 - 203.11568)/(std::pow(std::pow(10.0, -x237 - 9.3055640000000004)*x230, x240) + 1.0)); - Real x213 = -2.0563129999999998*x206 + 0.58640729999999996*x207 + 0.82274429999999998; + Real x242 = 4.6331670000000003*x217; - Real x214 = std::pow(10.0, 69.700860000000006*x138*x211 - 68.422430000000006*x139 + 43.20243*x140*x141 - x210 - 23705.700000000001*x52 - 2080.4099999999999*x52/(std::pow(std::pow(10.0, -x212 - 13.656822)*x198, x213) + 1.0) - 178.4239 - (69.700860000000006*x138*x211 + 19.734269999999999*x138*std::log(16780.950000000001*x52 + 1.0) - 14.509090000000008*x139 + 37.886913*x142 - x210 - 307.31920000000002)/(std::pow(std::pow(10.0, -x212 - 14.82123)*x198, x213) + 1.0)); + Real x243 = std::log(40870.379999999997*x10 + 1.0); - Real x215 = X(10)*X(4); + Real x244 = -133.82830000000001*x10 - 4.8909149999999997*x176*x68 + 0.47490300000000002*x180; - Real x216 = 1.5e-32*x51 + 5.0000000000000004e-32*x78 + 2.0e-99; + Real x245 = -2.0563129999999998*x238 + 0.58640729999999996*x239 + 0.82274429999999998; - Real x217 = X(10)*((X(4))*(X(4))); + Real x246 = std::pow(10.0, -23705.700000000001*x10 - 2080.4099999999999*x10/(std::pow(std::pow(10.0, -x244 - 13.656822)*x230, x245) + 1.0) + 69.700860000000006*x176*x243 - 68.422430000000006*x177 + 43.20243*x178*x179 - x242 - 178.4239 - (69.700860000000006*x176*x243 + 19.734269999999999*x176*std::log(16780.950000000001*x10 + 1.0) - 14.509090000000008*x177 + 37.886913*x180 - x242 - 307.31920000000002)/(std::pow(std::pow(10.0, -x244 - 14.82123)*x230, x245) + 1.0)); - Real x218 = X(4)*x198*x5/((1.0 + 10000.0*std::exp(-600.0/(Tdust + 9.9999999999999993e-41)))*(7.9999999999999996e-6*x2 + x3 + 0.040000000000000001*std::sqrt(Tdust + x1))); + Real x247 = X(10)*X(4); - Real x219 = 1.1194685000500219e-20*std::pow(x1, 2.7400000000000002)*std::exp(-4740.0*x52) + 1.0e-99; + Real x248 = 5.0000000000000004e-32*x44 + 1.5e-32*x96 + 2.0e-99; - Real x220 = X(10)*X(20); + Real x249 = X(10)*((X(4))*(X(4))); - Real x221 = 2.1555555555555562e-20*((x1)*(x1)*(x1))*std::exp(-4045.0*x52) + 1.0e-99; + Real x250 = X(4)*x230*x4/((1.0 + 10000.0*std::exp(-600.0/(Tdust + 9.9999999999999993e-41)))*(7.9999999999999996e-6*x1 + x2 + 0.040000000000000001*std::sqrt(Tdust + x0))); - Real x222 = X(23)*X(4); + Real x251 = 1.1194685000500219e-20*std::pow(x0, 2.7400000000000002)*std::exp(-4740.0*x10) + 1.0e-99; - Real x223 = x219*x220 - x221*x222; + Real x252 = X(10)*X(20); - Real x224 = 1.9643512449766835e-11*std::pow(x1, 0.17000000000000001)*std::exp(-6400.0*x52) + 1.0e-99; + Real x253 = 2.1555555555555562e-20*((x0)*(x0)*(x0))*std::exp(-4045.0*x10) + 1.0e-99; - Real x225 = X(10)*X(19); + Real x254 = X(23)*X(4); - Real x226 = 1.0e-99 + 1.0e-10*std::exp(-7600.0*x52); + Real x255 = x251*x252 - x253*x254; - Real x227 = X(20)*X(4); + Real x256 = 1.9643512449766835e-11*std::pow(x0, 0.17000000000000001)*std::exp(-6400.0*x10) + 1.0e-99; - Real x228 = x224*x225 - x226*x227; + Real x257 = X(10)*X(19); - Real x229 = 2.1999999999999999e-10*X(19)*X(4); + Real x258 = 1.0e-99 + 1.0e-10*std::exp(-7600.0*x10); - Real x230 = 1.0e-99 + 5.4599999999999998e-10*std::exp(-1943.0*x52); + Real x259 = X(20)*X(4); - Real x231 = X(10)*X(18); + Real x260 = x256*x257 - x258*x259; - Real x232 = -x229 + x230*x231; + Real x261 = 2.1999999999999999e-10*X(19)*X(4); - Real x233 = 1.69e-9*X(10)*X(21); + Real x262 = 1.0e-99 + 5.4599999999999998e-10*std::exp(-1943.0*x10); - Real x234 = x169*x170 + x170*x171 + x233; + Real x263 = X(10)*X(18); - Real x235 = 1.0e-13*X(25)*X(30); + Real x264 = -x261 + x262*x263; - Real x236 = 1.0e-99 + 2.54e-10*std::exp(-13300.0*x52); + Real x265 = 1.69e-9*X(10)*X(21); - Real x237 = X(33)*X(4); + Real x266 = x201*x202 + x202*x203 + x265; - Real x238 = x235 - x236*x237; + Real x267 = x20 - x21*x22; - Real x239 = 6.1739095063118665e-10*std::pow(x1, 0.40999999999999998) + 1.0e-99; + Real x268 = 6.1739095063118665e-10*std::pow(x0, 0.40999999999999998) + 1.0e-99; - Real x240 = -x239; + Real x269 = -x268; - Real x241 = -x105; + Real x270 = -x144; - Real x242 = x106*x240 + x106*x241 + x107*x239; + Real x271 = x145*x269 + x145*x270 + x146*x268; - Real x243 = 1.0e-99 + 6.6399999999999998e-10*std::exp(-11700.0*x52); + Real x272 = 1.0e-99 + 6.6399999999999998e-10*std::exp(-11700.0*x10); - Real x244 = X(10)*X(17); + Real x273 = X(10)*X(17); - Real x245 = 1.0e-99 + 1.3100000000000001e-10*std::exp(-80.0*x52); + Real x274 = 1.0e-99 + 1.3100000000000001e-10*std::exp(-80.0*x10); - Real x246 = X(18)*X(4); + Real x275 = X(18)*X(4); - Real x247 = x243*x244 - x245*x246; + Real x276 = x272*x273 - x274*x275; - Real x248 = 1.01e-9*X(10)*X(24); + Real x277 = 1.01e-9*X(10)*X(24); - Real x249 = x149*x150 + x150*x152 + x248; + Real x278 = x277 + x42*x43 + x43*x45; - Real x250 = x157*x162 + x158*x160 - x161*x162; + Real x279 = x189*x194 + x190*x192 - x193*x194; - Real x251 = -x178 + x182*x183; + Real x280 = -x210 + x214*x215; - Real x252 = -x173*x174 + x175*x176; + Real x281 = -x205*x206 + x207*x208; - Real x253 = -x130*x131 + x132*x133 + x133*x134; + Real x282 = -x168*x169 + x170*x171 + x171*x172; - Real x254 = ((x148) ? ( + Real x283 = ((x29) ? ( 1.371e-8 ) : ( 0 )) + 1.0e-99; - Real x255 = X(16)*X(25); + Real x284 = X(16)*X(25); - Real x256 = ((x151) ? ( - 4.1551773299893819e-8*x78 + 5.673e-10 + Real x285 = ((x32) ? ( + 4.1551773299893819e-8*x44 + 5.673e-10 ) : ( 0 )) + 1.0e-99; - Real x257 = -x129 + x254*x255 + x255*x256; + Real x286 = -x14 + x283*x284 + x284*x285; - Real x258 = ((x1 <= 1167.4796423742259) ? ( - std::pow(10, 5.8888600000000002*x139 + 7.1969200000000004*x142 + 2.2506900000000001*x185 - 56.473700000000001 - 2.1690299999999998*x180/((((x137)*(x137)))*(((x137)*(x137)))) + 0.31788699999999998*x181/((x137)*(x137)*(x137)*(x137)*(x137))) + Real x287 = ((x0 <= 1167.4796423742259) ? ( + std::pow(10, 5.8888600000000002*x177 + 7.1969200000000004*x180 + 2.2506900000000001*x217 - 56.473700000000001 - 2.1690299999999998*x212/((((x175)*(x175)))*(((x175)*(x175)))) + 0.31788699999999998*x213/((x175)*(x175)*(x175)*(x175)*(x175))) ) : ( 0 )) + 1.0e-99; - Real x259 = X(10)*X(7); + Real x288 = X(10)*X(7); - Real x260 = ((x1 > 1167.4796423742259) ? ( - 3.1699999999999999e-10*std::exp(-5207.0*x52) + Real x289 = ((x0 > 1167.4796423742259) ? ( + 3.1699999999999999e-10*std::exp(-5207.0*x10) ) : ( 0 )) + 1.0e-99; - Real x261 = ((x1 > 200.0) ? ( - 5.25e-11*std::exp(-4430.0*x52 + 173900.0/x2) + Real x290 = ((x0 > 200.0) ? ( + 5.25e-11*std::exp(-4430.0*x10 + 173900.0/x1) ) : ( 0 )) + 1.0e-99; - Real x262 = X(12)*X(4); - - Real x263 = x258*x259 + x259*x260 - x261*x262; - - Real x264 = 3.5999999999999998e-16*std::pow(T, 1.52)*std::exp(-1740.0*x52) + 1.0e-99; - - Real x265 = X(10)*X(25); - - Real x266 = 6.3600000000000002e-14*x1*std::exp(-9610.0*x52) + 1.0e-99; + Real x291 = X(12)*X(4); - Real x267 = X(27)*X(4); + Real x292 = x287*x288 + x288*x289 - x290*x291; - Real x268 = x264*x265 - x266*x267 + x80*x81; - - Real x269 = ((x1 <= 150.0) ? ( + Real x293 = ((x0 <= 150.0) ? ( 4.9970000000000003e-11 ) : ( 0 )) + 1.0e-99; - Real x270 = X(22)*X(25); + Real x294 = X(22)*X(25); - Real x271 = ((x1 > 150.0) ? ( - 2.4000000000000001e-11*std::exp(110.0*x52) + Real x295 = ((x0 > 150.0) ? ( + 2.4000000000000001e-11*std::exp(110.0*x10) ) : ( 0 )) + 1.0e-99; - Real x272 = 1.0e-99 + 3.1319999999999996e-10*std::exp(-8156.0*x52); + Real x296 = 1.0e-99 + 3.1319999999999996e-10*std::exp(-8156.0*x10); - Real x273 = X(32)*X(4); + Real x297 = X(32)*X(4); - Real x274 = x269*x270 + x270*x271 - x272*x273; + Real x298 = x293*x294 + x294*x295 - x296*x297; - Real x275 = x118*x66 - x119*x73 - x72*x73 - x73*x76; + Real x299 = x109*x157 - x115*x116 - x116*x119 - x116*x158; - Real x276 = -x164*x165 + x166*x167; + Real x300 = -x196*x197 + x198*x199; - Real x277 = ((x148) ? ( - 7.0509999999999997e-11 -) -: ( - 0 -)) + 1.0e-99; + Real x301 = 1.0e-25*X(4)*X(7); - Real x278 = X(17)*X(25); + Real x302 = x163*x164 - x301 + x62*x63; - Real x279 = ((x151) ? ( - 1.5557027645121538e-10*std::pow(x1, -0.33900000000000002)*std::exp(-0.108*x52) -) -: ( - 0 -)) + 1.0e-99; - - Real x280 = 1.0e-99 + 1.02e-10*std::exp(-914.0*x52); - - Real x281 = X(18)*X(22); - - Real x282 = x277*x278 + x278*x279 + x280*x281; + Real x303 = 1.0e-99 + 1.46e-12*std::exp(-9650.0*x10); - Real x283 = 1.0e-25*X(4)*X(7); + Real x304 = X(10)*X(22); - Real x284 = x124*x125 + x17*x18 - x283; - - Real x285 = 1.0e-99 + 1.46e-12*std::exp(-9650.0*x52); - - Real x286 = X(10)*X(22); - - Real x287 = ((x1 > 280.0) ? ( + Real x305 = ((x0 > 280.0) ? ( 5.4499999999999998e-17 ) : ( 0 )) + 1.0e-99; - Real x288 = X(25)*X(4); + Real x306 = X(25)*X(4); - Real x289 = ((x1 <= 280.0) ? ( - 8.1009819130659227e-21*std::pow(x1, 2.7999999999999998)*std::exp(-1950.0*x52) + Real x307 = ((x0 <= 280.0) ? ( + 8.1009819130659227e-21*std::pow(x0, 2.7999999999999998)*std::exp(-1950.0*x10) ) : ( 0 )) + 1.0e-99; - Real x290 = x285*x286 - x287*x288 - x288*x289 + x83*x84; - - Real x291 = 8.648509591313017e-18*std::pow(x1, -0.38) + 1.0e-99; + Real x308 = x122*x123 + x303*x304 - x305*x306 - x306*x307; - Real x292 = X(22)*X(4); + Real x309 = 8.648509591313017e-18*std::pow(x0, -0.38) + 1.0e-99; - Real x293 = x20*x21 - x291*x292; + Real x310 = X(22)*X(4); - Real x294 = 4.9999999999999996e-6*x10 + 1.0e-99; + Real x311 = -x309*x310 + x65*x66; - Real x295 = X(5)*X(8); + Real x312 = 4.9999999999999996e-6*x55 + 1.0e-99; - Real x296 = x202*(6.0000000000000001e-32*x51 + 2.0000000000000002e-31*x78 + 1.0e-99) + x215*(-x209 - x214 - 3.0000000000000001e-99) + x294*x295; + Real x313 = X(5)*X(8); - Real x297 = -x294*x295; + Real x314 = x234*(2.0000000000000002e-31*x44 + 6.0000000000000001e-32*x96 + 1.0e-99) + x247*(-x241 - x246 - 3.0000000000000001e-99) + x312*x313; - Real x298 = x106*x239 + x107*x240 + x107*x241; + Real x315 = -x312*x313; - Real x299 = X(6)*X(9); + Real x316 = x145*x268 + x146*x269 + x146*x270; - Real x300 = x299*(-9.8726896031426014e-7*x78 - 1.0e-99); + Real x317 = X(6)*X(9); - Real x301 = -x143*x144 + x145*x146; + Real x318 = x317*(-9.8726896031426014e-7*x44 - 1.0e-99); - Real x302 = -x258*x259 - x259*x260 + x261*x262; + Real x319 = -x181*x182 + x183*x184; - Real x303 = -x248; + Real x320 = -x287*x288 - x288*x289 + x290*x291; - Real x304 = 1.6974097914175e-15*x77*std::exp(-7550.0*x52) + 1.0e-99; + Real x321 = -x277; - Real x305 = X(10)*X(33); + Real x322 = -x11*x12; - Real x306 = -x304*x305; + Real x323 = -x223; - Real x307 = -x191; + Real x324 = -x265; - Real x308 = -x233; + Real x325 = -x251*x252 + x253*x254; - Real x309 = -x219*x220 + x221*x222; + Real x326 = -x256*x257 + x258*x259; - Real x310 = -x224*x225 + x226*x227; + Real x327 = x261 - x262*x263; - Real x311 = x229 - x230*x231; + Real x328 = X(10)*x225; - Real x312 = X(10)*x193; + Real x329 = -x272*x273 + x274*x275 - x328; - Real x313 = -x243*x244 + x245*x246 - x312; + Real x330 = -x47*x48 + x51*x52; - Real x314 = -x264*x265 + x266*x267; + Real x331 = x123*x124 - x303*x304 + x305*x306 + x306*x307; - Real x315 = -x285*x286 + x287*x288 + x288*x289 + x84*x85; + Real x332 = X(17)*X(31); - Real x316 = X(17)*X(31); + Real x333 = 5.2000000000000001e-11*x332; - Real x317 = 5.2000000000000001e-11*x316; + Real x334 = x333 - 7.9500000000000005e-10*x7; - Real x318 = X(16)*X(32); + Real x335 = -x283*x284 - x284*x285; - Real x319 = x317 - 7.9500000000000005e-10*x318; - - Real x320 = -x254*x255 - x255*x256; - - Real x321 = ((x1 <= 300.0) ? ( + Real x336 = ((x0 <= 300.0) ? ( 2.5000000000000002e-18 ) : ( 0 )) + 1.0e-99; - Real x322 = X(16)*X(22); - - Real x323 = ((x1 > 300.0) ? ( - 7.3875058242066447e-18*x188*std::exp(-68.0*x52) -) -: ( - 0 -)) + 1.0e-99; - - Real x324 = -x321*x322 - x322*x323; + Real x337 = X(16)*X(22); - Real x325 = -1.04e-10*x316; - - Real x326 = x40*x41; - - Real x327 = ((x1 <= 1052.0) ? ( - 3.2682787072746432e-10*x58 + Real x338 = ((x0 > 300.0) ? ( + 7.3875058242066447e-18*x220*std::exp(-68.0*x10) ) : ( 0 )) + 1.0e-99; - Real x328 = X(17)*X(32); - - Real x329 = ((x1 > 1052.0) ? ( - 3.7991369449339675e-16*std::pow(x1, 1.54)*std::exp(613.0*x52) -) -: ( - 0 -)) + 1.0e-99; - - Real x330 = -x327*x328 - x328*x329; - - Real x331 = -x277*x278 - x278*x279; - - Real x332 = -x280*x281; - - Real x333 = 4.5299999999999999e-10*x318; - - Real x334 = ((X(25))*(X(25))); - - Real x335 = std::pow(x1, 1.1399999999999999)*std::exp(-50.0*x52); - - Real x336 = x334*(2.4749459212614312e-15*x335 + 1.0e-99); - - Real x337 = ((X(22))*(X(22))); - - Real x338 = std::pow(x1, 1.5800000000000001); - - Real x339 = 1.0e-99 + 2.8e-11*std::exp(-26500.0*x52); + Real x339 = -x336*x337 - x337*x338; - Real x340 = X(22)*X(33); + Real x340 = -1.04e-10*x332; - Real x341 = 1.0e-99 + 4.1999999999999999e-12*std::exp(-24000.0*x52); + Real x341 = x85*x86; - Real x342 = X(30)*X(32); + Real x342 = -x25*x26 - x26*x27; - Real x343 = -x339*x340 + x341*x342; + Real x343 = -x30*x31 - x31*x33; - Real x344 = x317 + 3.4200000000000001e-10*x318; + Real x344 = -x34*x35; - Real x345 = x327*x328 + x328*x329; + Real x345 = ((X(22))*(X(22))); - Real x346 = -x269*x270 - x270*x271 + x272*x273; + Real x346 = std::pow(x0, 1.5800000000000001); - Real x347 = -x235 + x236*x237; + Real x347 = -x15*x16 + x17*x18; - Real x348 = x304*x305; + Real x348 = x333 + 3.4200000000000001e-10*x7; - Real x349 = x339*x340 - x341*x342; + Real x349 = -x293*x294 - x294*x295 + x296*x297; - ydot(1) = X(30)*x7 - x0; + ydot(1) = X(30)*x6 + x37; - ydot(2) = X(27)*x9 - x8; + ydot(2) = X(27)*x38 + x54; - ydot(3) = x104 + x108 + x113 + x117 + x120 + x123 + x13 + x16 + x19 + x22 + x34*(8.7048158229535104e-117*x32 + 2.0e-99) + x35 + x39 + x42 + x45*(7.5807998548789035e-18*x43 + 2.0e-99) + x46 + x49 + x56*(2.6047246489532126e-6*x54 + 2.0e-99) + x57 + x62*(4.2024322823973014e-6*x60 + 2.0e-99) + x63 + x65*x66 + x68*(7.7143746719362896e-209*x67 + 2.0e-99) + x70*(1.1816487727453014e-69*x69 + 2.0e-99) + x72*x73 + x73*x76 + x82 + x86 + x92 + x97; + ydot(3) = x100*(2.6047246489532126e-6*x98 + 2.0e-99) + x101 + x105*(4.2024322823973014e-6*x103 + 2.0e-99) + x106 + x108*x109 + x111*(7.7143746719362896e-209*x110 + 2.0e-99) + x113*(1.1816487727453014e-69*x112 + 2.0e-99) + x115*x116 + x116*x119 + x121 + x125 + x131 + x136 + x143 + x147 + x152 + x156 + x159 + x162 + x58 + x61 + x64 + x67 + x79*(8.7048158229535104e-117*x77 + 2.0e-99) + x80 + x84 + x87 + x90*(7.5807998548789035e-18*x88 + 2.0e-99) + x91 + x94; - ydot(4) = x117 + x121*x70 + x126 + x128 + x129 + x135 + x147 + x153 + x155 + x163 + x168 + x172 + x177 + x184 + x190; + ydot(4) = x113*x160 + x14 + x156 + x165 + x167 + x173 + x185 + x187 + x195 + x200 + x204 + x209 + x216 + x222 + x46; - ydot(5) = x114*x115 + x115*x116 + x123 + x127 + x190 + x191 + x192 + x195 + x196*x201 + x202*(-1.8e-31*x51 - 6.0000000000000005e-31*x78 - 3.0000000000000001e-99) + x215*(3*x209 + 3*x214 + 9.0000000000000006e-99) - x216*x217 - 5.9999999999999997e-18*x218 + x223 + x228 + x232 + x234 + x238 + x242 + x247 + x249 + x250 + x251 + x252 + x253 + x257 + x263 + x268 + x274 + x275 + x276 + x282 + x284 + x290 + x293 + x296 + x44*x45 + x66*(8.7599999999999997e-10*x64 + 2.0e-99) + x73*(2*x71 + 2.0e-99) + x73*(2*x75 + 2.0e-99) + x80*(8.9373821670554065e-6*x78 + 2.0e-99) + x89*(2*x88 + 2.0e-99) + x91*(2*x154 + 2.0e-99); + ydot(5) = x109*(8.7599999999999997e-10*x107 + 2.0e-99) + x116*(2*x114 + 2.0e-99) + x116*(2*x118 + 2.0e-99) + x128*(2*x127 + 2.0e-99) + x130*(2*x186 + 2.0e-99) + x153*x154 + x154*x155 + x162 + x166 + x222 + x223 + x224 + x227 + x228*x233 + x234*(-6.0000000000000005e-31*x44 - 1.8e-31*x96 - 3.0000000000000001e-99) + x247*(3*x241 + 3*x246 + 9.0000000000000006e-99) - x248*x249 - 5.9999999999999997e-18*x250 + x255 + x260 + x264 + x266 + x267 + x271 + x276 + x278 + x279 + x280 + x281 + x282 + x286 + x292 + x298 + x299 + x300 + x302 + x308 + x311 + x314 + x36 + x50*(8.9373821670554065e-6*x44 + 2.0e-99) + x53 + x89*x90; - ydot(6) = x122*x70 + x155 + x275 + x297 + x298 + x46; + ydot(6) = x113*x161 + x187 + x299 + x315 + x316 + x91; - ydot(7) = x250 + x300 + x301 + x39; + ydot(7) = x279 + x318 + x319 + x84; - ydot(8) = x163 + x284 + x298 + x299*(1.9745379206285203e-6*x78 + 2.0e-99) + x302 + x37*x38 + x49; + ydot(8) = x195 + x302 + x316 + x317*(1.9745379206285203e-6*x44 + 2.0e-99) + x320 + x82*x83 + x94; - ydot(9) = x186*x187 + x187*x189 + x251 + x297 + x92; + ydot(9) = x131 + x218*x219 + x219*x221 + x280 + x315; - ydot(10) = x126 + x242 + x300 + x47*x48; + ydot(10) = x165 + x271 + x318 + x92*x93; - ydot(11) = x120 + x184 - x196*x201 + x196*(x200 + 1.0e-99) + x216*x217 + x217*(-7.5000000000000001e-33*x51 - 2.5000000000000002e-32*x78 - 1.0e-99) + 2.9999999999999998e-18*x218 + x296 + x301 + x302 + x303 + x306 + x307 + x308 + x309 + x310 + x311 + x313 + x314 + x315; + ydot(11) = x159 + x216 - x228*x233 + x228*(x232 + 1.0e-99) + x248*x249 + x249*(-2.5000000000000002e-32*x44 - 7.5000000000000001e-33*x96 - 1.0e-99) + 2.9999999999999998e-18*x250 + x314 + x319 + x320 + x321 + x322 + x323 + x324 + x325 + x326 + x327 + x329 + x330 + x331; - ydot(12) = x156*x158 + x156*x162 + x19 + x192; + ydot(12) = x188*x190 + x188*x194 + x224 + x64; - ydot(13) = x108 + x147 + x263 + x283; + ydot(13) = x147 + x185 + x292 + x301; - ydot(14) = x103*x68 + x13; + ydot(14) = x111*x142 + x58; - ydot(15) = x104 + x11*x12 + x253 + x33*x34; + ydot(15) = x143 + x282 + x56*x57 + x78*x79; - ydot(16) = x100*x68 + x102*x68 + x135 + x35; + ydot(16) = x111*x139 + x111*x141 + x173 + x80; - ydot(17) = x113 + x276 + x319 + x320 + x324 + x55*x56; + ydot(17) = x100*x99 + x152 + x300 + x334 + x335 + x339; - ydot(18) = x109*x110 + x110*x111 + x110*x112 + x168 + x195 + x313 + x325 + x326 + x330 + x331 + x57; + ydot(18) = x101 + x148*x149 + x149*x150 + x149*x151 + x200 + x227 + x329 + x340 + x341 + x342 + x343; - ydot(19) = x194 + x247 + x311 + x332; + ydot(19) = x226 + x276 + x327 + x344; - ydot(20) = x232 + x310 + x312; + ydot(20) = x264 + x326 + x328; - ydot(21) = x228 + x309; + ydot(21) = x260 + x325; - ydot(22) = x252 + x308 + x333 + x61*x62 + x97; + ydot(22) = x104*x105 + x136 + x281 + x324 + x8; - ydot(23) = x15*(2.1137298733407573e-5*x14 + 2.0e-99) + x177 + x293 + x315 + x324 + x326 + x332 + x336 + x337*(-1.1950181274729681e-23*x338 - 2.0e-99) + x343 + x344 + x345 + x346 + x63 + x93*x94 + x94*x96; + ydot(23) = x106 + x132*x133 + x133*x135 + x209 + x28 + x311 + x331 + x339 + x341 + x344 + x345*(-1.1950181274729681e-23*x346 - 2.0e-99) + x347 + x348 + x349 + x41 + x60*(2.1137298733407573e-5*x59 + 2.0e-99); - ydot(24) = x223; + ydot(24) = x255; - ydot(25) = x22 + x234 + x303; + ydot(25) = x266 + x321 + x67; - ydot(26) = x172 + x290 + x291*x292 + x314 + x320 + x331 + x334*(-4.9498918425228625e-15*x335 - 2.0e-99) + x346 + x347 + x79*x80; + ydot(26) = x120*x50 + x204 + x23 + x308 + x309*x310 + x330 + x335 + x343 + x349 + x39*(-4.9498918425228625e-15*x40 - 2.0e-99); - ydot(27) = x249 + x307 + x86; + ydot(27) = x125 + x278 + x323; - ydot(28) = -X(27)*x9 + x153 + x268 + x336 + x348 + x8; + ydot(28) = 2.0e-99*X(1) + x54; - ydot(29) = x191 + x82; + ydot(29) = x121 + x223; - ydot(30) = x257 + x321*x322 + x322*x323 + x344 + x42; + ydot(30) = x286 + x336*x337 + x337*x338 + x348 + x87; - ydot(31) = -X(30)*x7 + x0 + x129 + x282 + x333 + x345 + x347 + x348 + x349; + ydot(31) = 2.0e-99*X(0) + x37; - ydot(32) = x127 + x16 + x325; + ydot(32) = x166 + x340 + x61; - ydot(33) = x128 + x274 + x319 + x330 + x337*(5.9750906373648406e-24*x338 + 1.0e-99) + x349; + ydot(33) = x167 + x19 + x298 + x334 + x342 + x345*(5.9750906373648406e-24*x346 + 1.0e-99); - ydot(34) = x238 + x306 + x343; + ydot(34) = x267 + x322 + x347; @@ -1103,7 +1103,7 @@ Real rhs_eint(const burn_t& state, Real x24 = X(11) + X(12) + X(18) + X(24) + X(25) + X(5); - Real x25 = 2.0*X(10) + 2.0*X(19) + 3.0*X(20) + 4.0*X(23) + 2.0*X(26) + 2.0*X(27) + 3.0*X(28) + X(3) + X(4) + 2.0*X(8) + x24; + Real x25 = 2.0*X(1) + 2.0*X(10) + 2.0*X(19) + 3.0*X(20) + 4.0*X(23) + 2.0*X(26) + 3.0*X(28) + X(3) + X(4) + 2.0*X(8) + x24; Real x26 = 1.0/x25; @@ -3962,7 +3962,7 @@ Real rhs_eint(const burn_t& state, ) : ( -9.9999999999999997e+98 - ))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))*Z)*((0.5*X(10) + 0.25*X(11) + 0.25*X(12) + 0.25*X(18) + 0.5*X(19) + 0.75*X(20) + X(23) + 0.25*X(24) + 0.25*X(25) + 0.5*X(26) + 0.5*X(27) + 0.75*X(28) + 0.25*X(3) + 0.25*X(4) + 0.25*X(5) + 0.5*X(8))*(0.5*X(10) + 0.25*X(11) + 0.25*X(12) + 0.25*X(18) + 0.5*X(19) + 0.75*X(20) + X(23) + 0.25*X(24) + 0.25*X(25) + 0.5*X(26) + 0.5*X(27) + 0.75*X(28) + 0.25*X(3) + 0.25*X(4) + 0.25*X(5) + 0.5*X(8))) + (1 - x176)*(((x451 || x453 || x455 < 9.0) ? ( + ))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))*Z)*((0.5*X(1) + 0.5*X(10) + 0.25*X(11) + 0.25*X(12) + 0.25*X(18) + 0.5*X(19) + 0.75*X(20) + X(23) + 0.25*X(24) + 0.25*X(25) + 0.5*X(26) + 0.75*X(28) + 0.25*X(3) + 0.25*X(4) + 0.25*X(5) + 0.5*X(8))*(0.5*X(1) + 0.5*X(10) + 0.25*X(11) + 0.25*X(12) + 0.25*X(18) + 0.5*X(19) + 0.75*X(20) + X(23) + 0.25*X(24) + 0.25*X(25) + 0.5*X(26) + 0.75*X(28) + 0.25*X(3) + 0.25*X(4) + 0.25*X(5) + 0.5*X(8))) + (1 - x176)*(((x451 || x453 || x455 < 9.0) ? ( 0.0 ) : ( @@ -4117,8 +4117,8 @@ std::pair compute_Semenov_Tdust(Real Tgas, const Array1D