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<Real, Real> compute_Semenov_Tdust(Real Tgas, const Array1D<Real, 0, Nu
     Real nH = 0.0;
     nH = composition(3) + composition(4) + composition(5) + composition(10) * 2.0 + composition(8) * 2.0 +
          composition(25) + composition(18) + composition(19) * 2.0 + composition(20) * 3.0 +
-         composition(27) * 2.0 + composition(23) * 4.0 + composition(24) + composition(26) * 2.0 +
-         composition(28) * 3.0 + composition(11) + composition(12);
+         composition(23) * 4.0 + composition(24) + composition(26) * 2.0 +
+         composition(28) * 3.0 + composition(11) + composition(12) + composition(1) * 2.0;
 
     //std::cout << "compute semenov tdust nH: " << nH << std::endl;
 
diff --git a/unit_test/burn_cell_metal_chem/reference_solution_1.out b/unit_test/burn_cell_metal_chem/reference_solution_1.out
index f81ada569..27a1ffb23 100644
--- a/unit_test/burn_cell_metal_chem/reference_solution_1.out
+++ b/unit_test/burn_cell_metal_chem/reference_solution_1.out
@@ -1,45 +1,45 @@
-T final =   2033.7719
+T final =   2033.663091
 Tdust initial = 2.73
-Tdust final =   2033.610864
+Tdust final =   2033.629887
 Eint initial = 3.061675486e+10
-Eint final =   1.764013907e+11
+Eint final =   1.763898902e+11
 rho initial = 2.203783613e-25
-rho final =   1.969175663e-06
-numdens final =   5.200225029e+17
+rho final =   1.969712207e-06
+numdens final =   5.200218306e+17
 ------------------------------------
 New number densities: 
 8.329619835e+13
-2.373088458e+14
-2.040789842
-6.529533105e-08
-1.944958889e+15
-2.554410631e-05
-6.055749336e-13
-9860014349
-1.817920945e-10
-1.199955649e-10
-4.481098213e+17
-2.444960273e-20
-8.975750696e+12
-2.630091902e-09
-8.82330545e-21
+2.373088567e+14
+2.040909947
+1.074117805e-06
+1.943614057e+15
+2.553528356e-05
+9.961696151e-12
+9853431998
+2.99089071e-09
+1.199571059e-10
+4.481104939e+17
+4.021635872e-19
+8.975754713e+12
+1.746215894e-09
+5.862794549e-21
 6.963813204e+16
-9.871889761e-21
-6.739096013e-14
-2.597905176e-13
-5.714209204e-11
-1.703374921e-08
-2.53711863e-40
-3.383022442e-15
-1.998329514e-07
-3.16125345e-35
-5.387126028e-11
-1.825600604e-31
-1.766483831e-07
-2.040750381
-1.559962404e-34
-8.563324351e-29
-6.493533745e-05
-2.06061654e-32
-2.973172266e-55
+1.092755289e-19
+4.6280631e-13
+1.784798056e-12
+3.928263544e-10
+1.171832259e-07
+5.089162611e-39
+4.125246338e-15
+1.375692779e-06
+5.196274139e-34
+5.382850015e-11
+3.0339037e-30
+1.766616999e-07
+2.040849206
+1.406529162e-33
+7.715200326e-28
+8.519569305e-05
+1.858338724e-31
+2.678201778e-54
 AMReX (22.07-589-g9b0d63e40412) finalized
diff --git a/unit_test/burn_cell_metal_chem/reference_solution_1_z10.out b/unit_test/burn_cell_metal_chem/reference_solution_1_z10.out
index e485c3b44..9158256e4 100644
--- a/unit_test/burn_cell_metal_chem/reference_solution_1_z10.out
+++ b/unit_test/burn_cell_metal_chem/reference_solution_1_z10.out
@@ -1,47 +1,47 @@
-Tgas = 31.53250536
-Tgas = 31.52704064
-T final =   2033.846636
+Tgas = 31.53251506
+Tgas = 31.5270498
+T final =   2033.744028
 Tdust initial = 30.03
-Tdust final =   2033.692617
+Tdust final =   2033.708884
 Eint initial = 3.061675486e+10
-Eint final =   1.764079063e+11
+Eint final =   1.763969463e+11
 rho initial = 2.203783613e-25
-rho final =   1.969175663e-06
-numdens final =   5.20022965e+17
+rho final =   1.969712208e-06
+numdens final =   5.200223311e+17
 ------------------------------------
-New number densities: 
+New number densities:
 8.329620029e+13
 2.373088525e+14
-2.226446088
-0.000297491388
-1.945883073e+15
-2.787567024e-05
-2.759062723e-09
-9864501378
-8.281880969e-07
-1.309455569e-10
-4.481093592e+17
-1.114010181e-16
-8.97571489e+12
-1.601389437e-09
-5.858032102e-21
+2.22753933
+0.001917923831
+1.944614328e+15
+2.787872492e-05
+1.778752712e-08
+9858318342
+5.339966185e-06
+1.309633359e-10
+4.481099941e+17
+7.18142559e-16
+8.975743083e+12
+6.009135232e-10
+2.200891197e-21
 6.963813204e+16
-1.415072594e-20
-1.996517884e-16
-7.694480337e-16
-1.691689432e-13
-5.040356928e-11
-5.482635383e-37
-1.604555857e-15
-5.910343021e-10
-1.715097579e-31
-6.415084542e-11
-1.001202566e-27
-2.10230157e-07
-2.226165405
-7.279438235e-38
-3.99809758e-32
-1.023291941e-05
-9.614519081e-36
-1.652332248e-58
+2.032448461e-18
+4.174154493e-16
+1.609285669e-15
+3.540277018e-13
+1.055530723e-10
+2.763957945e-36
+1.254725284e-15
+1.238523826e-09
+1.105069858e-30
+6.41122004e-11
+6.456231456e-27
+2.102757459e-07
+2.225642417
+2.245978095e-38
+1.232689556e-32
+1.509127756e-06
+2.967131355e-36
+5.094305713e-59
 AMReX (22.07-875-gb05f69fe3748) finalized
diff --git a/unit_test/burn_cell_metal_chem/reference_solution_1e-1.out b/unit_test/burn_cell_metal_chem/reference_solution_1e-1.out
index 5681bb98d..3a930088c 100644
--- a/unit_test/burn_cell_metal_chem/reference_solution_1e-1.out
+++ b/unit_test/burn_cell_metal_chem/reference_solution_1e-1.out
@@ -1,45 +1,45 @@
-T final =   2580.849045
+T final =   2580.84748
 Tdust initial = 2.73
-Tdust final =   2580.903736
+Tdust final =   2580.896055
 Eint initial = 3.074857008e+10
-Eint final =   2.252782311e+11
+Eint final =   2.25277789e+11
 rho initial = 2.193512909e-25
-rho final =   1.989589323e-06
-numdens final =   5.416474951e+17
+rho final =   1.98964343e-06
+numdens final =   5.416473705e+17
 ------------------------------------
 New number densities: 
 8.412917854e+12
-2.39681895e+13
-2.207979967
-6.56484998e-06
-3.496528463e+16
-0.0001134569489
-5.932464766e-11
-1.640674045e+11
-9.26261325e-09
-4.948497742e-10
-4.363062397e+17
-3.520507679e-18
-8.912257632e+12
-4.109464144e-09
-6.788802573e-22
+2.396818951e+13
+2.208070971
+4.917822618e-05
+3.496505964e+16
+0.0001134618309
+4.44409566e-10
+1.640663651e+11
+6.938765994e-08
+4.948711095e-10
+4.363063401e+17
+2.637260563e-17
+8.912258695e+12
+9.51889453e-09
+1.613163821e-21
 7.033451336e+16
-1.460565197e-19
-4.226130576e-11
-2.962304509e-11
-4.321097966e-10
-6.410593413e-09
-2.847622487e-35
-3.523568564e-12
-4.116188787e-09
-1.579040656e-30
-2.800477243e-08
-2.517542637e-28
-2.481297579e-06
-2.208086624
-1.820442242e-35
-2.493711987e-27
-2.172346472e-07
-1.03303842e-33
-1.308968162e-52
+4.313617275e-21
+5.114046531e-11
+3.584704032e-11
+5.229021966e-10
+7.757603344e-09
+2.133310512e-34
+3.523760341e-12
+4.981125586e-09
+1.182951876e-29
+2.800643178e-08
+1.886056751e-27
+2.481465418e-06
+2.208135
+1.690486453e-35
+2.334074296e-27
+1.667014038e-07
+9.596752642e-34
+1.22525698e-52
 AMReX (22.07-589-g9b0d63e40412) finalized
diff --git a/unit_test/burn_cell_metal_chem/reference_solution_1e-2.out b/unit_test/burn_cell_metal_chem/reference_solution_1e-2.out
index e6a2063ee..27b0abf28 100644
--- a/unit_test/burn_cell_metal_chem/reference_solution_1e-2.out
+++ b/unit_test/burn_cell_metal_chem/reference_solution_1e-2.out
@@ -1,45 +1,45 @@
-T final =   2811.207219
+T final =   2810.374696
 Tdust initial = 2.73
-Tdust final =   2811.228114
+Tdust final =   2811.399316
 Eint initial = 3.076181952e+10
-Eint final =   2.478723199e+11
+Eint final =   2.477868218e+11
 rho initial = 2.192485838e-25
-rho final =   1.989661731e-06
-numdens final =   5.670189018e+17
+rho final =   1.989668195e-06
+numdens final =   5.668962568e+17
 ------------------------------------
 New number densities: 
-8.412914755e+11
+8.412914754e+11
 2.396807474e+12
-2.236235976
-6.405982886e-05
-8.572287708e+16
-0.0001611585155
-5.332298561e-10
-4.037425745e+11
-6.855532471e-08
-7.073005674e-10
-4.109491978e+17
-3.761054318e-17
-8.671678263e+12
-1.456049659e-10
-1.650308794e-21
+2.237388464
+6.169419758e-07
+8.547697113e+16
+0.0001610647153
+5.137205748e-12
+4.02556912e+11
+6.608419452e-10
+7.068343694e-10
+4.110724588e+17
+3.621113956e-19
+8.672863935e+12
+6.173898911e-08
+3.888219294e-21
 7.033451336e+16
-2.473602851e-09
-0.000502828669
-0.000195826651
-0.001167253794
-0.006497881953
-1.865521147e-32
-2.352999816e-10
-0.001602637809
-5.537533357e-28
-9.719529679e-07
-2.822461066e-26
-2.69412668e-05
-2.131195609
-4.275851205e-23
-2.534923404e-13
-0.1051373945
-1.86318788e-21
-8.869168988e-21
+2.542238142e-09
+0.0004989372543
+0.0001946900676
+0.001163928336
+0.00650045374
+1.786797422e-34
+2.340063332e-10
+0.001608372808
+5.313130942e-30
+9.68531334e-07
+2.718197375e-28
+2.694990265e-05
+2.129003659
+4.392878947e-23
+-4.560006446e-12
+0.1085451266
+2.537886844e-21
+-4.677302387e-19
 AMReX (22.07-589-g9b0d63e40412) finalized
diff --git a/unit_test/burn_cell_metal_chem/reference_solution_1e-3.out b/unit_test/burn_cell_metal_chem/reference_solution_1e-3.out
index da53a37ac..5387e4fc2 100644
--- a/unit_test/burn_cell_metal_chem/reference_solution_1e-3.out
+++ b/unit_test/burn_cell_metal_chem/reference_solution_1e-3.out
@@ -1,45 +1,45 @@
-T final =   2831.090224
+T final =   2831.019478
 Tdust initial = 2.73
-Tdust final =   2831.289176
+Tdust final =   2831.686588
 Eint initial = 3.076314515e+10
-Eint final =   2.499224717e+11
+Eint final =   2.499151417e+11
 rho initial = 2.192383131e-25
-rho final =   1.989668738e-06
-numdens final =   5.700271327e+17
+rho final =   1.98966925e-06
+numdens final =   5.700162122e+17
 ------------------------------------
 New number densities: 
-8.412915268e+10
+8.412915269e+10
 2.396823174e+11
-4.033710506
-6.777594776e-05
-9.174097417e+16
-0.0002983374384
-5.593606493e-10
-4.328369028e+11
-7.097532349e-08
-1.311908161e-09
-4.079422459e+17
-4.006694765e-17
-8.64258412e+12
-2.245209778e-14
-2.215435849e-21
+4.033631317
+7.268898769e-05
+9.171911853e+16
+0.0002983039113
+5.999266876e-10
+4.327309304e+11
+7.612608637e-08
+1.311751089e-09
+4.07953181e+17
+4.297040753e-17
+8.642689996e+12
+2.354952341e-10
+2.398645645e-21
 7.033451336e+16
-1.304995286e-08
-8.203048039e-05
-3.050552029e-05
-0.0001694822058
-0.0008735603933
-5.43479155e-31
-6.483863045e-09
-0.0001998289953
-1.534470988e-26
-2.532177639e-05
-7.160676166e-25
-0.0006407136702
-0.4462317747
-2.224187897e-22
-8.810161969e-12
-3.587709209
-9.496345982e-21
-1.037846593e-34
+1.310908423e-08
+8.203187782e-05
+3.051098989e-05
+0.0001695545589
+0.0008741705667
+5.825557904e-31
+6.480280989e-09
+0.0002000217592
+1.645050698e-26
+2.53122014e-05
+7.67908607e-25
+0.0006406765401
+0.4462423571
+2.229330167e-22
+-6.520895357e-12
+3.587614486
+-2.492556219e-18
+-4.418583788e-19
 AMReX (22.07-589-g9b0d63e40412) finalized
diff --git a/unit_test/burn_cell_metal_chem/reference_solution_1e-4.out b/unit_test/burn_cell_metal_chem/reference_solution_1e-4.out
index 02948a346..9b01e0cfc 100644
--- a/unit_test/burn_cell_metal_chem/reference_solution_1e-4.out
+++ b/unit_test/burn_cell_metal_chem/reference_solution_1e-4.out
@@ -1,45 +1,45 @@
-T final =   2884.322233
+T final =   2884.408445
 Tdust initial = 2.73
-Tdust final =   2885.222529
+Tdust final =   2885.384517
 Eint initial = 3.076327772e+10
-Eint final =   2.555031196e+11
+Eint final =   2.555122564e+11
 rho initial = 2.19237286e-25
-rho final =   1.989670334e-06
-numdens final =   5.787722231e+17
+rho final =   1.98966997e-06
+numdens final =   5.787870732e+17
 ------------------------------------
 New number densities: 
-8412915325
+8412915330
 2.396818542e+10
-5.962439927
-0.0005137258492
-1.092306741e+17
-0.0004716327583
-4.142389741e-09
-5.183837544e+11
-5.083763652e-07
-2.087320476e-09
-3.991979278e+17
-3.093116661e-16
-8.557037312e+12
-4.381451187e-08
-3.7384033e-20
+5.962823224
+0.001071375051
+1.092606255e+17
+0.000471712832
+8.638627485e-09
+5.185317229e+11
+1.060123122e-06
+2.087700204e-09
+3.991828266e+17
+6.450887085e-16
+8.556889344e+12
+2.412076997e-07
+9.717358293e-20
 7.033451336e+16
-1.431469116e-07
-3.125096897e-05
-1.030321161e-05
-4.764610267e-05
-0.000200901188
-1.246503729e-28
-1.967791718e-07
-3.775902673e-05
-3.167115752e-24
-0.000678373265
-1.175791827e-22
-0.01353231898
-0.03355686142
-1.176079223e-22
-5.66867762e-11
-5.928840231
-4.761295975e-21
-1.381019684e-32
+1.420660126e-07
+3.101380105e-05
+1.022306037e-05
+4.726165595e-05
+0.0001992165272
+2.601230644e-28
+1.969050239e-07
+3.743069657e-05
+6.608038109e-24
+0.0006786635515
+2.452346741e-22
+0.01353299323
+0.03354856481
+1.166801e-22
+5.625847973e-11
+5.928673305
+4.620188648e-21
+1.370612271e-32
 AMReX (22.07-589-g9b0d63e40412) finalized
diff --git a/unit_test/burn_cell_metal_chem/reference_solution_1e-5.out b/unit_test/burn_cell_metal_chem/reference_solution_1e-5.out
index 69a2fca45..8304a7591 100644
--- a/unit_test/burn_cell_metal_chem/reference_solution_1e-5.out
+++ b/unit_test/burn_cell_metal_chem/reference_solution_1e-5.out
@@ -1,45 +1,45 @@
-T final =   2770.253479
+T final =   2770.260502
 Tdust initial = 2.73
-Tdust final =   2770.396907
+Tdust final =   2771.031412
 Eint initial = 3.076329098e+10
-Eint final =   2.437024468e+11
+Eint final =   2.437031597e+11
 rho initial = 2.192371833e-25
-rho final =   1.989670057e-06
-numdens final =   5.612534538e+17
+rho final =   1.989669807e-06
+numdens final =   5.61254438e+17
 ------------------------------------
 New number densities: 
-199898176.6
-3038212036
-2.259991418
-0.09252421839
-7.419371216e+16
-0.000154168177
-7.834842422e-07
-3.484514097e+11
-0.0001036533126
-6.744074981e-10
-4.167161489e+17
-5.354504797e-14
-8.726969668e+12
-3.849333801e-12
-1.421503406e-18
+199897170.8
+3038213042
+2.203629741
+0.0003491632575
+7.419566308e+16
+0.0001503269971
+2.956664728e-09
+3.484607261e+11
+3.91158392e-07
+6.576046293e-10
+4.167151823e+17
+2.020661145e-16
+8.726960346e+12
+5.643506785e-09
+5.771775307e-21
 7.033451336e+16
-0.00281312011
-21789165.94
-9352664.009
-64650614.55
-423270591.2
-1.537707812e-26
-1.341686703e-07
-122330252.1
-5.051058737e-22
-0.0006197903089
-3.102910261e-20
-0.02081899296
-2.164703811
-4.251138744e-23
-0.00765019574
-1.387355782e-34
-1.697403972e-27
-3.012881382e-24
+1.063497408e-05
+21790854.1
+9353212.954
+64652603.8
+423271205.5
+5.758672328e-29
+1.331455347e-07
+122327033.3
+1.891574487e-24
+0.0006150522586
+1.161974271e-22
+0.02065917213
+2.203419865
+1.59480483e-25
+0.007594234715
+-3.519279316e-31
+-3.283812312e-21
+2.967855969e-24
 AMReX (22.07-589-g9b0d63e40412) finalized
diff --git a/unit_test/burn_cell_metal_chem/reference_solution_1e-6.out b/unit_test/burn_cell_metal_chem/reference_solution_1e-6.out
index 372f204db..e9152d6c8 100644
--- a/unit_test/burn_cell_metal_chem/reference_solution_1e-6.out
+++ b/unit_test/burn_cell_metal_chem/reference_solution_1e-6.out
@@ -1,45 +1,45 @@
-T final =   2812.050392290817
+T final =   2813.763186
 Tdust initial = 2.73
-Tdust final =   2812.235837368595
-Eint initial = 30763292303.77594
-Eint final =   247958162899.1352
-rho initial = 2.192371730563381e-25
-rho final =   1.989669757431512e-06
-numdens final =   5.671429626664702e+17
+Tdust final =   2813.526575
+Eint initial = 3.07632923e+10
+Eint final =   2.481342021e+11
+rho initial = 2.192371731e-25
+rho final =   1.98967001e-06
+numdens final =   5.673973163e+17
 ------------------------------------
-New number densities: 
-6438308.410300142
-317322562.9538858
-2.486641265170912
-2.274499630590923e-09
-8.597266896217706e+16
-0.0001794039501637353
-1.8925886867936e-14
-404944748318.3216
-2.431869387308584e-12
-7.874340701282926e-10
-4.108267045200611e+17
-1.335778922651524e-21
-8670476376287.059
-3.501410855685503e-14
-6.130614268934865e-26
-7.033451336163859e+16
-2.268083235873295e-12
-3948863.726870825
-1534860.916182474
-9121355.520456472
-50610347.48632588
-4.088062980895315e-28
-0.1452240345710818
-12442512.37027289
-1.216206468005654e-23
-601.0948939083418
-6.175692426101071e-22
-16596.84230665624
-2.486820667629221
-2.85628902206954e-26
-32951.98449995967
-1.632493779270958e-29
-1.471009580362413e-15
-1.014011717966889e-11
+New number densities:
+6472049.128
+317338389.7
+2.473688688
+2.437815373e-09
+8.648123216e+16
+0.0001788715759
+2.026991226e-14
+4.073976485e+11
+2.601564064e-12
+7.85222469e-10
+4.10572495e+17
+1.432546932e-21
+8.668023467e+12
+5.262395419e-09
+3.603849018e-22
+7.033451336e+16
+2.490210408e-12
+4012171.391
+1553236.62
+9174440.11
+50566325.55
+4.2605529e-28
+0.1412494562
+12350866.07
+1.262582542e-23
+582.0669493
+6.362171187e-22
+15945.26769
+2.473867548
+3.01732562e-26
+32833.01421
+1.610264744e-29
+1.375039614e-15
+9.699753847e-12
 AMReX (22.07-589-g9b0d63e40412) finalized