program bb parameter (n0=1000) real*8 as(n0,n0),ap(n0,n0),ad(n0,n0) character*80 s80 open(9,file='G') nsc=0 npc=0 ndc=0 ns=0 np=0 nd=0 do i=1,n0 do j=1,n0 as(i,j)=0.0d0 ap(i,j)=0.0d0 ad(i,j)=0.0d0 enddo enddo 1 read(9,60,end=999,err=999)s80 60 format(a80) write(6,60)s80 if(s80(2:3).eq.'S ')then ns=ns+1 read(s80(4:6),*)ncs do 3 i=1,ncs nsc=nsc+1 3 read(9,*)as(nsc,1),as(nsc,1+ns) endif if(s80(2:3).eq.'P ')then np=np+1 read(s80(4:6),*)ncp do 4 i=1,ncp npc=npc+1 4 read(9,*)ap(npc,1),ap(npc,1+np) endif if(s80(2:3).eq.'SP')then np=np+1 ns=ns+1 read(s80(4:6),*)ncp do 5 i=1,ncp npc=npc+1 nsc=nsc+1 read(9,*)ap(npc,1),as(nsc,1+ns) ,ap(npc,1+np) 5 as(nsc,1)=ap(npc,1) endif if(s80(2:3).eq.'D ')then nd=nd+1 read(s80(4:6),*)ncd do 6 i=1,ncd ndc=ndc+1 6 read(9,*)ad(ndc,1),ad(ndc,1+nd) endif goto 1 999 close(9) open(19,file='GN') write(19,191),nsc,ns,0 191 format('$ S-TYPE FUNCTIONS',/,3I5) do 10 i=1,nsc 10 write(19,88)(as(i,ii),ii=1,1+ns) 88 format(f16.9,6f12.9) write(19,192),npc,np,0 192 format('$ P-TYPE FUNCTIONS',/,3I5) do 12 i=1,npc 12 write(19,88)(ap(i,ii),ii=1,1+np) write(19,193),ndc,nd,0 193 format('$ D-TYPE FUNCTIONS',/,3I5) do 13 i=1,ndc 13 write(19,88)(ad(i,ii),ii=1,1+nd) close(19) stop end