real a(50) sum=0. stdev=0. n=0 1 read(5,*,end=999,err=999)amag n=n+1 a(n)=amag sum=sum+amag goto 1 999 continue if (n.ge.1)then ave=sum/n if(n.ge.2)then do i=1,n stdev=stdev+(a(i)-ave)**2 enddo stdev=stdev/(n-1) stdev=sqrt(stdev) else stdev= 999 endif write(6,*)ave, stdev else write(6,*)'nodata' endif end