1subroutine averms(x,n,nskip,ave,rms)
2  real x(n)
3  integer ipk(1)
4
5  ns=0
6  s=0.
7  sq=0.
8  ipk=maxloc(x)
9  do i=1,n
10     if(abs(i-ipk(1)).gt.nskip) then
11        s=s + x(i)
12        sq=sq + x(i)**2
13        ns=ns+1
14     endif
15  enddo
16  ave=s/ns
17  rms=sqrt(sq/ns - ave*ave)
18
19  return
20end subroutine averms
21