subroutine timeav(seq,tim,av,sdv)
c
      include 'itg.par'
      real seq(nez),tim,av,sdv,deltm
      integer i1a,i2b,ii
c
      include 'itg.cmn'
c
      i1a=navgflag+1
      i2b=ne-2
c
      av=0.5*seq(i1a)*(timo(i1a)-avgflag)
      do 20 i=i1a,i2b
      ii=i+1
      av=av+seq(ii)*(timo(ii)-timo(i))
   20 continue
      av=av+0.5*seq(ne-1)*(timo(ne)-timo(ne-1))
      deltm=tim-avgflag
      if (deltm.gt.1.e-33) av=av/deltm
c
      sdv=0.5*(seq(i1a)-av)*(seq(i1a)-av)
     *  *(timo(i1a)-avgflag)
      do 30 i=i1a,i2b
      ii=i+1
      sdv=sdv+(seq(ii)-av)*(seq(ii)-av)*(timo(ii)-timo(i))
   30 continue
      sdv=sdv+0.5*(seq(ne-1)-av)*(seq(ne-1)-av)*
     *    (timo(ne)-timo(ne-1))
      if (deltm.gt.1.e-33) sdv=sdv/deltm
      sdv=sqrt(sdv)
c
      return
      end