[KUMAR7491]

I am trying to add bottom stress term in momentum equation of shallowwaterfoam.
The momentum equation is
ddt(hu) + h(phiv)u + f x hu = -magg*h*laplacian(h+h0)- bottom stress

here, h : mean surface height, u : velocity vector, h0 deviation from mean surface height, g : gravitational accl. and

Bottom stress along x direction =Cd*(sqrt (u^2+v^2))*u
Bottom stress along y direction =Cd*(sqrt (u^2+v^2))*v
here u : velocity along x direction,
y : velocity along y direction.
Cd: constant
Bottom stress is function of both u and v.
I gone through Patanlar book, It is necessary to linearise this source term using Tyler method. Is any one know how it is linearised? and how should i add these source tern in shallow water equation to solve implicitly.