Hello! I think it should be ustar = SQRT (wStress(i,j,bi,bj) / atmrho) rather than ustar = wStress(i,j,bi,bj)*recip_sqrtRhoA in exf package in S/R exf_winf.F (line 179)? In the current version wStress is not affected by the function SQRT, but should be, in my opinion. Best regards, Stanislav