pro stnd792,nopr,iout,$
    ave,ave2,aven,stnd

; ******************************
; Calculate the standard deviation
      s0=double(aven)
      s0m1=double(s0-1.00d0)
      stnd=-999.0d0
     if (s0m1 ge 2.0d0) then begin
      s1=ave
      s2=ave2
      f1=(s0*s2)-(s1*s1)
      f2=s0*s0m1
      if (f1 gt 1.0d-10) then begin
       stnd=sqrt(f1/f2)
      endif
     endif

     if (aven ge 0.5d0) then begin
      ave=ave/aven
     endif

; ******************************
    if (nopr eq 1) then begin
     printf,iout,'  '
     printf,iout,'  stnd792: ave,ave2,aven ',ave,ave2,aven
     printf,iout,'  stnd792: stnd ',stnd
    endif

; ******************************
   return
   end