Correction of the infrared flux in the 1D model when callslope

[llange]

When using the 1D model with a slope angle > 0, the atmosphere is in equilibrium with the cold sloped surface while it should be in equilibrium with the warm flat surface (Lange et al., GRL 2023). As a result, the infrared flux is underestimated and surface temperatures too low, especially during winter.
To correct this, I think a option should be done in the callphysdef (coldatmosphere / warmatmosphere) which will be consider only in the 1D and if callslope = .true., and the infrared flux / atmospheric temperature should be adapted if the user want a cold atmosphere or a warm atmosphere above his slope.