subroutine umasfl(flow,kflow,kstep,kinc,time,node,coords)
C     
c      include 'cdp.cmn'
      include 'ABA_PARAM.INC'
      dimension flow(kflow),time(2),coords(3)
C     
      h=sqrt(coords(1)**2+coords(3)**2)
      if(time(2).le.0.2)Then
         c=1.113E-6*h
      else if(time(2).gt.0.2.and.time(2).le.5.0)Then
         c=1.113E-6*h*(5.0-time(2))/4.8+1.113E-6*h*1.6E-5
      else
         c=1.113E-6*h*1.6E-5
      endif
      flow(1)=c*coords(3)/h
      flow(3)=-c*coords(1)/h
      flow(2)=0.0
      return
      end