The surface boundary layer option is not enabled in the Phillips two-layer example (BULKMIXEDLAYER = False). The wind stress is distributed over the top-most 20 m.
DIRECT_STRESS = True ! [Boolean] default = False
! If true, the wind stress is distributed over the topmost HMIX_STRESS of fluid
! (like in HYCOM), and KVML may be set to a very small value.
HMIX_STRESS = 20.0 ! [m] default = 20.0
! The depth over which the wind stress is applied if DIRECT_STRESS is true.