yifan zhuang
New member
Dear WRF developers,
I am running WRF v4.3.3 and using the following surface-layer scheme:
sf_sfclay_physics = 1 # Revised MM5 surface layer
I would like to modify the surface drag coefficient cd so that it no longer follows the Monin–Obukhov formulation. Instead, I would like Cd to follow a user-defined function of the 10-m wind speed (U10),
Any guidance on which routines to modify, or examples of similar modifications, would be greatly appreciated.
Thank you very much!
I am running WRF v4.3.3 and using the following surface-layer scheme:
sf_sfclay_physics = 1 # Revised MM5 surface layer
I would like to modify the surface drag coefficient cd so that it no longer follows the Monin–Obukhov formulation. Instead, I would like Cd to follow a user-defined function of the 10-m wind speed (U10),
My questions:
- Where in the code is the “active” Cd used by the surface-layer momentum flux actually computed?
I found the diagnostic Cd inside module_sf_sfclayrev.F:
Cd(I) = (karman/psix10)*(karman/psix10) - If I want WRF to truly use my custom Cd(U10), should I also modify ZNT, GZ10OZ0, so that the surface stress is consistent with my new drag law?
Any guidance on which routines to modify, or examples of similar modifications, would be greatly appreciated.
Thank you very much!