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