add checks at intermediate densities for high-z test

.github/workflows/burn_cell_metal_chem.yml

@@ -80,7 +80,7 @@ jobs:
         run: |
           set -e
           cd unit_test/burn_cell_metal_chem
-          line_numbers=(4461 4463 4465 4467 4468 4471 4472 4475 4481 4486)
+          line_numbers=(1612 1614 4479 4481 4483 4485 4486 4489 4490 4493 4499 4504)
           threshold=0.01
           error_found=false
           for line_number in "${line_numbers[@]}"; do

unit_test/burn_cell_metal_chem/burn_cell.H

@@ -344,6 +344,11 @@ auto burn_cell_c() -> int {
     state_over_time << std::endl;
     std::cout << "step " << n << " done with dt = " << dt << std::endl;
 
+    // print some output at intermediate densities for testing purposes if redshift > 0
+    if (state.rho > 1e-18 && state.rho < 1.2e-18 && redshift > 0) {
+      std::cout << "Tgas = " << state.T << std::endl;
+    }
+
   }
   state_over_time.close();
 

unit_test/burn_cell_metal_chem/reference_solution_1_z10.out

@@ -1585,23 +1585,41 @@ step 1537 done with dt = 2.137708209e+10
 step 1538 done with dt = 2.127099153e+10
 step 1539 done with dt = 2.116542749e+10
 step 1540 done with dt = 2.106038734e+10
+Tgas = 31.53250536
 step 1541 done with dt = 2.09558685e+10
+Tgas = 31.52704064
 step 1542 done with dt = 2.085186837e+10
+Tgas = 31.52147312
 step 1543 done with dt = 2.074838438e+10
+Tgas = 31.51580941
 step 1544 done with dt = 2.064541397e+10
+Tgas = 31.51005587
 step 1545 done with dt = 2.054295459e+10
+Tgas = 31.50421855
 step 1546 done with dt = 2.04410037e+10
+Tgas = 31.49830335
 step 1547 done with dt = 2.033955877e+10
+Tgas = 31.49231585
 step 1548 done with dt = 2.023861731e+10
+Tgas = 31.48626145
 step 1549 done with dt = 2.01381768e+10
+Tgas = 31.48014529
 step 1550 done with dt = 2.003823477e+10
+Tgas = 31.47397231
 step 1551 done with dt = 1.993878873e+10
+Tgas = 31.46774723
 step 1552 done with dt = 1.983983623e+10
+Tgas = 31.46147455
 step 1553 done with dt = 1.974137481e+10
+Tgas = 31.4551586
 step 1554 done with dt = 1.964340205e+10
+Tgas = 31.44880348
 step 1555 done with dt = 1.954591551e+10
+Tgas = 31.44241076
 step 1556 done with dt = 1.944891278e+10
+Tgas = 31.43598915
 step 1557 done with dt = 1.935239145e+10
+Tgas = 31.42953383
 step 1558 done with dt = 1.925634915e+10
 step 1559 done with dt = 1.91607835e+10
 step 1560 done with dt = 1.906569212e+10