Add more support to the Grackle interface

doc/wiki/Runtime-Parameters-related/Runtime-Parameters:-Chemistry-and-Radiation.md

@@ -9,6 +9,13 @@ Parameters described on this page:
 [GRACKLE_PE_HEATING](#GRACKLE_PE_HEATING),  
 [GRACKLE_PE_HEATING_RATE](#GRACKLE_PE_HEATING_RATE),  
 [GRACKLE_CLOUDY_TABLE](#GRACKLE_CLOUDY_TABLE)  
+[GRACKLE_THREE_BODY_RATE](#GRACKLE_THREE_BODY_RATE)  
+[GRACKLE_CIE_COOLING](#GRACKLE_CIE_COOLING)  
+[GRACKLE_H2_OPA_APPROX](#GRACKLE_H2_OPA_APPROX)  
+[GRACKLE_USE_V_HEATING_RATE](#GRACKLE_USE_V_HEATING_RATE)  
+[GRACKLE_USE_S_HEATING_RATE](#GRACKLE_USE_S_HEATING_RATE)  
+[GRACKLE_HYDROGEN_MFRAC](#GRACKLE_HYDROGEN_MFRAC)  
+[OPT__UNFREEZE_GRACKLE](#OPT__UNFREEZE_GRACKLE)  
 
 
 Parameters below are shown in the format:   **`Name`   (Valid Values)   [Default Value]**
@@ -24,62 +31,110 @@ Only applicable when enabling the compilation option
 <a name="GRACKLE_VERBOSE"></a>
 * #### `GRACKLE_VERBOSE` &ensp; (0=off, 1=on) &ensp; [1]
     * **Description:**
-Map to the "grackle_verbose" option in GRACKLE.
+Map to the ["grackle_verbose" option in GRACKLE](https://grackle.readthedocs.io/en/latest/Interaction.html#enabling-output).
     * **Restriction:**
 
 <a name="GRACKLE_COOLING"></a>
 * #### `GRACKLE_COOLING` &ensp; (0=off, 1=on) &ensp; [1]
     * **Description:**
-Map to the "with_radiative_cooling" runtime parameter in GRACKLE.
+Map to the ["with_radiative_cooling" runtime parameter in GRACKLE](https://grackle.readthedocs.io/en/latest/Parameters.html#c.with_radiative_cooling).
     * **Restriction:**
 
 <a name="GRACKLE_PRIMORDIAL"></a>
 * #### `GRACKLE_PRIMORDIAL` &ensp; (0=Cloudy, 1=6-species, 2=9-species, 3=12-species) &ensp; [0]
     * **Description:**
-Map to the "primordial_chemistry" runtime parameter in GRACKLE.
+Map to the ["primordial_chemistry" runtime parameter in GRACKLE](https://grackle.readthedocs.io/en/latest/Parameters.html#c.primordial_chemistry).
 One must increase
 [[--passive | Installation:-Option-List#--passive]]
-by 3, 6, or 9 for GRACKLE_PRIMORDIAL=1, 2, or 3, respectively.
+by 6, 9, or 12 for GRACKLE_PRIMORDIAL=1, 2, or 3, respectively.
     * **Restriction:**
 
 <a name="GRACKLE_METAL"></a>
 * #### `GRACKLE_METAL` &ensp; (0=off, 1=on) &ensp; [0]
     * **Description:**
-Map to the "metal_cooling" runtime parameter in GRACKLE. One must increase
+Map to the ["metal_cooling" runtime parameter in GRACKLE](https://grackle.readthedocs.io/en/latest/Parameters.html#c.metal_cooling). One must increase
 [[--passive | Installation:-Option-List#--passive]]
 by 1 and initialize the field `Metal` using the field index `Idx_Metal` properly.
     * **Restriction:**
 
 <a name="GRACKLE_UV"></a>
 * #### `GRACKLE_UV` &ensp; (0=off, 1=on) &ensp; [0]
     * **Description:**
-Map to the "UVbackground" runtime parameter in GRACKLE.
+Map to the ["UVbackground" runtime parameter in GRACKLE](https://grackle.readthedocs.io/en/latest/Parameters.html#c.UVbackground).
     * **Restriction:**
 
 <a name="GRACKLE_CMB_FLOOR"></a>
 * #### `GRACKLE_CMB_FLOOR` &ensp; (0=off, 1=on) &ensp; [1]
     * **Description:**
-Map to the "cmb_temperature_floor" runtime parameter in GRACKLE.
+Map to the ["cmb_temperature_floor" runtime parameter in GRACKLE](https://grackle.readthedocs.io/en/latest/Parameters.html#c.cmb_temperature_floor).
     * **Restriction:**
 
 <a name="GRACKLE_PE_HEATING"></a>
 * #### `GRACKLE_PE_HEATING` &ensp; (0=off, 1=on) &ensp; [0]
     * **Description:**
-Map to the "photoelectric_heating" runtime parameter in GRACKLE.
+Map to the ["photoelectric_heating" runtime parameter in GRACKLE](https://grackle.readthedocs.io/en/latest/Parameters.html#c.photoelectric_heating).
     * **Restriction:**
 
 <a name="GRACKLE_PE_HEATING_RATE"></a>
 * #### `GRACKLE_PE_HEATING_RATE` &ensp; (&#8805;0.0) &ensp; [8.5e-26]
     * **Description:**
-Map to the "photoelectric_heating_rate" runtime parameter in GRACKLE.
+Map to the ["photoelectric_heating_rate" runtime parameter in GRACKLE](https://grackle.readthedocs.io/en/latest/Parameters.html#c.photoelectric_heating_rate).
 Note that the input value should always be in units of
-<var>erg</var>&#8287;<var>cm</var><sup>-3</sup>&#8287;<var>s</var><sup>-1</sup>.
+<var>erg</var>&#8287;<var>cm</var><sup>-3</sup>&#8287;<var>s</var><sup>-1</sup>&#8287;<var>n</var><sub>H</sub><sup>-1</sup>.
     * **Restriction:**
 
 <a name="GRACKLE_CLOUDY_TABLE"></a>
 * #### `GRACKLE_CLOUDY_TABLE` &ensp; (string) &ensp; [none]
     * **Description:**
-Map to the "grackle_data_file" runtime parameter in GRACKLE.
+Map to the ["grackle_data_file" runtime parameter in GRACKLE](https://grackle.readthedocs.io/en/latest/Parameters.html#c.grackle_data_file).
+    * **Restriction:**
+
+<a name="GRACKLE_THREE_BODY_RATE"></a>
+* #### `GRACKLE_THREE_BODY_RATE` &ensp; (0=Abel+02, 1=Palla+83, 2=Cohen+83, 3=Flower+07, 4=Glover+08, 5=Forrey+13) &ensp; [0]
+    * **Description:**
+Map to the ["three_body_rate" runtime parameter in GRACKLE](https://grackle.readthedocs.io/en/latest/Parameters.html#c.three_body_rate).
+    * **Restriction:**
+
+<a name="GRACKLE_CIE_COOLING"></a>
+* #### `GRACKLE_CIE_COOLING` &ensp; (0=off, 1=on) &ensp; [0]
+    * **Description:**
+Map to the ["cie_cooling" runtime parameter in GRACKLE](https://grackle.readthedocs.io/en/latest/Parameters.html#c.cie_cooling).
+    * **Restriction:**
+
+<a name="GRACKLE_H2_OPA_APPROX"></a>
+* #### `GRACKLE_H2_OPA_APPROX` &ensp; (0=off, 1=Ripomonti+04) &ensp; [0]
+    * **Description:**
+Map to the ["h2_optical_depth_approximation" runtime parameter in GRACKLE](https://grackle.readthedocs.io/en/latest/Parameters.html#c.h2_optical_depth_approximation).
+    * **Restriction:**
+
+<a name="GRACKLE_USE_V_HEATING_RATE"></a>
+* #### `GRACKLE_USE_V_HEATING_RATE` &ensp; (0=off, 1=on) &ensp; [0]
+    * **Description:**
+Map to the ["use_volumetric_heating_rate" runtime parameter in GRACKLE](https://grackle.readthedocs.io/en/latest/Parameters.html#c.use_volumetric_heating_rate).
+    * **Restriction:**
+
+<a name="GRACKLE_USE_S_HEATING_RATE"></a>
+* #### `GRACKLE_USE_S_HEATING_RATE` &ensp; (0=off, 1=on) &ensp; [0]
+    * **Description:**
+Map to the ["use_specific_heating_rate" runtime parameter in GRACKLE](https://grackle.readthedocs.io/en/latest/Parameters.html#c.use_specific_heating_rate).
+    * **Restriction:**
+
+<a name="GRACKLE_HYDROGEN_MFRAC"></a>
+* #### `GRACKLE_PE_HEATING_RATE` &ensp; (&#8805;0.0) &ensp; [0.76]
+    * **Description:**
+Map to the ["HydrogenFractionByMass" runtime parameter in GRACKLE](https://grackle.readthedocs.io/en/latest/Parameters.html#c.HydrogenFractionByMass).
+Note that the input value will only set to Grackle when it is in non-equilibrium mode [GRACKLE_PRIMORDIAL](#GRACKLE_PRIMORDIAL)>0,
+because the tables for tabulated mode were created assuming hydrogen mass fraction of about 0.716.
+See the Grackle document for the details.
+When in tabulated mode ([GRACKLE_PRIMORDIAL](#GRACKLE_PRIMORDIAL)==0), this value may still be used elsewhere in GAMER and
+it may be different from HydrogenFractionByMass inside Grackle.
+    * **Restriction:**
+
+<a name="OPT__UNFREEZE_GRACKLE"></a>
+* #### `OPT__UNFREEZE_GRACKLE` &ensp; (0=off, 1=on) &ensp; [0]
+    * **Description:**
+Allow the evolution by Grackle solver when the fluid is frozen
+(for [[OPT__FREEZE_FLUID | Runtime-Parameters:-Hydro#OPT__FREEZE_FLUID]]==1 only).
     * **Restriction:**
 
 

doc/wiki/Runtime-Parameters-related/Runtime-Parameters:-Outputs.md

@@ -21,6 +21,9 @@ Parameters described on this page:
 [OPT__OUTPUT_LORENTZ](#OPT__OUTPUT_LORENTZ),  
 [OPT__OUTPUT_3VELOCITY](#OPT__OUTPUT_3VELOCITY),  
 [OPT__OUTPUT_ENTHALPY](#OPT__OUTPUT_ENTHALPY),  
+[OPT__OUTPUT_GRACKLE_TEMP](#OPT__OUTPUT_GRACKLE_TEMP),  
+[OPT__OUTPUT_GRACKLE_MU](#OPT__OUTPUT_GRACKLE_MU),  
+[OPT__OUTPUT_GRACKLE_TCOOL](#OPT__OUTPUT_GRACKLE_TCOOL),  
 [OPT__OUTPUT_USER_FIELD](#OPT__OUTPUT_USER_FIELD),  
 [OPT__OUTPUT_MODE](#OPT__OUTPUT_MODE),  
 [OPT__OUTPUT_RESTART](#OPT__OUTPUT_RESTART),  
@@ -211,6 +214,27 @@ Output the SRHD reduced enthalpy on grids.
     * **Restriction:**
 For [[--srhd | Installation:-Option-List#--srhd]] only.
 
+<a name="OPT__OUTPUT_GRACKLE_TEMP"></a>
+* #### `OPT__OUTPUT_GRACKLE_TEMP` &ensp; (0=off, 1=on) &ensp; [0]
+    * **Description:**
+Output the temperature calculated by Grackle.
+    * **Restriction:**
+For [[--grackle | Installation:-Option-List#--grackle]] only.
+
+<a name="OPT__OUTPUT_GRACKLE_MU"></a>
+* #### `OPT__OUTPUT_GRACKLE_MU` &ensp; (0=off, 1=on) &ensp; [0]
+    * **Description:**
+Output the mean molecular weight calculated by Grackle.
+    * **Restriction:**
+For [[--grackle | Installation:-Option-List#--grackle]] only.
+
+<a name="OPT__OUTPUT_GRACKLE_TCOOL"></a>
+* #### `OPT__OUTPUT_GRACKLE_TCOOL` &ensp; (0=off, 1=on) &ensp; [0]
+    * **Description:**
+Output the cooling time calculated by Grackle.
+    * **Restriction:**
+For [[--grackle | Installation:-Option-List#--grackle]] only.
+
 <a name="OPT__OUTPUT_MODE"></a>
 * #### `OPT__OUTPUT_MODE` &ensp; (1=const step, 2=const dt, 3=dump table) &ensp; [none]
     * **Description:**

doc/wiki/Runtime-Parameters-related/Runtime-Parameters:-Refinement.md

@@ -34,6 +34,7 @@ Parameters described on this page:
 [OPT__FLAG_LRTZ_GRADIENT](#OPT__FLAG_LRTZ_GRADIENT),  
 [OPT__FLAG_VORTICITY](#OPT__FLAG_VORTICITY),  
 [OPT__FLAG_JEANS](#OPT__FLAG_JEANS),  
+[OPT__FLAG_COOLING_LEN](#OPT__FLAG_COOLING_LEN),  
 [OPT__FLAG_CURRENT](#OPT__FLAG_CURRENT),  
 [OPT__FLAG_CRAY](#OPT__FLAG_CRAY),  
 [OPT__FLAG_LOHNER_DENS](#OPT__FLAG_LOHNER_DENS),  
@@ -225,6 +226,23 @@ An example file can be found at `example/input/Input__Flag_Jeans`.
 Recommended values: ≥4.
     * **Restriction:**
 
+<a name="OPT__FLAG_COOLING_LEN"></a>
+* #### `OPT__FLAG_COOLING_LEN` &ensp; (0=off, 1=on) &ensp; [0]
+    * **Description:**
+Refinement criterion: gas cooling length. It ensures that the cooling length
+is resolved by at least <var>N</var> cells. Specifically, a cell
+on level <var>l</var> will be flagged for refinement if its estimated
+cooling length <var>L</var><sub>cool</sub> satisfies
+<var>L</var><sub>cool</sub>&#8287;&#8801;<var>t</var><sub>cool</sub><var>c</var><sub>s</sub>&#8287;&#8287;<&#8287;<var>N</var><sub>l</sub>&Delta;&xi;<sub>l</sub>,
+where <var>t</var><sub>cool</sub> is the cooling time (currently calculated by Grackle), <var>c</var><sub>s</sub> is sound speed of gas, &Delta;&xi;<sub>l</sub> is the cell width along &xi; on level <var>l</var>, and <var>N</var><sub>l</sub> is the refinement threshold on level <var>l</var>.
+Specify the refinement
+thresholds on different levels in the input file `Input__Flag_CoolingLen`
+with the [[specific format | Runtime-Parameters:-Input__Flag_*]].
+An example file can be found at `example/input/Input__Flag_CoolingLen`.
+Recommended values: &#8805;1.
+    * **Restriction:**
+Must compile with [[--grackle | Installation:-Option-List#--grackle]].
+
 <a name="OPT__FLAG_CURRENT"></a>
 * #### `OPT__FLAG_CURRENT` &ensp; (0=off, 1=on) &ensp; [0]
     * **Description:**

doc/wiki/Runtime-Parameters-related/Runtime-Parameters:-Timestep.md

@@ -18,6 +18,7 @@ Parameters described on this page:
 [DT__PARVEL](#DT__PARVEL),  
 [DT__PARVEL_MAX](#DT__PARVEL_MAX),  
 [DT__PARACC](#DT__PARACC),  
+[DT__GRACKLE_COOLING](#DT__GRACKLE_COOLING),  
 [DT__CR_DIFFUSION](#DT__CR_DIFFUSION),  
 [DT__MAX_DELTA_A](#DT__MAX_DELTA_A),  
 [DT__SYNC_PARENT_LV](#DT__SYNC_PARENT_LV),  
@@ -108,6 +109,14 @@ for the exact formula.
 Only applicable when adopting the compilation option
 [[--store_par_acc | Installation:-Option-List#--store_par_acc]].
 
+<a name="DT__GRACKLE_COOLING"></a>
+* #### `DT__GRACKLE_COOLING` &ensp; (&#8805;0.0; <0.0 &#8594; off) &ensp; [-1.0]
+    * **Description:**
+Factor for Grackle cooling time.
+    * **Restriction:**
+Only applicable when adopting the compilation option
+[[--grackle | Installation:-Option-List#--grackle]].
+
 <a name="DT__CR_DIFFUSION"></a>
 * #### `DT__CR_DIFFUSION` &ensp; (&#8805;0.0) &ensp; [0.3]
     * **Description:**

example/input/Input__Flag_CoolingLen

@@ -0,0 +1,13 @@
+# Level        cooling_length/cell_size
+      0                             1.0
+      1                             1.0
+      2                             1.0
+      3                             1.0
+      4                             1.0
+      5                             1.0
+      6                             1.0
+      7                             1.0
+      8                             1.0
+      9                             1.0
+     10                             1.0
+     11                             1.0

example/test_problem/Hydro/AGORA_IsolatedGalaxy/Input_Options/Input__Parameter.high-res

@@ -131,8 +131,15 @@ GRACKLE_METAL                 1           # ...    "metal_cooling" (must increas
 GRACKLE_UV                    1           # ...    "UVbackground" [0]
 GRACKLE_CMB_FLOOR             1           # ...    "cmb_temperature_floor" [1]
 GRACKLE_PE_HEATING            1           # ...    "photoelectric_heating" [0]
-GRACKLE_PE_HEATING_RATE       8.5e-26     # ...    "photoelectric_heating_rate (in erg/cm^3/s)" [8.5e-26]
+GRACKLE_PE_HEATING_RATE       8.5e-26     # ...    "photoelectric_heating_rate" (in erg/cm^3/s/n_H) [8.5e-26]
 GRACKLE_CLOUDY_TABLE          CloudyData_UVB=HM2012.h5   # "grackle_data_file"
+GRACKLE_THREE_BODY_RATE       0           # map to "three_body_rate" (0=Abel+02, 1=Palla+83, 2=Cohen+83, 3=Flower+07, 4=Glover+08, 5=Forrey+13) [0]
+GRACKLE_CIE_COOLING           0           # map to "cie_cooling" [0]
+GRACKLE_H2_OPA_APPROX         0           # map to "h2_optical_depth_approximation" (0=off, 1=Ripomonti+04) [0]
+GRACKLE_USE_V_HEATING_RATE    0           # map to "use_volumetric_heating_rate" [0]
+GRACKLE_USE_S_HEATING_RATE    0           # map to "use_specific_heating_rate" [0]
+GRACKLE_HYDROGEN_MFRAC        0.76        # map to "HydrogenFractionByMass" (NOT set to Grackle when GRACKLE_PRIMORDIAL==0) [0.76]
+OPT__UNFREEZE_GRACKLE         0           # allow the evolution by Grackle solver when the fluid is frozen (for OPT__FREEZE_FLUID==1 only) [0]
 CHE_GPU_NPGROUP              -1           # number of patch groups sent into the CPU/GPU Grackle solver (<=0=auto) [-1]
 
 

example/test_problem/Hydro/AGORA_IsolatedGalaxy/Input_Options/Input__Parameter.low-res

@@ -131,8 +131,15 @@ GRACKLE_METAL                 1           # ...    "metal_cooling" (must increas
 GRACKLE_UV                    1           # ...    "UVbackground" [0]
 GRACKLE_CMB_FLOOR             1           # ...    "cmb_temperature_floor" [1]
 GRACKLE_PE_HEATING            1           # ...    "photoelectric_heating" [0]
-GRACKLE_PE_HEATING_RATE       8.5e-26     # ...    "photoelectric_heating_rate (in erg/cm^3/s)" [8.5e-26]
+GRACKLE_PE_HEATING_RATE       8.5e-26     # ...    "photoelectric_heating_rate" (in erg/cm^3/s/n_H) [8.5e-26]
 GRACKLE_CLOUDY_TABLE          CloudyData_UVB=HM2012.h5   # "grackle_data_file"
+GRACKLE_THREE_BODY_RATE       0           # map to "three_body_rate" (0=Abel+02, 1=Palla+83, 2=Cohen+83, 3=Flower+07, 4=Glover+08, 5=Forrey+13) [0]
+GRACKLE_CIE_COOLING           0           # map to "cie_cooling" [0]
+GRACKLE_H2_OPA_APPROX         0           # map to "h2_optical_depth_approximation" (0=off, 1=Ripomonti+04) [0]
+GRACKLE_USE_V_HEATING_RATE    0           # map to "use_volumetric_heating_rate" [0]
+GRACKLE_USE_S_HEATING_RATE    0           # map to "use_specific_heating_rate" [0]
+GRACKLE_HYDROGEN_MFRAC        0.76        # map to "HydrogenFractionByMass" (NOT set to Grackle when GRACKLE_PRIMORDIAL==0) [0.76]
+OPT__UNFREEZE_GRACKLE         0           # allow the evolution by Grackle solver when the fluid is frozen (for OPT__FREEZE_FLUID==1 only) [0]
 CHE_GPU_NPGROUP              -1           # number of patch groups sent into the CPU/GPU Grackle solver (<=0=auto) [-1]
 
 

example/test_problem/Hydro/CMZ/Input__Parameter

@@ -137,11 +137,15 @@ GRACKLE_METAL                 1           # ...    "metal_cooling" (must increas
 GRACKLE_UV                    1           # ...    "UVbackground" [0]
 GRACKLE_CMB_FLOOR             1           # ...    "cmb_temperature_floor" [1]
 GRACKLE_PE_HEATING            1           # ...    "photoelectric_heating" [0]
-GRACKLE_PE_HEATING_RATE       8.5e-26     # ...    "photoelectric_heating_rate (in erg/cm^3/s)" [8.5e-26]
+GRACKLE_PE_HEATING_RATE       8.5e-26     # ...    "photoelectric_heating_rate" (in erg/cm^3/s/n_H) [8.5e-26]
 GRACKLE_CLOUDY_TABLE          CloudyData_UVB=HM2012.h5   # "grackle_data_file"
-GRACKLE_THREE_BODY_RATE       0           # used Glover+08 rate
-GRACKLE_CIE_COOLING           0           # 0: off; 1:on
-GRACKLE_H2_OPA_APPROX         0           # H2 opacity from Ripamonti+04; 0:off, 1:Ripomonti+04
+GRACKLE_THREE_BODY_RATE       0           # map to "three_body_rate" (0=Abel+02, 1=Palla+83, 2=Cohen+83, 3=Flower+07, 4=Glover+08, 5=Forrey+13) [0]
+GRACKLE_CIE_COOLING           0           # map to "cie_cooling" [0]
+GRACKLE_H2_OPA_APPROX         0           # map to "h2_optical_depth_approximation" (0=off, 1=Ripomonti+04) [0]
+GRACKLE_USE_V_HEATING_RATE    0           # map to "use_volumetric_heating_rate" [0]
+GRACKLE_USE_S_HEATING_RATE    0           # map to "use_specific_heating_rate" [0]
+GRACKLE_HYDROGEN_MFRAC        0.76        # map to "HydrogenFractionByMass" (NOT set to Grackle when GRACKLE_PRIMORDIAL==0) [0.76]
+OPT__UNFREEZE_GRACKLE         0           # allow the evolution by Grackle solver when the fluid is frozen (for OPT__FREEZE_FLUID==1 only) [0]
 CHE_GPU_NPGROUP              -1           # number of patch groups sent into the CPU/GPU Grackle solver (<=0=auto) [-1]
 
 

example/test_problem/Hydro/CR_Diffusion/Input__Parameter

@@ -160,11 +160,15 @@ GRACKLE_METAL                 0           # ...    "metal_cooling" (must increas
 GRACKLE_UV                    0           # ...    "UVbackground" [0]
 GRACKLE_CMB_FLOOR             1           # ...    "cmb_temperature_floor" [1]
 GRACKLE_PE_HEATING            0           # ...    "photoelectric_heating" [0]
-GRACKLE_PE_HEATING_RATE       8.5e-26     # ...    "photoelectric_heating_rate (in erg/cm^3/s)" [8.5e-26]
+GRACKLE_PE_HEATING_RATE       8.5e-26     # ...    "photoelectric_heating_rate" (in erg/cm^3/s/n_H) [8.5e-26]
 GRACKLE_CLOUDY_TABLE          CloudyData_noUVB.h5  # "grackle_data_file"
-GRACKLE_THREE_BODY_RATE       4           # used Glover+08 rate
-GRACKLE_CIE_COOLING           1           # 0: off; 1:on
-GRACKLE_H2_OPA_APPROX         1           # H2 opacity from Ripamonti+04; 0:off, 1:Ripomonti+04
+GRACKLE_THREE_BODY_RATE       0           # map to "three_body_rate" (0=Abel+02, 1=Palla+83, 2=Cohen+83, 3=Flower+07, 4=Glover+08, 5=Forrey+13) [0]
+GRACKLE_CIE_COOLING           0           # map to "cie_cooling" [0]
+GRACKLE_H2_OPA_APPROX         0           # map to "h2_optical_depth_approximation" (0=off, 1=Ripomonti+04) [0]
+GRACKLE_USE_V_HEATING_RATE    0           # map to "use_volumetric_heating_rate" [0]
+GRACKLE_USE_S_HEATING_RATE    0           # map to "use_specific_heating_rate" [0]
+GRACKLE_HYDROGEN_MFRAC        0.76        # map to "HydrogenFractionByMass" (NOT set to Grackle when GRACKLE_PRIMORDIAL==0) [0.76]
+OPT__UNFREEZE_GRACKLE         0           # allow the evolution by Grackle solver when the fluid is frozen (for OPT__FREEZE_FLUID==1 only) [0]
 CHE_GPU_NPGROUP              -1           # number of patch groups sent into the CPU/GPU Grackle solver (<=0=auto) [-1]