c     Dump NetCDF SAR data: calculate statistics

      subroutine sub_stat(nx,ny,uv,obs,obs3,std,mean,min,max,missing)
      implicit none
      integer k,kk,kkk,nx,ny,j,i,obs,obs3
      real    uv(nx,ny,1),missing,min,max,xx
      double precision mean, std
c    ***************************************************************
      mean=0.
      obs=0
      std=0
c    ***************************************************************

      do j=1,ny
      do i=1,nx

      if(uv(i,j,1).ne.missing) then
      mean=mean+uv(i,j,1)
      obs=obs+1
      endif
      
      enddo 
      enddo 
      if(obs.eq.0) mean=missing
      if(obs.ne.0) mean=mean/obs
c    ***************************************************************

      do j=1,ny
      do i=1,nx

      if(uv(i,j,1).ne.missing) then
      xx=uv(i,j,1)-mean
      std=std+xx*xx
      endif

      enddo
      enddo
      if(obs.eq.0) std=missing
      if(obs.ne.0) std=SQRT(std/obs)

c    ***************************************************************
      max=mean
      min=mean

      do 120 j=1,ny
      do 124 i=1,nx
      if(uv(i,j,1).ne.missing) then

      if(uv(i,j,1).lt.min) min=uv(i,j,1)
      if(uv(i,j,1).gt.max) max=uv(i,j,1)
      endif
 124  continue
 120  continue
c    ***************************************************************
      obs3=0
      do 130 j=1,ny
      do 134 i=1,nx

      if(uv(i,j,1).ne.missing) then
      xx=ABS(uv(i,j,1)-mean)/std
      if(xx.gt.3.0) obs3=obs3+1
      endif

 134  continue
 130  continue

c    ***************************************************************
      return
      end