add cr test at low Z

.github/workflows/burn_cell_metal_chem.yml

@@ -48,6 +48,7 @@ jobs:
             ["Z=1e-5"]="4438 4440 4442 4444 4445 4448 4449 4452 4458 4463"
             ["Z=1e-6"]="4463 4465 4467 4469 4470 4473 4474 4477 4483 4488"
             ["Z=1_cr"]="652 654 656 658 659 662 663 666 672 677"
+            ["Z=1e-4_cr"]="652 654 656 658 659 662 663 666 672 677"
           )
 
           declare -A ref_map=(
@@ -60,6 +61,7 @@ jobs:
             ["Z=1e-5"]="reference_solution_1e-5.out"
             ["Z=1e-6"]="reference_solution_1e-6.out"
             ["Z=1_cr"]="reference_solution_1_cr.out"
+            ["Z=1e-4_cr"]="reference_solution_1e-4_cr.out"
           )
 
           ref_line_number_z10=(1 2 3 5 7 10 17 23 28)
@@ -68,25 +70,30 @@ jobs:
           original_input_file="inputs_metal_chem_1"
           modified_input_file="inputs_metal_chem_modified"
 
-          for Z in "Z=1" "Z=1_z10" "Z=1e-1" "Z=1e-2" "Z=1e-3" "Z=1e-4" "Z=1e-5" "Z=1e-6" "Z=1_cr"; do
+          for Z in "Z=1" "Z=1_z10" "Z=1e-1" "Z=1e-2" "Z=1e-3" "Z=1e-4" "Z=1e-5" "Z=1e-6" "Z=1_cr" "Z=1e-4_cr"; do
             cp $original_input_file $modified_input_file
             Z_val=${Z//Z=/}
 
             if [[ "$Z" == "Z=1_z10" ]]; then
               # Modify the redshift line for Z=1_z10
               sed -i 's/network.redshift = 0.0/network.redshift = 10.0/g' $modified_input_file
-            elif [[ "$Z" != "Z=1_cr" ]]; then
+            elif [[ "$Z" != "Z=1_cr" && "$Z" != "Z=1e-4_cr" ]]; then
               # Replace the metallicity and dust2gas_ratio values for other Z values
               sed -i 's/network.metallicity = .*/network.metallicity = '"$Z_val"'/g' $modified_input_file
               sed -i 's/network.dust2gas_ratio = .*/network.dust2gas_ratio = '"$Z_val"'/g' $modified_input_file
             fi
 
-            if [[ "$Z" == "Z=1_cr" ]]; then
+            if [[ "$Z" == "Z=1_cr" || "$Z" == "Z=1e-4_cr" ]]; then
               # Modify the number of steps, since CRs impact lowest densities
               sed -i 's/unit_test.nsteps = 100000/unit_test.nsteps = 600/g' $modified_input_file
               # Modify the crate line for Z=1_cr
               sed -i 's/network.crate = 0.0/network.crate = 3e-17/g' $modified_input_file
-            fi
+              if [[ "$Z" == "Z=1e-4_cr" ]]; then
+                 # Replace the metallicity and dust2gas_ratio values for other Z values
+                 sed -i 's/network.metallicity = 1/network.metallicity = 1e-4/g' $modified_input_file
+                 sed -i 's/network.dust2gas_ratio = 1/network.dust2gas_ratio = 1e-4/g' $modified_input_file 
+              fi
+           fi
 
             output_file="test_${Z_val}.out"
             ./main1d.gnu.DEBUG.ex $modified_input_file amrex.fpe_trap_{invalid,zero,overflow}=1 > $output_file
@@ -106,7 +113,7 @@ jobs:
                 reference_line_number=${ref_line_number_z10[$index]}
               elif [[ "$Z" == "Z=1e-6" ]]; then
                 reference_line_number=$((line_number - 4462))
-              elif [[ "$Z" == "Z=1_cr" ]]; then
+              elif [[ "$Z" == "Z=1_cr" || "$Z" == "Z=1e-4_cr" ]]; then
                 reference_line_number=$((line_number - 651))
               else
                 reference_line_number=$((line_number - 4437))