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.