Thanks for reply. I use the MM5 surface layer scheme, thus I need to modify the module_sf_sfclay.F. However, I found that the Ck and Cd is only diagnostic variable and is determined by ZNT in the source code. When I set Cd=0 and Ck=0.001, I wonder if these changes will be adopted for the calculation of surface heat and momentum fluxes?