Andreas (1987, BLM) is a well-known model for scalar roughness over snow and sea ice. While many of us utilize his quadratic approximation, is there anyone who has implemented the model using the theoretical equations? My calculation shows the slope of the normalized scalar roughness against the roughness Reynolds number becomes less steep than the results obtained from the quadratic approximation. The set of equations introduced in Andreas (1987) also seems to suggest that the relationship between the normalized roughness length and the Reynolds number (R∗) approaches a slope of -1 on a logarithmic scale as R∗→∞, which also indicates a less steep decline compared to the quadratic results. If zeta_hat​ (the nondimensional matching height over a rough surface) is assumed to be constant, the results align with the quadratic equation. However, I am uncertain whether this interpretation aligns with what Prof. Andreas intended.