In dust_related_props, we calculate the dust-to-gas ratio if photoelectric_heating is enabled. However, we either use the dust_density field if use_dust_density_field == 1 or we use the value of metallicity * local_dust_to_gas_ratio. Unless I am missing something, we are not doing the proper dust-to-gas calculation of summing all the Chiaki dust species when dust_chemistry == 2.
In
dust_related_props, we calculate the dust-to-gas ratio ifphotoelectric_heatingis enabled. However, we either use thedust_densityfield ifuse_dust_density_field == 1or we use the value of metallicity *local_dust_to_gas_ratio. Unless I am missing something, we are not doing the proper dust-to-gas calculation of summing all the Chiaki dust species whendust_chemistry == 2.