It looks like we compute the snow physics twice when running with groundwater? We call soil_snow in cbm here, which does all the snow physics here. But then the groundwater code additionally calls snow_processes_soil_thermal, which also does all the same things? @ccarouge @JhanSrbinovsky Is this duplication intentional, or is there some logic I'm missing that prevents this happening twice?