!
! File: max_speed.f
! Author: Christopher Moore
!
! Calculates the maximum speed (as max_slab calculates maximum amplitude)
! Called from main program, variable output in SIFT output file
!
      subroutine maxspeed(u,v,q,maxs,n1,n2,mm,nn,dry)
      integer*4 i,j,n1,n2,mm,nn
      real*8 u(mm,nn),v(mm,nn),q(mm,nn)
      real*4 maxs(mm,nn)
      real*8 holdit, dry

      common /nc/ ns1,ne1,ns2,ne2,istep,jstep,scl_m2cm

      do i = 1,n1
         do j = 1,n2
            if(q(i,j).GE.dry) then
              holdit = scl_m2cm*sqrt(u(i,j)**2+v(i,j)**2)
              if(holdit.GE.maxs(i,j)) then
                   maxs(i,j) = holdit
              endif
            endif
          end do
       end do

       return
       end