Analytical model for neutral river flow is started

source/model/turb_mod.f90

@@ -1836,5 +1836,19 @@ enddo
 
 END SUBROUTINE ED_TEMP_HOSTETLER2
 
+!> An analytical model for neutrally stratified river flow using Obukhov coordinate
+SUBROUTINE ZIRIVMODEL(M,h1,u_star_surf,gradp,z_half,eps_ziriv)
+use PHYS_CONSTANTS, only : kappa
+use NUMERIC_PARAMS, only : pi
+implicit none
+integer(kind=iintegers), intent(in) :: M      !> Number of layers
+real(kind=ireals), intent(in ) :: h1          !> Riverflow depth, m
+real(kind=ireals), intent(in ) :: u_star_surf !> Friction velocity in water at the surface, m/s
+real(kind=ireals), intent(in ) :: gradp       !> Longitudinal pressure gradient, Pa/m
+real(kind=ireals), intent(in ) :: z_half(1:M) !> A set of depths where analytical solution is to be evaluated
+real(kind=ireals), intent(out) :: eps_ziriv(1:M) !> TKE dissipation rate, m**2/s**3
+
+
+END SUBROUTINE ZIRIVMODEL
 
 END MODULE TURB