Change interpolation of near surface wind speed from linear to logarithmic
The surface wind speed is computed in file "module_sf_snowpacklsm.F":
!! GET U_SFC and V_SFC , i.e., getting surface layer wind speeds
do j=jts,jte
do i=its,ite
Z_ML = 0.5*dz8w(I,1,J)
do k=1,snc_N
U_SFC(i,k,j) = (U_PHY(I,1,J)/Z_ML) * snc_Z(I,k,J)
V_SFC(i,k,j) = (V_PHY(I,1,J)/Z_ML) * snc_Z(I,k,J)
enddo
enddo
enddo
A linear interpolation is used instead of a logarithmic one. The current option leads to an underestimation of the mass flux of blowing snow particles.