subroutine est_weirdim(depwid,desdis,wwidth,wdepth,cd) !! This program estimates rectangular weir dimensions based on !! width-depth ratio of wier at different stages implicit none real, intent(in) :: depwid ! | real, intent(in) :: cd ! | real, intent(in) :: desdis ! | real, intent(out) :: wdepth ! | real, intent(out) :: wwidth ! | real :: tempvar = 0. ! | tempvar=depwid**1.5 wwidth=(desdis*tempvar/(1.84*cd))**0.4 wdepth=wwidth/depwid return end subroutine est_weirdim