Program readff95 c double precicion to avoid round-off errors real*8 lat,long,popli,ff95,perc real*8 spop,sff95,sperc character*20 name open(10,file='ff95.dat',status='old') c prepare for summations: spop=0.d0 sff95=0.d0 sperc=0.d0 do 100 i=1,19969 read(10,10,end=911)lat,long,ff95,popli,perc,unid,name 10 format(f6.1,x,f6.1,2x,f8.2,2x,e10.0,2x,f8.4,2x,f3.1,2x,a20) write(12,10)lat,long,ff95,popli,perc,unid,name spop=spop+popli sperc=sperc+perc sff95=sff95+ff95 100 continue 911 continue i=i-1 c write summations to screen: write(*,*)'spop=',spop write(*,*) 'sperc=',sperc write(*,*) 'sff95=',sff95 write(*,*) 'lines read=',i stop end