git commit -m "changes in velocities"

setup/windforc_setup.dat

@@ -134,7 +134,7 @@ stabfunc 2
 kepsbc   1
 soiltype 5
 soil_depth  10.
-soilswitch  1
+soilswitch  0
 nsoilcols 5
 tricemethhydr 0.
 carbon_model 1

source/model/momentum_mod.f90

@@ -171,6 +171,7 @@ real(kind=4), allocatable :: Amatrix(:,:), rhs(:) !to be passed to lapack routin
 real(kind=ireals) :: tau_ix, tau_iy
 real(kind=ireals) :: Lxi, Lxj, Lyi, Lyj
 real(kind=ireals) :: tau_grx, tau_gry
+real(kind=ireals) :: velx_log, vely_log
 real(kind=ireals) :: time1, time2, totime
 real(kind=ireals) :: scross, hydr_radius, disch, speedav 
 real(kind=ireals), allocatable :: row0_(:)
@@ -334,8 +335,10 @@ endif
  !coef2 = g*Cz_gr**(-2)*sqrt(u1(M+1)**2+v1(M+1)**2) !*0.005 !*row0
 
  !Bottom momentum exchange coefficient from logarithmic profile
- coef2 = sqrt(u1(M+1)**2+v1(M+1)**2)*kappa**2/log(h1*0.25*ddz(M)/z0_bot)**2
- velfrict_bot = sqrt(u1(M+1)**2+v1(M+1)**2)*kappa/log(h1*0.25*ddz(M)/z0_bot)
+ velx_log = 0.5*(u1(M+1)+u1(M))
+ vely_log = 0.5*(v1(M+1)+v1(M))
+ coef2 = sqrt(velx_log**2+vely_log**2)*kappa**2/log(h1*0.5*ddz(M)/z0_bot)**2
+ velfrict_bot = sqrt(velx_log**2+vely_log**2)*kappa/log(h1*0.5*ddz(M)/z0_bot)
  tau_gr = row0*velfrict_bot**2
 
 ! Friction at the water - upper ice interface
@@ -446,7 +449,7 @@ SELECT CASE (Turbpar%par)
       endif
     enddo
     do i = 1, M
-      k2(i) = 0.00632!max(CE(i)*E1(i)**2/(eps1(i) + ACCk),min_visc) + niu_wat + knum(i)
+      k2(i) = 0.00632 !max(CE(i)*E1(i)**2/(eps1(i) + ACCk),min_visc) + niu_wat + knum(i) !!
     end do
 
 ! 3. Nickuradze (NICK) formulation: Rodi (1993)

source/model/turb_mod.f90

@@ -991,7 +991,7 @@ allocate (rhotemp(1:M), rhosal(1:M))
   endif
 
   do i = 1, M
-    KT(i) = 0.00632*cw_m_row0 !max(CEt(i)*E2(i)**2/(eps2(i) + ACCk),min_diff)*cw_m_row0 !E2, eps2
+    KT(i) = 0.00632*cw_m_row0 !max(CEt(i)*E2(i)**2/(eps2(i) + ACCk),min_diff)*cw_m_row0 !E2, eps2 ! !
   enddo
 
   ! Diagnostic calculation