Add Dust Cooling Support in Grackle

[Xiang-Zhe-Wei]

Summary

• Implement dust cooling with subcycling (modified solve_rate_cool_g.F)
• Enable switch via COOLING=0/1 in Input__Parameter
• Provide setup & analysis scripts for constant-density cooling test

Compilation

• MODEL=HYDRO
• SUPPORT_GRACKLE

Simulation Setup

1. Units: [L]=kpc, [M]=1e9 Msun, [T]=Myr
2. Gas Init: test radiative cooling at fixed density
   • COOLING=0 → no cooling
   • COOLING=1 → enable cooling, replace:
     grackle/src/clib/solve_rate_cool_g.F
     with dust version:
     https://github.com/Xiang-Zhe-Wei/grackle-dust/blob/tommy_dust/src/clib/solve_rate_cool_g.F
3. Metallicity & Dust: Gas_MetalMassFrac, dust/gas=0.01
4. Output & Time Step: fixed-step
5. Dust Cooling Mod: subcycling

Notes

• Mimics constant-density cooling test
• Scripts: dust_density.py, dust_temp.py, gas_temp.py, metal_density.py
• Confirm Grackle edot setting before running
• Install: see ./example/grackle/README

[technic960183]

@hyschive I noticed that the file src/TestProblem/Hydro/Dust/Init_TestProb_Hydro_Dust_Rho_Temp_copy_const_dt.cpp
was committed and then deleted in consecutive commits.
I suggest cleaning this up to improve the repository health.

[hyschive]

@Xiang-Zhe-Wei Just a reminder. Please check the above suggestion.

[Xiang-Zhe-Wei]

@hyschive I noticed that the file src/TestProblem/Hydro/Dust/Init_TestProb_Hydro_Dust_Rho_Temp_copy_const_dt.cpp was committed and then deleted in consecutive commits. I suggest cleaning this up to improve the repository health.

Thanks for the reminder. I’ve addressed the issue by cleaning up the commit history. Please take another look.

[hyschive]

@technic960183 Could you help confirm it? Thanks!

[technic960183]

@hyschive Looks good to me.
@Xiang-Zhe-Wei Thanks for your contribution.