program reconstruct implicit none integer*4 ic,n,ie,i,j,np,ee character*80 s80 real*8 x,y real*8,allocatable::c(:),sx(:),s(:) character*80,allocatable::sn(:) ee=2 open(9,file='COEF.TXT',status='old') do 1 ic=1,2 n=0 222 read(9,90,end=99,err=99)s80 90 format(a80) n=n+1 if(ic.eq.2)then read(s80(5:23),*)c(n) sn(n)=s80 endif goto 222 99 if(ic.eq.1)then write(6,*)n,' coefficients in COEF.TXT' allocate(c(n),sn(n)) rewind 9 endif 1 continue close(9) do 3 ie=len(sn(1)),1,-1 3 if(sn(1)(ie:ie).eq.' ')goto 33 33 open(9,file=sn(1)(25:ie)) np=0 111 read(9,*,end=88,err=88)x,y np=np+1 goto 111 88 write(6,*)np,' points in '//sn(1)(25:ie) close(9) allocate(sx(np),s(np)) s=0.0d0 do 2 i=1,n do 4 ie=len(sn(i)),1,-1 4 if(sn(i)(ie:ie).eq.' ')goto 44 44 open(9,file=sn(i)(25:ie)) do 22 j=1,np read(9,*)sx(j),y 22 s(j)=s(j)+c(i)**ee*y 2 close(9) open(9,file='SUM.PRN') do 5 i=1,np 5 write(9,900)sx(i),s(i) 900 format(f12.4,e12.4) close(9) end